تعداد نشریات | 43 |
تعداد شمارهها | 1,640 |
تعداد مقالات | 13,343 |
تعداد مشاهده مقاله | 29,979,776 |
تعداد دریافت فایل اصل مقاله | 12,001,483 |
الگوریتم بهینهسازی ژنتیک هیبرید تقویت شده با پایگاه داده مجازی جهت ارزیابی طرحهای توسعه تولید و شبکه انتقال | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
هوش محاسباتی در مهندسی برق | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقاله 6، دوره 5، شماره 2، تیر 1393، صفحه 79-90 اصل مقاله (298.63 K) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نوع مقاله: مقاله پژوهشی فارسی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نویسندگان | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
محمد صادق جوادی* 1؛ محسن صنیعی2؛ حبیب رجبی مشهدی3 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1دانشجوی دکتری، دانشکده مهندسی برق- دانشگاه شهید چمران- اهواز- ایران | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2استادیار، دانشکده مهندسی برق- دانشگاه شهید چمران- اهواز- ایران | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3استاد، دانشکده مهندسی برق- دانشگاه فردوسی- مشهد- ایران | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
چکیده | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
این مقاله، به معرفی یک روش بهینهسازی هیبرید به منظور تحلیل مسائل بهینهسازی مختلط با عدد صحیح میپردازد. مدل هیبرید یاد شده به صورت یک برنامه بهینهسازی دو سطحی، شامل مسئله اصلی و تعدادی زیرمسئله پیاده سازی شده است. در روش پیشنهادی از یک پایگاه داده مجازی به منظور ذخیرهسازی مجموعه حالات ارزیابی شده در حین اجرای الگوریتم بهینهسازی استفاده شده است. این پایگاه داده مجازی ضمن سرعت بخشیدن به روند شبیهسازی، امکان به اشتراک گذاشتن و استفاده از آن در برنامههای شبیهساز با بیش از یک پروسسور را نیز فراهم مینماید. مسئله برنامهریزی توسعه همزمان تولید و شبکه انتقال به عنوان یک مسئله بهینهسازی مختلط با عدد صحیح دینامیک جهت ارزیابی مدل بهینهسازی هیبرید پیشنهادی مورد توجه قرار گرفته است. نتایج شبیهسازی بر روی شبکههای تست نمونه نشان دهنده کارایی آن در حل مسئله برنامهریزی بلندمدت در سیستم قدرت میباشد. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
کلیدواژهها | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
الگوریتم ژنتیک؛ برنامهریزی توسعه؛ بهینهسازی مختلط با عدد صحیح؛ پایگاه داده مجازی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
اصل مقاله | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
بخش عمدهای از مسائل بهینهسازی در سیستمهای قدرت از نوع مسائل بهینهسازی مختلط با عدد صحیح هستند. برنامهریزی مشارکت واحدهای نیروگاهی، در افق کوتاه مدت، برنامهریزی تعمیرات واحدهای نیروگاهی و تجهیزات شبکه انتقال، در افق میانمدت و برنامهریزی توسعه تولید و شبکه انتقال در افق بلندمدت، از جمله این مسائل بهینهسازی بهشمار میروند [1]. الگوریتمهای بهینهسازی ریاضی و ابتکاری بیشماری به منظور تحلیل چنین مسائلی ارائه شدهاند و نرمافزارهای تجاری متعددی بدین منظور توسعه یافته و در دسترس هستند. اگرچه الگوریتمهای ریاضی در دستیابی به پاسخهای بهینه در چنین مسائلی به اندازه کافی توسعه یافتهاند، اما با این حال در تحلیل مسائل با افزایش ابعاد و غیرخطی شدن قیود و تابع هدف چندان کارآمد نیستند [2]. به علاوه، ماهیت دینامیک مسائل بهینهسازی یاد شده در سیستمهای قدرت، در افق برنامهریزی، خود بر پیچیدگی مسأله میافزاید [3-4]. در سالهای اخیر الگوریتمهای بهینهسازی ابتکاری متعددی به منظور تحلیل مسائل بهینهسازی مختلط با عدد صحیح ارائه شدهاند. الگوریتم ژنتیک [5]، جستجوی ممنوعه [6]، اجتماع ذرات [7]، کلونی مورچگان [8]، جستجوی هارمونی [9] و.... از جمله الگوریتمهای ابتکاری هستند که در مطالعات سیستم قدرت استفاده شدهاند. یکی از مهمترین معایب الگوریتمهای ابتکاری وابستگی شدید جوابها به انتخاب و تنظیم پارامترهای الگوریتم بهینهسازی مورد نظر است. در برخی موارد انتخاب پارامترها و تنظیم آنها به کمک روش سعی و خطا صورت میگیرد و در برخی دیگر نیز با بهرهگیری از آنالیز حساسیت میتوان پارامترهای الگوریتم بهینهسازی را تنظیم نمود. با این حال، برای دستیابی به پاسخ بهینه نهایی باید چندین مرتبه فرآیند بهینهسازی تکرار گردد [10]. در مسائل بهینهسازی با ابعاد بزرگ، نظیر برنامهریزی توسعه شبکه قدرت، میتوان مسأله مورد نظر را به چندین مرحله تقسیم نمود و سپس پیوستگی بین مراحل را به کمک یک مدل دینامیکی برقرار کرد. از این تکنیک در روشهای بهینهسازی مبتنی بر قاعده تفکیک نظیر روش Benders Decomposition یا Dantzing Wolf استفاده شده است [11-12]. در این مقاله یک مدل هیبرید شامل مدل بهینهسازی مختلط با عدد صحیح خطی، به منظور تحلیل زیر مسألههای برنامهریزی تولید، و الگوریتم بهینهسازی ابتکاری ژنتیک، به منظور تحلیل دینامیک مسأله برنامهریزی توسعه تولید و شبکه انتقال، ارائه شده است. به علاوه، یک پایگاه داده مجازی به منظور ذخیرهسازی حالات ارزیابی شده از فضای جستجوی دینامیک در کنار الگوریتم ژنتیک در نظر گرفته شده است تا از حجم محاسبات غیر ضروری کاسته شود. این موضوع در تحلیل برنامههای بهینهسازی مختلط با عدد صحیح که تعداد حالات ممکن آن قابل شمارش هستند بسیار کارا و مفید خواهد بود. از آنجا که تحلیل زیرمسألههای موجود در این مدل زمانبر هستند، در تکرار مجدد الگوریتم، به منظور تنظیم پارامترها و یا آنالیز حساسیت، سرعت محاسبات مجدد به شکل درخور توجهی بالا خواهد رفت. این موضوع، به طور خاص، در مسائل برنامهریزی که دارای زیرمسألههای زمانبر هستند، سودمند خواهد بود. در این مقاله، موضوع برنامهریزی توسعه تولید و شبکه انتقال، به عنوان یک مسأله بهینهسازی مختلط با عدد صحیح که ماهیتی دینامیکی در افق برنامهریزی نیز دارد، ارزیابی شده است. بخشهای مختلف مقاله به شرح زیر است: بخش بعد به معرفی برنامهریزی توسعه تولید و انتقال در سیستم قدرت میپردازد؛ بخش سوم به مدلسازی ریاضی مسأله برنامهریزی توسعه همزمان تولید و انتقال اختصاص یافته است؛ مبانی نظری الگوریتم بهینهسازی هیبرید و روندنمای حل مسأله برنامهریزی توسعه در بخش چهارم تشریح گردیده است؛ روش پیشنهادی در بخش پنجم با معرفی دو شبکه آزمون پذیرفته شده ارزیابی و نتایج آن با روش پیشنهادی در مقاله مرجع مقایسه شده است و بخش نهایی این مقاله به نتیجهگیری اختصاص یافته است.
1- برنامهریزی توسعه تولید و انتقالدر برنامهریزی بلندمدت، طرحهای توسعه واحدهای تولیدی و توسعه یا تقویت شبکه انتقال به منظور تأمین بار آینده سیستم قدرت تحلیل و ارزیابی میشوند. اهداف اصلی در این برنامهریزی تعیین زمانبندی احداث خطوط انتقال جدید، توسعه زیرساختهای لازم برای توسعه شبکه انتقال موجود، مکانیابی برای احداث واحدهای تولیدی جدید و در نهایت، برنامهریزی به مدار آمدن تجهیزات تولید و انتقال در افق برنامهریزی مورد نظر است [1]. از دیدگاه نهاد برنامهریز شبکه قدرت اضافه شدن تجهیزات مورد نیاز شبکه باید به نحوی متناسب با نیازمندیهای شبکه قدرت در آینده سیستم باشد و از نظر اقتصادی کمترین هزینه را به سیستم تحمیل نموده، امنیت و کفایت سیستم قدرت در حد مطلوب تضمین گردد [13]. ارزیابی طرحهای توسعه تولید و شبکه انتقال نیازمند ابزارهای قدرتمند برای بررسی و تحلیل راهبردهای مختلف توسعه شبکه از دیدگاه فنی، اقتصادی، قابلیت اطمینان، پایداری و ... است. در برنامهریزی بلندمدت، معمولاً، طرحهای توسعه و یا تقویت شبکه در فاز نخست، از دیدگاه اقتصادی و قابلیت اطمینان ارزیابی میشوند. در فازهای بعدی مطالعات سیستم معطوف به ارزیابی پایداری، سطوح حفاظت، تقویت یا تغییر آرایش پستها و سیستم کلیدزنی و ... میگردد. ارزیابی طرحهای توسعه در فاز نخست، خود شامل بخشهای متعدد و محاسباتی حجیم است. پیشبینی بار، پیشبینی قیمت سوخت، محدودیتهای زیست- محیطی، محدودیت منابع آب و ... به عنوان مهمترین و اصلیترین پارامترهای تصادفی در آینده سیستم مورد توجه نهاد برنامهریز سیستم قدرت قرار میگیرند. همان گونه که اشاره شد، اهداف اصلی در این برنامهریزی، تعیین زمانبندی احداث خطوط انتقال جدید، توسعه زیرساختهای لازم برای توسعه شبکه انتقال موجود، مکانیابی برای احداث واحدهای تولیدی جدید و در نهایت برنامهریزی به مدار آمدن تجهیزات تولید و انتقال در افق برنامهریزی مورد نظر است. از دیدگاه نهاد برنامهریز شبکه قدرت اضافه شدن تجهیزات مورد نیاز شبکه باید به نحوی با نیازمندیهای شبکه قدرت در آینده سیستم متناسب باشد و از نظر اقتصادی کمترین هزینه را به سیستم تحمیل کند و امنیت و کفایت سیستم قدرت در حد مطلوب تضمین گردد. برنامهریزی توسعه تولید به دنبال پاسخگویی به نیازمندیهای بار در افق مورد مطالعه بوده ، شامل مکانیابی، تخصیص منابع سوخت، انتخاب سایز، تکنولوژی، زمانبندی احداث و در نهایت، ورود نیروگاه به مدار است [16-14]. برنامهریزی توسعه شبکه انتقال از طریق احداث خطوط جدید یا ارتقای تجهیزات موجود، میزان توان قابل مبادله به طور ایمن بین مولدها و مصرفکنندگان شرکتکننده در بازار را افزایش میدهد [14-15]. لذا توسعه شبکه انتقال، رقابتی شدن بازار را افزایش میدهد [16]. از سوی دیگر، سرمایهگذاری در تجهیزات جدید انتقال، بسیار هزینهبر است؛ بههمین علت فقط در شرایطی باید انجام شود که توجیه اقتصادی داشته باشد. به منظور ارائه حداکثر رفاه اقتصادی به اجتماع، صنعت تأمین برق باید روند توسعه بلندمدت را درکمترین هزینه اتخاذ کند. این امر مستلزم روشی هماهنگ در بهینهسازی بهرهبرداری و توسعه تولید و انتقال است [17]. بهینهسازی مجزای شبکه انتقال و منابع تولید، به احتمال قریب به یقین، هدف مذکور را برآورده نمیسازد [22-21]. در یک سیستم قدرت متمرکز پتانسیلهای سرمایهگذاری در بخش تولید محدود و مشخص است. توسعه شبکه انتقال نیز عموماً تابعی از شرایط تولید و بار سیستم است و معمولاً زمانی که نقاط بار و تولید جدید به شبکه اضافه گردد، توسعه انتقال نیز در دستور کار قرار میگیرد. در این حالت نیز طرحهای توسعه شبکه انتقال و یا تقویت آن محدود بوده و سرمایهگذاری در این بخش ماهیتی استراتژیک نخواهد داشت. با این حال، برخی از راهکارهای پیشنهاد شده در سالهای اخیر بر این برنامهریزی اثرگذار خواهند بود. برای مثال، نشان داده شده است که توسعه تولید در محل با بهرهگیری از واحدهای تولید پراکنده[1]، DG [24-23]، و همچنین برنامههای مدیریت مصرف به حذف یا تعویق برخی طرحهای توسعه انتقال غیرضروری در شبکه قدرت منجر خواهد شد [18].
2- مدل سازی ریاضیدر این بخش از مقاله مدلسازی ریاضی مسأله برنامهریزی توسعه تولید و انتقال ارائه شده است. این مدل در قالب یک مسأله بهینهسازی مختلط با عدد صحیح به شکل (1) قابل توصیف است:
که در آن بخش اول و دوم به آنالیز اقتصادسنجی سرمایهگذاری نیروگاهها و خطوط انتقال جدید مربوط بوده و بخش سوم به هزینههای تأمین بار شبکه در سالهای تحت مطالعه اشاره دارد. قیود مسأله بهینهسازی به شرح زیر است:
قیود (2-4) به محدودیتهای تصمیمگیری در خصوص اضافه شدن واحدهای نیروگاهی اشاره دارد. قید (2) بیان میکند که تنها در صورتی یک واحد نیروگاهی امکان ورود به مدار را دارد که حداقل زمان مورد نیاز برای احداث، MTGIki، را پشت سر گذاشته باشد. به علاوه، در صورت ورود نیروگاه به مدار در کل دوره مطالعه این نیروگاه در دسترس خواهد بود؛ هر چند ممکن است بهرهبرداری نشود. قیود (3) و (4) به ترتیب به محدودیت بودجه و ظرفیت اضافه شونده تولید به شبکه در هر یک از سالهای مطالعه اشاره دارد.
مشابه واحدهای نیروگاهی کاندید، برای خطوط انتقال و تجهیزات آن نیز قیود زمانبندی دینامیکی در نظر گرفته میشود. به این ترتیب، قید (5) به حداقل زمان مورد نیاز برای احداث و به مدار آمدن تجهیزات انتقال توان کاندید به شبکه، قید (6) به محدودیت بودجه سرمایهگذاری سالیانه در بخش انتقال و قید (7) به محدودیت افزایش ظرفیت شبکه انتقال در سال تصمیم اشاره دارند. در مورد واحدهای تولیدی موجود و کاندید محدودیت میزان تولید متناسب با ظرفیت نامی نیروگاه در نظر گرفته میشود. این مورد توسط روابط (8) و (9) به صورت قیود مختلط با عدد صحیح ارائه شده است. در مورد واحدهای موجود، میزان تولید به وضعیت روشن و خاموش بودن نیروگاه بستگی داشته؛ حال آنکه برای واحدهای کاندید، در دسترس بودن نیروگاه نیز ملاک خواهد بود. به بیان دیگر، در قید (8) تنها یک متغیر باینری، IGEkiby، به چشم میخورد، ولی در مورد قید (9) دو متغیر باینری IGCkibyو Gnkiy در نظر گرفته شده است.
معادلات مربوط به توان عبوری از خطوط انتقال موجود، پخش توان DC، به شکل زیر هستند:
مدلسازی توان عبوری خطوط کاندید نیز مشابه با خطوط موجود در شبکه است؛ با این تفاوت که در روابط نظیر باید در دسترس بودن خط کاندید نیز مدلسازی گردد. بدین ترتیب روابط (12) تا (14) با لحاظ کردن متغیر باینری دوم؛ یعنی Lnjy، به صورت زیر بیان شده است:
مقدار MCj باید به میزان کافی بزرگ انتخاب گردد تا مدل برنامهریزی به شکل مختلط با عدد صحیح قابل پیادهسازی باشد. همچنین، قید (14) نیز به این موضوع اشاره دارد که توان عبوری از هر یک از کریدورهای موجود، در صورت در مدار بودن خط، باید بین مقادیر کمینه و بیشینه مجاز بارگیری از خط قرار داشته باشد. قید توازن تولید و مصرف در هر یک از نقاط مصرف به صورت (15) است. شایان ذکر است که در این رابطه PDlbyمیزان توان مورد نیاز در شین lام، به بلوک زمانی bام در سال yام مربوط است.
زاویه شین مرجع برای اجرای برنامه پخشبار ”0“ در نظر گرفته میشود که این موضوع در (16) مورد تأکید قرار گرفته است.
3- مدل هیبرید بهینهسازی مختلط با عدد صحیح دینامیکالگوریتم هیبرید پیشنهادی شامل دو مسأله بهینهسازی به صورت ترکیبی است. به کمک الگوریتم بهینهسازی ژنتیک، در هر تکرار یک رشته عدد صحیح تولید میگردد. این رشته عددی، بردار تصمیم تولیدی توسط الگوریتم ژنتیک است. ارزیابی این رشته عددی صحیح به منظور بازگرداندن مقدار عددی تابع هدف توسط الگوریتم ژنتیک صورت میگیرد. همچنین، بازتولید رشتههای عددی جدید در تکرارهای بعدی به کمک الگوریتم ژنتیک و بر اساس مقادیر عددی تابع هدف و تکرارهای قبلی صورت میگیرد. ارزیابی رشته عددی صحیح به منظور تعیین مقدار عددی تابع هدف، در یک الگوریتم بهینهسازی قطعی صورت میپذیرد. با توجه به این موضوع که تابع هدف در برنامهریزی بلندمدت معرفی شده بر کمینهسازی هزینه سرمایهگذاری و بهرهبرداری استوار شده است؛ لازم است در هر یک از سالهای افق برنامهریزی، هزینه بهرهبرداری از طریق حل معادله پخش بار بهینه محاسبه گردد و با هزینه سرمایهگذاری که نگاشتی از رشته عدد صحیح تولیدی توسط الگوریتم ژنتیک در مطالعه توسعه سیستم است، جمع گردد. برای محاسبه هزینههای بهرهبرداری، میتوان از روشهای بهینهسازی ریاضی به منظور محاسبات پخش بهینه توان در شبکه بهره جست. در این مقاله، محاسبات پخش توان با کمک نرمافزار GAMS صورت گرفته است. کارایی روش حل به کار گرفته شده در محیط GAMS نسبت به نرمافزار MATPOWER برای مسائل پخش بار بهینه در شبکههای قدرت در ابعاد بزرگ در مقاله [19] تأیید شده است. از این رو، در این مطالعه از تکنیک حل مسأله پخش بار یاد شده در محیط GAMS بهره گرفته شده است. فرض نمایید تعداد کاندیدهای پیشنهاد شده برای توسعه تولید G نیروگاه و تعداد خطوط کاندید T مورد باشد. با فرض در نظر گرفتن افق برنامهریزی برای Y سال، فضای جستجوی دارای ابعاد *Y (G+T) خواهد بود. هدف برنامهریزی توسعه شبکه، تعیین زمان ورود هر یک از تجهیزات تولید یا انتقال کاندید به مدار است. از این رو، یک متغیر باینری به منظور تعیین وضعیت هر یک از تجهیزات کاندید نیاز است تا مشخص گردد که کدامیک از تجهیزات کاندید در افق برنامهریزی در دسترس بهرهبردار قرار خواهند داشت. این متغیر باینری در حقیقت همان Gnkiy و Lnjy هستند که در فرمولاسیون مسأله به آنها اشاره شده است. به این ترتیب فضای جستجو در این حالت، *Y (G+T)2 خواهد بود. این فضای جستجو برای مسائل با ابعاد بزرگ غیرقابل تصور خواهد بود. از این رو، لازم است راهکاری برای حذف فضای جستجوی غیرمؤثر مورد توجه قرار گیرد. از آنجا که افق برنامهریزی نسبت به عمر مفید تجهیزات کمتر در نظر گرفته میشود، در صورت ورود تجهیزات کاندید به مدار تا آخر دوره برنامهریزی در دسترس خواهند بود. این موضوع به کاهش فضای جستجو به رقم (G+T)(Y+1) منجر خواهد شد که در نظر گرفتن حداقل زمان مورد نیاز برای ساخت هر یک از تجهیزات نیز به کاهش بیشتر فضای جستجوی مفید منجر خواهد شد. شکل (1) نگاشتی از یک رشته عددی صحیح را که توسط الگوریتم ژنتیک تولید شده است، به فضای تفکیک شده در افق برنامهریزی نشان میدهد. فرض کنید تعداد واحدهای تولیدی کاندید 4 مورد، {G4, G5, G6, G9} و تعداد خطوط انتقال کاندید نیز 4 مورد، {T.8, T.9, T.10, T.11} باشند. با فرض در نظر گرفتن افق برنامهریزی 10 ساله، فضای پاسخ یک ماتریس 8*10 خواهد بود که درایههای آن میتوانند مقادیر ”صفر“ یا ”یک“ را به خود اختصاص دهند. با توجه به تکنیک توصیف شده در شکل (1)، رشته عددی صحیح به صورت اعداد باینری رمزگشایی میگردد. نحوه رمزگشایی به این ترتیب است که عدد صحیح تصادفی متناظر با هر یک از تجهیزات در افق برنامه تعیین کننده وضعیت هر تجهیز تا انتهای افق برنامهریزی خواهد بود.
شکل (1): نحوه رمزگشایی رشته عدد صحیح تولیدی
برای مثال، در رشته عددی تصادفی تولید شده در برنامه ژنتیک، عدد صحیح متناظر با واحد G4 رقم ”0“ و این عدد برای واحد G5، برابر با ”6“ است که مفهوم آنها این است که واحد G4 در این طرح توسعه انتخاب نشده و واحد G5، از سال ”ششم“ تا انتهای افق برنامهریزی در دسترس خواهد بود؛ هر چند بهرهبرداری از آن تابع نظر بهرهبردار سیستم خواهد بود. به این ترتیب، برای هر یک از ستونها در افق برنامهریزی Y سال، (Y+1) حالت ممکن وجود خواهد داشت ({0, 1, 2, …Y}). نظر به اینکه (G+T) تجهیز به عنوان کاندید مورد توجه هستند و توسعه آنها مستقل از یکدیگر است، تعداد حالتها در این مورد از 280 به 118 کاهش خواهد یافت؛ ضمن آنکه محاسبات مربوط به هر یک از سالها، در پایگاه داده تعبیه شده ذخیره میگردد تا از محاسبات غیرضروری و اضافه در تکرارهای بعدی جلوگیری شود؛ برای مثال، در صورتی که آرایه تولیدی، نسبت به حالت نشان داده شده در شکل (1) در سال ورود واحد G5 تفاوت داشته باشد؛ مثلاً سال ورود آن در سال ”پنجم“ پیشنهاد شده باشد، نیازی به انجام کلیه محاسبات نخواهد بود، بلکه تنها کافی است از رشته جدید تولیدی رمزگذاری شده فقط محاسبات سال ”پنجم“ صورت پذیرد و به پایگاه داده اضافه گردد. شایان ذکر است که بقیه حالات در تکرار قبلی محاسبه و در پایگاه داده ذخیره شدهاند و تنها از آنجا فراخوانی میگردند. روندنمای اجرای این برنامه در شکل (2) به تصویر کشیده شده است. از آنجا که ارزیابی هر یک از طرحها فرآیندی زمانبر است و در برنامه ژنتیک احتمال آنکه یک رشته عددی چندین بار آنالیز شود وجود دارد؛ یک پایگاه داده مجازی در کنار الگوریتم ژنتیک در نظر گرفته شده است تا در صورتی که اطلاعات هزینه هر یک از طرحها در پایگاه موجود باشد، زمانی برای محاسبه مجدد آن طرح صرف نشود. این قابلیت خصوصاً در تکرارهای بعدی برای تنظیم بهینه پارامترهای الگوریتم ژنتیک، زمان محاسبات را به شکل درخور توجهی کاهش میدهد.
4- بررسی نتایج شبیهسازیبه منظور شبیهسازی الگوریتم پیشنهادی، دو مورد مطالعاتی که در مرجع [20] به منظور ارزیابی برنامهریزی همزمان توسعه تولید و انتقال مطرح شدهاند، مطالعه و بررسی شدهاند: مورد مطالعاتی نخست، یک شبکه تست با 6 شین و مورد دوم شبکه 118 شین اصلاح شده استاندارد است. شایان ذکر است که در هر یک از موردهای مطالعاتی، تعداد جمعیت اولیه 200، تعداد تکرارها 100، تعداد متغیرها برابر با تعداد تجهیزات کاندید انتخاب شدهاند. محدوده متغیرهای عدد صحیح نیز بین ”صفر“ تا ”افق برنامهریزی“ انتخاب شده است، که افق برنامهریزی در هر دو مورد مطالعاتی، 10 سال در نظر گرفته شده است.
شکل (2): روندنمای الگوریتم هیبرید پیشنهادی با بهرهگیری از پایگاه داده مجازی
5-1- مورد مطالعاتی 6 شین دیاگرام تکخطی شبکه مورد نظر در شکل (3) و اطلاعات مربوط به خطوط انتقال موجود و کاندید در جدول (1) ارائه شده است. اطلاعات مربوط به نیروگاههای موجود و کاندید در جدول (2) و پیشبینی بار پیک برای 10 سال آینده در جدول (3) ارائه شده است. جدول (4) ضرایب مورد نیاز برای استخراج منحنی تداوم بار سیستم تحت مطالعه را نشان میدهد. شایان ذکر است که این ضرایب برای سال مرجع داده شده و برای سالهای آتی نیز ضرایب متناسب با سال مرجع در نظر گرفته شده است. سهم بار شین شماره 3 معادل 40% کل بار و سهم دو شین دیگر؛ یعنی شین 4 و 5 برابر 30% کل بار در نظر گرفته شده است. حداقل زمان مورد نیاز برای به مدار آمدن نیروگاههای کاندید، سه سال در نظر گرفته شده، در حالی که این زمان برای خطوط انتقال کمتر از یک سال فرض شده است. نتایج عددی بهدست آمده با بهرهگیری از روش پیشنهادی به همراه نتایج ارائه شده در مرجع [20] به صورت مقایسهای در جدول (5) ارائه شده است. هزینه کل، شامل هزینه سرمایهگذاری و بهرهبرداری در افق مورد برنامهریزی در مقاله یاد شده 6/366 میلیون دلار بهدست آمده است، در حالی که هزینه کل با بهرهگیری از روش هیبرید به 06/342 میلیون دلار کاهش یافته است.
جدول (1): اطلاعات خطوط انتقال موجود و کاندید-6شین
شکل (3): دیاگرام تکخطی شبکه مورد مطالعه [20]
جدول (2): اطلاعات واحدهای نیروگاهی موجود و کاندید 6شین
جدول (3): پیک بار پیشبینی شده
جدول (4): بلوکهای منحنی تداومی بار (برای سال اول)
شایان ذکر است که در این مطالعه، به منظور فراهم آوردن شرایط مقایسه، نرخ تنزیل سالیانه ”صفر“ در نظر گرفته شده و هیچ محدودیت بودجه یا ظرفیت برای تجهیزات منظور نشده است. این موضوع به تولید چندین طرح بهینه با هزینههای مشابه منجر خواهد شد. تعیین نقاط بهین متعدد یا شبهبهین در مسائل بهینهسازی تنها با بهرهگیری از روشهای بهینهسازی ابتکاری میسر است؛ حال آنکه روشهای کلاسیک ریاضی قادر به شناسایی چنین ویژگی نیستند و در تکرارهای بعدی پاسخ بهدست آمده مجدداً تکرار خواهد شد. با توجه به این موضوع که نرخ تنزیل سالیانه صفر منظور شده است و اضافه شدن واحدهای نیروگاهی در سالهای نخست برنامهریزی عملاً در کاهش هزینههای بهرهبرداری مؤثر واقع نشدهاند، لذا تقدم یا تأخر زمانی آنها بر برنامهریزی بهینه اثربخش نبوده است. جدول (5): نتایج مقایسهای الگوریتم هیبرید و مرجع [20]
این موضوع بیانگر توانایی الگوریتم پیشنهادی در دستیابی به پاسخی بهینهتر نسبت به مقاله مرجع [20] است که در آن روش تجزیه بِندِر [ii] استفاده شده است. آنالیز مقایسهای دوم مربوط به ارزیابی تأثیر پایگاه داده مجازی بر زمان اجرای محاسبات در روش پیشنهادی است. شکل (4) زمان اجرای 15 مرتبه تکرار شبیهسازی با در نظر گرفتن و بدون در نظر گرفتن اثر پایگاه داده را نشان میدهد. همچنین روند تکمیل پایگاه داده برای 15 مرتبه تکرار در شکل (5) به تصویر کشیده شده است.
شکل (4): تأثیر پایگاه داده بر زمان اجرای شبیهسازی
شکل (5): روند تکمیل پایگاه داده در تکرارهای متوالی
5-2- مورد مطالعاتی 118 شین اصلاح شده به منظور ارزیابی نتایج عددی در یک شبکه قدرت بزرگ، سیستم 118 شین اصلاح شده معرفی شده در مرجع [20]، نیز بررسی شده است. اطلاعات خطوط انتقال کاندید و واحدهای نیروگاهی به ترتیب در جداول (6) و (7) ارائه شدهاند.
جدول (6): اطلاعات خطوط انتقال کاندید-118 شین
در این مطالعه نرخ رشد بار معادل با 7% در نظر گرفته شده است. در حالت پایه، مشابه با حالت پایه مرجع [20]، نرخ تنزیل سالیانه 10% لحاظ شده است تا شرایط برای مقایسه عددی هر دو طرح فراهم باشد. ظرفیت نصب شده نیروگاهی 5850 مگاوات است. حداقل زمان لازم برای احداث تجهیزات کاندید پیشنهادی سه سال در نظر گرفته شده است. از این رو، سیستم قدرت مورد نظر برای سه سال نخست برنامهریزی هیچ تغییری از نظر توپولوژی و ظرفیت نیروگاهی نخواهد داشت. نتایج مربوط به طرح نهایی توسعه به صورت مقایسهای در جدول (8) ارائه شده است.
جدول (7): اطلاعات واحدهای نیروگاهی کاندید- 118 شین
با توجه به پرشدگی خطوط 41، 128، 129، 153 و 159 در سه سال نخست و ادامه این روند با افزایش تقاضای بار در سیستم، امکان تأمین بار در سالهای آتی وجود نخواهد داشت. با ورود واحدهای 1 و 5 در سال ”چهارم“، بخشی از پرشدگی خطوط انتقال برطرف میگردد. در سال ”پنجم“، با ورود واحد کاندید شماره 2 در شین 12، ضمن افزایش ظرفیت نیروگاهی ارزان قیمت، پرشدگی مربوط به خطوط 54 و 129 به طور کامل برطرف میشود. با توجه به این موضوع که خط انتقال کاندید شماره 1، لینک بین شینهای 30 و 38، نیز میتوانست در طرح توسعه پرشدگی خط 54، که بین شینهای 30 و 38 از پیش نصب شده است را بر طرف نماید، اما توسعه تولید در شین 12 در آینده سیستم مؤثرتر خواهد بود. وضعیت تأمین بار سیستم تا سال ”هشتم“ با کمک واحدهای از پیش احداث شده بدون مشکل ادامه مییابد. با این حال، در این سال ورود واحد 100 مگاواتی شماره 11، پرشدگی خطوط 23 و 37 را تا انتهای سال ”نهم“ بر ظرف مینماید. ظرفیت نصب شده در انتهای سال ”هشتم“ 6550 مگاوات خواهد شد و ظرفیت رزرو بالقوه موجود در سیستم در این سال، در حدود 655 مگاوات است. با این حال، در دو سال انتهایی، نیاز به ظرفیت نیروگاهی جدید شدیداً احساس میشود. به این ترتیب، در ابتدای سال ”نهم“ چهار واحد نیروگاهی 100 مگاواتی به شبکه اضافه خواهند شد و در نهایت، در سال ”دهم“، علاوه بر واحدهای 6، 12، 13 و 14 نیاز است خط کاندید شماره 1 نیز به مدار اضافه گردد تا پرشدگی خط 37 نیز در این سال به طور کامل بر طرف گردد.
جدول (8): نتایج مقایسهای الگوریتم هیبرید و مرجع [20]
در مقایسه با نتایج ارائه شده در مقاله مرجع، هر چند هزینه سرمایهگذاری در آن طرح کمتر است، با این حال در طرح پیشنهادی با بهرهگیری از روش هیبرید، هزینه کل سیستم، نسبت به هزینهای که از طرح پیشنهادی در مقاله مرجع قابل محاسبه است، کمتر خواهد شد. از این رو، طرح پیشنهادی مقبولیت بیشتری نسبت به طرح ارائه شده در مقاله مرجع خواهد داشت. شایان ذکر است که با توجه به ارائه طرح نهایی در مقاله مرجع میتوان هزینه کل سیستم را برای افق 10 سال محاسبه نمود. این رقم برای طرح پیشنهادی در مقاله مرجع 4644099954 دلار و در روش هیبرید پیشنهادی، 4634617644 دلار بهدست آمده است.
5- نتیجهگیریدر این مقاله یک روش هیبرید به منظور تحلیل مسائل بهینهسازی مختلط با عدد صحیح ارائه گردید. اضافه شدن پایگاه داده در کنار مسأله بهینهسازی ابتکاری در سطح نخست، قابلیت این تکنیک شبیهسازی را در یافتن سریع نقاط بهینه به شکل قابل توجهی بالا میبرد. مهمترین اثر پایگاه داده سرعت بخشیدن به فرآیند شبیهسازی در حین اجرای تکرار نخست و همچنین، در اجراهای بعدی است. به علاوه، در برنامههای بهینهسازی در ابعاد بسیار بزرگ امکان پردازش به کمک چندین پروسسور را نیز فراهم میآورد. نتایج شبیهسازی بر روی مسأله برنامهریزی توسعه در شبکه قدرت، شامل برنامهریزی توسعه تولید و انتقال، که در اصل یک مسأله بهینهسازی مختط با عدد صحیح دینامیک است، نشان میدهد که الگوریتم مورد نظر نه تنها به پاسخهای بهینهتری نسبت به روشهای کلاسیک ریاضی منتهی شده است، بلکه از نظر زمان اجرای برنامه نیز کارایی بهتری داشته است. لیست نمادها و علائم
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مراجع | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
[1] Mazer, A. Electric Power Planning for Regulated and Deregulated Markets, New Jersey: John Wiley & Sons, 2007. [2] Goldberg, D. E., Genetics Algorithms in Search, Optimization and Machine Learning, Reading, MA: Addison-Wesley, 1989. [3] Binato, S., Oliveira, G. C., "Dynamic planning applied to transmission network expansion", 9th Automation Brazilian Congress, Vitória, Brazil, 1992. [4] Haffner, S., Garcia, A., Monticelli, A., Romero, R. A., "Dynamic power transmission expansion planning considering multiple stages," Third International Conference on Electric Utility Deregulation and Restructuring and Power Technologies, pp. 635-640, Thailand, 2008. [5] Sheble, G. B., Charles, W., Richter, J., "A Profit-Based Unit Commitment GA for the Competitive environment", IEEE Trans. Power Syst., Vol. 15, No. 2, pp. 715 - 721, May 2000 [6] Gallego, R. A., Romero, R. A., Monticelli, A., "Tabu search algorithm for network synthesis", IEEE Trans. Power Syst., Vol. 15, No. 2, pp. 490-495, May 2000. [7] Swarup, K. S., Kumar, P. R., "A new evolutionary computation technique for economic dispatch with security constraints", International Journal of Electrical Power and Energy Systems, Vol. 28, No. 4, pp. 273-283, 2006. [8] Song, Y. H., Chou, C. S., Stonham, T. J., "Combined heat and power dispatch by improved ant colony search algorithm.," Electric Power Systems Research, Vol. 52, pp. 115-121, 1999. [9] Afkousi-Paqaleh, M., Rashidinejad, M., Pourakbari-Kasmaei, M., "An implementation of harmony search algorithm to unit commitment problem", Electrical Engineering Journal, Vol. 92, No. 6, pp. 215-225, 2010. [10] Deb, K., Pratap, A., Agarwal, S., Meyarivan, T., "A fast and elitist multiobjective genetic algorithm: NSGA-II", IEEE Trans. Evolut. Comput., Vol. 6, No. 2, pp. 182-197, 2002. [11] Binato, S., Pereira, M. V., Granville, S., "A new benders decomposition approach to solve power transmission network design problems", IEEE Trans. Power Syst., Vol. 16, No. 2, pp. 235-240, May 2001. [12] Dantzing, A. Wolfe, P., "Decomposition principle for linear programs", Operation Research, Vol. 8, No. 1, pp. 101-111, 1960. [13] Seifi, H., Sepasian, M. S., Electric Power System Planning, Issues, Algorithms and Solutions. Verlag Berlin Heidelberg: Springer, 2011. [14] Fischer, R., Joo, S. K., "Economic Evaluation of Transmission Expansion for Investment Incentives in a Competitive Electricity Market", International Journal of Control, Automation, and Systems, Vol. 6, No. 5, pp. 627-638, Oct. 2008. [15] Rocha, M. C. D. Saraiva, J. T., "Discrete Evolutionary Particle Swarm Optimization for Multiyear Transmission Expansion Planning", 17th Power Systems Computation Conference, Stockholm Sweden, 2011. [16] Fang, R., Hill, D. J. "A New Strategy for Transmission Expansion in Competitive Electricity Markets", IEEE Trans. Power Syst.,Vol. 18, No. 1, pp. 374-380, Feb. 2003. [17] Saboori, H., Mohammadi, M., Taghe, R., "Composite Generation and Transmission Expansion Planning Considering the Impact of Wind Power Penetration", Power and Energy Engineering Conference (APPEEC), 2011. [18] Shyesteh, E., ParsaMoghadam, M., Yousefi, G. R., "An Economic Comparison between Incorporation of FACTS Devices and Demand Response Programs for ATC Enhancement", 16th Power Systems Computation Conf. (PSCC), Glasgow, U.K, 2008. [19] Javadi, M. S., Javadinasab, A., Shishebori, A., Basri, H., Azami, R., "New Approach for LMP Calculation in Large Scale Power System Based on Network Incident Matrix", International Review on Modelling and Simulations (IREMOS), Vol. 4, No. 5, pp. 2462-2469, 2011. [20] Khodaei, A., Shahidehpour, M., Kamalinia, S., "Transmission Switching in Expansion Planning", IEEE Transactions on Power Systems, Vol. 25, No. 3, pp. 1722-1733, 2010.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
آمار تعداد مشاهده مقاله: 789 تعداد دریافت فایل اصل مقاله: 350 |