УДК 519.832.3 ББК 22.18
ТЕОРЕТИКО-ИГРОВОЙ ПОДХОД К АНАЛИЗУ ПЕРСПЕКТИВ РАЗВИТИЯ СИСТЕМЫ ЭЛЕКТРОСНАБЖЕНИЯ РЕГИОНА
Беленький А. С.1, Болкунов Д. С.2
(НИУ «Высшая школа экономики», Москва, Россия)
Рассматривается проблема построения элементов инструментария для количественного анализа перспектив развития системы электроснабжения региона, исходя из оценки администрацией региона существующего уровня обеспеченности электроэнергией населения и предприятий региона и тенденций развития этого региона. Предлагается математическая модель, с использованием которой перспективы инвестирования в строительство новых электростанций и хранилищ электроэнергии в регионе оцениваются в рамках бескоалиционной игры двух лиц на выпуклых многогранниках. В этой игре платежная функция одного из игроков является суммой билинейной функции и линейной функции векторных аргументов, а платежная функция другого игрока является билинейной функцией тех же векторных аргументов. Доказывается возможность отыскания точек равновесия в рассматриваемой игре из решения задач линейного программирования, образующих двойственную пару. Эти точки определяют в частности, приемлемые для обеих сторон объемы инвестирования в систему электроснабжения региона и договорные цены на электроэнергию для потребителей.
Ключевые слова: выпуклые многогранники, игры двух лиц, линейное программирование, равновесие Нэша, седловые точки, управление большими системами.
1 Александр Соломонович Беленький, доктор технических наук, профессор ([email protected]).
2Дмитрий Сергеевич Болкунов, магистр прикладной математики и информатики ([email protected]). 76
Введение
Электроэнергетика является одной из основных отраслей экономики любого государства, и успешность функционирования и развития этой отрасли определяется прежде всего тем, какие меры предусматривает государство для обеспечения стабильности поставок электроэнергии всем категориям пользователей по приемлемым для них ценам и для создания условий, стимулирующих инвестиции в эту отрасль, в частности, в строительство новых мощностей по производству электроэнергии. В последние два десятилетия существенное внимание стало также уделяться инвестированию в строительство хранилищ электроэнергии, особенно вследствие появления новых идей в использовании физических и химических процессов в этих системах, а также в связи с бурным развитием систем производства электроэнергии из энергии солнца и ветра. Кроме того, хранилища электроэнергии оказываются конкурентоспособными с пиковыми генераторами (т.е. с генераторами, используемыми для подачи электроэнергии потребителям в часы пик) как с экономической точки зрения, так и с точки зрения обеспечения бесперебойных поставок электроэнергии всем пользователям, подключенным к сетям электроснабжения.
Интерес инвесторов к системам производства, хранения и передачи электроэнергии в условиях рынка определяется прежде всего экономическими соображениями, т.е. тем, насколько инвестиции в эти сферы деятельности выгодны с точки зрения умножения вложенного в них капитала как в краткосрочной, так и в долгосрочной перспективе, что, в свою очередь, зависит от правил, по которым функционирует этот рынок. Поскольку правила игры на рынке электроэнергии, устанавливаемые регулятором рынка (чаще всего государством), существенно влияют не только на инвестиции в строительство новых электростанций и хранилищ электроэнергии (или новых блоков действующих электростанций и хранилищ), но и на оптовые цены на электроэнергию, формирующиеся на краткосрочных и долгосрочных рынках про-
дажи электроэнергии, регулятору рынка необходим аппарат для количественного анализа рынка, который позволил бы оценить последствия введения регулятором тех или иных правил работы на рынке с точки зрения возможных действий участников рынка.
В настоящей работе рассматривается проблема оценки объемов электроэнергии, продаваемых по прямым двусторонним договорам между производителями и потребителями электроэнергии на оптовом рынке электроэнергии в регионе, оптовых цен на нее и ожидаемых объемов инвестиций производителей как в строительство новых мощностей, так и в строительство хранилищ электроэнергии (или их блоков) в условиях конкретного рынка электроэнергии. В ней предлагается теоретико-игровой подход к решению этой проблемы, в основе которого лежит идея моделирования взаимодействия индивидуальных производителей электроэнергии (или группы производителей электроэнергии) в регионе с индивидуальными потребителями (или с группой потребителей электроэнергии) в регионе на определенном временном промежутке (в несколько лет) в форме бескоалиционной игры двух лиц. В этой игре множества допустимых стратегий каждого (из двух) игроков являются выпуклыми многогранниками, а платежные функции игроков являются суммой билинейной и линейной функций векторных переменных (для одного из игроков) и билинейной функцией тех же векторных переменных для другого игрока. Компонентами векторов этих переменных являются объемы электроэнергии, производимые группой производителей (которая, в частности, может состоять из одного крупного производителя), рассматриваемой в качестве одного из игроков, и цены на поставляемую электроэнергию, приемлемые для группы потребителей (которая может состоять из одного крупного потребителя), рассматриваемой в качестве другого игрока в игре.
В модели учитываются технологические особенности производства предполагаемых (производителями электроэнергии в регионе) к строительству электростанций и хранилищ электроэнергии (или их блоков), а также существующие технические ограничения, связанные с производством, хранением и переда-
чей электроэнергии. Показано, что отыскание равновесных объемов производства и хранения электроэнергии и (оптовых) цен на нее для потребителей эквивалентно отысканию седловых точек в некоторой антагонистической игре на выпуклых многогранниках допустимых стратегий игроков, являющихся множествами допустимых стратегий игроков в исходной бескоалиционной игре двух лиц. Отыскание же седловых точек в указанной антагонистической игре осуществляется решением задач линейного программирования, образующих двойственную пару, что позволяет рассматривать взаимодействие любого конечного числа производителей и потребителей электроэнергии в регионе и реали-зовывать поиск равновесных стратегий в игре с использованием стандартного программного обеспечения (для решения задач линейного программирования), позволяющего за приемлемое время решать задачи с миллионами переменных и ограничений [13].
Следует отметить, что хотя, наряду с прямыми долгосрочными договорами, аукционы по продаже электроэнергии представляют возможный (альтернативный) механизм взаимодействия производителей и потребителей электроэнергии, на сегодняшний день во многих странах (в частности, в странах северной и западной Европы) только прямые договора являются действующим, общепринятым механизмом по продаже больших объемов электроэнергии на перспективу. Обзор различных типов аукционов по продаже электроэнергии на короткую перспективу, в частности, на сутки вперед, предлагается в [2], а обсуждение аукционов по продаже электроэнергии на долгосрочную перспективу, которые существуют на некоторых рынках (к примеру, в России и в некоторых штатах США), можно найти в [22]. Следует отметить, что хотя в России у участников рынка электроэнергии существует возможность заключать прямые долгосрочные договора на продажу электроэнергии, в настоящее время основным способом взаимодействия производителей и потребителей электроэнергии являются аукционы по продаже электроэнергии и мощности [3].
Статья состоит из введения, двух разделов и заключения. В первом разделе описывается содержательная постановка задачи
и предлагается краткий обзор подходов к моделированию проблем инвестирования в новые мощности электроэнергии (электростанции и их блоки). Эти подходы тесно связаны с рассматриваемой в статье задачей и дают определенные представления о направлениях разработки математического аппарата для анализа и исследования крупномасштабных систем электроснабжения региона. Во втором разделе устанавливается эквивалентность отыскания точек равновесия в исходной бескоалиционной игре двух лиц отысканию седловых точек в некоторой вспомогательной антагонистической игре двух лиц на тех же множествах допустимых стратегий игроков, что и в исходной бескоалиционной игре двух лиц (описывающей взаимодействие производителей электроэнергии и ее потребителей на оптовом рынке электроэнергии в регионе). Отыскание же седловых точек в указанной антагонистической игре может осуществляться решением задач линейного программирования, образующих двойственную пару [1].
1. Содержательная постановка задачи и краткий обзор публикаций, близких по содержанию к проблематике статьи
1.1. СОДЕРЖАТЕЛЬНАЯ ПОСТАНОВКА ЗАДАЧИ
Рассмотрим взаимодействие производителя (или группы производителей) электроэнергии с потребителем (или с группой потребителей) электроэнергии в некотором регионе в какой-либо момент времени в условиях, когда производители могут заключать с потребителями только двусторонние договора на поставку им электроэнергии на последующие Т лет (период планирования).
На основании этих двухсторонних договоров производители электроэнергии принимают решения об инвестировании в строительство новых мощностей и хранилищ электроэнергии и/или в строительство новых блоков действующих электростанций и хранилищ электроэнергии. Отыскание взаимоприемлемых условий поставки и покупки электроэнергии в рамках указанных двухсторонних договоров можно осуществлять, рассматривая взаи-80
модействие договаривающихся сторон в форме бескоалиционной игры двух лиц, в которой в качестве первого игрока выступает группа всех производителей электроэнергии в регионе (которая, в частности, может состоять из одного крупного производителя), а в качестве второго - группа всех потребителей электроэнергии (которая, в частности, может состоять из одного крупного потребителя). В этой игре множеством допустимых стратегий каждого (из двух) игроков является некоторый выпуклый многогранник, задаваемый совместной системой линейных неравенств балансового типа, а платежные функции игроков представляют собой сумму билинейной и линейной функций двух векторных переменных (для первого игрока) и билинейную функцию тех же векторных переменных (для второго игрока), причем каждая из векторных переменных принадлежит одному из этих двух выпуклых многогранников.
При моделировании взаимодействия производителей и потребителей электроэнергии региона в форме игры предполагается, что оптовый рынок электроэнергии имеет узловую структуру, т.е. разбит на несколько непересекающихся частей, называемых узлами. Если конкретный потребитель и конкретный производитель электроэнергии находятся в одном узле, то ограничения на пропускную способность по передаче электроэнергии между этим производителем и этим потребителем предполагаются несущественными и не учитываются в модели. Если же они находятся в разных узлах, то, напротив, ограничения на пропускную способность линий электропередач, соединяющих узлы, в одном из которых находится рассматриваемый производитель, а в другом рассматриваемый потребитель, являются частью модели. Предполагается, что линии электропередач принадлежат некоторой компании, которая устанавливает цены на передачу электроэнергии между узлами, и производители электроэнергии не могут влиять на эти цены.
В рамках рассматриваемого взаимодействия с потребителями электроэнергии предполагается, что группа производителей электроэнергии (первый игрок) в начале каждого года периода
планирования инвестирует в строительство а) новых или в модернизацию имеющихся электростанций различных типов (к примеру, газотурбинных или угольных) и в строительство различных типов новых хранилищ электроэнергии (к примеру, химических) с учетом ограничений на максимально возможные ежегодные инвестиционные расходы, а также с учетом стоимости строительства новых электростанций, и б) новых хранилищ электроэнергии и/или модернизацию имеющихся хранилищ.
Непостоянство спроса на электроэнергию и высокая стоимость её хранения принадлежат к числу основных рыночных проблем электроэнергии, которые вынуждают производителей электроэнергии использовать одновременно генераторы базовой и пиковой нагрузки, а также хранилища электроэнергии. Ясно, что рассмотрение взаимодействия производителей и потребителей электроэнергии в течение каждого часа каждого года периода планирования представляет значительный практический интерес, так как позволяет учесть как дневные, так и сезонные перепады спроса на электроэнергию, которые могут влиять на решения о строительстве новых электростанций и хранилищ электроэнергии различных типов. Однако использование часовой дискретности в течение всего периода планирования связано с существенным увеличением размера оптимизационных задач, в форме которых удается математически описывать взаимодействие производителей электроэнергии и ее потребителей в регионе. Кроме того, достаточно достоверная оценка ожидаемого спроса потребителей в каждый час в течение каждого года периода планирования представляет принципиальные трудности.
Для того чтобы учесть дневные и сезонные перепады спроса и при этом избежать указанных трудностей, предлагается рассматривать некоторую «усреднённую» структуру спроса на электроэнергию. Идея этого усреднения состоит в разбиении каждого года из периода планирования на «типичные» (непересекающиеся) периоды взаимодействия производителей и потребителей электроэнергии, каждый из которых повторяется несколько раз в году, и отражающих, например, различия в объемах потребления
электроэнергии в регионе в разное время года. Каждый «типичный» период может длиться как несколько дней, так и несколько недель, что позволяет учесть суточные колебания уровня спроса на электроэнергию в рамках каждого «типичного» периода.
Для упрощения понимания дальнейших рассуждений можно, например, считать, что в каждом году периода планирования существуют четыре «типичные» недели, каждая из которых повторяется тринадцать раз подряд в году в пределах одного и того же квартала. Таким образом, в каждом году периода планирования, в каждом квартале года выбирается некоторая «типичная неделя» спроса на электроэнергию. Тем самым, охватываются 364 (364 = 13 ■ 7 ■ 4) дня в году из 365 или 366 дней, а 365-й день обычного и 366-й день високосного года считаются особыми днями или особыми «типичными» периодами, первый из которых повторяется один раз в году, а второй один раз в четыре года. При этом, например, в каждой из 13 недель квартала, которому принадлежит выбранная «типичная» неделя, уровень спроса на электроэнергию в каждый конкретный день недели (один и тот же для всех 13 недель), например, в среду, считается одним и тем же для всех 13 недель и полагается равным уровню спроса в среду «типичной» недели. (Ясно, что этот уровень спроса может быть разным для разных дней «типичной» недели, т.е. уровень спроса, например, во вторник «типичной» недели и уровень спроса в среду той же самой «типичной» недели могут быть различны.)
Производители выбирают объемы электроэнергии, поставляемые потребителям исходя из заявок потребителей с учетом ограничений на мощность электростанций и хранилищ электроэнергии, которые могут использоваться в каждом году периода планирования, а также с учетом возможности изменения за короткий промежуток времени уровня производства электроэнергии на различных типах электростанций. В свою очередь потребители выбирают цены, которые они готовы заплатить за единицу объема поставляемой им электроэнергии при условии поставки им определенного объема электроэнергии ее производителями.
Эти объемы и цены фиксируются в договорах между производителями и потребителями электроэнергии с взаимоприемлемым для договаривающихся сторон уровнем детализации их численных значений в каждом конкретном году периода планирования.
Возможность расчета равновесных объемов поставок электроэнергии и цен на нее в каждый «типичный» период каждого года, входящего в период планирования, позволяет каждой из сторон договора на стадии заключения договора выбрать желаемый для нее уровень детализации объемов поставок электроэнергии и цен на нее, вплоть до часового, в рамках каждого «типичного» периода (повторяющегося соответствующее число раз в каждом году периода планирования). Это может быть сделано решением задач распределения ресурсов в какой-либо постановке (в частности, например, в форме задач линейного программирования) с учетом важных для каждой из сторон особенностей (соответственно) производства и потребления электроэнергии в любом временном промежутке внутри каждого из повторяющихся «типичных» периодов. Кроме того, указанная возможность расчета равновесных объемов поставок электроэнергии и цен на нее позволяет договаривающимся сторонам предусмотреть возможность их пересчета по взаимному согласию при возникновении непредвиденных, в частности, форс-мажорных обстоятельств у какой-либо или у обеих сторон.
Необходимо заметить, что, в отличие от предлагаемой в настоящей работе модели, в известных авторам существующих математических моделях инвестирования в строительство новых мощностей электроэнергии, которые оценивают будущие инвестиции производителей по ценам, формирующимся в условиях конкуренции производителей при меняющихся спросе и ценах на энергоносители, не учитывается влияние крупных потребителей на цены на электроэнергию и, как следствие, на инвестиции производителей в строительство и/или модернизацию имеющихся мощностей по производству и хранению электроэнергии.
1.2. КРАТКИЙ ОБЗОР ПУБЛИКАЦИЙ, БЛИЗКИХ
ПО СОДЕРЖАНИЮ К ПРОБЛЕМАТИКЕ СТАТЬИ
В [7] конкуренцию п производителей на рынке электроэнергии в течение некоторого периода времени Т (более одного года) предлагается описывать в виде динамической игры п лиц, в которой каждый из игроков принимает как инвестиционные, так и производственные решения. В рамках этой игры в начале каждого года в течение периода Т производители принимают решения об инвестировании в новые мощности, которые могут быть введены в эксплуатацию в течение некоторого времени (определяемого возможностями строительства), выбирая эти решения из некоторых конечных множеств допустимых решений. При этом выбор каждого из производителей осуществляется с целью максимизации ожидаемой полезности от своих решений в условиях неопределенности, порождаемой двумя основными факторами. Во-первых, предполагается, что всем производителям известен лишь ожидаемый спрос на электроэнергию, а, во-вторых, предполагается, что часть производителей более мелких, чем конкретный производитель, принимают решения об инвестициях после того как это решение принимается указанным конкретным производителем и становится известным этим более мелким производителям.
По окончании периода инвестиций все производители участвуют в аукционе Курно по продаже электроэнергии на спотовом рынке с известной для всех функцией рыночного спроса на электроэнергию и с учетом построенных к этому моменту новых мощностей. При этом предполагается, что аукцион проводится ежечасно с момента начала торговли электроэнергией и до начала следующего периода инвестирования. По окончании этапа торговли электроэнергией в последнем периоде Т оценивается ожидаемая прибыль от ежегодной продажи электроэнергии в будущем (т.е. за пределами периода планирования). Предполагается, что ежегодный доход каждого производителя в будущем будет таким же, каким он ожидается в течение последнего периода Т с учетом «стоимости» денег в будущем (отражаемой введением
некоторого коэффициента дисконтирования). При этом будущий период предполагается бесконечным, так что ожидаемая прибыль за пределами периода планирования (т.е. за пределами периода Т) оценивается суммой некоторой убывающей геометрической прогрессии.
Авторы работы [7] обсуждают в ней результаты численных экспериментов на построенном ими модельном примере, предлагая собственные выводы о работе моделируемого ими рынка электроэнергии. В частности, они полагают, что а) с целью увеличения собственной прибыли за счет повышения цены на электроэнергию доминирующий на рынке производитель может поставить электроэнергию в объеме меньшем, чем ему позволяют его технологические и технические возможности, б) в то время как контрактные обязательства участников рынка влияют на объемы инвестиций, указанные выше неопределенности, учитываемые в модели, могут влиять на динамику инвестиций (например, приводить к их отсрочке).
Рассмотренная в [7] модель имеет ряд недостатков, которые ограничивают её применимость на практике. Во-первых, в ней не учитываются ни ограничения на передачу электроэнергии, ни технологические ограничения при производстве электроэнергии. Авторы [7] предполагают, что в каждый час года производство электроэнергии ограничено только существующей мощностью электростанций, в то время как такое ограничение не является единственным существенным ограничением, которое следует учитывать. Например, уровень производства электроэнергии на базовых электростанциях (угольных, ядерных, гидроэлектростанциях) не может быть быстро изменен, поскольку уровни производства электроэнергии в соседние часы оказываются зависимыми между собой, чего предложенная в [7] модель не учитывает. Во-вторых, предложенный авторами [7] алгоритм поиска равновесия в игре перебирает всевозможные реализации спроса и инвестиций производителей в новые мощности, что приводит к вычислительным сложностям уже в случае 4 производителей и 10 инвестиционных планов у каждого из них. Действительно, даже
в этом случае необходимо рассмотреть 104 возможных комбинаций инвестиционных планов, для каждой из которых необходимо решить задачу нелинейного математического программирования.
В [12] рассматривается взаимодействие п производителей электроэнергии, конкурирующих на рынке в течение двух временных промежутков (периодов). В первом периоде каждый производитель инвестирует в строительство новых мощностей, которые будут построены к началу второго периода, а также выбирает объемы производства электроэнергии, которые будут поставлены на рынок в первом периоде. Во втором периоде производители выбирают только объемы производства электроэнергии с учетом новых мощностей, построенных к началу второго периода.
Моделирование указанного взаимодействия предлагается осуществлять с использованием двух игровых моделей, каждая их которых отражает схему принятия решений производителями электроэнергии как по объемам инвестиций в строительство новых и модернизацию имеющихся электростанций, так и по уровню производства электроэнергии. В соответствии с первой моделью, все производители одновременно принимают решения как об инвестициях, так и об уровне производства электроэнергии в обоих периодах, стремясь максимизировать свою ожидаемую прибыль. В соответствии со второй моделью, все производители сначала выбирают объемы инвестиций в строительство новых электростанций, стремясь максимизировать ожидаемую полезность от продажи электроэнергии во втором периоде, и лишь затем одновременно выбирают уровни производства электроэнергии в обоих периодах, исходя из известных им объемов инвестиций их конкурентов.
В [12] доказано, что если все производители планируют свои действия в соответствии со второй моделью, одному из производителей оказывается выгодно принимать решения об инвестициях в соответствии с первой моделью, если остальные производители предполагают, что он принимает решение об инвестициях в соответствии с его оптимальной стратегией в игре, формулируемой на основе второй модели.
Для произвольного числа периодов Т (т.е. при Т > 2) авторами [12] рассматривается модификация только первой модели, и в рамках этой модификации все производители одновременно выбирают объемы инвестиций и выпусков во всех Т периодах, что обусловлено меньшей вычислительной сложностью поиска равновесия в игре, формулируемой на основе этой модели. В этом случае неопределенность будущего спроса описывается деревом возможных уровней спроса, где в каждый момент времени £ £ 1,Т, множество возможных уровней спроса (представленных в виде £-го уровня этого дерева) зависит от реализации спроса в момент £ — 1. В рамках рассмотренной авторами модификации первой модели все производители, стремясь максимизировать свою ожидаемую полезность, в начале первого периода планирования одновременно принимают решения об инвестициях во всех периодах, кроме последнего, а также выбирают объемы производства электроэнергии при всех возможных реализациях спроса в каждом периоде с учетом ограничений на максимальную мощность электростанций, построенных к этому периоду.
Поиск равновесия в предложенной в [12] модели для произвольного числа периодов сводится к решению линейной задачи о дополнительности (LCP), что позволяет рассматривать задачи большего размера, по сравнению с размерностями задач, которые позволяет рассматривать использование модели из [7]. Однако в общем случае решение задачи LCP (так же как и при использовании модели из [7]) может представлять вычислительные трудности из-за отсутствия эффективных (т.е. с полиномиальными оценками) методов решения этой задачи. Кроме того, в [12] не учтены, в частности, ограничения на изменение уровня производства электроэнергии, а также не рассмотрены ограничения на передачу электроэнергии.
В [18] исследуется влияние ряда экономических параметров (например, периода планирования и эластичности спроса на электроэнергию) на уровень инвестиций на оптовом рынке электроэнергии. С этой целью рассматривается модель кон-
куренции п производителей электроэнергии, формулирующаяся в форме некооперативной статической игры, в которой, как и в [12], неопределенность описывается деревом возможных событий. Используемая в [18] модель близка к модели, рассматриваемой в [12], однако в ней учтены некоторые особенности рынка, неучтенные в [12]. В частности, авторы [18] различают два типа спроса на электроэнергию: пиковый и базовый, и учитывают устаревание существующих мощностей.
Для изучения влияния некоторых параметров рынка на инвестиции в электроэнергетику Финляндии авторы [18] провели серию численных экспериментов, в которых исследовали, как изменится общий уровень инвестиций в электроэнергетику при изменении первоначальных параметров модели, значения которых взяты из [17]. На основании полученных результатов авторы считают, что небольшая зависимость между пиковым и базовым спросом на электроэнергию не влияет на инвестиции, однако учет устаревания и более длинный период планирования увеличивают общий объем инвестиций.
Основные недостатки предложенной в [18] модели такие же, как и у модели, предложенной в [12]: а) не учтены важные технические ограничения производства и ограничения на передачу электроэнергии; б) равновесия Нэша не могут быть найдены эффективным методом (с полиномиальными оценками сложности и времени решения соответствующих игровых задач).
В [21] исследуется целесообразность инвестиций в тепловые электростанции с учетом существующей в Германии государственной поддержки развития возобновляемых источников электроэнергии. Одна из форм такой поддержки проявляется в том, что электроэнергия, генерируемая ветряными электростанциями, оплачивается потребителями по фиксированному тарифу вне зависимости от цены электроэнергии, сформировавшейся на рынке. Такая поддержка способствует вытеснению пиковых генераторов, производящих дорогую электроэнергию (в сравнении с базовыми электростанциями). Однако производительность ветряных электростанций существенно зависит от погодных условий,
и владельцы таких электростанций, как правило, не могут обеспечить стабильной подачи электроэнергии на рынок, что приводит к необходимости использования электроэнергии, производимой тепловыми электростанциями, которые могут полностью покрыть спрос при отсутствии поступления электроэнергии от ветряных электростанций (при соответствующих неблагоприятных погодных условиях). Для того чтобы оценить, насколько существующие правила рынка, поддерживающие развитие возобновляемых источников энергии, в то же время создают необходимые стимулы для инвестирования в тепловые электростанции, авторы рассмотрели модель конкуренции производителей на рынке электроэнергии в форме некооперативной статической игры n лиц.
Рассматриваемая ими модель описывает конкуренцию производителей на протяжении недели, в каждый час которой игроки (производители электроэнергии) участвуют в аукционе Кур-но по продаже электроэнергии. Все производители одновременно в начале периода планирования выбирают уровни производства электроэнергии в каждый час недели, при этом каждый производитель максимизирует свою полезность с учетом ограничений на мощность электростанций, а также возможности изменения объемов производства электроэнергии в течение часа. Поиск равновесия в рассматриваемой в [21] игре сводится к решению смешанной задаче о дополнительности (MCP) [8].
Авторы [21] количественно оценивают целесообразность инвестиций в тепловые электростанции на рынке электроэнергии Германии по значению разности между средним значением стоимости инвестиций на единицу объема произведенной электроэнергии по так называемым скрытым ценам (которым соответствуют множители Лагранжа к ограничениям на мощность электростанций в системе ограничений решаемой авторами оптимизационной задачи) и фактическим значением стоимости инвестиций в какой-то конкретный вид электростанций на единицу объема произведенной электроэнергии. Такой подход к оценке целесообразности инвестиций применим только к небольшим фирмам, потому что существенное изменение предложения элек-
троэнергии со стороны крупного производителя может заметно уменьшить рассматриваемое значение стоимости инвестиций (на единицу объема произведенной электроэнергии). В [21] сравниваются значения указанного показателя целесообразности инвестиций а) при отсутствии ветряных электростанций и б) при удвоении текущего количества ветряных электростанций (которое ожидается в будущем) с учетом существующей производительности таких электростанций в Германии. Проведенное авторами [21] исследование показало, что увеличение числа ветряных электростанций в два раза по сравнению с существующим числом электростанций снижает указанный выше показатель целесообразности инвестиций по сравнению с его значением при отсутствии таких электростанций, например, в газотурбинные генераторы практически на 50%, в парогазовые генераторы на 30% и в угольные электростанции на 25%. Такое снижение инвестиционной привлекательности тепловых электростанций и в особенности пиковых генераторов в перспективе может привести к нехватке электроэнергии во время безветрия. Для того чтобы избежать появления подобной ситуации, но при этом сохранить стимулирование инвестиций в развитие возобновляемых источников энергии, авторы предлагают, в частности, обязать производителей электроэнергии из возобновляемых источников компенсировать нехватку электроэнергии в периоды неблагоприятной погоды, что создаст для таких производителей стимулы к строительству пиковых генераторов. Для избежания нехватки электроэнергии предлагается также развивать хранилища электроэнергии и развивать линии электропередач, стимулируя таким образом торговлю с другими регионами.
Предложенный в [21] подход к инвестированию в новые мощности не лишен недостатков. Во-первых, в рамках рассматриваемого подхода оказывается, что стимулы к инвестированию в тепловые электростанции существенно уменьшатся при увеличении числа ветряных электростанций, однако невозможно оценить, будет ли этих стимулов всё же достаточно для инвестирования или нет. Во-вторых, в рамках такого подхода невозможно
утверждать, что инвестирование в строительство новых средних и крупных электростанций вообще может быть выгодным, потому что такое инвестирование может существенно уменьшить цены на электроэнергию, что сделает инвестиции убыточными. В-третьих, в данной модели не учитываются ограничения на передачу электроэнергии. В-четвертых, для того чтобы найти равновесие Нэша в рассматриваемой авторами игре, необходимо решать задачу математического программирования (MCP), для которой неизвестны эффективные вычислительные алгоритмы.
В [20] исследовалось влияние неопределенностей в ценах на топливо и в законах, регулирующих нормы выбросов в атмосферу, на инвестиции производителей в строительство новых электростанций, с использованием модификации игры из [21]. В этой модифицированной игре предполагается, что производители имеют возможность инвестировать в строительство новых мощностей на протяжении нескольких лет (в [20] приводятся результаты численных экспериментов авторов с периодом планирования в 25 лет). Как и в [21], оценка будущих прибылей производителей в течении одного года производится на участке в несколько часов (авторы рассматривают 120 часов), где часовой шаг по времени используется при моделировании технических ограничений на производство электроэнергии (на максимальную мощность и на максимально возможное изменение уровня производства электростанции).
Автор работы [20] рассматривает два подхода к моделированию неопределенности. В рамках первого подхода исследуются несколько возможных сценариев изменения цен на топливо и законов, ограничивающих выбросы в атмосферу, где для каждого сценария находится равновесие в детерминированной игре. В рамках второго подхода рассматривается стохастическая игра, в которой неопределенность задана деревом возможных событий (так же как ив [12,18]), и для этой стохастической игры находится равновесие. Автор приводит результаты численных экспериментов с данными по рынку Германии и сравнивает равновесие в стохастической игре (используемой в рамках второго под-
хода) с вектором, каждая компонента которого является математическим ожиданием соответствующих компонент равновесий в детерминированных играх (решаемых в рамках первого подхода), а каждая игра соответствует одному возможному сценарию развития цен на топливо. На основании проведенных численных экспериментов автор делает вывод о том, что при использовании стохастической модели оценка ожидаемой прибыли оказывается ниже, чем при использовании детерминированных моделей из-за того, что производители откладывают инвестиции, однако общий уровень инвестиций оказывается выше, чем при оценке их объема по детерминированным моделям. Рассмотренные автором модели (как детерминированные, так и стохастическая) учитывают как технические ограничения на производство, так и влияние инвестиций производителей на будущие цены. Тем не менее, в рамках рассмотренных автором работы [20] моделей не учтены ограничения на передачу электроэнергии, а поиск равновесия в рассмотренных ими играх сводится к сложной с вычислительной точки зрения задаче о дополнительности (MCP), что делает предложенные модели трудно применимыми при решении задач большого размера.
В [15] исследуется целесообразность использования хранилищ электроэнергии в рамках анализа работы оптового рынка PJM (Пенсильвании, Нью-Джерси и Мэриленда). Авторы предполагают, что с 2010 года условия для инвестирования в хранилища электроэнергии не сильно изменились, поэтому они изучают то, как изменился бы рынок PJM в 2010 году, если бы туда были добавлены хранилища электроэнергии, по сравнению с тем, как функционировал рынок в действительности. Авторами рассматривается оптимизационная модель, описывающая торговлю электроэнергией в каждый час 2010 года, в которой минимизируются суммарные издержки всех производителей. В модели учитывается множество технологических ограничений, возникающих при производстве электроэнергии, а также ограничения на её передачу. Сформулированная на её основе задача оптимизации является задачей целочисленного линейного программирования.
Используя данные спроса за 2010 год, авторы верифицируют модель, показывая, к примеру, что средняя часовая ошибка в значениях между ценами (множителями Лагранжа к ограничениям на максимальную мощность электростанций в решаемой ими оптимизационной задаче), предсказанными по модели, и ценами, которые были в 2010 году на рынке, равняется нулю. Авторы исследуют, как изменились бы цены на рынке Р.]Ш в 2010 году, если бы тогда существовали большие хранилища электроэнергии. С этой целью они модифицируют предложенную ими модель с учетом существования хранилищ и проводят численные эксперименты с новой моделью при различных параметрах хранения электроэнергии, соответствующих нескольким типам хранилищ, а также при различной общей мощности воображаемых хранилищ. Результаты численных экспериментов авторов работы [15] показывают, что использование хранилищ может уменьшить цены на электроэнергию, при этом уменьшается и прибыль производителей. Однако с учетом стоимости строительства хранилищ электроэнергии общественное благосостояние тоже уменьшается, а прибыль компаний, оперирующих построенными хранилищами, становится отрицательной.
Использование оптимизационной модели позволяет авторам учесть множество ограничений на производство и передачу электроэнергии. Однако, в отличие от игрового подхода, предложенный ими подход не учитывает влияния производителей и потребителей на рыночные цены. Следует также иметь в виду, что результаты были получены авторами [15] на данных рынка за 2010 год, которые нельзя непосредственно, без соответствующих изменений перенести на настоящее и будущее время, так как они не учитывают, например, возможного строительства новых электростанций, устаревания существующих электростанций, будущих изменений спроса и цен на энергоносители, что существенно ограничивает применимость полученных результатов.
Для исследования влияния инвестиций в новые электростанции на цены на электроэнергию на оптовом рынке и на объемы выбросов вредных веществ в атмосферу при производстве элек-
троэнергии в [14] предлагается моделировать взаимодействие производителей электроэнергии в условиях конкуренции в виде некооперативной динамической игры, которая учитывает ограничения на максимальную мощность электростанции, ограничения на передачу электроэнергии и устаревание мощностей. Результаты проведенных в [14] численных экспериментов показывают, что в случае совершенной конкуренции рыночные цены оказываются ниже, а уровень выбросов в атмосферу меньше, чем при монополии или олигополии. Инвестиции на рынках электроэнергии в условиях совершенной конкуренции также изучаются в [9], где рассматриваются модели для риск-нейтрального и риск-предпочитающего инвестора, учитывающие, в частности, неопределенность будущего спроса и цен на энергоносители.
Проблемы, тесно связанные с отысканием оптимальных объемов инвестирования, но не относящихся к проблеме инвестирования непосредственно, рассмотрены в [10, 11, 16, 19]. Так, обсуждение различных методов поддержки возобновляемых источников энергии можно найти в [11,16], влияние гидроаккумулиру-ющих хранилищ электроэнергии на цены и благосостояние производителей и потребителей в условиях рынка (Германии) исследуется в [19], а обзор теоретико-игровых моделей, посвященный проблемам «умных» (smart grid) электрических сетей, представлен в [10].
В [3] предлагается математическая модель для планирования работы энергосистемы в условиях конкурентного оптового рынка электроэнергии и мощности. В рамках задачи планирования загрузки действующих генераторов электроэнергии и дополнительных генерирующих мощностей, в частности, осуществляется расчёт планов как так называемого «балансирующего рынка» (т.е. рынка электроэнергии, на котором покупка и продажа электроэнергии осуществляется в режимах приближенных к реальным режимам производства и потребления), так и плана работы генерирующего оборудования на сутки вперед, наряду с расчетом которых рассматривается задача выбора генерирующего оборудования на некотором временном промежутке существенно боль-
шем суток. Этот выбор осуществляется с учетом как технических возможностей генераторов, так и ценовых заявок, отражающих стоимость производства заявленного объема электроэнергии производителями электроэнергии по каждому из генераторов, работой которых они управляют. Приводятся математические формулировки задач оптимизации, которые решаются при составлении планов работы подключенных в сеть генераторов электроэнергии на следующие сутки вперед и на «оставшиеся» часы текущих суток, а также задач, которые решаются на неделю вперед и на конец рабочей недели с учетом подаваемых в течение недели ценовых заявок на производство электроэнергии. Эти заявки используются для определения значения верхней границы объемов производства электроэнергии, зафиксированных в этих заявках, которая (верхняя граница) является параметром в указанных двух задачах планирования (на неделю вперед и на конец рабочей недели), наряду с подаваемым каждым поставщиком электроэнергии значения стоимости пуска генерирующего оборудования, подключаемого к функционирующей сети. Задача планирования «балансирующего рынка» электроэнергии формулируется в сетевой постановке как задача минимизации суммарной стоимости электроэнергии, поступающей в сеть в течение некоторого периода времени, начинающегося в интересующий диспетчера момент времени текущих суток с учетом ценовых заявок производителей, балансовых ограничений по активной и реактивной мощности в узле, физических закономерностей, определяющих перетоки мощности по ветвям сети, а также ограничений по перетокам активной мощности по так называемым контролируемым сечениям, производственных возможностей генераторов, включенных в сеть, и по скорости набора и сброса нагрузки по каждому генератору.
Задача планирования рынка электроэнергии «на сутки вперед» формулируется как задача максимизации некоторой функции благосостояния, в которой помимо ценовых заявок производителей электроэнергии учитываются ценовые заявки её покупателей при ограничениях указанной выше задачи оптимизации
планирования «балансирующего рынка» и дополнительных ограничениях по ценовым заявкам покупателей электроэнергии.
В задаче планирования рынка электроэнергии «на неделю вперед» помимо генераторов, включенных в сеть, рассматриваются все генерирующие мощности, не находящиеся в ремонте. Эта задача формулируется как задача минимизации суммарной стоимости электроэнергии, большая часть ограничений в которой близка по структуре к системе ограничений задачи планирования «балансирующего рынка» электроэнергии, а остальные ограничения отражают требования а) по обязательному присутствию в сети резервов мощностей; б) по числу включенных генераторов в каждом узле сети; в) по распределению включений в сеть генераторов, относящихся к одному узлу сети, во времени, а также ограничения на общее число включений небольших генераторов электроэнергии, работающих в режиме теплофикации, за всю неделю. Эта задача является частично целочисленной (с булевыми переменными) задачей нелинейного программирования, для которой известны подходы к построению методов ее решения ([4,5]).
Поскольку реальная стоимость электроэнергии для участников оптового рынка, как поставляющих, так и потребляющих эту электроэнергию, помимо цен на рынке на сутки вперед и на «балансирующем рынке» зависит также от стоимости мощности, необходимо учитывать особенности этого «товара» при расчете указанной реальной стоимости электроэнергии как в краткосрочной перспективе (на месяц вперед, в ходе биржевых торгов по свободным ценам, разрешенных в России с 2009 года), так и в долгосрочной (на год и на несколько лет). В частности, в [3] предлагается процедура отбора проектов по строительству новых мощностей, исходя из предполагаемых сценариев роста потребления электроэнергии (как во временном разрезе, так и по регионам) по критерию полной стоимости обеспечения потребления электроэнергии (определяемого указанными сценариями роста потребления) с учетом строительных и топливных затрат, который подлежит минимизации при ограничениях, аналогичных
по структуре ограничениям задачи планирования рынка на неделю вперед.
2. Математическая формулировка задачи
Пусть
К - число узлов (непересекающихся частей рынка электроэнергии региона, внутри которых практически отсутствуют ограничения на передачу электроэнергии);
I - число производителей электроэнергии, каждый из которых может иметь электростанции и хранилища электроэнергии в одном или нескольких узлах из множества 1, К;
.] - число потребителей электроэнергии, каждый из которых может иметь энергопринимающие устройства в одном или нескольких узлах из множества 1, К;
Т - продолжительность (в годах) периода планирования;
Mt - число «типичных» периодов взаимодействия производителей и потребителей электроэнергии в году Ь, где Ь € 1, Т;
- число раз, которое повторяется «типичный» период т в году Ь, где т € 1, М4, Ь € 1, Т;
Ньт - длина «типичного» периода т года Ь (в часах), где т € ТМ, Ь € 1, Т;
5 - число различных типов электростанций, которые могут уже иметь или построить производители в течение периода Т (при этом предполагается, что каждый из производителей электроэнергии 1,1 в каждом узле 1, К может иметь или построить не более одной электростанции каждого из 5 типов, однако каждый производитель может увеличивать мощность существующих электростанций, достраивая новые блоки);
К - число различных типов хранилищ электроэнергии, которые могут уже иметь или построить производители в течение периода Т (при этом предполагается, что каждый из производителей электроэнергии 1,1 в каждом узле 1, К может иметь или построить не более одного хранилища электроэнергии каждого из К типов, однако каждый производитель может увеличивать
мощность существующих хранилищ электроэнергии, достраивая новые блоки);
4hk2 (t) - максимальный объем электроэнергии, который может быть передан в течение одного часа между узлами k\ и k2 (ki = k2) в году t (индекс йг соответствует английскому термину throughput), где k1, k2 е 1, K, t е 1, T;
Lfc°fc2 (t) - стоимость передачи 1 МВт электроэнергии между узлами k1 и k2 в году t, где k1, k2 е 1, K, t е 1, T;
Lfc°Sfc2 - постоянная (в течение периода планирования 1, T) доля передаваемого объема электроэнергии, которая теряется в процессе её передачи между узлами k1 и k2, где k1, k2 е 1, K;
9isk ( ) - функция стоимости строительства генерирующей мощности (новой электростанции или нового блока существующей электростанции) типа s производителем i в узле k от её будущей мощности (в МВт), где i е 1,1, s е 1, S, k е 1, K;
CcpSk - приблизительная стоимость строительства в расчете (приходящаяся) на 1 МВт генерирующей мощности (новой электростанции или нового блока существующей электростанции) типа s производителем i в узле k (индекс p соответствует английскому термину power plant), где i е 1,1, s е 1, S, k е 1, K;
CliSk - максимальная мощность новой электростанции (или нового блока существующей электростанции) типа s, которая (который) может быть построена (построен) производителем i в узле k в течение одного года, где i е 1,1, s е 1, S, k е 1, K;
/irk ( ) - функция стоимости строительства емкости для хранения электроэнергии (нового хранилища или нового блока существующего хранилища) типа r производителем i в узле k от его будущей мощности (в МВт), где r е 1, R, i е 1,1, k е 1, K;
Ccbrk - приблизительная стоимость строительства в расчете (приходящаяся) на 1 МВт емкости для хранения электроэнергии (нового хранилища или нового блока существующего хранилища) типа r производителем i в узле k (индекс b соответствует английскому термину battery), где r е 1, R, i е 1,1, k е 1, K;
Clbirk - максимальная мощность нового хранилища электроэнергии (или нового блока существующего хранилища электро-
энергии) типа r, которое (который) может быть построено (построен) производителем i в узле k в течение одного года, где i G ÏJ, r G ITS, k G ;
zïsk (t) - мощность (в МВт) новой электростанции (или нового блока существующей электростанции) типа s в узле k, которая (который) будет построена (построен) производителем i в течение года и будет введена (введен) в эксплуатацию к началу года t, где s G IS, k G IK, i G Ï7, t G Ï7T;
zbrk (t) - мощность (в МВт) нового хранилища электроэнергии (или нового блока существующего хранилища электроэнергии) типа r в узле k, которое (который) будет построено (построен) производителем i в течение года и будет введено (введен) в эксплуатацию к началу года t, где r G I, R, k G Ï, K, i G I, I, t G I,T;
Bj(t) - максимальный объем инвестиций в строительство новых электростанций и хранилищ электроэнергии, который может сделать производитель i в году t, где i G I, I, t G I, T ;
WSk ( ) - функция оценочной стоимости электростанции типа s в узле k, принадлежащей производителю i, от её мощности (в МВт) на момент окончания периода T, где s G I, S, i G I, I, k G IK;
wbrk ( ) - функция оценочной стоимости хранилища электроэнергии типа r в узле k, принадлежащего производителю i, от его мощности (в МВт) на момент окончания периода T, где r G I, R, i G I7I, k g ;
QpSk (t) - мощность (в МВт) или существующей электростанции типа s в узле k, которую производитель i рассчитывает эксплуатировать в году t, или строящейся электростанции, которая будет введена в эксплуатацию к началу года t периода планирования ITT, где s G ITS, i G I,!, k G I^K, t G IT ;
cfSk ( ) - функция издержек годового содержания электростанции типа s в узле k производителем i от установленной мощности этой электростанции, где s G I, S, i G I, I, k G I, K ;
c/irk ( ) - функция издержек годового содержания хранилища электроэнергии типа r в узле k производителем i от установлен-
ной мощности этого хранилища, где r G 1, R, i G 1,1, k G 1, K;
costisk(■) - функция издержек производства электроэнергии на электростанции типа s в узле k, принадлежащей производителю i, от объема часового производства электроэнергии, где s G 17S, i G 17т, k G 1K;
as - отношение максимально возможного размера изменения (увеличения/снижения) объема производства электроэнергии на электростанции типа s к объему выработки электроэнергии в момент т «типичного» периода (в часах), вычисляемое за время А (не превосходящего одного часа), такое, что выполняется неравенство т + А ^ т + 1;
Qirfc (t) - мощность (в МВт) или существующего хранилища электроэнергии типа r в узле k, которое планирует эксплуатировать производитель i в году t, или строящегося хранилища, которое будет введено в эксплуатацию к началу года t периода планирования 1, T, где s G 1, S„ где r G 1, R, i G 1,1, k G 1, K, t G ITT;
pr - доля электроэнергии, которая в среднем теряется за один полный час хранения в хранилище типа r, где r G 1, R;
Djk(t, m, h) - минимальный объем электроэнергии (в МВт), который может понадобиться потребителю j в узле k в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где j G 1, J, k G 1, K, t G 1, T, m G 1, Mt, h G 1, Htm;
Djk(t, m, h) - максимальный объем электроэнергии (в МВт), который может понадобиться потребителю j в узле k в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где j G 1, J, k G 1, K, t G 1, T, m G 1, Mt, h G 1, Htm;
PjTx(t, m, h) - максимальная цена за 1 МВт электроэнергии, которую готов заплатить потребитель j в узле k в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где j G 1, J, k G 1, K, t G 1, T, m G 1, Mt, h G 1, Him;
Pisk (t,m,h) - минимальная цена, по которой производитель i готов продать 1 МВт электроэнергии с электростанции типа s в узле k в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где i G 1,1, s G 1, S, k G 1, K,
t е l,T, m е l,Mt, h е i,Him;
Pirk (t, m, h) - минимальная цена, по которой производитель i готов продать 1 МВт электроэнергии из хранилища электроэнергии типа r в узле k в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где i е 1,1, r е 17R, k е TK, t е lTT, m е 17Ж, h е l, Him;
xpSkljk2 (t,m,h) - объем электроэнергии (в МВт), который поставит производитель i с электростанции типа s в узле ki потребителю j в узле k2 в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где i е 1,1, ki е TK, j е TJ, k2 е T,K, s е TS, t е T^T, m е 1ГЖ,
h е T Htm;
xbrkljk2 (t, m,h) - объем электроэнергии (в МВт), который поставит производитель i из хранилища электроэнергии типа r в узле ki потребителю j в узле k2 в течение часа, начинающегося в момент времени h из «типичного» периода m года t, где i е 1,1,
ki е lTK, j е TTJ, k2 е lTK, r е lTR, t е TTT, m е 1ГЖ, h е l, Htm;
xbSr^^(t,m,h) - объем электроэнергии (в МВт), который производитель i поставит с электростанции типа s в узле ki в собственное хранилище типа r в узле k2 в течение часа, начинающегося в момент времени h из «типичного» периода m года t (индекс bc соответствует английскому термину battery charge), где i е TTT, ki е 17K, s е lTS, r е lTR, k2 е IK, t е TTT, m е l,Mt, h е l,Htm;
qirk(t, m, h) - объем электроэнергии (в МВт), который находится в хранилище типа r, принадлежащем производителю i, в узле k в течение часа, начинающегося в момент времени h из «типичного» периода m года t (в рамках упрощения, связанного с использованием «типичных» периодов, предполагается, что в начале каждого «типичного» периода m все хранилища не заполнены, то есть могут вместить электроэнергию в объеме равном емкости хранилищ), где i е 1,T, k е 1,K, r е 1,R, t е 1,T, m е l,Mt, h е l,Htm;
yfskijk2 (t, m,h) - цена, которую готов платить потребитель
3 в узле за 1 МВт электроэнергии, поставленный с электростанции типа в в узле производителем г в течение часа, начинающегося в момент времени Н из «типичного» периода т года г, где 3 е 177, к2 € 17К, в е 1,5, е 17К, г е 177, г е 17Т, т € Т7Ж, Н € 1,Я<ТО;
УЬгк^къ (г, т, Н) - цена, которую готов платить потребитель з в узле ^2 за 1 МВт электроэнергии, поставленный из хранилища электроэнергии типа г в узле производителем г в течение часа, начинающегося в момент времени Н из «типичного» периода т года г, где з € 1,7, к2 € 1,К, г € 1,Д, к1 € 1,К, г е 1,1,
г е 17Т, т е 1,м4, н е Т7нт.
Предположение 1. Издержки производства электроэнергии, стоимость строительства электростанций и новых хранилищ электроэнергии, а также оценочная стоимость электростанций и хранилищ электроэнергии описываются выпуклыми (или вогнутыми) функциями соответствующих переменных (в рамках математической формулировки задачи, эти функции аппроксимируются кусочно-линейными функциями).
Предположение 2. Как указано выше, в разделе 1.1, в каждом году г периода планирования 17 Т предлагается рассматривать М^ «типичных» периодов взаимодействия потребителей и производителей, отражающих, например, различие в объемном потреблении электроэнергии в регионе в разное время года. Каждый «типичный» период т (т € ^М^) состоит из Нт часов, для каждого из которых имеются прогнозные значения спроса каждого потребителя на электроэнергию; при этом Нт часов могут в сумме составлять как несколько дней, так и несколько недель, что позволяет учесть дневные колебания спроса на электроэнергию. Предполагается, что «типичный» период т (т € ^М^) в году г повторяется раз, так что год г разбивается на непересекающиеся «типичные» периоды т^ ■ ■ ■ 7 тм1, и каждый «типичный» период т^ (з € ^М^) имеет длину (в часах) Н^., 3 € 17 Mt и повторяется в году г ^т^ раз, откуда следует, что
364 * 24 = £^=1 Я*™. ■ 7<т..
Взаимодействие производителей и потребителей электроэнергии в регионе математически может быть описано в форме игры двух лиц, в которой множество допустимых стратегий каждого (из двух) игроков является выпуклым многогранником. В этой игре приведенные ниже неравенства (1)-(10) описывают множество допустимых стратегий первого игрока (группы всех производителей, снабжающих регион электроэнергией), а неравенства (11)-(14) - множество допустимых стратегий второго игрока (группы всех потребителей электроэнергии в регионе).
Неравенства (1) и (2) отражают ограничения, связанные с производством электроэнергии.
В частности, неравенства (1) показывают, что объем электроэнергии, продаваемый производителем электроэнергии г с электростанции типа находящейся в узле к, всем потребителям электроэнергии в регионе в течение часа, начинающегося в момент Н «типичного» периода т в году Ь, и объем электроэнергии, направляемый производителем электроэнергии г с электростанции типа находящейся в узле к, во все используемые им хранилища в регионе в течение часа, начинающегося в момент Н «типичного» периода т в году Ь в сумме не могут превышать мощности этой электростанции, доступной для использования в году Ь:
J к д к
(1) т>Н) + ^2 £ (Ь> т>Н) < ^=1 к'=1 г=1 к'=1
< с* (¿) + Е ^ (О,
г е 17Г, к е Т,К, 8 е Ь е 1Т, т е Н е 1,Я*ТО, причем предполагается, что энергия вырабатываемая электростанциями, строительство которых было начато более чем за год до начала года Ь = 1 периода планирования 1, Т и будет закончено к началу года Ь ^ Т, учитывается в слагаемом О?^ (Ь).
Неравенства (2) показывают, что отношение изменения (увеличения/снижения) объема производства электроэнергии на электростанции типа в производителя электроэнергии г в узле к в момент времени Н + 1 «типичного» периода т года г к объему этого производства в предыдущий момент времени не может превышать
.К ЯК
(2) 1 £ £ & т7 Н + 1) + £ £ хЬгС8ткк'& ^ Н + 1) -
^=1 к'=1 г=1 к'=1
.К ЯК
<
! £ £ xi>s kjk' (t, m, h) + ^ у ^ у xbsrkk' (t, m,
h)
j=1 k' = 1 r=1 k'=1
J K R K
^ as ■ ( ^ ^ хр
kjk' (t, m, h) + ^ у ^ у xbsrkk' (t, m, h) j , j=1 k' = 1 r=1 k'=1
i G 1,1, k G 1, K, s G 1, S, t G 1, T, m G 1, Mt, h G 1, Him - 1, (каждое из этих нелинейных неравенств эквивалентно двум линейным неравенствам).
Равенства и неравенства (3)-(5) описывают функционирование хранилищ электроэнергии.
Равенства (3) показывают, что объем электроэнергии в хранилище электроэнергии типа r производителя электроэнергии i в узле k в начале любого часа h + 1 «типичного» периода m года t складывается из а) объема электроэнергии, имевшегося в этом хранилище в начале предыдущего часа h за минусом потерь от хранения в течение часа; б) объема электроэнергии, поступающего в хранилище между моментами h и h + 1, за вычетом объема электроэнергии, направленного производителем электроэнергии i из этого хранилища всем потребителям электроэнергии в регионе в течение часа, начинающегося в момент h и заканчивающего
в момент h + 1:
(3) qirfc(t, m, h + 1) = (1 - pr) ■ qirfc(t, m, h)+
S к J к
+ (1 TT) ' E E ^srfc'fc(t, m, h) — E E Xb>rfcjfc'(t,m,h),
s=1 fc'=1 j=1 fc'=1
i e 1,1, r e 1, R, k e 1, K, t e 1, T, m e 1, Mt, h e 1, Him - 1, здесь предполагается, что (1 — Pp) отражает средние потери электроэнергии, поступающей в хранилище типа r в течение часа, начинающегося в момент h (т.е. до момента h + 1).
Равенства (4) отражают предположение о незаполненности (пустоте) хранилища электроэнергии типа r производителя электроэнергии i в узле k в начальный час «типичного» периода m года планирования t:
(4) (Ь т, 1) = 0,
г е 1,7, г е к е ^К, Ь е :Г7Г, т е , (где в указанных равенствах последний индекс равен 1 в силу включения Н е 1, Я*т, а также в силу предположения, сделанного при описании параметра (Ь, т, Н)).
Неравенства (5) отражают ограничения на объем электроэнергии, хранимой в хранилище типа г производителя электроэнергии г в узле к в каждый час Н «типичного» периода т года планирования Ь, который не может превосходить суммарной
мощности этого хранилища в году Ь:
*
(5) ^(Ь, т, Н) < (Ь) + Е 4* (Ь'),
г е 1,7, г е к е Ь е 1Т, т е Н е 1,Я*ТО, предполагается, что энергия доступная для использования из хранилищ, строительство которых было начато более чем за год до начала года Ь = 1 периода планирования 1, Т и будет закончено к началу года Ь ^ Т, учитывается в слагаемом (Ь).
Значения уровней (объемов) максимального и минимального часового спроса потребителя ] на электроэнергию в каждом
«типичном» периоде являются правыми и левыми частями двусторонних неравенств (6) соответственно. При этом учитывается, что потребители могут получать электроэнергию, как с электростанций, так и из хранилищ электроэнергии, с учетом её потерь при передаче по линиям электропередач:
(6)
/К , Б
,т7 Н) < £ £ (1 - ■ £х?,т7 Н) + ¿=1 к1 = 1 ^8=1
Я \ _
+ £ хЬгкик2 т7 НЛ < 7 т7 Н) 7
г=1 '
з е 1777 к2 е 17X7 г е 17^ т е 17Ж7 н е 17Н7т.
Неравенства (7) отражают ограничения на пропускную способность линий электропередач как при передаче электроэнергии от производителей к потребителям, так и при передаче электроэнергии с электростанций в хранилища электроэнергии в регионе:
/ . Б / Б Я
(7) £ £ £ <кик2т7 Н) + £ £ £ хЬсгк1к2& т Н)+
¿=1 ,7 = 1 8=1 ¿=1 8=1 Г=1
/ . Я
+£ £ £ хи,к2 (г т7 н) < (г)7
¿=17=1г=1
к17 к2 € 17К (к1 = к2>7 г € 17^ т € ИТМ Н € 17Н7т. Неравенства (8) отражают бюджетные ограничения на инвестиции в строительство новых электростанций (или в строительство новых блоков существующих электростанций) и в строительство новых хранилищ электроэнергии (или в строительство новых блоков существующих хранилищ электроэнергии) в течение периода планирования 17 Т для каждого производителя электроэнергии, снабжающего электроэнергией пользователей с электростанций, а также из хранилищ, находящихся на территории
региона: (8)
Е (Е ■ (Ь + 1) + Е Сс'гк ■ ^(Ь + 1)) < £*(Ь),
г е 1,7, ь е 1,Т - 1. Неравенства (9) и (10) отражают ограничения на максимальную мощность новых электростанций (или новых блоков существующих электростанций) и новых хранилищ электроэнергии (или новых блоков существующих хранилищ электроэнергии), которые могут быть построены за 1 год:
(9) ^(Ь) < , г е 1,7, к е 8 е 1,5, ь е 1,7,
(10) (ь) < С&д., г е 1,7, к е 1,К, г е 1,Д, Ь е 1,7.
Неравенства (11) и (12) отражают ограничения на максимально возможные цены на электроэнергию, которые готов заплатить потребитель за 1 МВт электроэнергии (получаемой как с электростанций, так и из хранилищ электроэнергии) в течение часа, начинающегося в момент Н «типичного» периода т года Ь:
(11) (Ь,т,Н) < p:;m2ax(í,m,н),
г е 1,7, 8 е 1,5, к1, к2 е 1,К, ; е 1,7,
Ь е 1,Т, т е 1,М*, Н е 1,Я*ТО,
(12) (Ь,т,Н) < (Ь,т,Н),
г е 1,7, г е к1, к2 е 1,К, ; е 1,7,
Ь е 1Т, т е 1,М*, Н е 1,Я*ТО.
Неравенства (13) и (14) отражают ограничения на минимально возможные цены, по которым производитель готов продать 1 МВт электроэнергии (как с электростанций, так и из хранилищ электроэнергии) в течение часа, начинающегося в момент Н «типичного» периода т года Ь: 108
(13) Ск17к2(^тН) > РЙТ^Н)^ _
г € 17/7 8 € 17 ^ к17 к2 € 17 X з € 17 77
г € 17^7 т е 17м*7 н € 17x^7
(14) уРгк17к2 (^ > рй?^ т7 Н)_ _
г € 17/7 г € к17 к2 € 17 X з € 17 77
г е 17^7 т е н е ^я^.
Предполагается, что система неравенств (1)-(14) совместна, что легко проверяется, например, по методике, предложенной в [1]. Суть этой методики применительно к исследованию совместности общей системы линейных неравенств с действительными коэффициентами - в построении и решении вспомогательной задачи линейного программирования с заведомо совместной системой ограничений, являющейся некоторой модификацией исследуемой на совместность системы неравенств. Значение этой вспомогательной задачи дает ответ на вопрос о том является ли исследуемая система неравенств (в рассматриваемом случае система (1)-(14)) совместной. Если исследуемая система неравенств несовместна, то вектор решения вспомогательной задачи определяет, как следует изменить правые части (каких-либо или всех) неравенств из исследуемой системы с тем, чтобы так «скорректированная» система неравенств была совместной. При этом корректировка правых частей неравенств исследуемой системы производится таким образом, чтобы «скорректированная» (совместная) система неравенств и исследуемая (несовместная) система неравенств были в определенном смысле «наиболее близкими» друг к другу. В случае системы линейных неравенств (1)-(14) указанная «близость» понимается в смысле незначительности (или приемлемости) структурных изменений в действующей системе электроснабжения региона (соответствующих вектору правых частей «скорректированной» системы неравенств
(1)-(14)) по сравнению с предложенными (и определенными значениями правых частей неравенств из несовместной системы (1)-(14)) изменениями.
Предполагается, что (возможно несовместная) система неравенств (1)-(14) в реальных задачах — отражающая основные количественные соотношения, определяющие структуру как существующей системы электроснабжения региона, так и планируемых к постройке новых систем электроснабжения и/или новых блоков существующей системы — не содержит неравенств, противоречащих здравому смыслу. В этом случае применение методики, предложенной в [1], позволяет найти приемлемые для инвесторов или, по крайней мере, интересные для них варианты развития действующей системы электроснабжения региона, рассматривая в качестве таких вариантов решения «скорректированной» системы ограничений.
Следует отметить, что выяснение факта совместности или несовместности системы линейных неравенств (1)-(14) принципиально может осуществляться классическими методами линейной алгебры, приведенными, например, в [6]. Однако смысл применения предложенной в [1] методики не только в установлении самого факта совместности или несовместности системы (1)-(14), но и в отыскании приемлемой для инвесторов «корректировки» исходной системы балансовых ограничений (если она несовместна), которой является система неравенств (1)-(14). Более того, как установление факта совместности или несовместности системы неравенств (1)-(14), так и отыскание «оптимальной» корректировки несовместной системы в соответствии с методикой легко реализуется алгоритмически и программно в системах поддержки принятия решений с использованием стандартного программного обеспечения для решения задач линейного программирования. Возможность применения этого программного обеспечения при исследовании на совместность систем ограничений вида (1)-(14) и при «корректировке» несовместных систем указанного вида в соответствии с предложенной в [1] методикой является ее существенным достоинством. Именно, применение
этой методики позволяет находить «приемлемые» (в указанном выше смысле) решения, в частности, в прикладных задачах развития систем электроснабжения региона с практически любыми интересными для практики размерами линейных моделей с (разумно составленными, но возможно несовместными) ограничениями вида (1)-(14). Пусть
xp =«111(1,1,1), ... , xpskjk (T, mt , Htmt )), = (x1nn(1,1,1), ••• , xjrkjk (t,mt,htmt )), xbc =( (x?cmi(1,1,1), ••• , xisrkk(t,mt,htmt)),
q = ((m(1,1,1), ■ ■ ■ , q/як)T, MT, htmt)),
zp =(z?n(1), ••• , zPsk (T)), =( zbii(1), ••• , zjrk (T)),
z = (zp, zb),
x = (xp, xb, xbc, q, zp, zb), УР = (yfxiii (1,1,1), ••• , y?SKJK (T,MT ,htmt )), = (y?iiii(1,1,1), ••• , y^RKJK (t,mt ,htmt )), У = (yp, yb),
cost(.) = (cosim(.), ••. , cost/sK(■); ■ ■ ■
••• ; cosi111(-), ••• , cost/sK (■)), где набор чисел (cosim(-), ■ ■ ■ , cost/sK(■)) повторяется ^^ 1 Him раз,
#(•) = (gmO^••• ,g/sK(0; ••• ; giiiO^••• ,g/sK(0), где набор чисел (g111( ), ■ ■ ■ ,g/sK( )) повторяется T раз,
/О = (/lll(•),••• ,f/RK(-); ••• ; /шО^••• ,f/RK(0), где набор чисел (/111( ), ■ ■ ■ , //RK( )) повторяется T раз,
/(■) = (c/fii(-),••• ,c//s^(-); ••• ; c/fn(.),■■■ ,c/psk(■)), где набор чисел (c/fn(-), ■ ■ ■ , c//SK(■)) повторяется T раз,
С/*(•) = (c/fii1-), ••• ,c/V(-); ... ; c/fn(-), ••• ,c/V (■)), где набор чисел (c/f11(.), ■ ■ ■ , c//RK(■)) повторяется T раз,
wp(0 = kii(-), ... , wpsk(-)),
wb(-) = (wbii(-), ••• , w/rk (■)),
Qp = (Qiii(T), ••• , qpsk (T)),
ep = (Qiii(1),••• ,qpsk(1); •••; QPii(T),••• ,qpsk(T)),
Qb = (Qiii(T), •••, q/rk(t)),
= (QinC1),••• (1); •••; QlnCO,••• (t)).
В векторно-матричной форме системы ограничений, описывающие множества допустимых значений векторов x и у, могут быть записаны в виде M = {x ^ 0 : Ax ^ b} и Q = {у ^ 0 : By ^ d}, где A, B - матрицы и b, d - векторы соответствующих размеров, формируемые из параметров системы ограничений (1)-(14). Нетрудно убедиться в том, что при естественных предположениях, отражающих природу векторов x и у, множества M и Q ограничены, т.е. являются выпуклыми многогранниками.
Пусть далее
- матрица линейного оператора, отображающего вектор x в вектор издержек производителей на передачу электроэнергии между производителями и потребителями, а также между электростанциями и хранилищами электроэнергии в каждый момент времени (число компонент вектора ALx равно
KK ' E=1 E^l i Htm);
Ac - матрица линейного оператора, отображающего вектор x в вектор объемов электроэнергии, которые были произведены электростанциями в каждый момент времени (число компонент вектора Acx равно ISK ■ ^^ 1 Him);
A - матрица линейного оператора, отображающего вектор zp (объемов (в МВт) строящихся новых электростанций или новых блоков действующих электростанций в каждом году) в вектор суммарных мощностей электростанций, которые будут находиться в эксплуатации к началу каждого периода t (число компонент вектора
AV равно ISKT), где t G 1, T;
Ap - матрица линейного оператора, отображающего вектор zp (объемов (в МВт) строящихся новых электростанций или новых блоков действующих электростанций в каждом году) в вектор суммарных мощностей электростанций, которые будут находиться в эксплуатации к окончанию периода T (число компонент вектора Apzp равно ISK);
A - матрица линейного оператора, отображающего вектор (объемов (в МВт) строящихся новых хранилищ электроэнер-
гии или новых блоков существующих хранилищ электроэнергии в каждом году) в вектор суммарных мощностей хранилищ электроэнергии, которые будут находиться в эксплуатации к началу каждого периода t (число компонент вектора A равно IRKT), где t е ITT;
- матрица линейного оператора, отображающего вектор (объемов (в МВт) строящихся новых хранилищ электроэнергии или новых блоков существующих хранилищ электроэнергии в каждом году) в вектор суммарных мощностей хранилищ электроэнергии, которые будут находиться в эксплуатации к окончанию периода T (число компонент вектора Abzb равно IRK);
Yi, Y2, 7з, 74 - матрицы соответствующего размера, ненулевыми элементами которых являются значения параметров Ytm;
£i, ■ ■ ■ , £g - единичные векторы, число компонент которых равно соответственно ISK ■ ^T=1 ^M= 1 Htm, KK ■ Et=1 EMl1 Htm, ISKT, IRKT. ISK, IRK, ISKT, IRKT.
С учетом введенных обозначений рассматриваемая игра может быть сформулирована в векторно-матричной форме следующим образом:
(15)
(xp,Yiyp) + (xb,Y2yb) - (ei,Y3Cost (Acx))-
-<£2,74ALx> -<£3,g (zp)>-
-<£4,/ (zb)> - <£7, c/p (AV + ep)> --<es, c/b (äV + e^ > + <£5, wp (Apzp + Qp)> + + <£6, (Abzb + Qb)> ^ шахуем
xp,Yiyp
+ <хь,72у^ ^ ттуеп .
В этой игре двух лиц все производители электроэнергии в регионе (первый игрок) стремятся максимизировать свою прибыль, выбирая объемы поставляемой электроэнергии в приемлемых для потребителей количествах, в то время как все потребители электроэнергии в регионе (второй игрок) стремятся получить
приемлемые для них объемы электроэнергии с наименьшими затратами.
3. Математические свойства задачи
Нетрудно убедиться в том, что введение новых переменных и ограничений-неравенств с целью перехода от кусочно-линейных функций (аппроксимирующих функции costisk(■),
wt-fc(-Х gisk/¿rfc(-Х c/iSfcOX c/tk())к линейным в
функции выигрыша первого игрока позволяет записать игру (15) в векторно-матричной форме следующим образом:
(16) i(x,Dy> + (v,x> ^ maxxeM
\<x,Dy> ^ minyen, где x = (ж, ж'), ж £ M, а ж' отражает дополнительные переменные, вводимые при аппроксимации выпуклых (или вогнутых) функций costisk(■), wpsk(■), wirk(■), giskО, /irk('), c/isk(), c/irk(') кусочно-линейными; M = |ж ^ 0 : AX ^ - выпуклый многогранник допустимых стратегий группы всех производителей электроэнергии в регионе (одного из игроков в рассматриваемой игре двух лиц); D - матрица коэффициентов линейного оператора, который отображает вектор y в вектор (Yiyp, Y2yb, 0, ■ ■ ■ , 0) того же размера, что и вектор ж; v - вектор параметров, а Q = {y ^ 0 : By ^ d} - выпуклый многогранник допустимых стратегий группы всех потребителей электроэнергии в регионе (второго из двух игроков в рассматриваемой игре двух лиц).
Утверждение 1. Пара векторов (ж*,у*) £ M х Q является равновесием Нэша в игре (16) тогда и только тогда, когда существуют векторы z* ^ 0 и t* ^ 0 такие, что (г*,ж*) и (t*,y*) являются решениями задач линейного программирования
(d, 2) + (v,ic) ^ max _ _ ,
(—b, t) ^ min ,
(t,y)e{(t,y)>0: tA<-v-Dy, By^d}
образующих двойственную пару.
Доказательство. 1. Поисольку из неравенства (xr*,Dy*) + (v, ж*) С + (v, ж*) Vy е Q следует неравенство
(ж*^у*) ^ (ж*^у) Vy е Q, пара векторов (ж*, у*) является точкой равновесия Нэша в игре (16) тогда и только тогда, когда эта точка является седловой точкой в игре двух лиц на выпуклых многогранниках М и Q с платежной функцией
(17) (x,DyC + (v,XC,
которая максимизируется по ж и минимизируется по у.
2. Игра (17) является частным случаем игры с платежной функцией
(18) (р',ж'С + (ж',£'у'С + (q' ,у'С
на выпуклых многогранниках M' = {ж' ^ 0 : А'ж' ^ b'} и Q' = {у' ^ 0 : B'y' ^ d'}, в которой платежная функция максимизируется по у' и минимизируется по ж'. В [1] доказано, что векторные компоненты седловой точки игры (18) находятся из решения задач линейного программирования
(b',2> + (q'^') ^ mt^/
(—d',t) + (р',ж') ^ min , N ' x ' (t,x')eP'
где Q' = {(2, у') ^ 0 : 2A' ^ p + D'y', B'y' ^ d'}, P' = {(t, ж') ^ 0 : tB' ^ —q' - ж'£', А'ж' ^ b'}, ж' е M' = {ж' ^ 0 : А'ж' ^ b'}, у' е Q' = {у' ^ 0 : B'у' ^ d'},
образующих двойственную пару, где М' и О' - выпуклые многогранники. Основываясь на этом результате и полагая, что х = у',
y = ж', v = q', 0 = p', b' = d, d' = b, B' = A, A' = B и D = (D')T
нетрудно показать, что задача отыскания седловой точки в игре (17) сводится к решению пары двойственных задач линейного программирования
Доказанное утверждение позволяет исследовать рынки электроэнергии практически любого размера, поскольку задачи линейного программирования (19) можно решать при любых интересных для практики числах ограничений и переменных.
Рассмотрим рынок электроэнергии в Москве и в Московской области. На этом рынке на четырех основных производителей электроэнергии (Мосэнерго, E-ON, ОГК-1 и РусГидро) приходится 96% производства электроэнергии в регионе (по данным годового отчета Мосэнерго за 2011 год). Указанным четырем компаниям принадлежит 17 электростанций (Мосэнерго принадлежит 15 электростанций, а E-ON и ОГК-1 принадлежит по одной электростанции) и одно хранилище электроэнергии (которое принадлежит РусГидро). В рамках предложенной модели имеет смысл рассматривать каждую из 15 электростанций, принадлежащих Мосэнерго, как отдельного производителя, так как это позволяет учесть технологические ограничения на производство электроэнергии на каждой электростанции. Таким образом, для рассматриваемого региона параметр модели I принимает значение 18, а поскольку все электростанции региона являются конденсационными, параметр модели S принимает значение 1. Далее, на сегодня в рассматриваемом регионе имеется одно крупное хра-
(19)
где, как и ранее, ж е M {y ^ 0 : By ^ d}.
|ж ^ 0 : Ax ^ bj и y е Ü
нилище электроэнергии (Загорская гидроаккумулирующая электростанция), однако в будущем возможно использование хранилищ электроэнергии других типов (например, химических), поэтому на ближайшее будущее можно положить значение параметра К равным 2. На оптовом рынке электроэнергии Москвы и Московской области крупнейшим потребителем электроэнергии, который закупает более 85% электроэнергии в регионе, является Мосэнергосбыт (по данным годового отчета компании за 2013 год), что позволяет считать значение параметра 3 модели равным 1. В рассматриваемой в работе задаче горизонт планирования имеет смысл выбирать так, чтобы дать производителям достаточно времени для построения новых электростанций малой или средней мощности и хранилищ электроэнергии, функционирование которых окупило бы их строительство. Например, можно положить значение параметра Т равным 15 и рассматривать взаимодействие производителей и потребителей электроэнергии в течение четырех «типичных» недель каждого года горизонта планирования.
Разделение сети электроснабжения на узлы (непересекающиеся части рынка электроэнергии региона, внутри которых ограничения на передачу электроэнергии незначительны) является сложной задачей. В настоящее время системный оператор рынка электроэнергии России использует так называемые зоны свободного перетока электроэнергии («...зоны оптового рынка электрической энергии (мощности), внутри которых в стандартных условиях функционирования сети отсутствуют существенные системные ограничения на объемы перетока электроэнергии в течение примерно 30% времени в каждом месяце...» см. Глоссарий системного оператора единой энергетической системы) при выборе им электростанций, которые будут эксплуатироваться в среднесрочной перспективе. Поскольку Москва и Московская область образуют единую зону свободного перетока, в рамках предложенной в работе модели можно положить значение параметра К равным 1. Нетрудно убедиться в том, что для указанных значений параметров модели поиск равновесия в игре, формулиру-
емой с использованием этой модели, сводится к решению пары двойственных задач линейного программирования с числом как прямых, так и двойственных переменных не превышающим трех миллионов.
4. Заключение
В статье рассмотрена проблема оценки будущих цен и объемов электроэнергии, продаваемых производителями электроэнергии её потребителям по прямым двусторонним договорам на оптовом рынке электроэнергии в регионе, и ожидаемых объемов инвестиций производителей как в строительство новых электростанций и новых блоков существующих электростанций, так и в строительство новых хранилищ электроэнергии и новые блоки существующих хранилищ электроэнергии в этом регионе.
Взаимодействие производителей и потребителей электроэнергии в регионе математически описывается бескоалиционной игрой двух лиц, в которой первым игроком является группа всех производителей электроэнергии, снабжающих электроэнергией потребителей региона или имеющих хранилища электроэнергии на территории региона, а вторым игроком является группа всех потребителей электроэнергии в регионе и учитываются технологические ограничения связанные с производством, хранением и передачей электроэнергии. Функции издержек производства, функции стоимости инвестирования в строительство новых электростанций или новых блоков существующих электростанций и в строительство новых хранилищ электроэнергии или новых блоков существующих хранилищ электроэнергии, а также функции оценочной стоимости электростанций и хранилищ электроэнергии, которые (вообще говоря) могут быть произвольными выпуклыми или вогнутыми функциями, с необходимой (для прикладных расчетов с использованием предложенной модели) точностью можно приблизить кусочно-линейными функциями. В построенной игре отыскание равновесия Нэша сводится к решению пары двойственных задач линейного программирования, что позволяет использовать предложенную модель для практических 118
расчетов при решении задач любых интересных для практики размеров.
Необходимо отметить, что в статье не приводится какого-либо численного примера отыскания равновесных стратегий в рассмотренной игре на модельных данных. По мнению авторов такой пример демонстрировал бы не более чем работу стандартных пакетов линейного программирования, что не представляет какого-либо научного и даже прикладного интереса, в то время как основным результатом статьи является продемонстрированная в ней возможность сведения задачи отыскания оптимальных объемов производства и потребления электроэнергии в регионе в рамках прямых двусторонних договоров к задачам линейного программирования, образующим двойственную пару. Ясно также, что решение задач с реальными данными, которые являются конфиденциальной коммерческой информацией энергетических компаний, не представляется возможным без согласия этих компаний предоставить такие данные (с учетом существующей конкуренции на рынке электроэнергии) и потому не являлось целью настоящей статьи.
Литература
1. БЕЛЕНЬКИЙ A.C. Минимаксные задачи планирования с линейными ограничениями и методы их решения // Автоматика и телемеханика. - 1981. - №10. - С. 157-170.
2. ВАСИН A.A., ДАЙЛОВА Е.А. Анализ краткосрочной эффективности механизмов оптового рынка электроэнергии // Журнал Новой экономической ассоциации. - 2013. -№2(18). - С. 35-60.
3. ДАВИДСОН M.P., ДОГАДУШКИНА Ю.В., КРЕЙ-НЕС Е.М. И ДР. Математическая модель управления энергосистемой в условиях конкурентного оптового рынка электроэнергии и мощности в России // Известия РАН. Теория и системы управления. - 2009. - №2. - С. 84-94.
4. ДАВИДСОН M.P., НОВИКОВА Н.М. Задача выбора состава генерирующего оборудования в рынке электроэнер-
119
гии // Тезисы докладов 5-й Московской междунар. конф. по исследованию операций. - М.: ВЦ РАН, 2007. -С. 34-36.
5. СТОФТ С. Экономика энергосистем. Введение в проектирование рынков электроэнергии - М.: Мир, 2006. - 623 с.
6. ЧЕРНИКОВ С.Н. Линейные неравенства. - Итоги науки. Сер. Мат. Алгебра. Топол. Геом. - М.: Наука, 1968. -С. 137-187.
7. BUSHNELL J., ISHII J. An Equilibrium Model of Investment in Restructured Electricity Markets: Working paper, University of California Energy Institute, 2007. - 37 p.
8. DIRKSE S.P., FERRIS M.C. Mcplib: a collection of nonlinear mixed complementarity problems // Optimization Methods and Software. - 1995. - №5(4). - P. 319-345.
9. EHRENMANN A., SMEERS Y. Generation Capacity Expansion in a Risky Environment: A Stochastic Equilibrium Analysis // Operations Research. - 2011. - №59(6). -
P. 1332-1346.
10. FADLULLAH Z.M., NOZAKI Y., TAKEUCHI A., KATO N. A survey of game theoretic approaches in smart grid // Wireless Communications and Signal Processing (WCSP), 2011 International Conference. - 2011. - P. 1-4.
11. GARCIA A., ALZATE J.M., BARRERA J. Regulatory design and incentives for renewable energy // Journal of Regulatory Economics. - 2012. - №41(3). - P. 315-336.
12. GENC T.S., REYNOLDS S.S., SEN S. Dynamic oligopolistic games under uncertainty: A stochastic programming // Journal of Economic Dynamics and Control. - 2007. -№31(1). - P. 55-80.
13. Gurobi Optimizer. - URL: http://www.gurobi.com/ products/gurobi-optimizer (дата обращения: 10.04.2015).
14. LISE W., KRUSEMAN G. Long-term price and environmental effects in a liberalised electricity market // Energy Economics. - 2008. - №30(2). - P. 230-248.
15. LUEKEN R., APT J. The effects of bulk electricity storage
on the PJM market // Energy Systems. - 2014. - №5(4). -P. 677-704.
16. NAGL S. Prices vs. Quantities: Incentives for Renewable Power Generation - Numerical Analysis for the European Power Market: Working paper, Energiewirtschaftliches Institut an der Universitaet zu Koeln. - 2013. - 41 p.
17. PINEAU P.-O., MURTO R. An Oligopolistic Investment Model of the Finnish Electricity Market // Annals of Operations Research. - 2003. - №121(1-4). - P. 123-148.
18. PINEAU P.-O., RASATAB H., ZACCOUR G. Impact of some parameters on investments in oligopolistic electricity markets // European Journal of Operational Research. -2011. -№213(1). - P. 180-195.
19. SCHILL W.-P., KEMFERT C. Modeling Strategic Electricity Storage: The Case of Pumped Hydro Storage in Germany // The Energy Journal. - 2011. - №3. - P. 59-88.
20. SCHRODER A. An electricity market model with generation capacity expansion under uncertainty // Energy Systems. -2014. - №5(2). - P. 253-267.
21. TRABER T., KEMFERT C. Gone with the wind? - Electricity market prices and incentives to invest in thermal power plants under increasing wind energy supply // Energy Economics -2011. - №33(2). - P. 249-256.
22. VASIN A., KARTUNOVA P., WEBER G.-W. Models for capacity and electricity market design // Central European Journal of Operations Research. - 2013. - №21(3). -
P. 651-661.
GAME-THEORETIC APPROACH TO ANALYSIS OF FUTURE REGIONAL POWER SUPPLY SYSTEM
Alexander S. Belenky, National Research University Higher School of Economics, Moscow, Russia, Doctor of Science, professor ([email protected]).
Dmitry Bolkunov, National Research University Higher School of Economics, Moscow, Russia, MSc in Applied Mathematics and Computational Science ([email protected]).
Abstract: We design the elements of an automated tool for the quantitative analysis of the future regional power supply system based on the opinion of regional authorities on the current and future availability of power supply services in the industry and in the domestic sector. A mathematical model is suggested to estimate the return ofinvestment into new power station and storage construction. The model includes a non-cooperative game on polyhedra, where the first player's payoff function adds up from the bilinear and the linear function of vector arguments, whereas the the payoff function of the second player is a bilinear function of the same arguments. We prove that equilibria in this game can be found by solving a dual pair of linear programming problems. These equilibria determine, in particular, the volume ofinvestment for generators and the electricity price for consumers, which are mutually acceptable to both parties.
Keywords: convex polyhedra, control in large systems, linear programming, Nash equilibrium, saddle points, two-person game.
Статья представлена к публикации членом редакционной коллегии В.Н. Бурковым
Поступила в редакцию 03.02.2015. Дата опубликования 30.09.2015.