ИНЖЕНЕРНАЯ ВЫСТАВКА ИЗОБРЕТАТЕЛЬ XXI Бирюков Георгий Владимирович, уч-ся 11 класса III курса физико-математического отделения лицея 1, г.Усолье-Сибирское. Научные руководители профессор ИрГТУ
Агарышев Анатолий Иванович, преподаватель физики Шулья Ирина Петровна, Лицей 1. г. Усолье-Сибирское СОДЕРЖАНИЕ 1. Введение 2. Методы расчета напряженности 3. Исходные данные 3.1 для ИОРТПЦ 11 3.2 для Усольского телецентра - 4. Расчет напряженности поля ИОРТПЦ 5. Расчет напряженности поля вблизи Усольского телецентра 13 1.
Сравнение результатов измерения и расчетов в зоне обслуживания Усольского телецентра 2. Заключение 3. Приложение 5.1 программа расчета напряженности поля 17 5.2 таблица измерения напряженности поля вблизи Усольского ретранслятора 4. Список литературы 24 ВВЕДЕНИЕ Для решения вопросов проектирования и эксплуатации радиотелепередающих цетров и других радиотехнических систем необходимо рассчитывать напряженности поля радиоволн
УКВ диапазона. На основе этизх расчетов устанавливаются санитарно-защитные зоны СЗЗ радиотехнических объектов, зоны ограничения застройки, а также зоны обслуживания объектов. Особый интерес вызывают вопросы электромагнитной экологии, что обусловлено резким увеличением числа передатчиков УКВ и СВЧ диапазонов, используемых в радио- и телевещании, для спутниковой, сотовой связи и т.д. источниками электромагнитных полей антропогенного происхождения являются также персональные компьютеры,
бытовые приборы, такие как СВЧ-печи, телевизоры. В результате возросли фоновые уровни электромагнитных полей, а также количество зон повышенной опасности, в которых напряженности поля существенно выше фоновых. Отрицательное влияние достаточно интенсивного электромагнитного поля на организмы людей в настоящее время доказано, на основе чего установлены санитарные нормы предельно допустимые уровни электромагнитного поля - ПДУ. В рамках нашей работы анализировались результаты измерений уровня электромагнитного поля
окрестностях иркутского областного и усольского радиотелепередающих центров для того, чтобы выяснить, не превышают ли значения напряженности поля ПДУ. Также было разработано программное обеспечение для расчетов напряженности поля, позволяющее учитывать диаграммы направленности антенн различного назначения. С помощью соответствующих программ можно определить зону обслуживания радиотелепередающего центра для заданной чувствительности приемников, а также санитарно-защитную зону объекта.
МЕТОДЫ РАСЧЕТА НАПРЯЖЕННОСТИ ПОЛЯ УКВ. Предположим, что в свободном пространстве т.е. в однородной непоглащающей среде, относительная диэлектрическая проницаемость которой равна единице помещен изотропный излучатель - воображаемый точечный излучатель, равномерно излучающий радиоволны во всех направлениях. Обозначая через Р1 излучаемую источником мощность, определим плотность потока энергии вектор Пойнтинга на расстоянии r от источника радио волн рис.1, основываясь на том, что излучаемая энергия
равномерно распределяется по поверхности сферы радиуса r. Выражая мощность излучателя в Вт, а линейные размеры - в м, получим для численного значения вектора Пойнтинга выражение , Втм2 1.1 Рис.1. К определению напряженности поля волны, создаваемой изотропным излучателем В принятой системе единиц среднее за период численное значение вектора Пойнтинга выражается формулой , Втм2 1.2 где напряженности электрического и магнитного полей связаны
между собой соотношением , ам. 1.3 Здесь величина 120 представляет собой волновое сопротивление свободного пространства и выражается в омах. Подставляя формулу 1.2 в 1.3, получаем , Втм1.4 Приравнивая выражения 1.1 и 1.4 и решая полученное уравнение относительно Ед, находим , Вм. 1.5 Рис. 2. Диаграммы направленности антенн, направленной А и изотропной В. В реальных условиях изотропные излучатели, конечно, не применяются, а используются
антенны, обладающие направленным действием. Предложим, что рядом расположены направленная А и изотропная В антенны. На рис. 2 схематически диаграммы направленности обеих антенн. Изотропная антенна, как и следовало ожидать, обладает круговой диаграммой направленности Если обе антенны излучают одинаковые мощности Р1, то ясно, что в пункте приема, который достаточно удален от антенн и на который ориентирована направленная антенна, большая напряженность поля создается от направленной
антенны, так как она концентрирует излучаемую энергию в желаемом направлении. Будем постепенно увеличивать подводимую к изотропной антенне мощность до тех пор, пока она не создаст такое же поле, что и направленная антенна. Множитель D1, показывающий, во сколько раз следует увеличить мощность, подводимую к изотропной антенне, чтобы она создавала такую же напряженность поле, что и направленная, носит название коэффициента направленности
или коэффициента усиления Коэффициент усиления несколько отличается от коэффициента направленности, так как он учитывает кпд антенны. Однако эту тонкость не будем принимать за внимание, потому что у многих направленных антенн кпд близок к 100 Таким образом, направленная антенна по создаваемой ею в месте приема напряженности поля эквивалентна изотропной антенне, которая излучает в D1 раз большую мощность. Это позволяет представить формулу для напряженности поля, создаваемой в свободном
пространстве направленной антенной, в следующем виде Вм. 1.6 Амплитудное значение напряженности поля выражается формулой , Вм. 1.7 Выражение для мгновенного значения напряженности электрического поля радио волны можно записать в форме Вм, 1.8 где - волновой множитель. Единицы измерения величин, входящих в формулы 1.5 - 1.7, не очень удобны для практического применения, так как напряженность поля выражена в
Вм, а расстояние - в м. Выражая мощность в кВт, расстояние - в км, а напряженность поля - в мВм, получаем , мВм 1.9 для действующего значения напряженности поля и , мВм 1.10 для амплитудного. В течение долгого времени условия распространения волн было принято оценивать напряженностью электрического поля, создаваемого передатчиком в месте приема. Такой критерий был более или менее оправдан в условиях, когда радиосвязь осуществлялась в диапазоне длинных, средних и, частично, коротких волн.
Степень направленности антенны характеризуется ее коэффициентом направленности D или усилением по отношению к изотропному излучателю коротких волн. В связи с широким применением в последние годы диапазона УКВ более рационально характеризовать условия приема мощностью, создаваемой на входе приемного устройства, ибо чувствительность современных приемных устройств принято выражать мощностью на входе, требуемой
для уверенного приема сигналов. Для этого необходимо знать направленной антенны D2. Однако это обстоятельство не ограничивает область применения такого метода, так как направленность передающей антенны D1 также должна быть известна. Наконец, чтобы исключить конкретные типы антенн, можно предположить, что обе антенны изотропны, т. е. D1 D2 1. z Pr y x Рис. 3. Сферические координаты точки наблюдения
Наглядное представление о распределении энергии волн дает амплитудная характеристика направленности, определяемая зависимостью амплитуды напряженности создаваемого антенной поля или величины, ей пропорциональной от направления в пространстве. Направление определяется азимутальным и меридиональным углами сферической системы координат, как это показано на рис. 3. При этом поле измеряется на одном и том же достаточно большом расстоянии r от антенны и предполагается, что потери в среде отсутствуют.
Графическое изображение характеристики направленности называют диаграммой направленности. Пространственная диаграмма направленности изображается в виде поверхности f Построение такой диаграммы неудобно. Поэтому на практике обычно строят диаграммы направленности в какой-нибудь одной плоскости, в которой она изображается плоской кривой f или f в полярной или декартовой системе координат. Пространственная диаграмма направленности, у которой максимальное значение равно
единице, называется нормированной диаграммой и обозначается как F Она легко получается из ненормированной диаграммы путем деления всех ее значений на максимальное F, f,fmax 1.12 Простейший излучатель в виде элементарного диполя имеет тороидальную диаграмму направленности, показанную на рис. 4 в полярных координатах и выражаемую уравнением Е Е0 sin , 1.13 где Е0 напряженность поля в направлении максимума т.е. при 90о угол, отсчитываемый от
оси диполя. На рис. 5, а показан пример игольчатой диаграммы. Основное излучение антенны с такой диаграммой направленности сконцентрировано в пределах небольшого телесного угла. На рис. 5, б показан пример диаграммы направленности специальной формы, определяемой в вертикальной плоскости уравнением Е Е0 cosec , 1.13 где Е0 коэффициент пропорциональности угол в вертикальной плоскости, отсчитываемый относительно горизонта.
Такие диаграммы желательно иметь в некоторых типах радиолокационных станций, например в самолетных радиолокаторах наземных объектов. При отражении от таких объектов, находящихся от поверхности земли на различных расстояниях от самолета в пределах радиуса действия радиолокатора, уровень отраженного сигнала на входе приемника будет сохраняться неизменным. Направленное действие антенны часто оценивают по углу раствора диаграммы направленности, который также
называют шириной диаграммы. Под шириной 20,5 диаграммы главного лепестка подразумевают угол между направлениями, вдоль которых напряженность поля уменьшается в раз, по сравнению с напряженностью поля в направлении максимума излучения, а поток мощности соответственно уменьшается вдвое. В некоторых случаях под шириной 20 подразумевают угол между направлениями ближайшими к направлению максимума, вдоль которых напряженность поля равна нулю.
Для сравнения между собой направленных антенн вводят параметр, называемый коэффициентом направленного действия КНД. Коэффициент направленного действия число, показывающее, во сколько раз пришлось бы увеличить мощность излучения антенны при переходе от направленной антенны к ненаправленной при условии сохранения одинаковой напряженности поля в месте приема при прочих равных условиях , 1.14 где мощность излучения ненаправленной антенны мощность излучения направленной антенны.
90о 180о 0о 360о а б Рис. 4. Диаграмма направленности элементарного диполя а проекция в плоскости, перпендикулярной оси диполя б проекция в плоскости, проходящей через ось диполя. Боковые лепестки Главный лепесток Направление главного излучения Задний лепесток а б Рис. 5. Коэффициент направленного действия в направлении максимального излучения для реальных антенн достигает значений от единиц до многих тысяч.
Он показывает тот выигрыш в мощности, который можно получить за счет использования направленного действия антенны, но он не учитывает возможных потерь в направленной антенне. Для суждения о выигрыше, даваемом антенной, при учете как ее направленного действия, так и потерь в ней служит параметр, называемый коэффициентом усиления антенны. Он равен произведению КНД на к.п.д 1.15 Учитывая 1.12, получаем .
1.16 Таким образом, коэффициент усиления показывает, во сколько раз нужно уменьшить или увеличить мощность, подводимую к направленной антенне, по сравнению с мощностью, подводимой к идеальной ненаправленной антенне без потерь, для того чтобы получить одинаковую напряженность поля в рассматриваемом направлении. Если не делается специальных оговорок, то под коэффициентом усиления так же, как и под коэффициентом направленного действия подразумевается его максимальное значение, соответствующее направлению максимума
диаграммы направленности. ДН антенны H R H h r Рис. 6. Расчеты действующих значений напряженности выполняются по методике 2 при задании излучаемых мощностей, КНД и нормированных ДН передающих антенн в вертикальной и азимутальной плоскостях. При этом учитываются уровни боковых лепестков ДН, а также рельеф местности и высоты зданий. Напряженность поля в нашей работе рассчитывалась по формуле 1.6, в которую были внесены выше сказанные
поправки 1. 17 где Р мощность, Вт R расстояние от фазового центра антенны до точки наблюдения с высотой h от основания опоры, м F нормированная диаграмма направленности ДН в вертикальной плоскости, arctgH hr угол места, r расстояния от основания опоры до проекции точки наблюдения на уровень h, так что R rsin F 1 нормированная ДН в горизонтальной плоскости. Рис. 6. Функция F для типовой передающей антенны антенная решетка задается
формулой 1.18 В этой формуле коэффициент b принимает значения 2 для антенной решетки всех радиотелепередающих центров до модернизации в 1998 г и 1,3 после модернизации согласно 3. Для антенн типа полуволновой вибратор функция F задается иначе 1.19 или, если угол отсчитывается от оси диполя рис. 7 1.20 При проектировании и эксплуатации современных радиотехнических объектов важно учитывать ПДУ воздействия электро-магнитных полей ЭМП на здоровье людей.
В таблице приведены значения ПДУ для некоторых частот. Частота МГц48,488,4192 300ПДУ, вм5,04,03,0 2,5В соответствии с официально утвержденной методикой 2 рассчитывались значения нормированной суммарной напряженности поля S в зависимости от расстояний от опоры 1.21 где индекс суммирования k соответствует номеру передатчика и меняется от 1 до 7, а санитарные нормы не нарушаются при
S 1. Направление излучения Рис. 7. Изложенные в этом разделе методы расчета напряженности поля использовались для расчетов напряженности поля в ряде конкретных ситуаций. ИСХОДНЫЕ ДАННЫЕ 1. Для ИОРТПЦ Обозначим действующие передатчики типов АРТС и Дождь- 2 номерами 1- 4. Эти передатчики имеют следующие параметры излучения мощности Р 1 Р 2 5 кВт и Р 3 Р 4 3 кВт рабочие частоты f 1 80
МГц, f 2 96 МГц, f 3 f4 70 МГц коэффициент усиления передающих антенн D 1 D 2 6.56, D 3 D 4 9.84 высоты фазовых центров антенн от основания опоры Н1 Н2 161 м, Н3 Н4 155 м Для планируемых к размещению передатчиков с номерами 5 - 7 мощности Р 5 Р 6 Р 7 1 кВт рабочие частоты f 5 474 МГц, f 6 506 МГц, f 7 570 МГц коэффициент усиления передающих антенн
D5 4, D6 D7 3.2 высоты фазовых центров антенн от основания опоры Н 5 Н 6 Н 7 180 м. Согласно СанПин 2.2.4 2.1.8.056 - 96 примем следующие предельно допустимые уровни ПДУ напряженности поля Е для f1 и f2 ПДУ1 ПДУ2 4 вм, для f3 и f4 ПДУ3 ПДУ4 5 вм, для f5 - f7 ПДУ5 ПДУ6 ПДУ7 6 вм. 2. Для Усольского телепередающего центра
Обозначим действующий передатчик номером 1. Этот передатчик имеет следующие параметры излучения мощность Р 0.1 кВт рабочая частота f 66 МГц коэффициент усиления передающей антенны D 8 дБ высота фазового центра антенны от основания опоры Н 127 м Согласно СанПин 2.2.4 2.1.8.056 - 96 примем следующий предельно допустимый уровень ПДУ напряженности поля Е ПДУ 4.5 вм. Расчет напряженности поля в окрестностях
ИОРТПЦ Напряженность поля для каждого из передатчиков ИОРТПЦ, всего их семь, рассчитывает компьютерная программа, которую мы написали в результате нашей научной работы на языке TURBO PASCAL 7.0. Программа разработана таким образом, что в нее входит четыре типа антенн первый тип антенная решетка с коэффициентом b, равным 2 второй тип антенная решетка с коэффициентом b, равным 1.3 третий тип антенны это полуволновой вибратор четвертый тип определяет сам пользователь
вводит формулу функции F для конкретной антенны. В программу вводятся исходные данные мощность P в кВт коэффициент усиления передающей антенны D высота фазового центра от основания опоры Н в м R- расстояние от фазового центра опоры до точки наблюдения с высотой h от основания опоры, в м и количество точек, где были проведены измерения напряженности поля. Затем программа производит расчет и выводит на экран систему координат, где строится график зависимости
напряженности поля, в мкВм, от расстояния, в км. Мы видим, что с увеличением расстояния от ретранслятора график убывает, а также на графике могут быть видны незначительные скачки излучения напряженности поля, это зависит от рельефа рассматриваемой местности. На некоторой высоте, где установлен передатчик, находится щит, который снижает излучение передатчика до некоторого расстояния r 300 м. Так как в этом радиусе расположен пункт слежения за радио и телевещанием.
На рисунке 8 построен примерный график, который может изобразить нам компьютерная программа. Расчет напряженности поля в близи зоны Усольского радиотелецентра Напряженность поля в близи Усольского района рассчитывает компьютерная программа, написанная на языке TURBO PASCAL 7.0. В программу входит четыре типа антенн первый тип - антенная решетка с коэффициентом b, равным 2 второй тип - антенная решетка с коэффициентом b, равным 1.3 третий тип - это полуволновой
вибратор четвертый тип определяет сам пользователь - вводит функцию F для конкретной антенны. В данном случае используется антенна типа 3-х элементный волновой канал с круговой поляризацией, направленная на город. Рис 9 На рис.9 представлена примерная карта города Усолья - Сибирского трансляции телепередач на город. Данные измерения напряженности поля указаны в таблице.
График зависимости напряженности поля Ед, мВм от расстояния R, км, построен на рис. 10. Как мы видим, с увеличением расстояния от ретранслятора напряженность поля убывает. На графике максимальное расстояние 18 км. Рис.10 Сравнение результатов измерения и расчетов Усольского радиотелецентра. Таблица 1. Расстояние от ретранслятора, км Напряженность поля, мкВмРасчетные данныеЭкспериментальные
данные24741,544662,43209,2199522,5299212 59041237,1120455767,931835,564839806537, 820897397,513518305,619959242,2233911162 ,9229,51593,6890 В таблице указаны расчетные данные, которые рассчитала компьютерная программа и экспериментальные, которые были измерены специальным прибором. Если сравнить данные полученные в результате расчета и экспериментальные, то они несколько отличаются друг от друга. Экспериментальные данные больше, чем расчетные, это может зависеть от рельефа рассматриваемой местности.
Также оказало влияние то, что в расчетах не учитывалась ДН передающей антенны в азимутальной плоскости. ЗАКЛЮЧЕНИЕ При выполнении данной работы были получены следующие результаты 1 были изучены методы расчета напряженности поля 2 была разработана программа, рассчитывающая напряженность электромагнитных волн, излучаемых телерадиопередатчиками, в зависимости от расстояния до опоры передающей антенны 3 были рассчитаны значения напряженности поля
вблизи ИОРТПЦ, также были рассчитаны значения нормированной суммарной напряженности, где санитарные нормы не нарушаются 4 были рассчитаны значения напряженности поля вблизи Усольского телерадиопередающего центра и сделаны сравнения с экспериментальными данными. ПРИЛОЖЕНИЕ Приложение 1 программа расчета напряженности поля. uses crt,graph,omenu const ffi 1 NBg blue1 NFg white15 HBg white15 HFg black0 BC black0
SC lightcyan11 col 200 deltarm 90 var vf text VMenu OVMenu HMenu OHMenu HVMenu OHVMenu p,d,hb,em real i,j,choice,errc, a,x,Hmenuchoice,len integer rm longint ord array1 col of real del array1 10 of real delstr,si,AStr,vstr string ch,rk char inputis boolean Процедуры ввода данных procedure inputvaluexi,yiinteger var znreal begin vstr while rk 13 DO begin rkreadkey if rk 47andrk 58orrk46andlen 10 then begin vstrvstrrk
lenlengthvstr gwritexyxilen,yi1,rk,3,2 end end valvstr,zn,errc end procedure input begin gwritexy1,5,Мощность ,3,2 inputvalue11,4,p readln gwritexy1,6,К. у. антенны ,3,2 inputvalue1,6,d readln gwritexy1,7,Высота передающей антенны ,3,2 inputvalue1,7,hb readln end Функция выводит осн. меню на экран и возвращает номер выбранного пункта меню Function ddtinteger begin HVMenu.init gwritexy0,1 0,0
HVMenu.SetHorItems00,00,80,01,NBg, NFg,HBg,HFg,BC,SC,1,1,BorderOn,ShadowOff , File Антенна HVMenu.SetVerItems01,00,01,10,03,NBg,NFg ,HBg,HFg,BC,SC,4,1,BorderOn,ShadowOff, Данные Выход HVMenu.SetVerItems2,6,01,29,04,NBg,NFg,H Bg,HFg,BC,Sc, 4,1,BorderOn,ShadowOff, Ант. решетка 1 - 1,3 Ант. решетка 2 -
2 Диполь HMenu.EraseOKFalse XHVMenu.MenuResultfalse,true ddtx end Функции расчета напряженности function falfareal begin case choice of 1 falfa12cos1.3pisinarctanhbrm3 2 falfa12cos2pisinarctanhbrm3 3 falfacospi2sinarctanhbrmcosarctanhbrm end end function Rbreal begin rbrmsinarctanhbrm end function E2real begin E230pdsqrfalfasqrffisqrRb end Заполнение массива ординат procedure ordinates begin rm1 for i1 to col
do begin rmrmdeltarm ordi1000SQRTE2 х1000, т.к. ед. изм мВм end end Максимальное значение напряженности procedure Emaximum var iinteger maxreal begin Maxord1 if col 1 then for i2 to col do if ordi Max then Maxordi if max0 then max1 Emmax end Сохранение результатов расчета в файл results.txt procedure ToFile begin assignvf,results.txt rewritevf rm0 for i1 to col do begin rmrmdeltarm writelnvf,rm, m ,
ordi05, mVm end end Инициализация графики procedure grinit var grDriver Integer grMode Integer ErrCode Integer begin grDriver Detect InitGraphgrDriver, grMode,cbpbgi ErrCode GraphResult if ErrCode 0 then WritelnGraphics error, GraphErrorMsgErrCode end procedure drawcoords Оси координат begin setcolordarkgray Oy line100,445,100,30 line99,445,99,30 line99,30,96,35 line100,30,103,35
outtextxy25,23, Е, мВм Ox line95,440,515,440 line95,441,515,441 line515,440,510,437line515,441,510,444 outtextxy525,445,R, м end procedure drawgridСетка begin setcolorlightgray Горизонтальная j40 for i1 to 10 do begin line100,440-j,500,440-j jj40 end Вертикальная jround80ln1.91 for i1 to 6 do begin line100roundj,440,100roundj,40 jjround80lni1.8 end end procedure valuesРазметка сетки begin По вертикали del1em10
Цена деления for i2 to 10 do delidel1deli-1 setcolordarkgray outtextxy90,445,0 j40 for i1 to 10 do begin strdeli01,delstr outtextxy90-lengthdelstr8,438-ji,delstr end По горизонтали j95round80ln1.91 outtextxyj,445,3 jjround80ln2.8 outtextxyj,445,6 jjround80ln3.8 outtextxyj,445,9 jjround80ln4.8 outtextxyj,445,12 jjround80ln5.8 outtextxyj,445,15 jjround80ln6.8 outtextxyj,445,18 end Построение графика procedure drawgrafic var dltinteger x1,x2,y1,y2integer begin setcolorchoice1 x1100-round2ln1.91
for i1 to col do begin y1440-round400ordiem y2440-round400ordi1em if i 1andi 40 then begin x1x1round2ln1.91 x2x1round2ln1.91 end if i 40andi 80 then begin x1x1round2ln3.71 x2x1round2ln3.71 end if i 80andi 120 then begin x1x1round2ln5.51 x2x1round2ln5.51 end if i 120andi 160 then begin x1x1round3ln7.31 x2x1round3ln7.31 end if i 160andi 200 then begin x1x1round4ln9.11 x2x1round4ln9.11 end linex1,y1,x2,y2 linex1,y1-1,x2,y2-1 linex1,y1-2,x2,y2-2 delay20 end end Графические процедуры procedure drawing1st
Инициализирует графику, подготавливает экран begin grinit setbkcolor15 cleardevice setcolordarkgray rectangle10,10,getmaxx-10,getmaxy-10 drawgrid drawcoords end procedure drawing2nd Выводит график на экран begin drawgrafic readln closegraph end begin ClrScr Inputp100 d8 hb127 grinit repeat cleardevice i2 repeat addt until a 0 Hmenuchoicea div 100 Case Hmenuchoice of 1 begin choicea mod 100 if choice2 then break else begin input
inputistrue end end 2 if notinputis then begin gwritexy17,10, Сначала необходимо ввести даннные ,5,1 chreadkey continue end else begin choicea mod 100 Drawing1st Ordinates EMaximum ToFile Values Drawing2nd end end until false HVMenu.Done cleardevice closegraph writep12, ,d12, ,hb12 end. Приложение 2 таблица измерения напряженности поля вблизи
Усольского ретранслятора. Места проведения измерений Расстояние Направленная антенна Круговая антенна кмЕ видео, мкВм Е звук, мкВмЕ видео, мкВмЕ звук, мкВм1.Тайтурка111781123162002.Средний112 801267083543.Мальта663035414127084.Белор еченск770744612587085.Тракт2446619952511 870806.Зеленый гор. ул.Энергетиков4177805010251214127.Зелены й гор. ул.Фурманова2,5125904466501211228.У-
С. Горбольница 5,539801258223812589.У-С. ул. Горького635481122177889010. У- С. ул.Крупской5354814122623123811. У- С. жд. переход2,41995270805623223812. У- С. площадь 4631019952512112213. У- С. Комсомольский пр-т 52818708177889014. У- С. мр-н Привокзальный7199570870840015. У- С. Ленинский пр-т8199579489044616. У- С. Восточ. окраина9233963070850017. Тельма, Зап. окраина15890354446224 г.