1- مقدمه^[1]

در دهه‏های اخیر یکی از مهم‌ترین منابع تولید برق، مربوط به واحدهای حرارتی بوده است. منابع اصلی تأمین این نیروگاه‏ها استفاده از سوخت‏های فسیلی است؛ اما با توجه به فناپذیری و همچنین میزان آلایندگی و هزینۀ زیاد استفاده از آنها پژوهشگران را در پی کشف روش‏هایی برای بهینه‌کردن استفاده از این سوخت‏های فسیلی وا‌داشت ]1[. در ابتدا هزینۀ سوخت، متغیر اصلی در توزیع اقتصادی در نظر گرفته می‌شد و با گسترش قوانین زیست‌محیطی، انتشار آلاینده‏ها بخشی از هزینه برای توزیع اقتصادی شد. سپس توزیع اقتصادی، مشکل چندمنظوره برای به حداقل‌رساندن هزینه و انتشار آلاینده‏ها شد ]2[.

هزینه تولید، به‌ویژه در نیروگاه‌های حرارتی، بیش‌ازحد است؛ بنابراین برنامه‌ریزی مناسب خروجی واحدهای نیروگاهی به صرفه‌جویی در هزینه‌های عملیاتی کمک چشمگیری می‌کند ]3[. بهره‏برداری مناسب از سیستم قدرت، برای برگشت سود سرمایه‏گذاری، بسیار مهم است و هزینه‏هایی که ارگان‏های دولتی تعیین می‏کنند و همچنین اهمیت صرفه‌جویی در سوخت و هزینه آن، شرکت‏های برق را برای به دست آوردن حداکثر بازدهی ممکن، وادار می‏کند ]4[.

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

در دهه گذشته متدهای ابتکاری بسیاری با استفاده از هوش مصنوعی برای حل مشکل پیچیده توزیع نیرو به کار رفته است؛ ازجمله الگوریتم‏های بهینه‌سازی‏ ژنتیک[4] (GA) ]7[، انتشار ذرات[5] ]2[، کلونی زنبورعسل[6] (ABC) ]8،9[، ازدحام ذرات[7](PSO)  ]10[، کرم شب‏تاب[8] (FA) ]11[، کلونی مورچه‏ها[9] (ACO) ]12[، خفاش بی‌نظم[10] ]13[. برای بهبود کارایی جستجوی این متدها، تعدادی الگوریتم تغییر داده‌ شده که از مدل اصلی آنها مشتق شده‏اند نیز پیشنهاد شده است که الگوریتم‌های بهینه‌سازی ژنتیک بهبودیافته[11] ]14[، روش ‏ تکامل‌یافتۀ ازدحام ذرات[12] ( ) ]3[، کلونی زنبورعسل مصنوعی بهبودیافته[13] ]4[ را شامل می‌شوند. این تکنیک‌های مدرن با محدودیت‌هایی در فرم منحنی هزینه مواجه نیستند؛ اما پارامترهای وابسته زیادی دارند که باید به‌درستی تعیین شوند.

الگوریتم جستجوی عقب‌گرد یا ردیابی بازگشتی[14] (BSA) ]15[، الگوریتم تکاملی[15] (EA) است که برای حل مسائل بهینه‏سازی استفاده می‏شود. ساختار این الگوریتم ساده است و تنها یک پارامتر کنترلی دارد که باید تعیین شود. برای بهبود عملکرد همگرایی و گسترش دامنۀ استفاده از آن، با الگوریتم آموزش و یادگیری[16] (TLBO) ]16[ ترکیب شد و الگوریتم جدیدی به نام الگوریتم بهینه‏سازی ‏یادگیری ردیابی بازگشتی[17] (LBSA) ]17[ به دست آمد که اساس بهینه‌سازی در این مقاله است.

در دنیای امروز و در مراکز تولید، پخش بار اقتصادی، مهم‌ترین هدف توزیع انرژی بوده است و نقش مهمی در عملکرد اقتصادی سیستم قدرت دارد. در صنعت، زمانی که ژنراتورهای دارای بارهای مختلف، به‌هم‌پیوسته شوند، ظرفیت تولیدی آنها بسیار بزرگ‌تر از بارها می‏شود؛ به همین دلیل، اختصاص‌دادن بار برای ژنراتورها می‏تواند متنوع شود و به سبب اینکه کاهش هزینۀ تولید برق مهم است، تقسیم بار اقتصادی، مدنظر قرار می‏گیرد.

هزینۀ سوخت، بیشتر هزینۀ تولید در نیروگاه‌ها را به خود اختصاص می‏دهد. سایر هزینه‌ها مانند هزینۀ کار، تعمیر و نگهداری، عوامل اقتصادی‌اند؛ درنتیجه، متخصصان توزیع انرژی باید کنترل تولید نیروگاه‏ها را برای کم‌کردن هزینه به دست گیرند که در این صورت مسئلۀ پخش بار اقتصادی به‌عنوان مسئلۀ بهینه‏سازی، شامل تابع هدف و قیود، بیان می‌شود ]3،2[.

در این مقاله برای حل مسئلۀ پخش بار اقتصادی برای اولین بار از الگوریتم بهینه‏سازی یادگیری ردیابی بازگشتی (LBSA) استفاده ‌شد. ابتدا به بیان مسئلۀ پخش بار اقتصادی و فرمول‌بندی آن پرداخته و در ادامه، الگوریتم بهینه‏سازی یادگیری ردیابی بازگشتی (LBSA) و شیوۀ کار آن بررسی شده است. سپس با استفاده از این الگوریتم مسئلۀ پخش بار اقتصادی، برای سیستم‏های آزمایشیِ ذکرشده، حل شده است و نتایج حاصل، تحلیل و بررسی و با نتایج حاصل از سایر الگوریتم‏ها، مقایسه و نتیجه‌گیری شده‌‌اند.

2- فرمول‌بندی مسئله

2-1- توصیف مسئلۀ پخش بار

پخش بار اقتصادی در مراکز مدرن کنترل انرژی و به‌منظور تعیین توان بهینه واحدهای تولیدی، یکی از مسائل مهم در بهره‌برداری از سیستم‌های قدرت است. با توجه به بازۀ زمانی، دو نوع پخش بار اقتصادی معرفی می‌شود: پخش بار اقتصادی استاتیکی که در آن برای یک میزان بار مشخص و در یک بازۀ زمانی، مسئله را به‌صورت بهینه حل می‌کند و پخش بار اقتصادی دینامیکی که مسئله را برای میزان بارهای مختلف و در چندین بازۀ زمانی بررسی می‌کند. بدیهی است فرایند محاسبات در پخش بار دینامیکی، پیچیده‌تر از حالت استاتیکی است؛ اما نتایج محاسبات آن تطابق بیشتری با نیازهای واقعی سیستم قدرت دارد ]18[. پخش بار اقتصادی دینامیکی، مسئلۀ بهینه‏سازی برای تعیین برنامۀ زمان‏بندی توان خروجی واقعی واحدهای تولید، با در نظر گرفتن تعادل توان واقعی با بار مصرفی و همچنین محدودیت‏های خروجی ژنراتورها است.

تورم سالانه و افزایش قیمت مواد سوختی باعث می‏شود همواره به بهره‏برداری مناسب از سیستم‏های تولید انرژی الکتریکی توجه شود ]1[. هدف از توزیع اقتصادی، به حداقل رساندن مجموع هزینۀ تولید سیستم قدرت در یک بازۀ زمانی تعریف‌شده (معمولاً 1 ساعت) است؛ درحالی‌که محدودیت‌های عملیاتی سیستم قدرت رعایت شود؛ بنابراین، می‌توان این هزینه را به‌صورت ریاضی، مانند یک مسئلۀ بهینه‌سازی با تابع هدف و دو محدودیت (معادله و نامعادله) فرمول‌بندی کرد ]3[.

2-2- تابع هدف بدون اثر شیر بخار

عامل اصلی در توزیع اقتصادی بار، تابع هزینۀ ژنراتورها است و هزینه‏های دیگر به‌صورت درصد ثابتی از هزینۀ سوخت در نظر گرفته می‏شوند. بدیهی است تابع هزینه، تنها برای نیروگاه‏های حرارتی مفهوم دارد.

معمولاً تابع هزینۀ سوخت برای مشخص‌کردن توان تولیدی در یک نیروگاه را با تقریب مناسب، با تابع درجه 2 مشخص می‏شود؛ درنتیجه، برای تولید انرژی الکتریکی ( ) در نیروگاه  ام، هزینه‏ای به‏دست می‏آید که با  مشخص می‏کنند. واضح است که هزینۀ کل سیستم، برابر با مجموع هزینه‏های همه واحدها است. در این صورت، مسئلۀ توزیع اقتصادی بار به‌صورت رابطه (1) تعریف خواهد شد:

در این رابطه  هزینۀ سوخت کل نیروگاه، ،  و  ضرایب مربوط به هزینۀ واحد نیروگاهی  ام و  تعداد نیروگاه‏های سوخت فسیلی موجود در شبکه است.

در تحلیل پخش بار اقتصادی، با ایجاد تغییر در توان تولیدی هر یک از نیروگاه‏ها، باید روشی را انتخاب کرد که تابع هزینۀ ارائه‌شده به‌صورت رابطه (2) مینیمم شود ]19[.

2-3- تابع هدف با اثر شیر بخار

مجموع هزینۀ تولید معمولاً با همان تابع درجه دوم از توان خروجی واحد تولیدی تقریب زده می‌شود. ازطرفی معمولاً در پخش بار اقتصادی تابع هزینه دارای نقاط زانویی مشتق‌ناپذیر است ]19[ که به دلیل در نظر گرفتن بارگذاری نقطۀ دریچه (تخلیه) ]3[ و یا در نظر گرفتن اثرات موقعیت شیر ورودی ]4[، تابع هزینۀ نیروگاه را به‏صورت حاصل جمع تابع هزینه درجه 2 و قدر مطلق تابع سینوسی در نظر می‏گیرند.

در این مقاله، مدل والتر- شب ]20[ برای نشان‌دادن این خاصیت‌ها استفاده می‌شود؛ بنابراین، تابع هزینه اصلاح و تابع سینوسی تصحیح‌شده در تابع درجه دوم گنجانده می‌شوند. با توجه به این اثر، تابع هدف توزیع اقتصادی به‌عنوان یک مسئله بهینه‌سازی محدود به‌صورت رابطه ریاضی (3) توصیف می‌شود.

که  و  ضرایب مربوط به اثر شیر بخار ژنراتورها هستند. اگر بارگذاری نقطۀ دریچه یا شیر ورودی در نظر گرفته شود و به دنبال آن، واحد تولیدی تابع هزینه‌ای غیرصاف (غیریکنواخت) داشته باشد؛  برابر با یک و در غیر این صورت  برابر صفر و تابع هزینۀ این واحد صاف (یکنواخت) است ]3[.

2-4- محدودیت‏های معادله (تساوی)

محدودیت‏هایی مانند تعادل توان یا محدودیت تقاضا را با عنوان محدودیت‏های معادله (تساوی) مشخص می‏کنند؛ به این صورت که مجموع توان تولیدی سیستم ( ) باید مانند رابطه (5) با حاصل جمع بارهای کل سیستم ( ) و تلفات ( ) برابر باشد ]3[.

در این رابطه، تلفات خطوط به‏وسیلۀ تابعی از توان حقیقی و ماتریس ضرایب B که فرمول تلفات Korn هستند، با استفاده از معادله (6) به دست می‌آید.

که در آن پارامترهای  ضرایب اتلاف نام دارند. این رابطه را به‏صورت برداری با رابطه (7) بیان می‌شود.

ضرایب اتلاف را با ماتریس متقارن B نشان داده و  بردار ستونی و  ثابت اسکالر این ضرایب است ]21[.

شکل (1) تابع هزینۀ واحد حرارتی را با اثر شیر بخار و بدون اثر آن نشان می‌دهد.

شکل (1): تابع هزینۀ واحد حرارتی ]4[

2-5- محدودیت‌های نامعادله (نامساوی)

محدودیت‏های نامعادله یا نامساوی شامل محدودیت ظرفیت تولید، محدودیت سرعت تولید، میزان‌های افزایشی کاهشی، منطقۀ عملیاتی ممنوع، محدودیت‏های قابلیت اطمینان و امنیت شبکه و ... است.

محدودیت ظرفیت تولید: قدرت خروجی هر ژنراتور باید بین محدودۀ خودش باشد؛ یعنی ازطرفی نباید بیشتر از مقدار نامی و از طرف دیگر، نباید کمتر از مقداری باشد که برای بهره‏برداری مناسب دیگ بخار لازم است؛ بنابراین، خروجی توان حقیقی هر ژنراتور، باید بین حد بالا و پایین مانند رابطه (8) باشد.

که  و  کمترین و بیشترین تولید واحد i ام هستند. این حدود، علاوه بر اینکه ناشی از محدودیت‏های فنی هر واحد است، باعث می‏شوند واحدهای با هزینه کمتر، بیش از حداکثر توان مجاز خود و واحدهای با هزینه بیشتر، کمتر از حد مجاز خود تولید نداشته باشند ]18[.

محدودیت سرعت تولید: یکی از فرض‌های غیرعملی ساده‌سازی مسئله در بسیاری از پژوهش‌های قبلی، تنظیمات توان خروجی هم‌زمان است؛ اما در شرایط عملی، در خروجی واحد تولید، محدودیت دامنۀ (اندازه) تغییرات وجود دارد؛ بنابراین، خروجی ژنراتور باید بین حد بالا و پایین اندازه تغییرات باشد. اعمال این کران‌ها محدودیت عملکرد ژنراتور را به‌صورت رابطه (9) اصلاح می‌کند.

که  توان خروجی قبلی واحد i ام و  میزان‏ کاهشی[18] واحد است ]3[.

میزان‌های افزایشی کاهشی: نیروگاه‏های حرارتی به دلایل فنی نباید به‌صورت آنی توان خود را تغییر دهند و این افزایش یا کاهش باید با شیب ملایمی همراه باشد. به عبارت دیگر، هر نیروگاه محدودیت‏هایی در شیب تغییرات توان تولیدی خود دارد که تجاوز از این محدودیت‏ها ایجاد خسارت به روتور و افزایش هزینۀ عملکرد را باعث می‏شود. در شرایط واقعی تغییر تولید در هر ساعت از شبانه روز نسبت به ساعت قبل محدودیت‏هایی دارد؛ یعنی محدودیت حداکثر افزایش و یا حداقل تولید در تحلیل سیستم، وارد و به شرایط واقعی نزدیک‌تر خواهد ‏شد. این شرایط با رابطه (10) نشان داده شده است.

که در آن  و  حداکثر افزایش و کاهش تولید نسبت به ساعت قبل نیروگاه i ام،  تعداد نیروگاه‏ها و  ساعتی از شبانه‌روز هستند ]19[.

منطقه عملیاتی ممنوعه: مناطق عملیاتی (کاری) ممنوعه ژنراتور مناطقی از توان خروجی‌اند که در آن، کارکرد ژنراتور با توجه به محدودیت‌های فیزیکی و یا ناپایداری محدود شده است. به سبب دشواربودن تعیین مناطق دقیق ممنوعه، معمولاً از کار در چنین مناطقی ممانعت به عمل می‌آید ]22[. مناطق عملیاتی مناسب واحد تولید در توزیع اقتصادی با رابطه‏های (11) به دست می‏آیند.

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

محدودیت‏های قابلیت اطمینان و امنیت شبکه: معیار قابلیت اطمینان و امنیت شبکه که محدودیت‏های فنی مسئلۀ توزیع بار اقتصادی‌اند، در حل مسئلۀ پخش بار اقتصادی گنجانده نمی‏شوند و بیشتر موارد، این محدودیت‏ها در سایر مطالعات یا برنامه‌ریزی‌ها در نظر گرفته و توزیع اقتصادی بار بدون این محدودیت‏ها حل می‏شود.

3- الگوریتم بهینه‏سازی یادگیری ردیابی بازگشتی (LBSA)

دبائوچن (Debaochen) و همکارانش، اولین بار در سال 2017 الگوریتم LBSA را برای حل مسائل ریاضی معرفی کردند که یک الگوریتم بهینه‌سازی تصادفی بر اساس جمعیت است و با استفاده از ترکیب دو روش بهینه‌سازی ردیابی بازگشتی (BSA) و آموزش و یادگیری (TLBO) ساخته شده است.

3-1- BSA بنیادی

BSA الگوریتم تکاملی مبتنی بر جمعیت است. پنج مرحله از BSA به شرح زیرند و جزئیات بیشتر آن در مرجع ]15[ بیان شده است.

مقداردهی اولیه: جمعیت اولیه ( ) و تاریخچه جمعیت ( ) با توجه به رابطه‏های (12) مقداردهی می‏شوند که U توزیع یکنواخت است و  و  مرزهای پایین و بالای متغیرها بوده‌اند و i فرد i ام جمعیت است.

گزینش 1: در هر تکرار، تاریخچه جمعیت ( ) با توجه به رابطه (13) معرفی می‏شود.

جهش: فرم اولیه جمعیت آزمایشی با عملیات جهش از رابطه (14) و با استفاده از مزیت تجربیات نسل‌های قبلی و استفاده از ماتریس جستجوی جهت‌دار F تولید می‏شود. مقدار مشترک F برابر با 3•randn است.

ترکیب: ماتریسی با داده‌های صحیح باینری (نگاشت)، مسیرهای ترکیب‌شده الگوریتم را با رابطۀ (15) هدایت می‏کند؛ به‌طوری‌که افرادی با برازش بهتر برای تکامل افراد جامعه هدف استفاده می‏شوند.

گزینش 2: جمعیت نسل بعدی با توجه به مکانیسم گزینش حریصانه تولید می‏شود. طبق رابطه (16) اگر برازش  نسبت به  کوچک‌تر باشد،  جایگزین  می‏شود.

3-2- مراحل اصلی TLBO

الگوریتم TLBO نیز مبتنی بر جمعیت است که آموزش و یادگیری را تقلید می‏کند. دو مرحله اصلی آن عبارت‌اند از: مراحل آموزش‌دهنده و یادگیرنده. در الگوریتم پیشنهادی LBSA از این دو مرحله استفاده شده است. جزئیات بیشتر در مرجع ]16[ بیان شده است.

مرحله آموزش‌دهنده: آموزش‌دهنده دانش خود را برای همۀ یادگیرنده‌ها مانند رابطه (17) در گروه توزیع می‏کند.

که  و  موقعیت‌های جدید و قبلی یادگیرنده i و rand(•) عددی تصادفی در بازۀ [0,1] است. ضریب آموزش ، مقدار میانگین را برای تغییر داده‌شدن تعیین می‌کند و به‌صورت ابتکاری روی 1 یا 2 تنظیم می‌شود که با احتمال برابر با رابطه (18) بیان می‌شود.

مرحله یادگیرنده: در هر تکرار، یادگیرنده k ام به‌صورت تصادفی، هدف یادگیری، یادگیرنده i ام انتخاب می‏شود. روش یادگیری، یادگیرنده i ام به‌صورت ریاضی با رابطه (19) بیان می‌شود.

که ، موقعیت جدید فرد i ام است و  و ، موقعیت‌های پیشین i ام و k ام یادگیرنده‌ها هستند. علاوه بر این، اگر برازش  نسبت به  بهتر باشد،  پذیرفته می‏شود.

3-3- عملیات جهش در LBSA

برای بهبود توانایی الگوریتم LBSA، هدایت‌کننده‏ای، یادگیری بهترین فرد را در فرآیند جهش الگوریتم با رابطه (20) معرفی می‏کند.

عملیات جهش اصلاح‌شده به‌صورت رابطۀ (21) خلاصه می‌شود:

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

4- روش حل مسئله پخش بار اقتصادی با استفاده از LBSA

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

مقداردهی اولیه الگوریتم: در این مرحله تعداد جمعیت (PopSize) برابر با 50 نفر در نظر گرفته می‏شود. بُعد مسئله (DimSize)، تعداد ژنراتورهای سیستم مطالعه‌شده است و همچنین شمارندۀ تکرار الگوریتم معرفی می‌شود.

مقداردهی اولیه مسئله: حداقل و حداکثر توان ژنراتورها ( )، ضرایب تابع هزینه ( ،  و ) و ضرایب اثر شیر بخار (  و )، مقداردهی می‌شوند.

تولید تصادفی جمعیت اولیه (Pop) و تاریخچه جمعیت (OldPop): به‌صورت تصادفی با اعدادی بین حداقل و حداکثر توان ژنراتورها ساخته می‌شود.

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

همان‌طور که در فلوچارت مشخص است مراحل انتخاب 1، جهش، ترکیب و انتخاب 2 با توجه به روابط ارائه‌شده انجام می‏شوند. پارامترهای کنترل‌کننده مربوط به الگوریتم در جدول (1) آورده شده‌‌اند.

شکل (2): فلوچارت اجرای LBSAبرای پخش بار اقتصادی

5- نتایج شبیه‌سازی

برای بررسی امکان‌پذیری و کارایی الگوریتم بهینه‏سازی پیشنهادی، حل مسئلۀ پخش بار اقتصادی، روی سیستم‌ آزمایشی 6 واحده و سیستم 30-BUS IEEE، اصلاح و با در نظر گرفتن محدودیت‌های مختلف در MATLAB R2014a اجرا شده است.

5-1- سیستم مطالعه‌شدۀ اول

سیستم شامل 30 شین، 46 خط انتقال و 6 واحد تولید حرارتی است که کل تقاضای بار سیستم برابر با 1263 مگاوات است. ضرایب تابع هزینه، حدود واحدهای تولیدی، میزان‌های افزایشی و کاهشی و نواحی ممنوعه سیستم در جداول (2،3) و شیوۀ اتصال شبکه و دیاگرام تک‌خطی برای این سیستم در شکل (3) نشان داده شده‌‌اند ]13،22[.

شکل (3): دیاگرام تک‌خطی سیستم 6 نیروگاهی ]4[

جدول (2): ضرایب هزینه و حدود واحدها ^*]13،22[

جدول (3): میزان افزایشی، کاهشی و نواحی ممنوع ^* ]13،22[

شبیه‌سازی برای سیستم 6 واحده با استفاده از الگوریتم پیشنهادی LBSA انجام شد. نتایج عددی این شبیه‌سازی‌ها در جدول (4) دیده می‏شود. این نتایج با روش بهینه‏سازی -PSO ]3[، الگوریتم بهبودیافته کلونی زنبورعسل (IABC) ]4[ و الگوریتم خفاش بی‌نظم (Bat.A) ]13[ مقایسه شد که در جدول (4) آورده شده‌اند. در شکل (4) نیز توان تلف‌شده و هزینۀ نهایی سیستم مطالعه‌شده در این الگوریتم‌ها مقایسه شده است.

جدول (4): مقایسه نتایج به‌دست‌آمده از الگوریتم‌های مختلف روی سیستم موردمطالعه اول (6 واحده) ^*

^* در این جدول هزینه‏ها بر حسب ($/h)، همه توان‏ها بر حسب MW)) و زمان بر حسب Second است.

با توجه به شکل (4) و جدول (4)، هزینۀ نهایی در روش LBSA نسبت به روش Bat Algorithm در حدود 174 واحد، کاهش و نسبت به الگوریتم  در حدود 162 واحد کاهش داشته است؛ اما نسبت به روش IABC حدود 22 واحد افزایش داشته است. همچنین مقدار تلفات نهایی برای الگوریتم LBSA در مقایسه با بهترین جواب از سایر الگوریتم‌های مربوط به الگوریتم IABC و مقدار آن 18/0 مگاوات است؛ حدود 0144/0 واحد افزایش داشته است. گفتنی است در تعداد دفعات شبیه‌سازی، همواره تلفات در روش LBSA، کمتر از یک واحد به‌دست‌آمده است که یکی از مزیت‌های الگوریتم پیشنهادی خواهد بود. با توجه به جدول (4)، زمان لازم برای به نتیجه رسیدن نهایی در الگوریتم پیشنهادی، کمتر از یک ثانیه است که بعد از روش IABC در رتبه دوم قرار دارد. روند بهینه‌سازی این سیستم در شکل (5) نشان داده شده است که روند صحیح بهینه‌سازی را نشان می دهد. نتایج نشان می‏دهند روش پیشنهادی برای این مطالعه‌شده، روش کارآمد و سریعی است و جواب‌هایی درخور ‌مقایسه با سایر الگوریتم‌های مشابه به دست می‌آورد.

5-2- سیستم مطالعه‌شدۀ دوم

سیستم مطالعه‌شدۀ دوم، سیستم 30-BUS IEEE اصلاح‌شده است که ژنراتورهای آن شامل سه واحد سوخت زغالی، دو واحد سوخت گازی، یک واحد سوخت نفتی و دو واحد توربین بادی است. ضرایب تابع هزینه و حدود واحدهای تولیدی سیستم در جدول (5) مشخص است ]2،23[. شیوۀ اتصال و دیاگرام تک‌خطی سیستم مطالعه‌شدۀ دوم نیز در شکل (6) نشان داده شده است ]23[.

شبیه‌سازی این سیستم بدون در نظر گرفتن ژنراتورهای بادی، برای تقاضای بارهای 1200، 1400 و 1600 مگاوات با محدودیت شیر بخار انجام شده است که بهترین جواب‌های به‌دست‌آمده در جدول (6) آورده شده‌‌اند.

در شکل (7) روند بهینه‌سازی در تکرارهای مختلف برای این مورد مطالعه نشان داده شده است که با توجه به جواب به‌دست‌آمده روند صحیح بهینه‌سازی را نشان می‏دهد.

نتایج به‌دست‌آمده از این شبیه‌سازی با روش‌های بهینه‏سازی انتشار ذرات^[19] (DPO)، الگوریتم کلونی زنبور مصنوعی^[20] (GABC) و الگوریتم ذرات کوانتومی الهام‌بخش^[21] (QPSO) که همین سیستم را مطالعه کرده ]2[، مقایسه شده‌اند. روش DPO بدون اثر شیر بخار و روش‌های GABC و QPSO با در نظر گرفتن اثر شیر بخار، شبیه‌سازی شده‌اند.

شکل (6): دیاگرام تک‏خطی و شیوۀ اتصال سیستم 30-BUS IEEE اصلاح‌شده ]23[

(الف)

(ب)

(ج)

شکل (7): روند بهینه‌سازی هزینه در تکرارهای مختلف برای سیستم مطالعه‌شدۀ دوم در دو حالت با و بدون اثر شیر بخار در بارهای الف) 1200 مگا وات ب) 1400 مگاوات ج) 1600 مگا وات

بار 1200 مگاوات: جدول (7) نتایج شبیه‏سازی را در بار درخواستی 1200 مگاوات نشان می‌دهد. با در نظر گرفتن اثر شیر بخار،  هزینه در الگوریتم پیشنهادی، نسبت به الگوریتم GABC و QPSO در حدود 52  و 460 واحد کاهش داشته است. بدون اثر شیر بخار نیز در حدود 62 واحد کاهش نسبت به الگوریتم DPO وجود دارد.

در این شبیه‌سازی توان تلف‌شده در حالتی که اثر شیر بخار در نظر گرفته شده، برای هر دو الگوریتم GABC و QPSO مقدار صفر به دست آمده است و در حالتی که اثر شیر بخار در نظر گرفته نشده، در الگوریتم DPO مقدار 35/0 مگاوات است که در هر دو حالت در روش LBSA نسبت به الگوریتم‌های مقایسه‌شده با اینکه افزایش را نشان می‌دهد، این افزایش کمتر از 4 واحد است.

بار 1400 مگاوات: جدول (8) مقایسه نتایج حاصل از شبیه‏سازی را در بار درخواستی 1400 مگاوات نشان می‌دهد. در اینجا نیز الگوریتم پیشنهادی نتایج بهتری نسبت به دیگر الگوریتم‌ها به دست آورده است.

بار 1600 مگاوات: با توجه به جدول (9)، نتیجه می‌گیریم در بار درخواستی 1600 مگاوات، نیز الگوریتم پیشنهادی عملکرد بهتری نسبت به دیگر الگوریتم‌ها به دست آورده است.

بررسی شبیه‌سازی ازنظر زمان محاسبات فقط با روش DPO مقایسه می‌شود؛ زیرا اطلاعاتی دربارۀ زمان شبیه‌سازی برای الگوریتم GABC و QPSO در منبع ]2[ که از اطلاعات آن استفاده شده، موجود نبوده است.

با توجه به زمان انجام شبیه‏سازی در همه حالت‏ها، با درنظر گرفتن اثر شیر بخار یا بدون این اثر و در بارهای 1200، 1400 و 1600 مگاوات با توجه به جدول‏های (5، 6 و 7) مشخص است که عملکرد الگوریتم پیشنهادی ازنظر زمانی در حدود 1 ثانیه بهتر از روش DPO بوده است.