Россия
Рассмотрены ветро-солнечные энергоустановки для автономного электроснабжения небольших нефтегазопромысловых объектов мощностью не более 20 кВт, ветроэнергоустановок в их составе – 10 кВт с коэффициентом мощности (энергетическим КПД) не более 0,35–0,45. Разработана усовершенствованная конструкция ортогональной ветроэнергоустановки на основе комбинированного ротора Н-Дарье-Савониуса. Эксперименты на лабораторной установке показали, что максимальные значения коэффициента мощности достигали 0,60, что гораздо выше, чем у роторов известных видов ветроэнергоустановок. Рассчитаны значения коэффициента мощности КРДС при различных геометрических размерах входящих в его состав роторов Н-Дарье и Савониуса. С увеличением скорости ветра до номинального для ветроэнергоустановок малой и средней мощности значения 9,0 м/c и диаметре ротора 15,0 м мощность комбинированного ротора достигает практически 60 кВт. Предлагаемые ортогональные ветроэнергоустановки с номинальной мощностью до 60 кВт смогут совместно с другими энергоустановками применяться для автономного энергоснабжения на нефтегазопромысловых объектах на суше и на море, а в дальнейшем и на интеллектуальных месторождениях.
ветро-солнечные энергоустановки, аксиальные ветроэнергоустановки, ортогональные ветро-энергоустановки, комбинированный ротор Н-Дарье-Савониуса, коэффициент мощности, энергоэффективность
Введение
Для электроснабжения нефтегазодобывающих промыслов, в частности в ООО «ЛУКОЙЛ», ПАО «Газпром нефть» и других инфраструктурах отдаленных месторождений на суше и на море, чаще всего применяются электростанции, работающие на природном и попутном газе. Однако стали использоваться и экологически чистые возобновляемые источники энергии (ВИЭ), такие как солнечная, ветровая и др., применяться комбинированные ветро-солнечные электростанции (КВСЭС), т. к. они более надежны и стабильны, чем ветровые или солнечные электростанции в отдельности [1, 2]. Успешно прошли испытания автономные блочно-комплектные энергетические установки серии БКЭУ-ВСМ на базе ветро-солнечных модулей (БКЭУ-ВСМ-ДГА) [3, 4]. Все эти ветро-солнечные энергоустановки используются в основном для автономного электроснабжения небольших линейных объектов малой мощности: площадки электроприводных узлов нефте- и газопроводов, станций электрохимической защиты от коррозии, систем обнаружения утечек и др. Как правило, мощность этих установок составляет не более 20 кВт, а мощность ветроэнергоустановок (ВЭУ) в их составе – 10 кВт, причем их коэффициент мощности (энергетический КПД) не превышает 0,45.
Однако на практике уже в настоящее время возникает потребность в применении более мощных и энергетически эффективных ВЭУ, которые совместно с другими энергоустановками смогли бы использоваться не только в инфраструктуре на небольших линейных объектах, но и для привода скважинных штанговых, центробежных, винтовых, винтовых штанговых и других установок для добычи нефти. В частности, ими могут быть оборудованы полностью автоматизированные, так называемые интеллектуальные скважины [5, 6]. Поэтому достаточно актуальной представляется разработка более мощных и энергетически эффективных ортогональных ветроэнергоустановок (ОВЭУ) и создание на их основе энергокомплексов с использованием ветро-солнечных энергоустановок для нефтегазопромысловых объектов на суше и на море [7, 8].
Целью работы является разработка и оценка параметров модернизированных ОВЭУ для автономного электроснабжения различных нефтегазопромысловых объектов на суше и на море, анализ существующих ВЭУ в блочных ветро-солнечных модулях, экспериментальное исследование основных параметров модернизированного ортогонального ротора и оценка параметров ОВЭУ для добычи нефти.
Краткая характеристика ВЭУ
В настоящее время наиболее широко применяются горизонтально-осевые или аксиальные ветроэнергоустаноки (АВЭУ) на основе роторов пропеллерного типа [7, 8]. Такие АВЭУ стали использоваться и в БКЭУ-ВСМ-ДГА производства ОАО «НИПОМ» в шкафном исполнении (рис. 1) [3, 4].
Рис. 1. БКЭУ-ВСМ-ДГА в шкафном исполнении
Fig. 1. BKEU-VSM-DGA in cabinet design
В последние десятилетия наряду с АВЭУ стали активно развиваться и вертикально-осевые или ОВЭУ [7–10]. Несомненным преимуществом всех многочисленных видов ОВЭУ является отсутствие необходимости использовать в их конструкции направляющие вдоль направления ветра устройства, т. к. работа этих энергоустановок не зависит от направления ветра. ОВЭУ имеются 2 основных типа ветроколеса (ротора): быстроходных типа ротора Дарье и тихоходных типа ротора Савониуса. Быстроходный ротор Дарье с криволинейными или вертикальными лопастями крылового профиля, так называемый ротор Дарье, вращается в основном за счет подъемных аэродинамических сил. В качестве примера на рис. 2 показана ОВЭУ на основе ротора Н-Дарье с вертикальными лопастями крылового профиля в составе БКЭУ-ВСМ-ДГА в блок-контейнере производства ОАО «НИПОМ» [3, 4].
Рис. 2. БКЭУ-ВСМ-ДГА в блок-контейнере
Fig. 2. BKEU-VSM-DGA in a block container
Тихоходные роторы Савониуса с лопастями полуцилиндрической формы работают за счет разности сил лобового сопротивления ветровому потоку воздуха. Преимуществами ОВЭУ с роторами Дарье различных конструкций являются их быстроходность с коэффициентом быстроходности Z = 4–5, показывающим соотношение между окружной скоростью ротора и скоростью ветра, а также достаточно высокий коэффициент мощности или энергетический КПД, показывающий соотношение между полезной мощностью ротора и мощностью ветрового потока, который достигает значения η = 0,45. Недостатком ОВЭУ с роторами Дарье является невозможность самозапуска при малых скоростях ветра, что является преимуществом ротора Савониуса. Недостатками ротора Савониуса являются низкий энергетический КПД, не превышающий величины η = 0,20, и низкий коэффициент быстроходности, не превышающий Z < 1.
Наряду с этими группами ОВЭУ, выпускаемых серийно, имеются комбинированные ортогональные ветроэнергоустановки (КОВЭУ), в которых сочетаются два описанных принципа действия. КОВЭУ типа Маглев с классическим ротором Дарье с криволинейными лопастями и ротора Савониуса гораздо меньшего диаметра с номинальной мощностью от 300 Вт до 3,0 кВт при номинальной скорости ветра v = 9,0 м/с уже выпускается серийно в ряде стран, в т. ч. и в России, и достаточно широко используется для энергоснабжения небольших объектов [11]. КОВЭУ типа Маглев в составе автономного энергокомплекса в сочетании с солнечный энергостанцией (СЭС) показана на рис. 3.
Рис. 3. КОВЭУ типа Маглев в сочетании с СЭС
Fig. 3. Maglev-type COVEU in combination with SPS
Однако эти установки имеют коэффициент мощности СР не более 0,30–0,40, т. е. они также недостаточно эффективны. Поэтому актуальным представляется разработка более энергетически эффективных ОВЭУ, которые могли бы быть использованы для автономного электроснабжения нефтегазопромысловых объектов как на суше, так и на море. Таким образом, целесообразно в нефтегазопромысловых автономных энергокомплексах использовать разработанные в Лаборатории нетрадиционной энергетики Астраханского государственного технического университета (АГТУ) ОВЭУ на основе модернизированных комбинированных роторов Н-Дарье-Савониуса (КРДС).
Экспериментальные исследования
В Лаборатории нетрадиционной энергетики АГТУ была разработана лабораторная установка для экспериментального исследования эффективности модернизированных КРДС, имеющих лопасти с зигзагообразными закрылками [12, 13]. Общий вид этой лабораторной экспериментальной установки показан на рис. 4.
|
|
|
|
а
|
б
|
Рис. 4. Лабораторная экспериментальная установка для исследования модернизированных КРДС
Fig. 4. Laboratory experimental setup for research of modernized CRDS
Как видно из рис. 4, ветровой поток в лабораторной экспериментальной установке имитировался вентилятором с диаметром рабочего колеса 0,45 м, который был заключен в квадратный кожух с предохранительной сеткой. Для определения скорости воздушного потока в 18 точках в каждом из сечений использовались цифровой анемометр с погрешностью не более 0,1 м/с и координатное устройство в виде рамки с сеткой из синтетических нитей на 3 × 6 ячеек размерами 0,10 × 0,10 м. На расстоянии х от вентилятора определялась средняя скорость воздушного потока V в каждом из сечений на основе осреднения скоростей потока в 18 ячейках. За счет изменения расстояния от вентилятора до оси КРДС осуществлялось изменение средней по сечению скорости воздушного потока. Вращающий момент КРДС уравновешивался моментом силы упругости пружины динамометра. Левая неподвижная часть была присоединена к синтетической нити, которая проходила по тормозному шкиву и крепилась к неподвижному стержню. Деформация предварительно протарированной пружины измерялась по соответствующей шкале с погрешностью не более 1 мм. КРДС монтировался на едином валу с тормозным шкивом и шкивом мультипликатора. При помощи электрогенератора постоянного тока с клиноременным мультипликатором для повышения частоты вращения вала электрогенератора и цифрового вольтметра определялась частота вращения КРДС. Трехлопастной ротор Н-Дарье диаметром DD = 0,320 м имеет лопасти высотой НD = 0,300 м, шириной bD = 0,100 м и толщиной δD = 0,012 м. Были использованы 2 двухступенчатых ротора Савониуса оптимальной геометрии различных размеров: диаметром DS = 0,185 м и DS = 0,0925 м, четырьмя лопастями с углом раскрытия θ = 135º (см. рис. 4, а), а также диаметром DPS = 0,100 м и DPS = 0,050 м, высотой НPS = 0,140 м и НPS = 0,070 м (см. рис 4, б). Все лопасти ротора Н-Дарье и Савониуса имели зигзагообразные закрылки с треугольными элементами.
По показаниям цифрового вольтметра, измеряющего электродвижущую силу электрогенератора E, об/мин, определялась частота вращения КРДС:
(1)
где Kn – тарировочный коэффициент, об/мин·В.
Полезная мощность экспериментального КРДС определялась по формуле:
(2)
где Kf – тарировочный коэффициент динамометра, н/м; Δx – деформация пружины, м; d – диаметр шкива для измерения силы упругости, м.
Результаты определения основных параметров КРДС показаны на рис. 5 [12, 13].
а б
Рис. 5. Основные параметры КРДС:
а – зависимость скорости воздушного потока от расстояния от корпуса вентилятора до оси КРДС;
б – зависимость мощности КРДС от скорости воздушного потока
Fig. 5. Main parameters of the CRDS:
а – dependence of the air flow speed on the distance from the fan housing to the axis of the fan;
б – dependence of the CRDS power on the air flow speed
Как видно из рис. 5, а, с увеличением расстояния от 0,45 до 3,0 м средняя скорость воздушного потока уменьшается с 7,2 до 2,7 м/с, т. е. в 2,7 раза, что приводит к изменению мощности воздушного потока, пропорциональной скорости в третьей степени в 19 раз. Такой диапазон скоростей воздушного потока практически соответствует реальному диапазону изменения скорости ветра, при которой могут эксплуатироваться КРДС в ОВЭУ малой и средней мощности. Зависимость мощности экспериментального КРДС от скорости воздушного потока приведена на рис. 5, б. С повышением скорости воздушного потока с 2,7 до 7,2 м/c мощность ротора увеличивалась с 0,60 до 9,91 Вт, т. е. в 17 раз. Несколько меньшая степень повышения полезной мощности КРДС с повышением скорости воздушного потока, по сравнению со степенью повышения затраченной мощности воздушного потока от скорости, объясняется снижением коэффициента мощности (энергетического КПД) при более высоких скоростях воздушного потока. Было установлено, что при малых скоростях воздушного потока ротор Савониуса способствует первоначальному страгиванию и раскрутке комбинированного ротора, а при больших скоростях начинает тормозить ротор Н-Дарье. Вопрос об оптимизации соотношения геометрических параметров ротора Н-Дарье и Савониуса в составе комбинированного ротора требует дальнейших исследований. Таким образом, КРДС обладает быстроходностью ротора Н-Дарье и возможностью самозапуска и устойчивой работой из-за ротора Савониуса даже при малых скоростях ветра 2–3 м/с, что обеспечит повышение выработки электроэнергии на 5–15 % при существующих для каждой местности времени работы ВЭУ при каждой градации скоростей.
Были также экспериментально определены коэффициенты мощности (энергетические КПД) нескольких роторов для ОВЭУ, имеющих лопасти с закрылками и коэффициенты быстроходности Z [12, 13]. На рис. 6 показаны зависимости коэффициентов мощности этих роторов от коэффициента быстроходности Z.
Рис. 6. Зависимости энергетического КПД ротора Н-Дарье (1), КРДС (2, 3)
и ротора Савониуса (4) от коэффициента быстроходности Z
Fig. 6. Dependences of the energy efficiency of the H-Darier rotor (1), KRDS (2, 3)
and Savonius rotor (4) on the speed coefficient Z
Как видно из рис. 6, все кривые зависимости коэффициентов мощности роторов от коэффициентов быстроходности Z имеют ярко выраженный максимум при различных их значениях, зависящих от конструкции ротора. Полученные максимальные значения коэффициента СР.
КРДС (кривая 2) составляет: СРDS = 0,60 при меньших размерах ротора Савониуса DS = 0,0925 м, DPS = 0,050 м, НS = 0,070 м и СРDS = 0,40 (кривая 3) при больших размерах ротора Савониуса DS = 0,185 м, DPS = 0,100 м, НS = 0,140 м. При малых скоростях воздушного потока ротор Савониуса способствует страгиванию и раскрутке комбинированного ротора, а при больших скоростях ветра начинает тормозить ротор Н-Дарье. Это приводит к значительно большим значениям коэффициента мощности КРДС с меньшими размерами входящего в его состав ротора Савониуса.
Для сравнения показаны также зависимости коэффициентов мощности ротора Н-Дарье (кривая 1) и ротора Савониуса (кривая 4), имеющие лопасти
с закрылками от коэффициента быстроходности Z. Максимальное значение коэффициентов мощности составляет для ротора Н-Дарье СР = 0,72 и ротора Савониуса СР = 0,32, что больше, чем полученного для роторов известных типов [7–10, 14, 15]. Максимумы кривых 1–4 соединены кривой 5 и зависят от коэффициентов быстроходности Z, изменяющихся от 0,71 для ротора Савониуса до 1,23 для ротора Н-Дарье. Следует отметить, что при уменьшении относительной ширины лопастей ротора Н-Дарье коэффициент быстроходности может быть увеличен до значения Z = 4–5 [7, 8, 14].
Таким образом, разработанный авторами КРДС обладает быстроходностью, характерной для ротора Н-Дарье, и возможностью самозапуска из-за большого пускового момента ротора Савониуса даже при небольших скоростях ветра 2–3 м/с.
В отличие от АВЭУ, работа ОВЭУ на основе КРДС не зависит от направления ветра, и поэтому удельная выработка энергии аналогичных по мощности конструкций у этих ВЭУ будет больше. Энергетический КПД КРДС может достигать 0,6, т. е. в 1,3–1,7 раза выше, чем у известных видов роторов ВЭУ.
Оценка параметров ОВЭУ
Величина полезной мощности КРДС для ОВЭУ при диаметре ротора, равной его высоте (D = H), определяется по формуле:
(3)
где ρ – плотность воздуха.
Коэффициент мощности КРДС при вынесении ротора Савониуса за пределы ротора Н-Дарье может быть определен по ранее полученной в [13] формуле:
(4)
где η1 – коэффициент, учитывающий снижение эффективности комбинированного ротора по сравнению с ротором Н-Дарье; CPD, CPS – коэффициенты мощности роторов Н-Дарье и Савониуса; ND, NS – количество ступеней роторов Н-Дарье и Савониуса; DD, DS – диаметры роторов Н-Дарье и Савониуса; HD, HS – высоты роторов Н-Дарье и Савониуса.
Эмпирический коэффициент η1 учитывает снижение эффективности КРДС по сравнению с ротором Н-Дарье за счет того, что максимальная эффективность роторов Савониуса и Н-Дарье достигается при разных значениях коэффициента быстроходности Z, среднее значение этого коэффициента на основе экспериментальных данных составило η1 = 0,78.
Результаты расчета мощности N по формуле (3) с использованием результатов расчета коэффициента мощности по формуле (4) в виде зависимости от скорости ветра V и диаметра КРДС DDS показаны в табл.
Результаты расчета мощности N в зависимости от скорости ветра V и диаметра ротора КРДС DPD
Results of calculation of power N depending on wind speed V and diameter of CRDS rotor DPD
|
DDS, м |
V, м/с |
||||||
|
3,0 |
4,0 |
5,0 |
6,0 |
7,0 |
8,0 |
9,0 |
|
|
5,0 |
0,245 |
0,581 |
1,13 |
1,96 |
3,11 |
4,65 |
6,62 |
|
10,0 |
0,980 |
2,33 |
4,54 |
7,84 |
12,5 |
18,6 |
26,5 |
|
15,0 |
2,20 |
5,24 |
10,2 |
17,64 |
28,1 |
41,9 |
59,6 |
Анализ таблицы показывает, что с увеличением скорости ветра от 3,0 до 9,0 м/c и диаметра ротора от 5,0 до 15,0 м мощность КРДС увеличивается
с 0,245 до 59,6 кВт, т. е. практически пропорционально квадрату диаметра КРДС DDS, а фактически ротора Н-Дарье DD и практически пропорционально кубу скорости ветра V.
Оценка технико-экономической эффективности оригинальных ОВЭУ на основе КРДС с лопастями с зигзагообразными закрылками показывает, что они при серийном производстве будут иметь удельные капитальные затраты не более 160 тыс. руб./кВт (1,68 тыс.$/кВт, по современному курсу доллара). Срок их окупаемости составит не более 6,6 лет. Поэтому эти ОВЭУ будут вполне конкурентоспособными с другими типами серийно выпускаемых ВЭУ. Кроме того, предлагаемые ОВЭУ будут оказывать минимальный вред окружаюшей среде, т. е. быть практически экологически чистыми.
Эти ОВЭУ с КРДС с номинальной мощностью до 60 кВт смогут использоваться совместно с другими энергоустановками, такими как СЭС, дизель-генераторными агрегатами (ДГА) и газопоршневыми электрогенераторами (ГПЭГ), работающими на природном и попутном газе. ДГА и ГПЭГ в настоящее время являются самыми доступными, имеют лучшее соотношение цены, надежности и качества, обладают высоким КПД, достигающим 40 %, и ресурсом около 40 тыс. моточасов до капитального ремонта [3, 4]. Эти предлагаемые автономные энергокомплексы с ВИЭ могли бы использоваться не только в инфраструктуре на небольших линейных объектах, но и для привода скважинных центробежных, винтовых, винтовых штанговых и других установок для добычи нефти на суше и на море. Номинальная мощность таких более энергетически эффективных энергокомплексов уже в ближайшие годы могла бы достигнуть мощности 120 кВт, что существенно выше, чем у известных современных энергокомплексов с ВИЭ, используемых на нефтегазопромысловых объектах. В частности, возможно использование таких энергокомплексов на базе ВИЭ для энергоснабжения бытовых модулей морских нефтедобывающих платформ [16]. В перспективе энергокомплексы с использованием ОВЭУ с КРДС смогли бы обеспечить работу интеллектуальных скважин, а затем и полностью автономное энергоснабжение интеллектуальных месторождений на суше и на море.
В дальнейшем будут более детально исследоваться вопросы агрегирования предлагаемых ОВЭУ с КРДС с другими энергоисточниками на основе возобновляемых и топливных энергоресурсов.
Заключение
1. Предложены усовершенствованные ОВЭУ на основе КРДС с зигзагообразными закрылками.
В них запуск осуществлялся с помощью роторов Савониуса, работа КРДС – с помощью ротора Н-Дарье, имеющего диаметр и высоту в 3–4 раза больше, чем у ротора Савониуса. Лабораторные эксперименты показали, что с увеличением средней скорости воздушного потока с 2,7 до 7,2 м/с мощность КРДС увеличивалась в 19 раз. Максимальные значения коэффициента мощности (энергетического КПД) КРДС достигали 0,6, что практически в 1,3–1,7 раза выше, чем у серийно выпускаемых ВЭУ, имеющих коэффициент мощности
не более 0,35–0,45.
2. Определены значения коэффициента мощности КРДС при различных геометрических размерах входящих в его состав роторов Н-Дарье и Савониуса
с различным количеством ступеней роторов. Расчеты показали, что с увеличением скорости ветра от 3,0 до 9,0 м/c и диаметра ротора от 5,0 до 15,0 м мощность КРДС увеличивается с 0,245 до 59,6 кВт.
3. Энергетически, экономически и экологически эффективные ОВЭУ с КРДС будут вполне конкурентоспособными с другими типами ВЭУ. Эти ОВЭУ с КРДС с номинальной мощностью до 60 кВт смогут совместно с другими энергоустановками, работающими на природном и попутном газе, использоваться не только в инфраструктуре нефтегазопромысловых объектов, но и для добычи нефти на суше и на море. В перспективе они могут обеспечить работу интеллектуальных скважин, а в дальнейшем и полностью автономное энергоснабжение интеллектуальных месторождений.
1. Туровин О. А., Огнев Е. Н., Кочнев А. Е. Применимость ветро-солнечной энергетики в качестве альтернативного источника электроснабжения нефтяных объектов компании // PROНЕФТЬ. Профессионально о нефти. 2017. № 2 (4). С. 69-74.
2. Голубев С. В. Возобновляемые источники энергии в энергетике газовой отрасли. Перспективы и аспекты применения ВИЭ на объектах ПАО «Газпром» // Газ. пром-сть. 2016. № 12/746. С. 72-76.
3. Голубев С. В. Применение энергоустановок и энергокомплексов на базе возобновляемых и альтернативных источников источников энергии для объектов газового комплекса // Газ. пром-сть. 2017. № 5. С. 2-6.
4. Автономные блочно-комплектные энергетические установки. URL: https://nipom.ru/blog/news/massmedia/avtonomnye-blochno-komplektnye-energeticheskie-ustanovki-bkeu-vsm (дата обращения: 25.09.2023).
5. Закиров Э. С., Закиров С. Н., Индрупский И. М., Аникеев Д. П. Интеллектуальные скважины: преимущества и проблемы // Актульные проблемы нефти и газа. 2018. Вып. 2 (21). URL: https://oilgasjournal.ru/issue_21/zakirov-zakirov-indrupskiy-anikeev-1.html (дата обращения: 02.10.2023).
6. Умные скважины и интеллектуальные месторождения. Нефть в цифре. URL: https://dprom.online/oilngas/umnye-skvazhiny-i-intellektualnye-mestorozhdeniya-neft-v-tsifre (дата обращения: 02.10.2023).
7. Ляхтер В. М., Шполянский Ю. Б. Аэродинамика ортогональных ветроагрегатов // Сб. науч. тр. Гидропроекта. М., 1988. Вып. 129: Ветроэнергетические станции. С. 113-127.
8. Горелов Д. Н. Аэродинамика ветроколес с вертикальной осью вращения. Омск: Полиграф. центр КАН, 2012. 66 с.
9. Моди Форандо. Характеристики ветродвигателя Савониуса // Современное машиностроение. Серия А. 1988. № 10. С. 139-148.
10. Соломин Е. В. Методология разработки и создания вертикально-осевых ветроэнергетических установок: монография. Челябинск: Изд-во ЮУрГУ, 2011. 254 с.
11. Вертикальные ветрогенераторы Маглев. URL: https://solar.manblan.ru/projects/14060/ (дата обращения: 05.12.2019).
12. Ильин Р. А., Шишкин Н. Д. Экспериментальное исследование параметров вертикально-осевых ветро-энергоустановок для привода гребных винтов маломерных судов // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технологии. 2019. № 2. С. 93-100.
13. Shishkin N. D., Ilyin R. A. The energy efficiency research of combined H-Darier-Savonius rotors for autonomous power supply of objects on land and at sea // Journal of Physics: Conference Series Volume 1652 15th International Scientific and Technical Conference (PESPC) 2020 6-9 October 2020. URL: https://iopscience.iop.org/issue/1742-6596/1652/1/01 2028 (дата обращения: 12.11.2020).
14. Баклушин П. Г., Вашкевич К. П., Самсонов В. В. Экспериментальное исследование аэродинамических характеристик ортогональных крыльчатых ветроколес // Сб. науч. тр. Гидропроекта. 1988. С. 98-105.
15. Горелов Д. Н. Энергетические характеристики ротора Дарье (обзор) // Теплофизика и аэромеханика. 2010. Т. 17, № 3. С. 325-333.
16. Шишкин Н. Д., Балтаньязов И. В., Герлов В. Н. Использование возобновляемых источников энергии для энергоснабжения нефтедобывающих морских платформ // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2009. № 2. С. 193-197.
