Разработка математической модели общей стоимости комбинированной установки на основе возобновляемых источников энергии Текст научной статьи по специальности «Энергетика и рациональное природопользование»

Разработка математической модели общей стоимости комбинированной установки на основе возобновляемых источников энергии

Д.В. Гринько, аспирант, Оренбургский ГАУ

Высокая стоимость оборудования является основным препятствием использования возобновляемых источников энергии (ВИЭ), что обусловливает необходимость поиска оптимальной конфигурации системы с максимальной производительностью при минимальных затратах. Нами была поставлена цель — разработать новый подход, позволяющий обеспечить высокое качество решения задачи нахождения общей стоимости комбинированных установок малой и средней мощности на основе ВИЭ, с последующей оптимизацией за счёт многовариантного анализа полученных результатов.

Общая годовая стоимость компонентов комбинированной установки на основе ВИЭ рассчитывается следующим образом [1]:

С

= 1 (С

+ с + С + С

acap,c arep,c aop,c -

emission

), (1)

c=1

где СасЩс1, — годовая стоимость капитальных вложений;

СаГеР,с — годовая стоимость замены компонентов; Саор,с — годовая стоимость эксплуатации.

Из уравнения 1, используя аддитивный закон суммирования затрат, автором выведена математическая модель расчёта общей годовой стоимости комбинированной установки на основе ВИЭ, состоящей из ветротурбины, солнечной батареи, дизельного генератора, аккумуляторной батареи, централизованной электрической сети, следующим образом:

Ccom I (С1,1

+ С2, w + С3, w ) +

Ns

IС i

s=1

w=1

+ С 2,s +С3,s

Nb

Ng

)+IС

+ С2, g + С3, g +

+C4,g ) + I (Clb

+ C2,b + C3,b ) +

(2)

b=1

+ С + Cst) + Cd + Cn + (Q, gr + С2, gr + С3, gr +

+C4,gr ) - С

5,gr>

C3, s — стоимость эксплуатации солнечных батарей ($/год);

Ng — количество дизельных генераторов (шт.); C1, — вложения в покупку дизельного генератора ($/год);

C2 — стоимость замены дизельного генератора ($/год);

C3, — стоимость эксплуатации дизельного генератора ($/год);

C4, — стоимость топлива для дизельного генератора ($/год);

Nb — количество аккумуляторных батарей (шт.); C1, b — вложения в покупку аккумуляторных батарей ($/год);

C2, — стоимость замены аккумуляторных батарей ($/год);

C3, — стоимость эксплуатации аккумуляторных батарей ($/год);

Cf — стоимость доставки (завоза) топлива ($/год); Cst — стоимость хранения топлива (в нашем случае коэффициент обнуляется из-за отсутствия необходимости хранения топлива) ($/год); Cd — стоимость доставки оборудования ($/год); Cn — суммарный налог на выбросы вредных веществ (плата за негативное воздействие на окружающую среду (НВОС) для Оренбургской области) ($/год);

C1, gr — капитальные вложения в строительство линии электропередач (ЛЭП) за 1 км. ($/км/ год);

C2, gr — стоимость эксплуатации ЛЭП ($/км/ год);

C3, gr — стоимость подключения к сети в случае наличия ЛЭП ($/год);

C4, gr — стоимость потребляемой энергии от ЛЭ П ($/кВт- ч/год);

C5, gr — стоимость продажи электроэнергии, выработанной от ВИЭ в сеть (в случае принятия в РФ соответствующих законов), ($/кВт ■ ч/год). Учёт коэффициента Cf особенно актуален при завозе топлива в отдалённые районы Севера, что значительно влияет на общую стоимость комбинированных установок, он рассчитывается в виде:

где N — количество ветротурбин (шт.);

С1, № — вложения в покупку ветротурбин ($/год); С2, » — стоимость замены ветротурбин ($/год); С3 № — стоимость эксплуатации ветротурбин

($/год);

N — количество солнечных батарей (шт.); С1 . — вложения в покупку солнечных батарей ($/год);

С2 . — стоимость замены солнечных батарей ($/год);

С/ = СР V,

(3)

где Ср — цена перевозки одной единицы объёма от поставщика к потребителю ($); V — объём перевозки (куб. ед.).

Коэффициент Сл применяется при доставке оборудования в отдалённые районы либо при доставке отсутствующего оборудования из других стран или регионов. В данном случае необходима доставка ветротурбины из Китая, всё остальное оборудование имеется в локальной доступности

и данный коэффициент по отношению к нему не учитывается. Коэффициент С рассчитывается следующим образом:

С, = (Dc + Dr) • M,

(4)

где Dc — цена доставки одной единицы массы по Китаю ($);

Dr — цена доставки одной единицы массы по России ($);

M — масса перевозки (ед.).

Недостатком централизованного электроснабжения с экологической точки зрения является наличие вредных выбросов при сжигании условного топлива в месте генерации электроэнергии. Поскольку одной из целей статьи является моделирование децентрализованной, экологически чистой комбинированной установки с практически 100-процентной выработкой энергии от ВИЭ, интеграция с сетью не рассматривается и все коэффициенты, связанные с ней, обнуляются.

Группа учёных под руководством Д. С. Стреб-кова проанализировала технико-экономические параметры комбинированной солнечно-ветробен-зиновой установки с использованием программного обеспечения (ПО) «HOMER» для условий Египта. Результаты показали, что себестоимость 1 кВтч электроэнергии, выработанной комбинированной установкой, равна $0,189, выработка электроэнергии от ВИЭ составляла 18 и 81% от бензогенератора [2].

Нами смоделирована комбинированная система электроснабжения части жилого комплекса «Экодолье», расположенного в районе с. Ивановка Оренбургской обл. Концепция создания посёлка, базирующаяся на принципах экономичности,

экологичности и энергоэффективности, делает его идеально подходящим для дальнейшего исследования.

Электрическая нагрузка оказывает самое большое влияние на стоимость и размер узлов установки. Рассматриваемая часть жилого комплекса состоит из 10 домов, предполагаемое потребление электроэнергии приведено в таблице 1.

На основе определённых почасовых значений энергопотребления за 24 часа программное обеспечение «HOMER 2», версия 2.81 [3], сгенерировало суточный профиль нагрузки, представленный на рисунке 1. Симулируя изменчивость нагрузки, чтобы график для каждого дня был уникален, добавлена хаотичность в виде 10% ежедневной и почасовой изменчивости нагрузки.

Определена средняя допускающая задержку нагрузка (требование, которое может быть удовлетворено в определённый период времени, при этом точная синхронизация не важна и существует гибкий график, типичный пример — водяной насос при наполненном резервуаре воды и заряженные аккумуляторные батареи). Как показано на рисунке 2, предполагается, что объём энергопотребления возрастёт с июня по август из-за полива в летний период.

Характер ВИЭ воздействует на режим и экономику энергосистем, так как ресурсы определяют объёмы и синхронизацию выработки энергии. Поэтому детальное и точное моделирование возобновляемых ресурсов является важным элементом проектирования устройства.

Солнечный ресурс. ПО «HOMER 2» генерирует синтезированные почасовые данные солнечного излучения на основе входных данных о ежемесячных

1. Предполагаемое среднесуточное потребление электроэнергии

Наименование Мощность, Вт Время работы (усреднённое), час. Энергия (из расчёта на 10 домов), Вт • ч/сут

Энергосберегающая лампа 10 Вт, 10 шт. 100 5 5000

Телевизор 100 1 1000

Компьютер 400 1 4000

Глубинный электронасос 1500 2 30000

Холодильник 150 24 36000

Итого 2250 - 76000

Рис. 1 - Среднесуточный профиль нагрузки, сгенерированный «HOMER 2»

Рис. 2 - Профиль среднегодовой допускающей задержку нагрузки, сгенерированный программным обеспечением «HOMER 2»

солнечных ресурсах с использованием алгоритма Грэхема [4]. Суммы ежемесячной солнечной радиации на горизонтальную поверхность (кВт-ч/ м2/сут) внесены из базы данных NASA SSE [5]. Индекс яркости определён ПО «HOMER 2» на основе указанного излучения, часового пояса (+5 GMT), координат (51.71 с.ш.; 55,21 в.д.), результаты приведены в таблице 2. Выходные данные — на основе значений ПО «HOMER 2» строит 8760-часовую базу данных солнечных излучений со статистическими характеристиками, подобными результатам реальных измерений с учётом изменчивости и автокорреляции.

Ветровой ресурс. Среднемесячная скорость ветра на высоте 10 м была взята из характеристики природных условий района изысканий [6]. Если измерение скорости ветра не было произведено на высоте втулки воздушной турбины (в нашем случае 25 м), необходимо установить зависимость скорости ветра от высоты втулки. Это может быть сделано с использованием логарифмического (Eurocode 1, EN 1991-1-4, ДБН В.1.2-2 2006 Нагрузки и воздействия) или степенного закона (СНиП 2.01.07-85, СП 20.13330.2011 Нагрузки и воздействия):

ЧZ) = Va

f \a z

z

V анем у

(5)

где vaHeM — скорость ветра на уровне расположения анемометра (zaHJ;

а — показатель степени зависит от шероховатости подстилающей поверхности и от самой скорости ветра.

В СНиП 2.01.07-85 используется следующая классификация поверхностей: А — открытые побережья морей, озёр и водохранилищ, пустыни, степи, лесостепи, тундра; В — городские территории, лесные массивы и др. местности, равномерно покрытые препятствиями высотой до 10 м; С — городские районы с застройкой зданиями высотой свыше 25 м. Показатели степени для указанных выше типов местности равны: аА= 0,16; аВ= 0,22 и аС= 0,33. С использованием уравнения 5 определена скорость ветра на высоте 25 м, результаты представлены в таблице 3.

Для обеспечения электроснабжения согласно расчётным нагрузкам, охарактеризованным выше, предлагается использовать комбинированную электростанцию, которая состоит из компонентов, представленных в таблице 4.

Схема моделируемого комбинированного устройства, сгенерированная при помощи программного обеспечения «HOMER 2» (v2.81), показана на рисунке 3. Результаты моделирования комбинаций узлов устройства представлены в таблице 5.

Оценка исследуемой области показала значительный потенциал ВИЭ. Поскольку для 3-й категории надёжности допустимое число часов отключения в год составляет 72 час. [7], а в си-

2. Среднегодовое солнечное излучение

Показатель Месяц

Янв. Февр. Март Апр. Май Июнь Июль Авг. Сент. Окт. Нояб. Дек. Ср.- год.

Индекс яркости Излучение, кВт • ч/м2/сут 0,448 0,990 0,514 1,880 0,569 3,390 0,545 4,660 0,565 5,990 0,559 6,440 0,563 6,220 0,559 5,200 0,555 3,790 0,494 2,140 0,467 1,180 0,443 0,800 0,545 3,565

3. Среднемесячная скорость ветра, м/сек

Месяц

Показатель Янв. Февр. Март Апр. Май Июнь Июль Авг. Сент. Окт. Нояб. Дек. Ср.- год.

10 м 4,5 4,6 4,2 4,3 4,3 3,7 3,4 3,4 3,7 4,0 4,0 4,1 4,0

25 м 5,2 5,3 4,9 5,0 5,0 4,2 4,0 4,0 4,2 4,6 4,6 4,7 4,6

4. Компоненты комбинированной электростанции

Наименование Модель Мощность, кВт Цена, $ Комментарий

Ветротурбина Windelectric-Europe 5000 5 кВт 4765 Генерирует переменный ток, стартовая скорость ветра 2,0 м/сек, мачта 25 м

Солнечная фотоэлектрическая батарея На основе аморфного кремния 5-22 кВт 4600 за 1 кВт С учётом тенденции снижения цены стоимость замены берётся в количестве 80% от начальных капитальных затрат

Инвертор HYUNDAI 15 кВт 1600 Эффективность инвертора принимается равной 90%

HYUNDAI 22 кВт 1800

Группа аккумуляторных батарей SURRETTE 6CS25P 6V, 1156 Ah, 6,94 kWh 950 Выбор основан на ключевых материальных свойствах -номинальное напряжение, ёмкость, срок службы, минимальный заряд и эффективность зарядки-разрядки

Дизельный генератор «Азимут» АД 15-Т400 15 кВт 5250 Для исследования рассматривают дизельное топливо с ценами: 0.8, 0.9 и 1 доллар США за литр. Текущая цена составляет 0.9 долл. за л

5. Возможные комбинации узлов устройства

%

m

И В

ак т о

б

а р

ы В

на

и

б р

урт

о рот

е В

«

е аре

ата б

к

S

чне

н л о С

р

ротра р

е н е г

«

ь л е

изе «

fr & J &

< ю

р

рот р

И

ь л а чаа На

у

л с б

,о

т S

£

«

,о се К

" S ^

S *

а

з з « е о

2 И а ч

i

о

s

о

л

m я s

об та

Й V H Л

св ,и

§

s Sa

О & с^

и к

в £ ^

с,

0 g ди

й ® § §

1 £ ир оо тб Со

оп топ

О g

а

а

2 ^ ыр Н о ее

аб

! s

s «

ов та

H чс

уо д

О X

я

s

К <4 Щ

о

,т о X

H

н Оа Вов

« g

н

а,

о U б

& I г

д 2 н

в лВ

§ «1 й «

ё ft и К m О

л п опт

изт

д д

I s

CS

СМ ^

о н ь

S «

изего./ д .с

тыча

I â

« ей S ft

I s

"2

В

52 9Q 1QQ

5 2x5 2x5

5

1Q 22

15 15 15

7

3Q 43

15 15 22

41,265 9Q,88Q 158,63Q

144,2Q7/Q.4Q8 146,142/Q.413 212,Q38/Q.599

515 1Q3Q 1Q3Q

9Q,4 2Q,6 Q,9

174,432 153,937 213,365

Q,469 Q,424 Q,6Q1

6779 1547

69,5

2142 451 16

6. Количество вредных выбросов и стоимость НВОС

Выбросы вредных веществ (кг/год)

Выработка НВОС

от ВИЭ, % двуокись окись несгоревшие дисперсные диоксид оксид ($/год)

углерода углерода углеводороды частицы серы азота

52 17,851 44,1 4,88 3,32 35,8 393 1Q98

9Q 4,Q73 1Q,1 1,11 Q,758 8,18 89,7 25Q

1QQ 183 Q,452 Q,Q5 Q,Q34 Q,367 4,Q3 11

АС ОС

Рис. 3 - Схема моделируемого устройства

стеме с практически 100-процентной (выработка от дизельного генератора 0,153%) выработкой от ВИЭ резервный дизельный генератор используется 16 час. в год, возможно его полное исключение из системы, что удешевит стоимость установки. Количество вредных выбросов и стоимость НВОС, рассчитанная от объёма вредных выбросов, приведены в таблице 6.

Таким образом, на основе разработанной математической модели представлен расчёт общей стоимости комбинированной установки. Достовер-

ность данного метода подтверждается совпадением результатов расчёта с ПО «HOMER 2», решающим аналогичные задачи. Математическая модель отличается близостью к реальным условиям (адекватностью) и простотой, что обеспечивает возможность и удобство использования модели по её прямому назначению. Модель, в отличие от уравнения 1 и аппарата ПО «HOMER 2», является более точной за счёт учёта ряда новых коэффициентов суммарной стоимости. Данная математическая модель может быть взята за основу написания ПО для расчёта стоимости комбинированных электроустановок.

Литература

1. M. Bailey, M. Eastwood, T Grieser, L. Borovick, V. Turner, and R.C. Gray, «Special Study: Data Center of the Future», New York, NY: IDC 06C4799, April 2007.

2. Стребков Д.С., Сокольский А.К., Ахмед Джайлани Т.А. Комбинированные электростанции для автономных сельскохозяйственных потребителей в Египте // Вестник ФГОУ ВПО «Московский государственный агроинженерный университет им. В.П. Горячкина». 2009. № 3. С. 42—44.

3. HOMER 2. URL: http://homerenergy.com/software.html. (Дата обращения 03.06.2012).

4. HOMER, the micropower optimization model. URL: http:// www.nrel.gov/homer, свободный. (Дата обращения 03.06.2012).

5. The NASA Surface Meteorology and Solar Energy Data Set // 2009. URL: http:// eosweb.larc.nasa.gov/sse/.

6. Климатическая характеристика. URL:http://www.ecoindustry. ru/phorum/viewtopic.html?geo=56&f=14&t=2402, свободный. (Дата обращения 01.10.2013).

7. Об утверждении Правил функционирования розничных рынков электрической энергии в переходный период реформирования электроэнергетики. Постановление Правительства РФ от 31 августа 2006 г. № 530. П. 113 // СЗ РФ. 2006. № 37. Ст. 3876.