Document Type : Research Article
Authors
1 Ph.Dd student, Faculty member of the University of Science and Technology, Tehran, Iran
2 Associate professor, Faculty member of the University of Science and Technology, Tehran Iran
3 professor, Faculty member of the University of Science and Technology, Tehran Iran
4 The member of Dept. of Power Systems Planning and Operation, Niroo Research Institute, Tehran, Iran
Abstract
Keywords
یکی از مهم‌ترین تغییرات درحال انجام در شبکه‌های قدرت، انتقال از منابع انرژی متمرکز سنتی به منابع انرژی پراکنده است. منابع انرژی پراکنده با توجه به فواید زیست‌محیطی بسیار آنها نقش کلیدی در تولید انرژی پاک و پایدار دارند [1]. منابع تولید پراکنده می‌توانند انتشار کربن، تلفات ارسال توان و هزینه ساخت زیربناها را کاهش دهند [2]. ریزشبکه‌ها شبکه‌های قدرت مقیاس کوچک و خودحمایت‌گر هستند که از منابع انرژی پراکنده استفاده می‌کنند. منابع انرژی پراکنده می‌تواند شامل منابع انرژی تجدیدپذیر مانند توربین‌های بادی و پنل‌های خورشیدی، منابع انرژی تجدیدناپذیر ازجمله دیزل ژنراتورها و سیستم‌های ذخیره انرژی باتری باشد. ریزشبکه‌ها می‌توانند در دو مد متصل به شبکه اصلی یا غیرمتصل (جزیره‌ای) کار کنند [3].
ریزشبکه‌ها قادرند قابلیت اطمینان شبکه را بهبود ببخشند و انرژی با کیفیت و پایدار را تأمین کنند؛ اما برنامه‌ریزی و عملکرد آنها به دلیل عدم‌قطعیت ناشی از پیش‌بینی بار مصرفی و توان خروجی انرژیهای تجدیدپذیر با چالش‌های زیادی مواجه است. اگرچه ریزشبکه‌ها با مهیاکردن یک مسیر انعطاف‌پذیر برای ترکیب منابع انرژی پراکنده تجدیدپذیر به شبکه‌های قدرت، نقش مهمی در مسیر مدرنیته‌کردن آنها دارند، منابع تجدیدپذیر با توجه به شرایط مختلف جوی، متغیر و پیش‌بینی‌ناپذیر است؛ برای مثال، ماژول‌های فتوولتائیک، تنها در حضور تابش خورشید می‌توانند برق تولید کنند. همچنین، مزارع بادی نیز در حضور باد کافی قادر به تولید برق هستند. علاوه بر مشکلات ناشی از وجود عدم‌قطعیت در ریزشبکه‌ها، لازم است از وابستگی ریزشبکه‌ها به شبکه اصلی نیز کاسته شود؛ به نحوی که سود کلیه واحدهای تولید انرژی افزایش و هزینه مصرف‌کنندگان داخلی در ریزشبکه‌ها نیز کاهش یاید. سیستم‌های ذخیره انرژی باتری به‌منظور ذخیره و تأمین برق در زمانهای مختلف در ریزشبکه‌ها استفاده می‌شوند. حداکثر میزان انرژی قابل شارژ یا تخلیه در یک زمان مشخص در یک باتری وابسته به قابلیت ذخیره انرژی، حالت شارژ فعلی (SOC) و مشخصات ذاتی آن است. همچنین، SOC فعلی باتری با رفتار شارژ/ تخلیه قبلی مشخص میشود؛ بنابراین، مدیریت انرژی باتری به یک مسئله تصمیم‌گیری ترتیبی برای یک سیستم دینامیکی تبدیل میشود که تصمیمات قبلی بر گزینه‌های موجود آتی تأثیرگذار است. باتری‌ها دارای خواص متغیربازمان و غیرخطی هستند؛ از این رو برنامه‌ریزی شارژ و تخلیه آنها چالش‌های جدیدی به سیستم مدیریت انرژی اضافه می‌کند و نیازمند الگوریتم‌های کنترلی پیشرفته‌تری است [4]. روش‌های مدیریت انرژی هوشمند به‌تازگی مورد توجه بسیاری قرار گرفته است. این روش‌ها برای حل مشکلات ناشی از عدم‌قطعیت و تغییرپذیری در ریزشبکه‌ها و ایجاد توان قابل اعتماد برای شبکه‌های قدرت استفاده میشود.
روش‌های مختلفی در سال‌های اخیر برای حل مشکل عدم‌قطعیت در مدیریت انرژی ریزشبکه‌ها توسعه داده شده است [5]. در [6] یک استراتژی کنترل برای عملکرد هماهنگ ریزشبکه‌ها در یک سیستم توزیع ارائه می‌دهد. اپراتور شبکه توزیع و هر ریزشبکه به‌عنوان یک واحد مجزا با توابع هدف منحصربه‌فرد در نظر گرفته میشود؛ به طوری که هزینههای عملیاتی کمینه شود. مسئله به‌عنوان یک مسئله دوسطحی تصادفی فرموله شده است. در سطح بالا، اپراتور شبکه و در سطح پایین، ریزشبکه‌ها در نظر گرفته شده است. در [7] یک مسئله تصادفی برای برنامه‌ریزی انرژی ریزشبکه‌ها و حل چالش‌های عملیاتی بارهای کنترل‌پذیر و منابع انرژی تنظیم می‌کند؛ به طوری که هزینههای عملیاتی مورد انتظار ریزشبکه و تلفات توان را به حداقل رساند و با طبیعت متناوب انرژیهای تجدیدپذیر سازگار است. در [8] یک مدل بهینه‌سازی محدب مقاوم برای مدیریت انرژی ریزشبکه‌ها ارائه شده است. در حالت متصل به شبکه اصلی، هزینه انرژیهای وارده به شبکه، بارهای تأمین‌شده با منابع انرژی پراکنده و باتری کمینه میشود. در حالت جزیره‌ای بارهای تأمین‌نشده، با در نظر گرفتن اولویت مصرف به حداقل می‌رسد. در روش‌های بهینه‌سازی لیاپانف، کنترل بهینه یک سیستم پویا ازطریق تابع لیاپانف به دست می‌آید. در [9] به توسعه یک روش مدیریت انرژی برخط برای کارکرد زمان حقیقی ریزشبکه‌ها با در نظر گرفتن پخش بار و قیود عملکردی سیستم پرداخته شده است. مدیریت انرژی برخط به‌عنوان مسئله تصادفی پخش بار بهینه، مدلسازی و از بهینه‌سازی لیاپانف استفاده شده است. در مقاله [10]، روش بهینه‌سازی دو مرحله‌ای ارائه شده است. در مرحله اول، برنامه ساعتی با استفاده از برنامه‌ریزی روز آینده انجام می‌گیرد. در مرحله دوم، پخش بار اقتصادی و مبادله انرژی با استفاده از بهینه‌سازی لیاپانوف به‌صورت زمان حقیقی برنامه‌ریزی میشود. در [11] در ساختار ریزشبکه از دو کنترل‌کننده مرکزی برای ریزشبکه و شبکه گاز استفاده شده است. برای کاهش هزینههای بهره‌برداری، مسئله خرید و فروش انرژی در ریزشبکه با استفاده از یک مدل خطی عدد صحیح‌آمیخته، مدلسازی و به کمک نرم‌افزار GAMS مسئله فوق حل شده است.
در روش‌های فوق فرض میشود مشخصات پروسه‌ها و متغیرهای تصادفی با مدل‌های پیش‌بینی‌کننده یا مقادیر مورد انتظار مربوط به آنها موجود است. وابستگی این روش‌ها به یک مدل تخمین‌گر منجر شده است تا دقت این روش‌ها به دقت مدل تخمین زده شده وابسته باشد. یادگیری تقویتی روش‌های مستقل از مدل برای حل مسائل کنترل بهینه سیستم‌های دینامیکی ارائه می‌دهد که نیازی به وجود خواص تصادفی پروسه‌های تصادفی نیست.
در [12]، سیستم مدیریت انرژی یک ساختمان برای کاهش پیک مصرف انرژی با استفاده الگوریتم فراابتکاری دومرحله‌ای پیاده‌سازی شده است. هزینه و مصرف انرژی با حضور انرژیهای تجدیدپذیر با استفاده از الگوریتم ژنتیک در [13] کمینه شده است. در [14]، از الگوریتم بهینه‌سازی بویایی کوسه و گرگ خاکستری برای ارزیابی تأثیر پاسخگویی بار در ریزشبکه‌ها استفاده شده و هزینه تولید و تلفات شبکه کمینه شده است. گرچه روش‌هایی که از الگوریتم‌های ابتکاری استفاده می‌کنند به مدل ریاضی نیازی ندارند و برای بهینه‌سازی جهانی برنامه‌ریزی غیرخطی و غیرهموار، به‌طور منعطف‌تر و موثرتری عمل می‌کنند، بیشتر این روش‌ها قادر به یادگیری تقلیدی و ذخیره دانش قبلی نیستند و در هر مرحله یک جمعیت جدید به‌صورت تصادفی انتخاب میشود؛ به همین دلیل، زمان اجرا برای محاسبه نقطه بهینه زیاد میشود [15].
از نظریه بازی‌ها نیز در طراحی سیستم مدیریت انرژی ریزشبکه‌ها استفاده میشود [16]. در مقاله [17]، برنامه‌ریزی روز آینده ریزشبکه‌ها و شرکت توزیع با استفاده از نظریه بازی‌ها پیاده‌سازی شده است. برنامه‌ریزی تقاضا و سیستم ذخیره انرژی به‌عنوان یک مسئله بهینه‌سازی چندتابع هدفه مدلسازی شده است. خروجی منابع تجدیدپذیر با استفاده از [18] تخمین زده و به کمک تکنیک نظریه بازی در مقاله [19] مسئله مدیریت انرژی حل شده است. در مقاله [20]، با استفاده از تئوری بازی‌های غیر مشارکتی، مدیریت انرژی چندعامله یک ریز شبکه با حضور منابع انرژی تجدیدپذیر و بارهای فصلی پیاده‌سازی شده است. در مقاله [21]، برنامه‌ریزی یک شبکه قدرت در حالت متصل به شبکه با استفاده از مدل نظریه بازی‌های مشارکتی و غیرمشارکتی به دست آمده است. توربین‌های بادی، پنل خورشیدی و باتری‌ها بازیکن در مسئله در نظر گرفته شده است. وجود نقطه تعادل نش با تحلیل تقعر توابع بازده و مدل عدم قطعیت اثبات شده است. برای یافتن نقطه تعادل در مقالات [20, 21]، از یک روش جستجو تکراری [22] استفاده شده است. همچنین، از پیش‌بینی سرعت باد، شدت تابش خورشید و تقاضای بار و الگوریتم ازدحام ذرات نیز استفاده شده است. در مسائل فوق برای حل مسئله به روش نظریه بازی‌ها، سیستم تنها شامل یک حالت در نظر گرفته میشود. در نظریه بازی‌ها یک فرض اساسی وجود دارد که بازیکن‌ها دارای توانایی یادگیری و سازگاری یکسانی هستند. بیشتر الگوریتم‌های نظریه بازی به مدلی از محیط نیاز دارند؛ زیرا از انتقال، T و پاداش، r استفاده می‌کنند. هدف این الگوریتم‌ها نیز محاسبه ارزش تعادلی بازی (یعنی پاداش‌های وزن‌دار مورد انتظار برای هر یک از عوامل)، به جای یافتن سیاست‌های تعادلی است. این بدان معنی است که آنها اغلب الزامات قوی برای رفتار همه عوامل ایجاد می‌کنند. درمقابل، الگوریتم‌های یادگیری تقویتی فرض می‌کنند جهان ناشناخته است و فقط مشاهدات و پاداش اعمال آنها از محیط موجود است. هدف در روش‌های یادگیری تقویتی این است که عامل خط‌مشی خود را در راه‌حل تعادل بازی بیابد و معمولاً الزامات کمتری برای رفتار عاملهای دیگر ایجاد می‌کند [23]؛ به همین دلیل به‌تازگی به استفاده از روش‌های یادگیری تقویتی توجه زیادی شده است. روش‌های یادگیری تقویتی در مقایسه با روش‌های یادگیری نظارتی و غیرنظارتی در حیطه انرژی، قابلیت‌های جالبی را در زمینه کاربردهای کنترلی ارائه می‌دهد [24]. در سیستم‌های کنترلی، با توجه به اینکه به دست آوردن اطلاعات اولیه از محیط و مدلسازی سیستم در عمل بسیار سخت است، یادگیری تقویتی قادر است روش‌هایی مستقل از مدل برای حل مسائل دارای عدم‌قطعیت ارائه دهد. این روش‌ها با استفاده از ذخیره اطلاعات تجربه‌شده قبلی، عملکرد عامل یادگیرنده را بهبود می‌بخشند. در [25]، برای مدیریت انرژی ساختمان از یک روش داده‌محور مبتنی بر شبکه عصبی و یادگیری Q استفاده شده است. یک تکنیک برنامه‌ریزی پویا تطبیقی براساس یادگیری تقویتی برای کنترل یک ریزشبکه هوشمند در [26] توسعه داده شده است. در [27]، مدیریت انرژی چندعامله غیرمتمرکز برای بارهای الکتریکی در یک ریزشبکه با استفاده از یادگیری Q پیاده‌سازی شده است. در [28] یادگیری تقویتی سلسله‌مراتبی برای محاسبه سیاست بهینه توسعه داده شده است. اگرچه روش ارائه‌شده برای حل معضل ابعاد در روش‌های یادگیری تقویتی سودمند است، به یک نوع سیاست بهینه محلی به نام سیاست بهینه بازگشتی همگرا میشود [29]. پخش بار اقتصادی توزیع‌شده در [30] با استفاده از الگوریتم یادگیری تقویتی مشارکتی براساس اطلاعات دریافتی از عاملهای همسایه و با استفاده از استراتژی Diffusion به دست آمده است.
با استفاده از الگوریتم یادگیری Nash Q بار درخواستی واحدهای ریزشبکه تخصیص داده شده و سود هرکدام بیشینه شده است [31]. الگوریتم یادگیری Nash Q یک توسعه از الگوریتم یادگیری Q معمولی برای سیستم‌های چندعامله غیرمشارکتی است [32]. در الگوریتم یادگیری Nash Q، یک عامل نه‌تنها پاداش خود، پاداش و اعمال سایر عامل‌ها را نیز دریافت می‌کند. در واقعیت، موجودبودن اطلاعات مربوط به اعمال و پاداش‌های سایر عامل‌ها برای همه عاملهای مصرف‌کننده و تولیدکننده یا حتی برای یک سیستم مرکزی نیز به‌راحتی امکان‌پذیر نیست. همچنین، با افزایش تعداد عاملهای یادگیرنده، سایز تابع Nash Q زیاد میشود؛ درنتیجه، زمان اجرا بسیار زیاد و انجام محاسبات پیچیده میشود.
روش یادگیری تقویتی عمیق از تلفیق روش‌های یادگیری عمیق با روش یادگیری تقویتی به دست می‌آید. در این روش‌ها، با تخمین توابع ارزش و توابع سیاست با استفاده از شبکه‌های عصبی عمیق مشکل ابعاد توابع Q حل میشود. در [33]، یک شبکه عمیق Q برای حل مسائل با تعداد زیاد سنسورهای ورودی توسعه داده شده است. برنامه‌ریزی زمان حقیقی یک ریزشبکه در [34] با استفاده از شبکه‌های عصبی عمیق برای تخمین توابع ارائه شده است. در [35]، برای افزایش انعطاف‌پذیری و قابلیت اطمینان یک ریزشبکه باحضور منابع تجدیدپذیر، از الگوریتم بهینه‌سازی Proximal Policy براساس یادگیری تقویتی عمیق و شبکه‌های عصبی مرکزی critic بهره گرفته شده است. برای کمینه‌کردن هزینههای عملیاتی ژنراتورها و کاهش هزینههای خرید انرژی از شبکه اصلی در [36]، توان اکتیو و راکتیو ژنراتورهای معمولی و مقدار توان شارژ/تخلیه باتری با استفاده از روش یادگیری تقویتی عمیق محاسبه شده است. در مقاله [37]، با استفاده از یک روش یادگیری تقویتی عمیق مبتنی بر سیاست با فضای تصمیم‌گیری و حالت پیوسته، هزینههای عملیاتی ریزشبکه شامل هزینه تبادل انرژی با شبکه توزیع‌شده، هزینه عملیاتی ریزتوربین و هزینه عملیاتی باتری کمینه شده است. به‌تازگی توجهی زیادی به ترکیب روش‌های یادگیری تقویتی و یادگیری عمیق شده است. روش‌های عمیق به کمک روش‌های یادگیری تقویتی می‌آید تا مشکل ابعاد برای محاسبه تابع Q با تعداد عاملهای زیاد را برطرف کند؛ اما همچنان سایر مشکلات موجود در روش‌های متمرکز در این روش‌ها برطرف نشده است. در روش‌های یادگیری عمیق لازم است اطلاعات مربوط به همه عامل‌ها شامل اعمال و پاداش‌ها برای یک واحد کنترل‌کننده مرکزی موجود باشد. بیشتر روش‌های ارائه‌شده برای مدیریت انرژی ریزشبکه‌ها از یک ساختار کنترل مرکزی استفاده می‌کنند. در کنترل‌کننده‌های متمرکز، یک واحد به‌عنوان کنترل‌کننده انتخاب میشود و مسئولیت مدیریت سایر واحدها را برعهده دارد. در ساختار متمرکز، واحد کنترل‌کننده با همه عامل‌ها در ارتباط است؛ در حالی که در ساختار کنترلی توزیع‌شده، کنترل‌کننده‌ها با همسایگی‌های خود در ارتباط هستند؛ اما در ساختار کنترل‌کننده غیرمتمرکز ارتباطی بین کنترل‌کننده‌ها به‌صورت جداگانه نیست [38, 39]. هر دو ساختار توزیع‌شده و غیرمتمرکز به‌عنوان سیستم کنترلی چندعامله در نظر گرفته میشود. در سیستم‌های قدرت ابعاد وسیع، هنگامی که واحدهای تولید برق در شبکه پراکنده شدهاند و ارتباطات، محدود و غیرقابل اعتماد هستند، استفاده از ساختارهای کنترلی متمرکز با چالش‌های جدی روبه‌رو میشود [40]. در ساختارهای کنترلی متمرکز، هر خطایی در کنترل‌کننده باعث خرابی کل سیستم میشود. برای کاهش سایز و پیچیدگی کنترل متمرکز از روش‌های کنترل غیرمتمرکز استفاده میشود. این ساختار کنترلی، یک معماری کنترلی ساده با توانایی بالا برای تضغیف اغتشاشات را فراهم می‌کند.
در جدول 1، به‌طور خلاصه روش‌های مدیریت انرژی با استفاده از یادگیری تقویتی مقایسه شده است. در بیشتر مقالات به مدیریت انرژی در زمینه بارهای الکتریکی پرداخته شده و به بارهای گرمایی توجه نشده است. برنامه‌ریزی انرژی در بسیاری از این روش‌ها به‌صورت متمرکز یا توزیع‌شده انجام می‌گیرد. همچنین، همان‌طور که مشاهده میشود در بیشتر مقالات روشی برای مدیریت مصرف و بهینه‌سازی سود مصرف‌کنندگان و تولیدکنندگان با ارائه پیشنهاد قیمت به‌صورت هم‌زمان انجام نشده است. در برخی مقالات، طول عمر باتری در نظر گرفته شده است؛ اما در هیچ‌کدام تأثیر مدل باتری بر محاسبه تعداد تعویض باتری لحاظ نشده است.
بنابراین، در این مقاله، با ارائه یک ساختار کامل غیرمتمرکز به طراحی سیستم مدیریت انرژی ریزشبکه‌ها پرداخته میشود. عامل‌ها بدون دردسترس بودن اطلاعات عاملهای همسایه و تنها با دریافت حالات محیط به یادگیری و آموزش می‌پردازند. طراحی سیستم به‌صورت مستقل از مدل و با استفاده از یادگیری تقویتی انجام میشود. در این پژوهش، مدیریت انرژی هر دو نوع بار الکتریکی و گرمایی در نظر گرفته میشود. فرض میشود ریزشبکه شامل منابع انرژی پراکنده گرمایی و الکتریکی، سیستم ذخیره انرژی باتری و بارهای مصرفی الکتریکی و گرمایی است. سیستم به گونه‌ای طراحی میشود که سود کلیه منابع افزایش یابد، هزینه مشتریان کمینه شود و از وابستگی ریزشبکه به شبکه اصلی کاسته شود. همچنین، با در نظر گرفتن طول عمر باتری، هزینه ناشی از تخریب باتری کمینه میشود.
به‌طور خلاصه نوآوری‌های اصلی مقاله به‌صورت زیر خلاصه شدهاند:
در همین راستا در بخش دوم، ساختار ریزشبکه معرفی میشود. در بخش سوم، طراحی سیستم مدیریت انرژی با استفاده از یادگیری تقویتی ارائه میشود. در بخش چهارم، شبیه‌سازی و نتایج بررسی می‌شوند. در قسمت آخر نتیجه‌گیری ارائه میشود.
جدول (1): مقایسه تحقیقات انجام‌شده در زمینه مدیریت انرژی ریزشبکه‌ها
|
مقاله |
بار الکتریکی |
بار گرمایی |
ساختار غیر متمرکز |
طول عمر باتری |
مدیریت بار |
ارائه پیشنهاد قیمت |
مدیریت تولیدکنندگان (هزینه تولید و سود فروش) |
مقایسه |
|
[25] |
ü |
× |
ü |
× |
ü |
× |
× |
ü |
|
[26] |
ü |
× |
× |
ü |
ü |
× |
× |
ü |
|
[27] |
ü |
× |
ü |
× |
ü |
ü |
ü |
× |
|
[28] |
ü |
× |
ü |
× |
ü |
× |
× فقط هزینه تولید |
ü |
|
[30] |
ü |
× |
× |
ü |
× |
× |
× فقط هزینه تولید |
ü |
|
[34] |
ü |
× |
× |
× |
× |
× |
× فقط هزینه تولید |
ü |
|
[31] |
ü |
ü |
× |
× |
× |
× |
× |
× |
|
[16] |
ü |
× |
ü |
× |
ü |
× |
× |
× |
|
[35] |
ü |
× |
× |
× |
ü |
× |
× |
ü |
|
[36] |
ü |
× |
× |
× |
× |
× |
× فقط هزینه تولید |
ü |
|
روش پیشنهادی |
ü |
ü |
ü |
ü |
ü |
ü |
ü |
ü |
به یک شبکه قدرت مقیاس‌کوچک، ولتاژ پایین و خودکار که منابع انرژی پراکنده و بارها را به هم متصل می‌کند، ریزشبکه گفته میشود. منابع انرژی پراکنده شامل انرژیهای تجدیدپذیر، انرژیهای تجدیدناپذیر و باتری است. ریزشبکه‌ها در دو حالت متصل به شبکه و جزیره‌ای کار می‌کنند [3]. در حالت کلی فرض میشود ریزشبکه‌ها به شبکه اصلی متصل‌اند [41, 42]. ریزشبکه‌ها ازطریق نقطه اتصال مشترک (Point of Common Coupling) به شبکه اصلی متصل می‌شوند. در مد متصل به شبکه، ریزشبکه‌ها تعادل تقاضا و تأمین را با فروش انرژی اضافی به شبکه اصلی و خریدن کسری انرژی از آن تأمین می‌کنند
یکی از اهداف مهم در مسئله مدیریت انرژی ریزشبکه‌ها، کاهش وابستگی آنها به شبکه اصلی است؛ بنابراین، سیستم مدیریت انرژی ریزشبکه‌ها باید به گونه‌ای طراحی شود که علاوه بر افزایش سود عامل‌ها، وابستگی ریزشبکه به شبکه اصلی نیز کاهش یابد.
بارها در ریزشبکه‌ها به دو دسته کنترل‌پذیر و غیرکنترل‌پذیر تقسیم می‌شوند. بارهای غیرکنترل‌پذیر مانند سیستم‌های مراکز درمانی و وظایف ضروری در صنعت هستند که می‌باید در زمان تقاضا تأمین شوند. درواقع، بارهای غیرکنترل‌پذیر انعطاف‌ناپذیر نسبت به زمان بوده‌اند و نمی‌توان آنها را در طول زمان جابه‌جا کرد؛ اما بارهای کنترل‌پذیر دارای قابلیت حذف یا انتقال به زمانهای کم‌باری هستند؛ حتی می‌توان برخی از بارهای کنترل‌پذیر را کاهش داد. شکل 1 ساختار یک ریزشبکه متشکل از پنل‌های خورشیدی، توربین بادی، دیزل ژنراتور، سلول سوختی الکتریکی و گرمایی، میکروتوربین گرمایی و الکتریکی، باتری و تعدادی بار محلی الکتریکی و گرمایی را نشان می‌دهد. اپراتور ریزشبکه (Microgrid Operator) یک عامل کنترل‌کننده سطح بالا در ریزشبکه‌های قدرت در نظر گرفته میشود.
|
شکل (1): ساختار ریزشبکه
2-2- توابع هدف:
هدف سیستم مدیریت انرژی برای یک ریزشبکه، بیشینه‌کردن سود کلیه عاملهای داخل شبکه در یک مدت زمان طولانی است [8, 27]؛ به همین دلیل، سود کلی عامل تولیدکننده iام برای مدت زمان طولانی به‌صورت زیر تعریف میشود:
|
(1) |
|
به طوری که t بازه زمانی است که در آن رابطه فوق محاسبه شده است. پارامتر نرخ تخفیف گفته میشود و مقداری بین صفر و یک دارد. این پارامتر بیان‌کنندۀ ارزش فعلی پاداش‌های آتی است. زمانی که به یک نزدیک میشود، به تولیدکنندۀ سود‌های آینده به‌شدت توجه می‌کند. و به‌ترتیب توان فروخته‌شده از ژنراتور iام به ریزشبکه و شبکه اصلی در بازه زمانی t است. قیمت پیشنهادی فروش انرژی از ژنراتور iام به ریزشبکه است. حالت محیط است که قیمت خرید انرژی توسط شبکه اصلی از ریزشبکه است. تابع هزینههای عملیاتی ژنراتور iام است. تابع هزینه به‌صورت عملی محاسبه میشود [43].
تابع هدف سیستم‌های ذخیره انرژی باتری در ادامه آورده شده است:
|
(2) |
|
به طوری که‌ترم اول سود حاصل از فروش انرژی، ترم دوم هزینه ناشی از خرید انرژی و ترم سوم هزینه ناشی از طول عمرباتری است. در هر بازه زمانی، باتری می‌تواند خریدار یا فروشنده انرژی باشد. قیمت پیشنهادی فروش انرژی از باتری است. و به‌ترتیب توان فروخته‌شده از باتری به ریزشبکه و شبکه اصلی در بازه زمانی t است و مقدار توان خریداری‌شده از باتری و قیمت بازار برق است. سیاست سیستم‌های ذخیره انرژی باید به گونه‌ای باشد که انرژی را در زمانهای با قیمت پایین خریداری و در پیک مصرف به شبکه بازگرداند. هزینه ناشی از کاهش طول عمر و تخریب باتری است که در اثر شارژ و تخلیه باتری رخ می‌دهد [30]. اهداف عاملهای مصرفی کمینه‌کردن هزینه‌ها است که به‌صورت زیر محاسبه میشود:
|
(3) |
|
و به‌ترتیب بارهای غیرکنترل‌پذیر و بارهای کنترل‌پذیر در بازه زمانی t است. نسبت بار خاموش‌شده به بارهای کنترل‌پذیر است. ضریب میزان نارضایتی در مصرف‌کنندگان به‌ازای حذف بار است. مقدار آن بستگی به نوع مصرف‌کننده و میزان اشتیاق آنها در مدیریت و بهینه‌سازی مصرف و هزینههای خود دارد.
برای تضمین قابلیت اطمینان و امنیت یک شبکه می‌باید توان مورد تقاضا توسط تولیدکننده‌ها در تمامی زمان‌ها تأمین شود. به‌منظور دستیابی به تعادل توان الکتریکی و تنظیم فرکانس شبکه به‌طور پیوسته از روش‌های مدیریت ذخیره استفاده میشود. در حالتی که ریزشبکه متصل به شبکه اصلی باشد، تأمین فرکانس و ذخیره عملیاتی توسط شبکه اصلی به‌عنوان یک ژنراتور اصلی تأمین میشود. برای حل مسئله مدیریت ذخیره در حالت جزیره‌ای، روش‌های کنترل مشارکتی برای تولیدکننده‌های پراکنده توسعه داده شده است [44-46]. در حالت متصل به شبکه اصلی، قید تعادل توان به معنای برابری میزان توان تولیدی با میزان بارهای مصرفی است [30]؛ بنابراین، برای بارهای الکتریکی قید تعادل توان به‌صورت زیر تعریف میشود:
|
(4) |
|
میزان تقاضای بار الکتریکی عامل مصرفی iام و n تعداد عاملهای مصرفی الکتریکی است. ، ، ، ، و به‌ترتیب توان تولیدی الکتریکی توربین بادی، پنل‌های خورشیدی، دیزل ژنراتور، باتری، ریزتوربین، سلول سوختی و شبکه اصلی است. قید تعادل توان برای بارهای گرمایی نیز به‌صورت زیر تعریف میشود:
|
(5) |
|
میزان تقاضای بار گرمایی عامل مصرفی iام و m تعداد عاملهای مصرفی گرمایی است. ، و به‌ترتیب توان تولیدی گرمایی ریز توربین، سلول سوختی و شبکه اصلی است.
قیود ظرفیت بیان‌کنندۀ بازه عملیاتی ژنراتورهای پراکنده است و دارای محدوده زیر است [31]:
|
(6) |
|
به طوری که توان خروجی ژنراتور پراکنده i در بازه زمانی t با بیان میشود. و به‌ترتیب کمینه و بیشینه توان خروجی ژنراتور i است.
SOC بیان‌کنندۀ سطح شارژ باتری نسبت به ظرفیت آن است. قید فنی زیر به‌منظور جلوگیری از شارژ و تخلیه بیش از حد سیستم ذخیره انرژی باتری اعمال میشود [30]:
|
(7) |
|
به طوری که و کمینه و بیشینه حالت شارژ باتری است. در این پژوهش SOC در بازه [0.2, 0.8] محدود شده است تا از آسیب به باتری جلوگیری شود.
اصلاح فعالیت‌ها براساس عمل و عکس‌العمل در تعامل با محیط، یادگیری تقویتی نامیده میشود. به یادگیری تقویتی، یادگیری براساس اعمال نیز گفته میشود. در این روش، یک عامل با اثرگذاشتن بر محیط و گرفتن پاسخ از آن، تلاش می‌کند اعمال و سیاست‌های کنترلی خود را بهبود ببخشد؛ درنتیجه، می‌تواند بازخورد بهتری از محیط دریافت کند. در یادگیری تقویتی، عامل با دریافت پاداش و عدم پاداش یا تنبیه در تعامل با محیط، سعی در اصلاح اعمال خود دارد.
الگوریتم‌های یادگیری تقویتی براساس این ایده ساخته شدهاند که تصمیمات کنترلی درست می‌بایست در حافظه سیستم توسط سیگنال تقویتی باقی بماند بطوری که در دفعات بعدی احتمال استفاده از آنها بیشتر باشد.
یادگیری تقویتی از ساختار رسمی پروسه‌های تصمیم‌گیری مارکوف استفاده می‌کند و ارتباط بین یک عامل یادگیرنده و محیط را با استفاده از حالات، اعمال و پاداش توصیف می‌کند [47]. در هر بازه زمانی t، عامل یادگیری تقویتی قادر است تا حالات محیط، St، را مشاهده نماید و براساس حالات مشاهده شده اعمال، At، را انجام دهد. در یک بازه زمانی بعد، به‌عنوان نتیجه عمل خود، عامل یک پاداش عددی، Rt+1، را دریافت می‌کند و به حالت جدید، St+1، می‌رود. بنابراین با استفاده عمل و عکس العمل با محیط، یک عامل یاد می‌گیرد اعمالی را انتخاب کند که پاداش خود را بیشینه کند. پاداش یک عدد است که با استفاده از تابع پاداش محاسبه میشود و مطابق با هدف مسئله یادگیری تقویتی تعریف میشود. هدف عامل هوشمند، بیشینه‌کردن کلیه پاداش‌های دریافتی در یک زمان طولانی است؛ حتی اگر پاداش لحظه‌ای منفی باشد؛ مانند رفتار بشر، عامل یادگیری تقویتی ممکن است پاداش لحظه‌ای را فدا کند تا در درازمدت پاداش تجمعی بیشتری به دست آورد.
در پروسه‌های تصمیم‌گیری مارکوف، احتمال هر ارزش ممکن برای حالت St و پاداش Rt، فقط به حالت و عمل لحظه قبل، St-1 و Rt-1، بستگی دارد. به عبارت دیگر، به حالات و اعمال اولیه وابسته نیست. حالت سیستم باید شامل کلیه اطلاعات مربوط به همه جنبه‌های عمل و عکس‌العمل محیط و عامل در گذشته باشد؛ بنابراین، یک حالت با این مشخصات دارای خاصیت مارکوف است.
یک پروسه تصمیم‌گیری مارکوف محدود شامل مجموعه اعمال، حالات و پاداش (S, A, R) با تعداد المان‌های محدود است. به‌منظور مدلسازی مسئله یادگیری تقویتی، دینامیک پروسه‌های تصمیم‌گیری با استفاده از تابع توزیع احتمال شرطی p به‌صورت زیر مشخص میشود [47]:
|
(8) |
|
برای همه . تابع شرطی احتمال ارزش متغیرهای تصادفی و در لحظه t است که فقط به حالت و عمل لحظه قبل وابسته است. به عبارت دیگر، زمانی که حالت و عمل قبلی داده میشود، مدل، حالت و پاداش بعدی را محاسبه می‌کند.
ارزش مورد انتظار جمع وزن دار پاداش‌ها در حالت s با سیاست π، تابع ارزش گفته میشود و به‌صورت زیر تعریف میشود [47]:
|
(9) |
|
پارامتر ، نرخ تخفیف است. پروسه‌های تصمیم‌گیری مارکوف باید بین پاداش‌های لحظه‌ای و آینده مبادله کند. نرخ تخفیف ارزش فعلی پاداش‌های آتی را مشخص می‌کند. این پارامتر مشخص می‌کند عامل یادگیری تقویتی چقدر به پاداش‌های آینده دور نسبت به پاداش‌های لحظه‌ای اهمیت دهد و زمانی که پارامتر به یک نزدیک میشود، عامل یادگیرنده به پاداش‌های آینده به‌شدت توجه می‌کند.
سیاست هر عامل مشخص می‌کند چگونه یک عامل در یک حالت مشخص عمل کند. درواقع یک نگاشت از حالات مشاهده‌شده به اعمال انجام‌شده در آن حالات است.
تابع ارزش - عمل، ، بیان‌کنندۀ ارزش عمل انجام‌شده در حالت s تحت سیاست π است و به‌صورت زیر بیان میشود [47]:
|
(10) |
|
به طوری که ارزش مورد انتظار جمع پاداش‌های وزن‌دار (با ضریب تخفیف) در حالت s، عمل a و تحت سیاست π است.
برای حل مسئله یادگیری تقویتی، عامل یادگیرنده باید سیاستی را محاسبه کند که مقدار زیادی پاداش در طول زمان زیادی به دست آورد. یک سیاست، سیاست بهینه گفته میشود، اگر جمع مورد انتظار پاداش‌ها بزرگ‌تر یا مساوی سایر سیاست‌ها باشد. سیاست‌های بهینه، تابع ارزش - عمل بهینه یکسانی را به اشتراک می‌گذارند و به‌صورت نمایش داده میشود. اگر یک عامل، تابع را برای هر حالت s داشته باشد، می‌تواند به‌راحتی هر عملی را که را بیشینه کند، پیدا کند.
حتی اگر یک مدل دقیق و کامل از دینامیک محیط موجود باشد، محاسبه سیاست بهینه با حل معادله بهینگی بلمن به‌راحتی امکان‌پذیر نیست. برای محاسبه سیاست بهینه از قانون بروزرسانی یادگیری Q استفاده میشود [48, 49]. یادگیری Q اساساً یک پروسه تصمیم‌گیری مارکوف است. فرض کنید از حالت S با اعمال A به حالت بعدی S انتقال داده شود و پاداش R دریافت شود؛ بنابراین، قانون به‌روزرسانی یادگیری Q به‌صورت زیر تعریف میشود [48]:
|
(11) |
|
به طوری که نرخ تخفیف و نرخ یادگیری است. تابع بهینه Q با استفاده از تابع Q یاد گرفته شده در رابطه 11 به‌صورت مستقیم تخمین زده میشود. یادگیری Q یک روش یادگیری تقویتی مستقل از مدل است. در صورتی که همه جفت حالت - عمل به‌صورت پیوسته به‌روزرسانی و گسسته‌سازی شود، با احتمال یک به مقدار ارزش - عمل بهینه همگرا میشود. هر عامل در تمامی حالات می‌باید به‌صورت تکرارپذیری آزمایش شود تا یک تخمین معتبر از پاداش مورد انتظار به دست آید.
3-3- مدیریت انرژی ریزشبکۀ پیشنهادی
منابع انرژی پراکنده و مشترکین برق، عاملهای مستقل و هوشمند در نظر گرفته می‌شوند. عامل‌ها دارای توانایی یادگیری هستند و می‌توانند با انتخاب تصمیمات درست، سود خود را بیشینه کنند. عاملهای یادگیری تقویتی با استفاده از بازخوردهای اعمال و تجربیات خود، سیاست بهینه را کشف می‌کنند. با توجه به خاصیت متغیر با زمان خروجی منابع انرژی تجدیدپذیر و تصادفی‌بودن مقدار بار مصرفی، از پروسه‌های تصمیم‌گیری مارکوف برای مدلسازی رفتار تصادفی عامل‌ها در ریزشبکه استفاده شده و سیاست بهینه عامل‌ها با الگوریتم مستقل از مدل یادگیری Q به دست آمده است. فرض شده است شبکه در مد متصل به شبکه اصلی کار می‌کند.
هدف مسئله یادگیری تقویتی، بهینه‌کردن توابع هدف (1)- (3) است؛ به طوری که قیود (4)-(7) نیز برآورده شود. حالات، اعمال و پاداش در ادامه به‌تفصیل توضیح داده می‌شوند.
حالات: برای کلیه عامل‌ها به‌جز باتری، حالات شامل است. t بازه زمانی، قیمت برق به هنگام فروش برق به شبکه اصلی و قیمت خرید برق از شبکه اصلی است. عامل باتری علاوه بر حالات فوق دارای یک حالت اضافی شامل سطح شارژ باتری است که مقدار آن بین 0 تا 100 درصد تغییر می‌کند.
اعمال: مجموعه اعمال برای ریزتوربین و سلول سوختی شامل مقدار توان الکتریکی تولیدشده، مقدار توان گرمایی تولیدشده، قیمت پیشنهادی فروش انرژی الکتریکی و گرمایی به ریزشبکه است. دیزل ژنراتور نیز دارای قدرت تصمیم‌گیری بر میزان توان تولیدی الکتریکی و قیمت پیشنهادی فروش انرژی است. مجموعه اعمال منابع انرژیهای انرژی تجدیدپذیر شامل توربین بادی و پنل خورشیدی فقط شامل قیمت پیشنهادی فروش انرژی است؛ زیرا توان خروجی منابع انرژیهای تجدیدپذیر کاملاً وابسته به شرایط آب‌وهوایی است و متغیرهای تصادفی هستند. برای بیشینه‌کردن سود منابع تجدیدپذیر فرض میشود آنها در حالت ردیابی نقطه حداکثر توان هستند. در صورتی که قیمت پیشنهادی منابع انرژی تجدیدپذیر با قیمت پیشنهادی منابع غیرتجدیدپذیر برابر باشد، اولویت فروش با منابع انرژی تجدیدپذیر است. ریزشبکه ابتدا توان تولیدشده با توربین بادی و پنل خورشیدی را خریداری می‌کند و در صورت نیاز از انرژی سایر منابع استفاده می‌کند. اگر مقدار توان تولیدی منابع انرژی بیش از نیاز ریزشبکه باشد، عامل‌ها می‌توانند انرژی مازاد را مستقیماً به شبکه اصلی بفروشند؛ اما با توجه به اینکه میزان قیمت خریداری‌شده توسط شبکه اصلی بسیار کمتر از قیمت فروش انرژی است، کلیه عاملهای تولیدی می‌باید به نحوی آموزش ببینند که با ارائه قیمت صحیح در بازار رقابتی برق، بتوانند توان تولیدشده خود را در داخل ریزشبکه بفروشند. درواقع، این موضوع کمک می‌کند تا ریزشبکه‌ها توان مورد نیاز خود را به جای خرید از شبکه اصلی از تولیدکننده‌های داخلی خریداری کنند؛ درنتیجه، ضمن افزایش سود تولیدکنندگان داخلی، از وابستگی ریزشبکه به شبکه اصلی نیز کاسته میشود. مجموعه اعمال باتری شامل حالت شارژ یا تخلیه، میزان توان مبادله‌شده و قیمت پیشنهادی است. در حالت شارژ توان باتری، منفی و در حالت تخلیه، مثبت است. میزان تقاضا فرض میشود یک متغیر تصادفی با تابع توزیع نمایی است. می‌توان تقاضا را به دسته بارهای غیرکنترل‌پذیر و کنترل‌پذیر تقسیم کرد. بروی دسته اول کنترلی نیست و در زمان تقاضا می‌باید برآورده شوند؛ اما میزان حذف بارهای کنترل‌پذیر با توجه به علاقه‌مندی آنها برای شرکت در مدیریت هزینه‌ها و الگوی مصرف قابل کنترل بوده و جزء مجموعه اعمال عاملهای مصرف‌کننده است.
پاداش‌ها: با توجه به اینکه هدف مسئله یادگیری تقویتی بیشینه‌کردن توابع هدف (1)-(3) است، پاداش لحظه‌ای به گونه‌ای تعریف میشود که توابع فوق را بیشینه کند؛ به همین دلیل، برای منابع انرژی پراکنده پاداش میزان سود لحظه‌ای حاصل از فروش انرژی است. پاداش مصرف‌کنندگان منفی صورتحساب برق مصرفی است. پاداش لحظه‌ای باتری در حالت شارژ منفی میشود و در این حالت باتری ممکن است برای جلوگیری از دریافت پاداش منفی هیچ‌گونه فعالیتی را انجام ندهد و پاداش دریافتی‌اش صفر شود. به‌منظور جلوگیری از تنبل‌شدن باتری، توابع پاداش لحظه‌ای باتری در حالت شارژ و تخلیه به‌صورت زیر تعریف و از یک ضریب اصلاح نیز استفاده میشود:
|
(12) |
|
ضریب مقداری بین صفر و یک دارد. به نحوی تنظیم میشود که سود باتری بیشینه شود. اگر یک باتری به‌خوبی آموزش ببیند، سود آن مقداری مثبت است. اگر سود باتری منفی شود، درواقع باتری انرژی را با قیمت بالا خریداری و با قیمت پایین فروخته است؛ بنابراین، باید آموزش باتری به نحوی صورت پذیرد که درنهایت سود باتری مثبت شود و میزان پاداش حاصل از فروش انرژی از هزینه خرید انرژی کمتر نشود. در هر لحظه باتری می‌تواند مصرف‌کننده یا تولیدکننده باشد و در آن واحد نمی‌تواند دارای هر دو حالت باشد. بیان‌کنندۀ هزینههای ناشی از کاهش طول عمر باتری است. برای محاسبه ، ابتدا هزینه تخریب و کاهش طول عمر باتری به‌صورت زیر محاسبه میشود [30]:
|
(13) |
|
ضریب تخریب و متناسب با هزینه باتری است. SOH (State Of Health) حالت سلامت باتری است و مقداری بین صفر و یک دارد [30]:
|
(14) |
|
برای دو حالت شارژ و تخلیه باتری به‌صورت زیر محاسبه میشود [50, 51]:
|
(15) |
|
به طوری که پارامترهای ثابت ، و وابسته به مشخصه تخریب باتری هستند. ضریب براساس مشخصه شارژ باتری محاسبه میشود. همه پارامترهای فوق با استفاده از تست‌های عملی به دست می‌آیند. مقدار توان شارژ و تخلیه باتری است. سایر پارامترهای و نیز با استفاده از تست‌های عملی محاسبه میشود [30].
الگوریتم مدیریت انرژی ریزشبکه شامل دو مرحله است. در مرحله اول، عامل‌ها تابع Q را یاد می‌گیرند. توان خروجی توربین بادی و پنل خورشیدی با استفاده از توابع توزیع Weibullو Beta به‌ترتیب مدل می‌شوند [52, 53]. پارامترهای توزیع Weibullو Beta از داده‌های خروجی منابع انرژی تجدیدپذیر شبکه قدرت تخمین زده میشود. از مدل خروجی توربین بادی و پنل خورشیدی صرفاً برای تولید داده عملی بیشتر برای استفاده در مراحل یادگیری و ارزیابی الگوریتم استفاده میشود. در مرحله دوم توابع Q تخمین زده شده ارزیابی می‌شوند. اپراتور شبکه، یک عامل حراج‌کننده در نظر گرفته میشود که قیمت بازار و میزان توان خریداری‌شده از ریزشبکه را مشخص می‌کند. به‌طور خلاصه، روش مدیریت انرژی ریزشبکه در الگوریتم 1 آورده شده است.
در این بخش، سیستم مدیریت انرژی پیشنهادی برای ریزشبکه هوشمند با استفاده از الگوریتم یادگیری تقویتی و داده‌های خروجی انرژیهای تجدیدپذیر و داده‌های دریافتی از بازار برق ایران شبیه‌سازی شده است. داده‌های واقعی شامل توان خروجی توربین بادی و پنل خورشیدی با همکاری پژوهشگاه هوا خورشید دانشگاه فردوسی مشهد و شرکت برق منطقه‌ای در بهار و تابستان 1399 به‌صورت ساعتی جمع‌آوری شدهاند. در جدول 2، پارامترهای مجموعه داده‌ها برای 24 ساعت شبانه روز ارائه شدهاند. مقدار خروجی در جدول 2 نرمالیزه شده است. ریزشبکه پیشنهادی مطابق شکل 1 متشکل از منابع انرژی گرمایی و الکتریکی، باتری و بارهای مصرفی الکتریکی و گرمایی است. مشخصات منابع پراکنده مطابق جدول 3 است. 4 عامل مصرف‌کننده بار الکتریکی، 3 عامل مصرف‌کننده بار گرمایی و یک عامل مصرف‌کننده بار الکتریکی و گرمایی به‌ترتیب به ظرفیت 8، 4 و 8 کیلووات در ریزشبکه در نظر گرفته شدهاند. ظرفیت دیزل ژنراتور به نسبت ظرفیت منابع تجدیدپذیر کمتر در نظر گرفته شده است تا استفاده از منابع غیرتجدیدپذیر در شبکه‌های قدرت محدودتر شود. مجموع میزان توان تولیدی از مجموع میزان توان مصرفی نیز به دلیل مدیریت مصرف کمتر در نظر گرفته شده است.
|
الگوریتم 1: // مقداردهی اولیه مقداردهی اولیه پارامترهای یادگیری و K1 مقداردهی Q0=0 مقداردهی k1 // یادگیری برای k=1:k1 هر عامل حالت‌های محیط را مشاهده می‌کند. میزان تقاضا و خروجی منابع تجدیدپذیر با استفاده توابع توزیع احتمال انتخاب میشود. هر عامل یک عمل را به‌صورت تصادفی انتخاب می‌کند. هر عامل پاداش خود را مشاهده می‌کند. جدول Q هر عامل به‌روزرسانی میشود. پایان // ارزیابی مقداردهی اولیه k2 برای k=1:k2 هر عامل حالات محیط را مشاهده می‌کند. مقدار تقاضا و خروجی منابع انرژی تجدیدپذیر مشخص میشود. هر عامل بهترین عمل را از جدول Q خود انتخاب می‌کند. هر عامل میزان پاداش خود را مشاهده می‌کند. میزان سود هر عامل، محاسبه و عملکرد سیستم ارزیابی میشود. پایان |
مصرف‌کنندگان قادرند به میزان حداکثر 70 درصد مصرف خود را مدیریت کنند. 30 درصد باقیمانده به‌عنوان بار ضروری در نظر گرفته میشود که در زمان تقاضا حتماً تأمین می‌شوند.
یک روز به 24 بازه زمانی یک‌ساعته تقسیم شده است. در هر بازه، نرخ خرید و فروش از شبکه اصلی در بازه 150-1200 ریال بر کیلو وات ساعت قرار دارد. با توجه به میزان خرید و فروش انرژی در بازار برق ایران در سایت IREMA [45] میزان قیمت پیشنهادی توسط تولیدکنندگان محدوده بین 200-1300 ریال بر کیلووات ساعت تعیین شده است.
روش ارائه‌شده تحت چهار سناریو: بدون یادگیری - یادگیری تولیدکنندگان - یادگیری مصرف‌کنندگان - یادگیری همه عامل‌ها ارزیابی شده است. عملکرد هر عامل طی چهار سناریو، هر سناریو به‌مدت 80 روز و جمعاً به‌مدت 320 روز شبیه‌سازی شده است. در 80 روز اول هیچ‌گونه یادگیری وجود ندارد و کلیه بارهای درخواستی در همان زمان برآورده میشود و منابع انرژی پراکنده به‌صورت تصادفی یک عمل را انتخاب می‌کنند. در 80 روز دوم، تنها منابع پراکنده آموزش دیده‌اند و دارای توانایی تصمیم‌گیری هوشمند هستند. در 80 روز سوم، تنها عاملهای مصرفی توانایی یادگیری دارند و در 80 روز آخر کلیه عامل‌ها آموزش دیده‌اند. فاز آموزش به‌مدت 10،000 روز اجرا شده و فاز ارزیابی برای هر سناریو به مدت 10 روز و جمعاً به‌مدت 800 روز شبیه‌سازی شده است. میانگین نتایج ارزیابی الگوریتم یادگیری تقویتی بروی سیستم مدیریت انرژی فوق در شکل‌های 2 الی 10 نمایش داده شده است. در این قسمت، مدل تخریب باتری در نظر گرفته نشده است. میانگین مقدار سود و توان کلیه عامل‌ها برای چهار سناریوها در جدول 4 نمایش داده شده است. در جدول 4، هزینه برای یک عامل مصرف‌کننده و توان شامل کل بار درخواستی در ریزشبکه است.
با توجه به شکل 2 و3، اگرچه میانگین خروجی توربین بادی و پنل خورشیدی در سناریو دوم (80 روز دوم) و در سناریو چهارم (80 روز چهارم) تقریباً تغییری نداشته، سود آنها به‌طور درخور توجهی افزایش یافته است؛ زیرا در این سناریوها منابع تولیدکننده دارای توانایی تصمیم‌گیری هستند و می‌توانند تصمیمات هوشمندانه‌تری اتخاذ کنند. در شکل 4 الی 6، میانگین روزانه سود و توان دیزل ژنراتور، سلول سوختی و ریزتوربین نشان داده شده است. با توجه به آموزش عاملهای تولیدکننده در سناریو دوم و چهارم، سود عاملهای دیزل ژنراتور، سلول سوختی و ریزتوربین در این سناریوها نیز افزایش یافته است. نسبت سود به تولید در دیزل ژنراتور در سناریو اول و چهارم به‌ترتیب 239 و 254 است؛ بنابراین، اگرچه تولید در سناریو چهارم افزایش پیدا کرده، نسبت سود به تولید (طبق جدول 4) برای دیزل‌ژنراتور نیز افزایش یافته است. درواقع دیزل ژنراتور توانسته است به‌طور هوشمندانه تولید خود را به ساعتی که درخواست و هزینه خرید بالا هست، منتقل کند. همچنین، این عامل با ارائه قیمت معقول برای پیشنهاد فروش انرژی توانسته است انرژی بیشتری در داخل ریزشبکه به فروش برساند و سود خود را افزایش دهد. نسبت سود به تولید در سناریو اول و چهارم به‌ترتیب برای سلول سوختی 225 و 291 و برای ریزتوربین 173 و 248 است؛ بنابراین، همانند دیزل ژنراتور نیز این عامل‌ها توانسته‌اند با اکتشاف و بهره‌برداری از محیط در حین آموزش تصمیمات بهینه‌تری اتخاذ کنند. شکل 7 نتایج شبیه‌سازی باتری را نشان می‌دهد. همان‌طور که مشاهده میشود در سناریوهایی که باتری آموزش داده شده است، سود آن مثبت و در سایر زمان‌ها منفی است. سود منفی به این معنا است که باتری در بیشتر زمان‌ها، انرژی را با هزینه بالا خریداری و هنگامی که قیمت برق پایین بوده، اقدام به فروش کرده است.
جدول (2): پارامترهای توابع توزیع احتمال Beta و Weibull برای مدلسازی توان خروجی پنل خورشیدی و توربین بادی.
|
Hour |
Weibull parameter (a,b) |
Beta parameter (a,b) |
||
|
1 |
0.565282 |
2.032585 |
0.019076 |
10000000 |
|
2 |
0.568418 |
2.290494 |
0.019076 |
10000000 |
|
3 |
0.5185 |
1.998599 |
0.019076 |
10000000 |
|
4 |
0.549884 |
1.692578 |
0.019076 |
10000000 |
|
5 |
0.577571 |
1.628994 |
0.019076 |
10000000 |
|
6 |
0.58976 |
1.899791 |
0.019076 |
10000000 |
|
7 |
0.605378 |
2.135819 |
0.009781 |
0.293334 |
|
8 |
0.597392 |
1.44144 |
7.895576 |
18.8233 |
|
9 |
0.561618 |
1.165105 |
5.990502 |
6.199672 |
|
10 |
0.474976 |
0.840914 |
9.22076 |
4.262388 |
|
11 |
0.492669 |
1.00863 |
26.1694 |
5.562986 |
|
12 |
0.476329 |
1.134804 |
2.451834 |
0.310795 |
|
13 |
0.445124 |
1.005857 |
1.932275 |
0.066813 |
|
14 |
0.419577 |
0.88521 |
7.61E+13 |
1.27E+13 |
|
15 |
0.419628 |
0.997423 |
25.38013 |
6.917985 |
|
16 |
0.460638 |
1.097642 |
95.42691 |
68.56888 |
|
17 |
0.459141 |
1.060481 |
3.88E+14 |
5.17E+14 |
|
18 |
0.53346 |
1.170101 |
6.438968 |
32.54034 |
|
19 |
0.64936 |
1.376157 |
0.000455 |
0.107348 |
|
20 |
0.716691 |
2.589421 |
0.000455 |
0.107348 |
|
21 |
0.70835 |
2.524027 |
0.019076 |
10000000 |
|
22 |
0.697165 |
2.909614 |
0.019076 |
10000000 |
|
23 |
0.671788 |
2.851862 |
0.019076 |
10000000 |
|
24 |
0.604807 |
1.991805 |
0.019076 |
10000000 |
جدول (3): ظرفیت واحدهای تولید انرژی
|
DER |
Wind |
PV |
BESS |
MT |
FC |
Diesel |
|
Prated(kW) |
10 |
10 |
5 |
6 |
6 |
5 |
|
(ب) |
|
(الف) |
شکل (2): (الف) میانگین سود، (ب)میانگین توان تولیدی روزانه پنل خورشیدی
|
(ب) |
|
(الف) |
شکل (3): (الف) میانگین سود، (ب) میانگین توان تولیدی روزانه توربین بادی
|
(ب) |
|
(الف) |
شکل (4): (الف) میانگین سود، (ب) میانگین توان تولیدی روزانه دیزل ژنراتور
در سناریو چهارم نسبت به سناریوی دوم سود عاملهای تولیدی کاهش یافته است؛ زیرا در حالت چهارم مصرف‌کنندگان نیز دارای قابلیت تصمیم‌گیری هوشمند هستند. در شکل 8 و 9 نتایج مصرف‌کنندگان الکتریکی و گرمایی به ترتیب نمایش داده شدهاند. در سناریو سوم و چهارم مصرف‌کنندگان آموزش دیده‌اند.
|
(ب) |
|
(الف) |
شکل (5): (الف) میانگین سود، (ب) میانگین توان تولیدی روزانه سلول سوختی
|
(الف) |
|
(ب) |
شکل (6): (الف) میانگین سود، (ب) میانگین توان تولیدی روزانه ریزتوربین
|
(ب) |
|
(الف) |
شکل (7): (الف) میانگین سود، (ب) توان تولیدی روزانه باتری
برای مقایسۀ صحیح بین سناریوها از مقایسه نسبت هزینه به مصرف در سناریو اول و دوم با سناریو سوم و چهارم استفاده میشود. برای مصرف‌کنندۀ الکتریکی این نسبت در سناریو اول تا چهارم به‌ترتیب 132، 137، 110 و 124 است. کاهش این نسبت‌ها به این معنا هست که عامل مصرف‌کننده توانسته است به مدیریت مصرف بپردازد و مصرف خود را در زمان‌هایی که قیمت برق بالا هست، کاهش دهد و در زمان‌هایی که قیمت پایین است، بیشتر کند و نیازهای خود را برآورده کند.
|
(ب) |
|
(الف) |
شکل (8): (الف) میانگین هزینه، (ب) توان الکتریکی مصرفی روزانه مصرف‌کننده الکتریکی
|
(الف) |
|
(ب) |
شکل (9): (الف) میانگین هزینه، (ب) میانگین توان گرمایی مصرفی روزانه مصرف‌کننده گرمایی
|
(ب) |
|
(الف) |
شکل (10): (الف) میانگین سود روزانه شبکه اصلی، (ب) میانگین توان روزانه تحویل داده شده به ریزشبکه
با انتخاب عدد 10 به‌عنوان ضریب نارضایتی ( ) عامل‌ها توانسته‌اند بین کاهش هزینه‌ها و به تبعِ کاهش مصرف و نیز ایجاد نارضایتی و عدم راحتی خود مصالحه کنند. با مقایسه این نسبت برای عاملهای گرمایی در سناریوهای مختلف مطابق جدول 4، مشاهده میشود توضیحات فوق برای مصرف‌کننده گرمایی نیز صادق است (شکل 9). اگرچه هزینه مصرف‌کنندگان در سناریو چهارم نسبت به سناریو سوم اندکی افزایش یافته است، سود تولیدکنندگان در سناریو چهارم به طور درخور توجهی نیز رشد کرده است. با توجه به اینکه در یک ریزشبکه هدف این است که هم سود تولیدکنندگان افزایش یابد و هم هزینه مصرف‌کنندگان کاهش یابد، این اختلاف هزینه قابل اغماض است. همچنین، در سناریو چهارم از وابستگی ریزشبکه به شبکه اصلی کاسته شده است (رجوع شود به شکل 10). مطابق شکل 10، به‌مجرد اینکه تعداد عاملهای بیشتری در ریزشبکه آموزش داده میشود، سود شبکه اصلی کمتر میشود. در سناریو چهارم حتی سود منفی شده است. سود منفی به معنای این است که سود حاصل از فروش انرژی به ریزشبکه از هزینه انرژی خریداری‌شده از ریزشبکه کمتر است. همچنین، مشاهده میشود در سناریو آخر توان خریداری‌شده از شبکه اصلی نیز منفی است؛ یعنی مجموع توان دریافتی از شبکه اصلی از مجموع توان داده شده به شبکه اصلی کمتر شده است؛ درنتیجه، وابستگی ریزشبکه به شبکه اصلی نیز به‌طور چشمگیری کاهش یافته است.
در شکل 11 مقدار ساعتی سود/ هزینه و توان مصرفی/تولیدی عاملهای ریزشبکه نشان داده شده است. پنل خورشیدی فقط در طول روز از ساعت 7 الی 18 در تابستان قادر است انرژی تولید کند. در سایر زمان‌ها، توان خروجی و سود پنل خورشیدی صفر است. خروجی توان توربین بادی تقریباً در ساعات شبانه‌روز یکسان است؛ زیرا این نمودار میانگین خروجی یک توربین بادی طی 800 روز را نشان می‌دهد. در ساعات اوج مصرف بین ساعت 12 الی 20 به دلیل درخواست بیشتر، قیمت انرژی نیز زیادتر شده و بنابراین سود توربین بادی و سایر تولیدکننده‌ها شامل دیزل ژنراتور، ریزتوربین و سلول سوختی نیز افزایش یافته است. طبق انتظار در ساعات اوج مصرف، هزینه و توان مصرفی عاملهای مصرفی نیز افزایش یافته است. روش پیشنهادی یک روش مدیریت انرژی ساعتی است. در مقاله [54] برای محاسبه مصرف انرژی در آینده، با کمک الگوریتم لونبرگ -مارکوارت شبکه‌های عصبی، مصرف انرژی به‌صورت کوتاه‌مدت پیش‌بینی شده است.
در این قسمت، هزینه ناشی از کاهش طول عمر باتری و میزان تخریب آن پس از هر استفاده، محاسبه و نتایج با حالت بدون مدل تخریب مقایسه میشود. قبل از در نظر گرفتن مدل تخریب باتری، تعداد دفعات تعویض باتری در مدت 800 روز به‌طور میانگین 23.1 بار است. به دلیل تخریب حاصل از شارژ و تخلیه بیش از حد و استفاده نادرست از باتری، تعداد تعویض باتری زیاد شده است؛ بنابراین، با توجه به قیمت اولیه باتری و تعداد دفعات زیاد تعویض باتری در قسمت قبل، لازم است مدل تخریب باتری در نظر گرفته شود. پس از اضافه‌کردن هزینه ناشی از تخریب باتری در تابع پاداش (رابطه 12)، تعداد تعویض باتری برای 800 روز به‌طور میانگین به مقدار 0.8 کاهش یافته است. اگرچه سود باتری کمتر شده و نزدیک به صفر است، با توجه به کاهش تعداد دفعات باتری قابل اغماض است؛ زیرا هزینه خرید باتری بسیار زیاد است و نیز وجود باتری‌ها برای تأمین بارهای ضروری در هنگام قطع برق ضروری است.
روش پیشنهادی با الگوریتم مونت‌کارلو [55] مقایسه شده است. روش مونت‌کارلو بر مبنای کسب تجارب بسیار و شبیه‌سازی زیاد انجام میشود و درنتیجه، تخمینی که از تابع Q به دست می‌آورد، ادعا میشود به مقدار بهینه خیلی نزدیک میشود [47]؛ به همین دلیل، روش مناسبی برای مقایسه است و این قابلیت را دارند که در سیستم‌های با ساختار غیرمتمرکز پیاده‌سازی شوند. جدول 5 نتایج شبیه‌سازی را نشان می‌دهد. با توجه به جدول 5، سود باتری در این روش منفی شده و باتری نتوانسته است به‌خوبی آموزش ببیند. همچنین، سود دیزل ژنراتور و سلول سوختی نیز نسبت به حالت قبل کاهش یافته است.
جدول (4): میانگین نتایج الگوریتم مدیریت انرژی ریزشبکه (بعد از 800 روز اجرا)
|
انرژی تولیدی/ مصرفی روزانه (کیلووات ساعت) |
سود/ هزینه روزانه (ریال) |
|
||||||
|
چهارم |
سوم |
دوم |
اول |
چهارم |
سوم |
دوم |
اول |
سناریو |
|
70 |
70 |
70 |
70 |
37318 |
27318 |
52206 |
31284 |
PV |
|
124 |
124 |
125 |
124 |
48966 |
36919 |
55431 |
40117 |
Wind |
|
108 |
63 |
108 |
63 |
27448 |
13086 |
40462 |
15042 |
Diesel |
|
130 |
76 |
130 |
76 |
37795 |
15378 |
54000 |
17092 |
FC |
|
132 |
76 |
131 |
74 |
32737 |
11885 |
41595 |
12766 |
MT |
|
221 |
226 |
625 |
621 |
27321 |
24893 |
85706 |
81896 |
Elec. load |
|
99 |
100 |
257 |
258 |
21751 |
21582 |
60308 |
60608 |
Heat load |
|
2.0 |
10.7 |
1.7 |
10.7 |
142 |
-1746 |
324 |
-2333 |
Battery |
|
-125 |
-15 |
435 |
539 |
-14606 |
64976 |
335280 |
457950 |
Maingrid |
جدول (5): میانگین نتایج شبیه‌سازی مدیریت انرژی ریزشبکه با استفاده از روش [55] (بعد از 800 روز اجرا)
|
انرژی تولیدی/ مصرفی روزانه (کیلووات ساعت) |
سود/ هزینه روزانه (ریال) |
|
||||||
|
چهارم |
سوم |
دوم |
اول |
چهارم |
سوم |
دوم |
اول |
سناریو |
|
70 |
70 |
70 |
70 |
39581 |
26766 |
55515 |
31366 |
PV |
|
124 |
124 |
124 |
124 |
57912 |
36618 |
64538 |
39913 |
Wind |
|
118 |
64 |
118 |
63 |
27127 |
12969 |
47637 |
14959 |
Diesel |
|
142.4 |
77 |
142 |
76 |
29429 |
16249 |
48588 |
17109 |
FC |
|
143 |
76 |
143 |
75 |
40907 |
11799 |
54168 |
13097 |
MT |
|
197 |
201 |
619 |
623 |
26492 |
21640 |
86893 |
81897 |
Elec. load |
|
84.8 |
85 |
256 |
257 |
19474 |
17850 |
61025 |
60276 |
Heat load |
|
0.16 |
10.6 |
0.18 |
10.7 |
-14.2 |
-1607 |
-9.5 |
-2399 |
Battery |
|
-241 |
-108 |
351 |
485 |
-60870 |
23115 |
280610 |
441970 |
Maingrid |
|
(ب) |
|
(الف) |
شکل (11): (الف) میانگین سود/هزینه ساعتی، (ب) میانگین توان تولیدی/مصرفی ساعتی پنل خورشیدی، توربین بادی، دیزل، سلول سوختی، ریز توربین و بارهای گرمایی و الکتریکی
سود سایر عامل‌ها افزایش یافته است. برای مقایسه عادلانه دو روش از شاخص مقایسه Fairness Factor (FF) در مقاله [27] استفاده شده است. در این شاخص، سود ریزشبکه با توجه به سود کلیه عاملهای تولیدی و مصرفی محاسبه میشود. مقدار شاخص FF در سناریوی چهارم برای روش مونت‌کارلو 1.63 است و برای روش ارائه‌شده در این مقاله 1.87 است. شاخصFF برای روش مونت‌کارلو به‌طور درخور توجهی از روش پیشنهادی کوچک‌تر است؛ بنابراین، از مقایسه مقدار فاکتور FF در این دو روش، می‌توان نتیجه گرفت سود ریزشبکه در روش پیشنهادی با در نظر گرفتن سود کلیه عامل‌ها بهبود یافته است.
در این مقاله، یک روش جدید غیرمتمرکز برای مدیریت انرژی الکتریکی و گرمایی ساعتی یک ریزشبکه پیشنهاد شد. در این روش، با در نظر گرفتن عدم‌قطعیت در تقاضای بارهای الکتریکی و گرمایی، انرژی تجدیدپذیر و قیمت برق، یک سیستم مدیریت انرژی مستقل از مدل با استفاده از یادگیری تقویتی ارائه شد. برخلاف روش‌های سنتی مبتنی بر مدل که نیازمند یک تخمین‌گر عدم قطعیت است، این روش براساس یادگیری است و نیاز به یک مدل صریح از عدم قطعیت ندارد. دردسترس بودن اطلاعات، برای یک واحد کنترل‌کننده مرکزی یا حتی برای عاملهای همسایه، در عمل سخت است. با افزایش ابعاد شبکه‌های قدرت این مشکل بیشتر میشود؛ بنابراین، با استفاده از روش غیرمتمرکز پیشنهادی، مشکلات ناشی از پیچیدگی ارتباطات و محاسبات برطرف شد. عملکرد روش ارائه‌شده تحت چهار سناریو: بدون یادگیری، یادگیری تولیدکنندگان، یادگیری مصرف‌کنندگان و یادگیری همه عامل‌ها شبیه‌سازی شد. برای ارزیابی مدل پیشنهادی از داده‌های واقعی توربین بادی و پنل خورشیدی و داده‌های بازار برق ایران استفاده شد. در قسمت شبیه‌سازی مقاله نشان داده شد سود کلیه واحدهای تولیدی افزایش، هزینه مصرف‌کنندگان کاهش و رضایت‌مندی آنها افزایش پیدا کرده است. همچنین، روش ارائه‌شده، از وابستگی ریزشبکه به شبکه اصلی نیز کاسته است.
علاوه بر آن، قابلیت پیاده‌سازی روش پیشنهادی به‌صورت ساعتی برای مدیریت انرژی ریزشبکه‌ها نشان داده شده است. در پایان، اثبات همگرایی روش پیشنهادی به جواب بهینه یا نزدیک به بهینه به‌عنوان کارهای آتی پیشنهاد میشود.
سپاسگزاری
بدین‌وسیله از حمایت‌های پژوهشی پژوهشگاه نیرو و همچنین، از حمایت‌های پژوهشگاه هوا خورشید دانشگاه فردوسی مشهد برای جمع‌آوری داده‌های عملی، صمیمانه تشکر و قدردانی می‌کنم.
[1] تاریخ ارسال مقاله: 20/01/1401
تاریخ پذیرش مقاله: 27/09/1401
نام نویسندۀ مسئول: سعید شمقدری
نشانی نویسندۀ مسئول: ایران، دانشگاه علم و صنعت، تهران، دانشکده مهندسی برق
)