СПОСОБ МОНИТОРИНГА ГОЛОЛЕДНО-ВЕТРОВЫХ НАГРУЗОК ВОЗДУШНЫХ ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧ Российский патент 2020 года по МПК G01R31/08 

Описание патента на изобретение RU2732037C1

Изобретение относится к области электроизмерительной техники и может быть использовано для обнаружения гололедно-ветровых нагрузок и прогнозирования гололедообразования посредством выявления коронирования воздушных высоковольтных линий.

Прототипом предлагаемого изобретения является патент №2612742 - Способ мониторинга воздушных линий электропередач. Здесь, в выражении (11), при расчетах принято заполнение междуфазного пространства в объема атмосферными осадками.

Однако, как показывают исследования (О.Г. Богаткин, Г.Г. Тараканов Основы метеорологии. - СПб, изд. РГГМУ 2006 - 232 с. - с. 85) «…Скорость падения капель дождя составляет 4-8 м/с, а самые крупные капли имеют скорость падения 10 м/с... Наиболее сильный дождь (дождь с наибольшей интенсивностью) выпадает из кучево-дождевых облаков. В этом случае интенсивность дождя может превышать 10 мм/ч...».

Для куба с сторонами 1 м, заполненного дождевыми осадками, при их выпадании верхний слой достигнет нижней поверхности за время 1 м/(4 м/с)=0,25 с. (скорость принята для наихудшего случая). Например, при интенсивности дождя 20 мм/ч=5,6*10-6 м/с за 0,25 с можно получить уровень воды - 5,6*10-6 м/с*0,25 с=1,4*10-6 м с объемом воды в кубе 1,4*10-6 м*1 м2=1,4*10-6 м3. От объема в 1 м3, данное количество воды составит 1,4*10-4%. Такое относительно малое заполнение практически не изменит электрическую емкость между фазным проводом и землей.

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

Далее, ветровые нагрузки опираются на данных (настройки значений частоты), полученных по значениям (теоретическим) атмосферных осадков, что практически осуществить нереально.

Задача изобретения - определение гололедно-ветровых нагрузок и начала гололедообразования с учетом температуры проводов воздушных линий.

Поставленная цель достигается тем, что генерируют в высоковольтную линию высокочастотный сигнал с определением размаха колебаний проводов при ветровой нагрузке по максимальному изменению взаимного сопротивления, определяемого первым отношением напряжения в конце и током в начале участка и измерении температуры провода и согласно изобретению определяют второе отношение тока в конце и напряжения в начале участка и изменением частоты высокочастотного сигнала добиваются равенства максимальных относительно быстрых и медленных изменений первого и второго отношений, характеризующее значение гололедно-ветровой нагрузки и интенсивности коронирования воздушной линии соответственно.

Вариант устройства, реализующий предлагаемый способ, изображен на фиг. 1, где к одной фазе линии приложен высокочастотный синусоидальный сигнал через генератор Г1 конденсатор связи С2.

В начале и конце контролируемого участка линии установлены датчики напряжения, тока и температуры ДН3-ДН4; ДТ5-ДТ6, ДТЕМ7-ДТЕМ8 соответственно. Они измеряют высокочастотные напряжения, токи и температуру провода U1 - U2; I1 - I2, τ1 - τ2 соответственно и передают информацию о них к приемному устройству ПрУ9.

Известные уравнения для линий с распределенными параметрами при пренебрежении активными потерями в линии с использованием высокочастотного сигнала (см. также прототип).

где, U1, U2, I1, I2 _ комплексные напряжения и токи в начале и конце рассматриваемого участка соответственно; Zc, β, х - комплексное характеристическое сопротивление, коэффициент фазы и длина линии соответственно.

Zc и β определяются известными соотношениями

где, f - частота генератора Г1, L0, Co - погонные индуктивность и емкость линии соответственно.

Для погонных параметров можно записать

где, Dcp, rпр, ε - среднее геометрическое расстояние между проводами; радиус провода; относительная диэлектрическая проницаемость среды соответственно.

С учетом (2) и (3) система (1) представится в виде

где, к=0,208*10-4; к1=138,44; α=Dcp/rпр.

Можно получить значение аргумента кfx1=(2n+1)π/2 (n-целое число),

откуда

Подставляя вышеприведенные значения, в частности при x1=100 км, n=1000 получим -

f=(2n+1)π/(2*100*0,208*10-4)=302,227 кГц.

При значении частоты f=302,227 кГц, посылаемого в линию сигнала генератора из (4) имеем

U 1=jI21lgα*sin (кf x1),

I 1=j(U2/(к1lgα)sin(кf x1)

или

Полученные 2 уравнения (5) позволяют определить значение α=Dcp/rпр, пропорциональное среднему геометрическому расстоянию между проводами 3-х фазной линии. Изменение Dcp через напряжения и токи датчиков, устанавливаемых по краям контролируемого участка характеризует ветровую нагрузку на линию.

Согласно (5) значение α можно получить по 2-м соотношениям и этим можно увеличить достоверность результатов измерений сравнением указанных значений.

Для линий 6-10 кВ

По данным источников:

Кабашов В.Ю. Исследование возможных сближений проводов сельских ВЛ 6-10 кВ при их маятниковых колебаниях под действием ветра // Научно-методический электронный журнал «Концепт». - 2016. - Т. 15. - С. 1331-1335.

Кабашов В.Ю. Повышение надежности сельских воздушных линий электропередачи 10 (6) кВ в условиях воздействия ветровых и гололедных нагрузок. Автореферат диссертации на соискание ученой степени доктора технических наук. Москва 2011. проведены нижеследующие расчеты.

Расчеты при ветровых нагрузках

Для простоты расчетов, во всех нижеприведенных случаях, средний провод принят неподвижным.

При горизонтальном расположении проводов при расстоянии между крайними и средним проводами 1,2 м и радиусе - rпр=0,48 см по известному выражению можно получить Соответственно lgα0=lg(Dcp0/rпр)=lg(151/0,48)=2,5.

Возьмем известное среднее значение минимального расстояния между проводами 0,9 м (с отклонением от нормального положения в 0,3 м) при скорости ветра 12 м/с и получим lgα1=lg(Dcp1/rпр)=lg(113/0,48)=2,37.

При отдалении друг от друга среднего и крайних проводов на 1,5 м lg(188/0,48)=2,59.

Относительное полное изменение отношений напряжения и тока будут: ((U1/I2=U2/I1)1 - (U1/I2=U2/I1)2/(U1/I2=U2/I1)0=((lg(188/0,48)-lg(113/0,48))/lg(151/0,48))*100%=8,8%.

Аналогично, имеем среднее значение минимального расстояния между проводами 0,7 м (с отклонением от нормального положения в 0,5 м) при скорости ветра 16-18 м/с и получим lgα1=lg(Dcp1/rпр)=lg(88/0,48)=2,26.

При отдалении друг от друга на 1,7 м, lg(214/0,48)=2,65.

Относительное полное изменение отношений напряжения и тока будут: ((U1/I2=U2/I1)1-(U1/I2=U2/I1)2)/(U1/I2=U2/I1)0=((2,65-2,26)/2,5)*100%=15,6%.

Расчеты при гололедно-ветровых нагрузках

При гололеде и скорости ветра 12 м/с сближение увеличивается в среднем на 28% и тогда (аналогично, как и выше) получим 0.9-(0,9/100)*28≈0,65 м и lgα1=lg(Dcp1/rпр)=lg(82/0,48)=2,23.

При отдалении друг от друга на 1,75 м, lg(220/0,48)=2,66.

Относительное полное изменение отношений напряжения и тока будут: ((U1/I2=U2/I1)1-(U1/I2=U2/I1)2)/(U1/I2=U2/I1)0=((2,66-2,23)/2,5)*100%=17,2%.

Аналогично, при скорости ветра 16-18 м/с (сближение проводов сохраняем прежним), получим 0.7-(0,7/100)*28≈0,5 м и lgα1=lg(Dcp1/rпр)=lg(63/0,48)=2,12.

При отдалении друг от друга на 1,9 м, lg(239/0,48)=2,7.

Относительное полное изменение отношений напряжения и тока будут: ((U1/I2=U2/I1)1-(U1/I2=U2/I1)2)/(U1/I2=U2/I1)0=((2,7-2,12)/2,5)*100%=23,2%.

Частота маятниковых колебаний проводов, составляющая доли одного Гц, снижающаяся с повышением массы гололедной нагрузки можно прослеживать по изменениям во времени отношений напряжений и токов датчиков контролируемого участка.

Линии 110 кВ и выше

Филиал ОАО «ИНЖЕНЕРНЫЙ ЦЕНТР ЕЭС» «ФИРМА ОРГРЭС» Методы борьбы в России с вибрацией, гололедом и пляской проводов и грозозащитных тросов на воздушных линиях. - презентация.

Согласно имеющимся размерам, например для напряжения 110 кВ, при горизонтальном расположении проводов при минимальном расстоянии между проводами 3,5 м и минимальном диаметре 11,4 мм=0,114 см (радиус - rпр=0,114/2=0,057 см), по известному выражению можно получить Соответственно lgα0=lg(Dcp0/rпр)=lg(441/0,57)=2,89. Поскольку размах колебаний при пляске проводов может достигать размера стрелы провеса проводов, то при рассматриваемой стреле провеса 5 м можно оценить в некоторый момент Dcp и lgα при неизменном rпр

Возьмем случай отклонения проводов в некоторый момент на половину стрелы провеса с отдалением друг от друга на 2,5 м от нормального положения (5/2=2,5 м) Можно получить см, lgα1=lg(Dcp1/rпр)=lg(877/0,57)=3,18.

При сближении всех фаз в момент отдаления друг от друга на 2 м lg(200/0,57)=2,54.

Для данного случая относительное полное изменение Dcp будет: (Dcp2-Dcp1)/Dcp0=(200-559 (877))/441*100%=81 (153)%.

Согласно (5) можно получить = U1/I21*lgα. Для указанных Dcp получим: (U1/I2=U2/I1)01*lgα0=138,44*2,89=400,09 Ом; (U1/I2=U2/I1)1=к1*lgα1=138,44*3,18=440,24 Ом; (U1/I2=U2/I1)21*lga2=138,44*2,54=351,63 Ом.

Относительное полное изменение отношений напряжения и тока будут: ((U1/I2=U2/I1)1-(U1/I2=U2/I1)2)/(U1/I2=U2/I1)0=((440,24-351,63)/400,09)*100%=22%.

Практическая реализация соотношения (5)

Устанавливается значение частоты сигнала генератора, посылаемого в линию согласно (4а).

При отсутствии ветровых нагрузок (например, по отсутствию изменений напряжений и токов, поступающих на центральный сервер с датчиков контролируемого участка) определяются отношения напряжений и токов согласно (5), равные характеристическому сопротивлению к1*lgα. Корректировкой значения частоты генератора добиваются равенства значений отношений напряжений в обоих уравнениях (5). Фиксируется установленное значение частоты fуст.

При наличии изменений отношений напряжений согласно (5), корректировкой значений частоты генератора относительно fуст добиваются максимальное совпадение значений по обоим выражениям (5) для одних и тех же моментов времени.

С использованием (5) можно определить скорость изменения (увеличение или уменьшение) амплитуды и периода колебаний во времени.

Увеличение скорости изменения амплитуды и периода колебаний (с снижением частоты внутри диапазона 2-0,2 Гц) свидетельствует о нарастании массы гололеда (снижение частоты колебаний обусловлено нарастанием массы гололеда).

Полученные данные можно использовать и при проведении мероприятий по повышению надежности линий по ветровым нагрузкам, поскольку выявление фактических таких нагрузок в районах прохождения воздушных линий является актуальной задачей.

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

Прогнозирование гололеда.

Для линий 110 кВ и выше (Техника высоких напряжений. Учебник для студентов электротехнических и электроэнергетических специальностей вузов. Под общей ред. Д.В. Разевига. Изд. 2-е, перераб. и доп. М., «Энергия», 1976.488 с. с ил. - с. 51; 57)

Объемный заряд короны, образовавшийся в один из полупериодов переменного напряжения, за время до изменения полярности провода может переместиться на несколько десятков сантиметров. Вследствие этого объемные заряды обоих знаков совершают возвратно-поступательное движение вблизи провода.

Туман способствует интенсивному образованию капелек воды на провод. Иней, гололед, изморозь вызывают формирование на поверхности провода кристаллов льда, которые значительно усиливают электрическое пол. Дождь и мокрый снег вызывают существенное увеличение неровностей на поверхности провода.

При отсутствии ветра, выше было получено значение (U1/I2=U2/I1)0=lgα0=lg(Dcp0/rпр)=lg(441/0,57)=2,89.

Объемный заряд короны равносильно увеличению диаметра провода (за счет увеличения проводимости объемного заряда) и при учете увеличения rпр всего на 10 см (вместо возможных нескольких десятков см), получим (U1/I2=U2/I1)к=lgαк=lg(Dcpк/rпр)=lg(441/10,57)=1,62,

с относительным изменением (2,89-1,62)/2,89=0,43=43%.

При непогоде интенсивность коронирования линий 110 кВ и выше значительно повышается, что приводит к изменениям отношений напряжений к токам датчиков контролируемого участка не ниже 43% и при минусовой температуре провода можно провести профилактический подогрев проводов.

Вариант реализации ПрУ9 на подстанции и передающего устройства в конце линии приведен на фиг. 2 с графиками прохождения сигналов на фиг. 3.

Генератор высокочастотных синусоидальных сигналов Г1 посылает сигнал в высоковольтную линию через конденсатор С2 и между фазой и землей выделяются напряжения U1, U2, измеряемое датчиками напряжения ДН3, ДН4. На выходах датчиков тока ДТ5, ДТ6 выделяются сигналы действующих значений высокочастотных токов фазы линии I1 и I2. Температуру проводов измеряются датчиками температур ДТЕМ7, ДТЕМ8.

На входы ППУ10 поступают сигналы с выходов ДН4, ДТ5, ДТ6, ДТЕМ7, ДТЕМ8. При этом сигналы с выходов ДН4, ДТ6, ДТЕМ8 поступают на вход ППУ10 посредством ППУ11, установленного в конце линии. На выходе приемо-передающего устройства ППУ10 непрерывно вырабатываются тактовые импульсы U3, поступающие на вход счетчика импульсов Сч12. Сигналы с выходов ДН3 и ППУ10 поступают на устройства деления УД13 и УД14 и отношения сигналов К1 и К2 соответственно приведены на фиг. 2 и фиг. 3. Выходные сигналы К1 и К2 УД13 и УД14 соответственно суммируются арифметическим сумматором СУМ15. Как показано на фиг. 2, например при частоте f1=100 кГц разность К1 и К2 не равна 0 и условие (5) не соблюдается. При f1=100 кГц за указанное время [6 периодов (6Т=6*1/100кГц=60 мкс)] значения К1 и К2 приняты постоянными.

В данном случае на выходе СУМ15 вырабатывается сигнал логической 1 - U4, благодаря которому частота Г1 повышается на заданную ступень f2, что видно на фиг. 2. Здесь же видно, что при этом разность К1 и К2 увеличилась и СУМ15 производит снижение частоты генератора до значения f3 через U5, при котором разность К1 и К2 равняется 0 (с заданной погрешностью). Через элемент ИЛИ16, логическими 1 U4 и U5 производится сброс Сч12.

На фиг. 3 показаны изменения К1 и К2 (случай, когда они равны), характеризующие колебания фазных проводов линии согласно (5) при изменении Dcp в известном диапазоне частот (0,2-2) Гц при ветровой нагрузке. Здесь приведен один период колебания и при частоте 0,2 Гц с значением 1/0,2=5 сек. В данном случае переменная составляющая (для удобства анализа) имеет синусоидальную форму.

Также здесь приведены практически постоянные составляющие К1ик, К2ик с равными значениями и обусловленные изменением эффективного радиуса (повышается проводимость пространства вокруг провода) фазного провода rпр (см. (5)) благодаря его коронированию.

На интервале Δt, согласно фиг. 3, в амплитудном детекторе АД17 и устройстве выделения постоянного напряжения УВПН18 записываются амплитудное значение переменной составляющей ветровой нагрузки и значение постоянной составляющей, обусловленное коронированием соответственно. На этом интервале от заднего фронта U6 срабатывает ждущий мультивибратор ЖМ19 и одновременно происходит сброс запоминающих устройств ЗУ20 и ЗУ21 с записью в них новых значений с выходов АД17 и УВПН18 соответственно.

Расчет значений U8 и U9, для приведенного на фиг. 3 сигнала в виде суммы постоянной и синусоидальной составляющих, производится по известным соотношениям (|К1maxвн|+|К1minвн|)/2 и (|К1maxвн|-|К1minвн|)/2.

При логическом 0 на выходе ИЛИ16, в конце интервала Δt через элементы НЕ22 и ЖМ19 на входах элемента И23 появляются логические 1, открывающая ключи К24 и К25 с индикацией значений ЗУ20 и ЗУ21 в индикаторах ветровой нагрузки ИВН26 и интенсивности коронирования ИИК27.

На интервалах фиг. 3 с частотами f1 и f2 на обоих входах И23 одновременно логические 1 отсутствуют, поэтому ключи закрыты и в индикаторах не будут зафиксированы недостоверные данные.

В индикаторах температур провода ИТЕМ28 и ИТЕМ29 на границах контролируемого участка отображаются соответствующие значения, поступающие от датчиков через ППУ10 и ППУ11.

Похожие патенты RU2732037C1

название год авторы номер документа
Способ мониторинга воздушных линий электропередач 2016
  • Гаджибабаев Гаджибуба Ражидинович
  • Гаджибабаев Эльдар Гаджибубаевич
RU2612742C1
СИСТЕМА ТЕЛЕИЗМЕРЕНИЯ ГОЛОЛЕДНО-ВЕТРОВОЙ НАГРУЗКИ 2006
  • Гаджибабаев Гаджибуба Ражидинович
RU2332765C1
Способ измерения расстояния до места замыкания на землю 2015
  • Гаджибабаев Гаджибуба Ражидинович
  • Гаджибабаев Эльдар Гаджибубаевич
RU2638088C2
Способ поиска земли в электрической сети постоянного оперативного тока 2016
  • Гаджибабаев Гаджибуба Ражидинович
  • Гаджибабаев Эльдар Гаджибубаевич
RU2653162C2
ФИКСАТОР НАПРАВЛЕНИЯ ЗАМЫКАНИЯ НА ЗЕМЛЮ 2010
  • Гаджибабаев Гаджибуба Ражидинович
  • Гаджибабаев Эльдар Гаджибубаевич
RU2468377C2
СПОСОБ ИЗМЕРЕНИЯ РАССТОЯНИЯ ДО МЕСТА КОРОТКОГО ЗАМЫКАНИЯ 2010
  • Гаджибабаев Гаджибуба Ражидинович
  • Гаджибабаев Эльдар Гаджибубаевич
RU2468378C2
ФИКСАТОР НАПРАВЛЕНИЯ КОРОТКОГО ЗАМЫКАНИЯ 2005
  • Гаджибабаев Гаджибуба Ражидинович
  • Гайдаров Рамазан Магомедович
RU2328752C2
Способ выделения ортогональных составляющих гармонического напряжения известной частоты 1990
  • Гаджибабаев Гаджибуба Ражидинович
  • Овчаренко Николай Ильич
SU1775681A1
СПОСОБ ОПРЕДЕЛЕНИЯ ПОСТОЯННОЙ СОСТАВЛЯЮЩЕЙ НАПРЯЖЕНИЯ 1990
  • Гаджибабаев Г.Р.
RU2028623C1
Устройство для автоматического управления компенсирующей установкой 1979
  • Тарикулиев Измир Ярахмедович
  • Кудинов Анатолий Иванович
  • Онищенко Александр Андреевич
  • Гаджибабаев Гаджибуба Ражидинович
SU941969A1

Иллюстрации к изобретению RU 2 732 037 C1

Реферат патента 2020 года СПОСОБ МОНИТОРИНГА ГОЛОЛЕДНО-ВЕТРОВЫХ НАГРУЗОК ВОЗДУШНЫХ ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧ

Изобретение относится к области электроизмерительной техники. Сущность: генерируют в высоковольтную линию высокочастотный сигнал с измерением температуры провода. Определяют размах колебаний проводов при ветровой нагрузке по максимальному изменению взаимного сопротивления, определяемого первым отношением напряжения в конце и током в начале участка. Определяют второе отношение тока в конце и напряжения в начале участка, и изменением частоты высокочастотного сигнала добиваются равенства максимальных относительно быстрых и медленных изменений первого и второго отношений, характеризующего значение гололедно-ветровой нагрузки и интенсивности коронирования воздушной линии соответственно. Технический результат: определение гололедно-ветровых нагрузок воздушных линий выше 6 кВ и начала гололедообразования на основе коронирования воздушных линий выше 110 кВ с учетом температуры провода. 3 ил.

Формула изобретения RU 2 732 037 C1

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

Документы, цитированные в отчете о поиске Патент 2020 года RU2732037C1

СПОСОБ ОБНАРУЖЕНИЯ ГОЛОЛЕДНЫХ, ИЗМОРОЗЕВЫХ И СЛОЖНЫХ ОТЛОЖЕНИЙ НА ПРОВОДЕ И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ 2013
  • Титов Дмитрий Евгеньевич
  • Петренко Станислав Александрович
  • Угаров Геннадий Григорьевич
RU2554718C2
СПОСОБ КОНТРОЛЯ ТЕМПЕРАТУРЫ ПРОВОДОВ ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧИ 2013
  • Мустафин Рамиль Гамилович
  • Мохаммед Кхалиль Сультан Абдулла
RU2547837C1
СПОСОБ ОБНАРУЖЕНИЯ ГОЛОЛЕДА НА ПРОВОДАХ ВОЗДУШНЫХ ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧИ 2013
  • Мустафин Рамиль Гамилович
  • Писковацкий Юрий Валерьевич
  • Хакимзянов Эльмир Фердинатович
RU2537380C1
СПОСОБ ИЗМЕРЕНИЯ ГОЛОЛЕДНОЙ И ВЕТРОВОЙ НАГРУЗОК НА ВОЗДУШНЫХ ЛИНИЯХ ЭЛЕКТРОПЕРЕДАЧИ 2011
  • Мустафин Рамиль Гамилович
  • Котельникова Елена Евгеньевна
RU2461941C1
СПОСОБ ОБНАРУЖЕНИЯ ГОЛОЛЕДНЫХ ОБРАЗОВАНИЙ НА ПРОВОДАХ И ГРОЗОЗАЩИТНЫХ ТРОСАХ ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧИ 2010
  • Куликов Александр Леонидович
RU2409882C1
СПОСОБ ОБНАРУЖЕНИЯ ПОЯВЛЕНИЯ ГОЛОЛЕДА НА ПРОВОДАХ ЛИНИИ ЭЛЕКТРОПЕРЕДАЧИ 2009
  • Мустафин Рамиль Гамилович
RU2399133C1
CN 106646115 A, 10.05.2017
CN 105182119 A, 23.12.2015
US 20130249321 A1, 26.09.2013
US 20130092678 A1, 18.04.2013
DE 102016222208 B3, 22.02.2018.

RU 2 732 037 C1

Авторы

Гаджибабаев Гаджибуба Ражидинович

Гаджибабаев Эльдар Гаджибубаевич

Даты

2020-09-10Публикация

2019-08-12Подача