Труды МАИ. Выпуск № 100
http://trudymai.ru/
УДК 621.396.967
Особенности реализации многолучевых режимов съемки с частотным разделением лучей в космических радиолокаторах синтезирования апертуры на базе активных фазированных антенных решеток
Булыгин М.Л
Научно-исследовательский институт точных приборов, ул. Декабристов, 51, Москва, 127490, Россия e-mail: Bulygin04@gmail.com
Аннотация
Одновременная регистрация аналоговым трактом АФАР космического РСА эхо-сигналов на разных несущих частотах при реализации многолучевых режимов съемки с частотным разделением приемных лучей приводит к искажению геометрии радиолокационного визирования. В статье проанализировано влияние данного эффекта на реализуемые характеристики многолучевого прожекторного и многолучевого сканирующего визирования, а также предложены алгоритмы компенсации данного эффекта за счет специального чередования номиналов несущих частот или оперативного управления положением многолучевой ДН в горизонтальной плоскости.
Ключевые слова: радиолокатор с синтезированной апертурой, активная фазированная антенная решетка, многолучевые режимы съемки.
Введение
Реализация в системах радиолокационного синтезирования апертуры (РСА) космического базирования режимов многолучевого визирования позволяет обходить фундаментальные ограничения, накладываемые эффектом неоднозначности на создание широкозахватных режимов съемки высокого разрешения [Ошибка! Источник ссылки не найден., 1, 6].
Одна из концепций реализации многолучевого визирования предусматривает пространственно-частотное кодирование формы сигнала на излучение с азимутальным цифровым формированием лучей диаграммы направленности (ДН) на прием [3, 4, 17]. Примером реализации являются многолучевой прожекторный и многолучевой сканирующий режимы съёмки, позволяющие существенно улучшить характеристики радиолокационного визирования по сравнению с однолучевыми режимами съемки [Ошибка! Источник ссылки не найден., 8, 9].
Применение многолучевых режимов съемки с азимутальной многолучевой ДН на прием возможно в РСА с сегментированными активными фазированными антенными решетками (АФАР), обеспечивающими многоканальную цифровую обработку сигналов в каждом из сегментов [10, 11, 19]. При этом антенные системы космических РСА реализуются на базе АФАР с аналоговыми фазовращателями приемо-передающиих каналов [12, 13, 19].
В многолучевых режимах съемки с азимутальной цифровой ДН на прием при использовании частотного кодирования формы зондирующего импульса реализуется следующий алгоритм формирования многолучевого визирования. На
излучение формируется импульсная последовательность сигналов, представляющих собой пачку из нескольких суб-импульсов, где каждому суб-импульсу соответствует своя несущая частота ^ и свое азимутальное положение максимума ДН в
горизонтальной плоскости ак [4, 7]. При этом в аналоговой АФАР по вертикали, для излучения каждого суб-импульса в одинаковом вертикальном направлении (соответствующему визированию полосы съемки под углом е), реализуется линейное фазовое распределение, которое определяется приращением фаз между элементами решетки. Данное приращение между соседними элементами решетки, определяется следующим соотношением из теории антенн [14]:
2л - /к . К . ), (1)
у С
где / - центральная частота излучаемого зондирующего суб-импульса, ку - шаг антенной решетки по вертикали, с - скорость распространения радиоволн.
Следовательно, на излучение в многолучевом режиме съёмки положения фазовращателей АФАР должны меняться между излучениями каждого суб-импульса в соответствии с несущей частотой последующего излучаемого суб-импульса.
На прием эхо-сигналы со всех азимутальных направлений, облученные зондирующими суб-импульсами, достигают АФАР одновременно, что делает невозможным выполнение своевременной подстройки фазовращателей по вертикали для приема сигналов каждой частоты аналогично излучению.
Таким образом, формулируется задача исследования - провести анализ одновременного приема многочастотного сигнала в аналоговой АФАР при
реализации многолучевых режимов съемки с многолучевой азимутальной ДН на
прием.
Искажение геометрии многолучевого визирования
На прием, приращение фаз между элементами решетки по вертикали определяется аналогично расчету на излучение:
2ж ■ Ы ■ ¡гу ■ ьш{вгОу)
Лфгу=-у-^, (2)
с
где ¥0 - центральная частота всего рабочего диапазона частот (в данном случае -частота на которую настроена АФАР по вертикали на прием), Иу - шаг антенной решетки по вертикали, вг0у - угол отклонения луча ДН в вертикальной плоскости на прием.
В случае приема сигнала с частотой / к при неизменном приращении фаз между каналами Лфгу и шаге антенной решетки по вертикали Иу соотношение (2) преобразуется к следующему виду:
2п ■ / ■ Иу ■ вт(в^)
Лфгу =-у-. (3)
' с
Из представленных соотношений, с учетом корректного нацеливания луча ДН при излучении суб-импульсов (в соответствии с соотношением 1) следует следующая зависимость углового положения максимума приемо-передающей ДН для сигналов с частотой ^ при их регистрации в АФАР, приемные каналы которой настроены на регистрацию сигнала частотой ¥0 с направления в0у:
вк—--
У
2
агоБт
/
Бт ( воу |
+ во.
(4)
Полная ширина многолучевой ДН в вертикальной плоскости (за счет данного
эффекта):
1
ЛБполн — — -2
агсБт
п ( / )
бШ
шт
( в0У )
Л /
- агсБт -
у
Ео
\
шах
(/)
бШ
(в«>)
+Ае.
(5)
где Ае - номинальная ширина одного луча приемо-передающей ДН по вертикали.
Таким образом, частотная дисперсия АФАР на прием ведет к отклонению результирующих приемо-передающих лучей ДН (каждый из которых реализует работу в своем азимутальном направлении на своей частоте /) от заданного угла визирования е - как в антеннах с частотным сканированием [17]. На рисунке 1 представлена конфигурация приёмо-передающей многолучевой ДН в вертикальной плоскости при номинальном отклонении луча ДН на угол -15°, полученная по
результатам моделирования.
Амплитудная приемо-передающая ДН
О.Э
о.а 0.7 0.6 0.5 0.4 0.3 0.2 0.1
О
-1
Угол в вертикальной плоскости, град
Рисунок 1 - Многолучевая приемо-передающая ДН в вертикальной плоскости при номинальном отклонении лучей на угол -15° в вертикальной плоскости
Результирующая приемо-передающая ДН каждого луча получает дополнительное отклонение по вертикали и дополнительные амплитудные потери, вызванные несоответствием частоты настройки АФАР частотам принимаемых эхо-сигналов. Рисунок 2 иллюстрирует зависимость ошибки нацеливания приемопередающих лучей ДН, вызванных частотной дисперсией АФАР на прием, как отношение отклонения максимума приемо-передающей ДН в вертикальной плоскости по отношению к ширине ДН в вертикальной плоскости (ву - ву0)/Ае.
9*
Номинальное значение отклонения луча, град.
Рисунок 2 - Ошибка нацеливания (ву - ву0)/Ае приёмо-передающей ДН по вертикали, вызванного приемом в АФАР сигнала с частотой, отличной от частоты настройки фазовращателей АФАР ^ = 9600 МГц
Ошибка нацеливания лучей в многолучевых режимах съемки приводит к существенному искажению геометрии визирования при реализации многолучевой радиолокационной съемки. Проиллюстрируем влияние данного эффекта с использованием имитационной модели РСА, реализующей комплексное суммирование монохроматического эхо-сигнала по апертуре АФАР РСА порожденного отражением каждой точкой на поверхности Земли (в пределах границ моделирования данной поверхности) зондирующего суб-импульса [15].
На рисунке 3 представлено геометрическое расположение кадров на поверхности Земли (в координатах путевая-горизонтальная дальность) при реализации многолучевого прожекторного режима съемки для системы РСА с параметрами, представленными в таблице 1 [Ошибка! Источник ссылки не
найден.]. Кадры расположены на ближнем краю полосы обзора (отклонение луча ДН по вертикали -15°).
Таблица 1 - Характеристики расчетной системы РСА
Параметры Значение
Центральная частота, ^ Х-диапазон 9600 МГц
Максимальное число формируемых каналов обработки на прием, К 4
Номиналы несущих частот суб-импульсов, 9367,5 МГц, 9522,5 МГц, 9677,5 МГц, 9832,5 МГц
Угол крена антенны, екр 35°
Высота орбиты, Н 650 км
Плоскость антенны Горизонтальная Вертикальная
Размер апертуры АФАР, Бх х Оу 4,032 м 1,6 м
Количество элементов решетки Ых х Ыу 9 64
Шаг между элементами решетки Их х Иу 0,448 м 0,025 м
Сектор сканирования лучом ДН, Аасшн х Лвскан ±1,0° ±15°
Ширина двусторонней ДН, Ла х Лв 0,30° 0,77°
Рисунок 3 - Расположение кадров на поверхности Земли, образованные К = 4 антенными лучами в режиме многолучевой прожекторной съемки (Синий - луч 1 (9367,5 МГц); оранжевый - луч 2 (9522,5 МГц); желтый - луч 3 (9677,5 МГц);
фиолетовый - луч 4 (9832,5 МГц)).
Эффективная ширина многолучевой ДН в вертикальной плоскости, в пределах которой формируется прямоугольный участок местности:
Л&эфф — Лвном —
arcsm
Fo
n ( f )
sin
min
( V )| —
arcsin
Fo
max
(f)
sin
(в°. v)
(6)
В случае, представленном на рисунке 3 поперечная протяженность полосы съемки сократится с Lcz = 9,6 км до Lcz = 4,5 км (диапазон горизонтальных
дальностей 236,0 км ... 240,5 км). В этом случае эффективная полоса съемки определяется следующим соотношением:
Г , , \
^сэфф z (е) Ь
8+
. 2 у
Аеэфф \ ( Аеэфф
- Ь
8+-
V
(7)
где Ь - горизонтальная дальность в соответствии с выражением [15]:
Ь(е) = (Я3 + Акпс ) • (агсзт(К3 • зт(е)) - е), (8)
где Яя - радиус Земли, ЛИПС - высота рельефа визируемого участка местности, К - коэффициент сферичности с учетом рельефа и изменений высоты КА:
+ Н + АН
3 = ^ +АНПС ' (9)
в - угол визирования, Лвэфф - эффективная ширина ДН в соответствии с выражением 6 вычисляемое для угла отклонения 00у, реализующего нацеливание ДН на угол визирования в:
в0у = 8 - 8кр , (10)
где вкр - угол крена, под которым ориентирована апертура АФАР РСА.
Для компенсации указанного сокращения полосы съемки возможно снизить уровень, по которому отсчитывается ширина ДН антенны по вертикали - то есть допустить дополнительные энергетические потери на краю полосы съемки в пользу расширения ширин луча ДН. Данная мера применима при реализации конкретного режима съемки с одновременным контролем энергетических и геометрических характеристик.
Амплитудные потери, вызванные данным частотной дисперсией АФАР на прием, могут быть оценены по форме приемо-передающей ДН АФАР по вертикали:
« Пу (вгку - воу ) • Fry (&ку - воу ), (11)
где ¥1у(6) и ¥ту(6) - ДН АФАР, соответственно, на излучение и прием в вертикальной плоскости (без отклонения), 6ук - угловое положение максимума приемо-передающей ДН, определяемое по соотношению (4), в0у - номинальное значение отклонения ДН по вертикали.
Таким образом, частотная дисперсия АФАР на прием, при использовании многолучевых режимов съемки, приводит к сокращению визируемого участка местности по горизонтальной дальности в многолучевом прожекторном режиме съемки и дополнительным энергетическим потерям.
Эффект частотной дисперсии в режиме многолучевой сканирующей съемки
Многолучевая сканирующая съемка основана на комбинации технологии многолучевого визирования и классического сканирующего режима БсаиБАК, алгоритмы реализации которой подробно рассмотрен в работе [4]. Согласно данным алгоритмам при визировании К = 4 антенными лучами М = 3 парциальных полос в Щ = 2 циклах визирования реализуется геометрия визирования, представленная в левой части рисунка 4, полученная по результатам моделирования. Центр полосы съемки визируется под углом в = 35° - что соответствует отклонению луча ДН в вертикальной плоскости 60у = 0). В правой части данного рисунка представлена геометрия при аналогичных параметрах съемки, реализуемая на ближнем краю полосы обзора при визировании центра полосы съемки под углом вц = 23° - что соответствует углу отклонения приемо-передающей ДН 60у = -12°.
Рисунок 4 - Моделирование расположения парциальных кадров при многолучевой сканирующей съемке. Слева - в середине полосы обзора (визирование без отклонения многолучевой ДН по вертикали, ец = 35°), справа - на ближнем краю полосы обзора (при отклонении многолучевой ДН по вертикали на угол 90у = 12°, ец = 23°) при М = 3, К = 4.
Частотная дисперсия АФАР на прием оказывает существенное влияние на геометрию радиолокационного визирования, и приводит к появлению пропусков при многолучевой сканирующей съемки. Компенсация данных пропусков возможно двумя способами.
Первый способ предполагает ведение специального чередования номиналов
несущих частот / всех К лучей при переходе на новую парциальную полосу,
реализующегося следующим образом. k-ый луч (где k = 1...K) при визировании m-ой парциальной полосы (где m = 1...M) для компенсации частотного эффекта частотной дисперсии на прием будет иметь номинал несущей частоты fKOMn = fn(k), где f = f, f, . f - массив номиналов частот данного многолучевого режима, n(k) -порядковый номер номинала несущей частоты в массиве fk, определяемый соотношением:
n(k)=(k -1 + m) mod K, (12)
где k - номер антенного луча, для которого определяется номинал несущей частоты, m - порядковый номер парциальной полосы съемки в которой осуществляется визирование, a mod b - операция взятия остатка от деления a на b, K - количество антенных лучей многолучевого визирования.
Положение лучей ДН во всех парциальных полосах определяется соотношением:
ak=Ja.M.Lt-iKlH. (13)
2 M + K
Полный азимутальные сектор, в котором осуществляется многолучевое визирование:
Аовиз=Аа •
2 • M • K M
( K +1)| ■ (14)
V М + К М + К у
На рисунке 2 представлена геометрия визирования при применении алгоритма
частотного чередования.
Рисунок 5 - Моделирование геометрии визирования М = 3 парциальных полос в многолучевом сканирующем режиме съемки с применением частотного чередования. Слева - визирование в К = 3 антенными лучами, справа - К = 3
антенными лучами.
Данный способ применим для визирования без пропусков только при соблюдении равенства М = К - что существенно ограничивает свободу применения многолучевой сканирующей съемки для повышения характеристик визирования в РСА. Однако не требует введения дополнительных отклонений лучей ДН по горизонтали.
Второй способ компенсации пропусков предполагает оперативное управление
положением многолучевой ДН в горизонтальной плоскости при перенацеливании на
новую парциальную полосу и их неизменном угловом положении в процессе синтеза апертуры. Для этого при визировании антенные лучи ориентируются в горизонтальных направлениях, определяемых следующим соотношением:
ак (ш)=Аа •
М 2• к-(К + 1) К 2 ^ М + К М + К
( М -1 т--
I 2
л Л
у
(15)
где Аа - ширина одного луча ДН в горизонтальной плоскости, к - номер антенного луча, для которого определяется угловое положение, М - количество парциальных полос визирования, К - количество антенных лучей многолучевого режима, т -номер парциальной полосы в которой осуществляется визирование.
Полный азимутальный сектора Аавиз, в котором осуществляется перенацеливание антенных лучей:
Аавиз=Аа ■
2 • М • К М + К
(16)
В результате компенсации пропусков при съемке на ближнем краю полосы обзора за счет азимутального нацеливания антенных лучей реализуется геометрия
визирования, представленная на рисунке 6.
Рисунок 6 - Геометрии визирования в многолучевом сканирующем режиме съемки с применением дополнительного отклонения антенных лучей по горизонтали. Слева -визирование в К = 4 антенными лучами в М = 3 парциальных полосах, справа - К =
3 антенными лучами в М = 5 парциальных полосах.
Увеличение числа визируемых полос М снижает влияние частотной дисперсии антенной системы при многолучевом визировании на поперечную
протяжённость реализуемого кадра, поскольку данный эффект сказывается только на первой и последней парциальных полосах визирования.
Таким образом, предложенные подходы позволяют снизить влияние эффекта частотной дисперсии АФАР на прием на реализуемые характеристики многолучевого визирования. Причем алгоритм частотного чередования задействует значительно меньший азимутальный сектор Лавиз, чем алгоритм оперативного управления положением многолучевой ДН по азимуту (рисунок 7), что может быть существенным, для систем РСА с ограниченными возможностями по отклонению луча ДН в горизонтальной плоскости.
Ширина задействуем ого углового секюра по горизонтали
1.5т
о
N 01
: 3 4 5 6 7 3 9 10
Количество парциальных полос
-Азим. нацелив.. К = 2
— Чередов. част.. К = 2 Азим. нацепив.. К = 3 Чередов. част., К = 3
-Азим. нацелив., К = 4
----Чередов. част.: К = 4
Рисунок 7 - Угловая ширина сектора в горизонтальной плоскости в котором осуществляется визирование в многолучевом сканирующем режиме при различных способах компенсации пропусков при радиолокационной съемке, вызванных частотной дисперсией аналоговой АФАР на прием, в зависимости от количества
парциальных полос визирования М.
Заключение
Применение многолучевых режимов съемки с частотным разделением приемных ДН в РСА с аналоговой АФАР сопряжено с проявлением сверхширокополосных эффектов, проявляющихся в дополнительных отклонениях лучей ДН в вертикальной плоскости. Данный эффект приводит к искажению геометрии визирования в многолучевых режимах съемки, выливающиеся в сокращение реализуемой протяженности визируемого участка местности по горизонтальной дальности в многолучевом прожекторном режиме съемки и появлению пропусков при съемке в многолучевом сканирующем съемке, а также к дополнительным энергетическим потерям.
Рассмотренные алгоритмы оценки искажений и методики компенсации данного эффекта позволяют частично нивелировать влияние данного эффекта на качество радиолокационной съемки. В режиме многолучевой сканирующей съемки применение алгоритма частотного чередования позволяет вести радиолокационную съемку без пропусков при соответствии числа формируемых антенных лучей количеству парциальных полос визирования без дополнительного расширения азимутального сектора визирования. Алгоритм оперативного управления положением многолучевой ДН по азимуту позволяет реализовать многолучевую сканирующую съемку без пропусков при произвольном количестве используемых антенных лучей и парциальных полос визирования за счет расширения азимутального сектора в котором осуществляется визирование.
Библиографический список
1. Krieger G., Gebert N. and Moreira A. Digital beamforming techniques for spaceborne radar remote sensing in Proc. EUSAR, Dresden, Germany, 2006, URL: http ://elib.dlr.de/43801/1/KrGeMo_EUSAR06_ 13Mar0 6.pdf
2. Gerhard Krieger, Nicolas Gebert, Alberto Moreira: Multidimensional Waveform Encoding: A New Digital Beamforming Technique for Synthetic Aperture Radar Remote Sensing // IEEE Transactions on Geoscience and Remote Sensing, January 2008, vol. 46, no. 1, pp. 31 - 46.
3. Marwan Younis, Felipe Queiroz de Almeida, Federica Bordoni, Paco Lopez-Dekker, Gerhard Krieger. Digital beamforming techniques for multi-channel synthetic aperture radar // IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 2016. URL: http://elib.dlr.de/104069
4. Булыгин М.Л., Внотченко С.Л., Коваленко А.И., Риман В.В., Чечина И.Н. Режим многоканальной съемки в многоапертурном космическом радиолокаторе с синтезированной апертурой // Успехи современной радиоэлектроники. 2015. № 5. C. 20 - 26.
5. Bulygin M., Kovalenko A., Riman V., Vnotchenko S. Multi-channel modes implementation in spaceborne SAR with digital active electronically scanned array, EUSAR 2016 // 11th European Conference on Synthetic Aperture Radar. Electronic Proceedings, Hamburg, Germany, 6-9 June 2016, pp. 315 - 318.
6. Pingping Huang, Wei Xu. A New Spaceborne Burst Synthetic Aperture Radar Imaging Mode for Wide Swath Coverage // Remote Sensing, 2014, vol. 6, no.1, pp. 801 -814.
7. Булыгин М.Л., Маркова А.С., Муллов К.Д. Реализация режима многоканальной сканирующей съемки в радиолокаторе космического базирования, // Труды МАИ. 2018. № 98. URL: http://trudymai.ru/published.php?ID=9043 8
8. Внотченко С.Л., Коваленко А.И., Риман В.В., Шишанов А.В. Технический облик многоапертурного космического радиолокатора с синтезированной апертурой на основе АФАР Х-диапазона // Всероссийские радиофизические научные чтения-конференции памяти Н.А. Арманда. Сборник докладов. (Муром, 28 июня - 1 июля 2010). - Муром: Полиграфический центр МИ ВлГУ, 2010, С. 91 - 95.
9. Moreira A., Prats-Iraola P., Younis M., Krieger G., Hajnsek I., Papathanassiou K.P. A tutorial on synthetic aperture radar // IEEE Geoscience and Remote Sensing Magazine, 2013, vol. 1, no. 1. pp. 6 - 43.
10. Kovalenko A., Riman V., Shishanov A., Vnotchenko S. Design of Prospective Spaceborne Multi-Aperture UWB Polarimetric High Perfor-mance SAR System // 4th Microwave and Radar Week MRW-2010, 11th International Radar Symposium, Vilnius, Lithuania, June 16-18, 2010, pp. 490 - 492.
11. Kovalenko A., Riman V., Shishanov A., Vnotchenko S. Architecture and Perfomance of the Spaceborne Multi-Aperture High-Resolution SAR System Based On Analog-Digital Active Array Antenna, EUSAR 2012 // 9th European Conference on Synthetic Aperture Radar. Electronic Proceedings, 23-26 April 2012, Nurnberg, Germany, pp. 429 - 432.
12. Нониашвили М.И., Крючков И.В., Лесников Г.А. и др. Обзор современных радиолокаторов с синтезированной апертурой космического базирования и анализ тенденций их развития // Инженерный журнал: наука и инновации. 2012. № 8. С. 10.
13. Potin P., Rosich B., Grimont P., Miranda N., Shurmer I., O'Connel A., Torres T., Krassenburg M. Sentinel-1 Mission Status. Proceedings of EUSAR 2016 // 11th European Conference on Synthetic Aperture Radar, June 6 - 9, 2016, Hamburg, Germany, pp. 59 -64.
14. Воскресенский Д.И., Гостюхин В.Л., Максимов В.М., Пономарев Л.И. Устройства СВЧ и антенны. - М.: Радиотехника, 2006. - 376 с.
15. Булыгин М.Л., Внотченко С.Л. Построение диаграмм слепых дальностей и надирных отражений радиолокатора с синтезированной апертурой в MATLAB // Труды МАИ. 2015. № 83. URL: http://trudymai.ru/published.php?ID=62290
16. Брызгалов А.П., Ковальчук И.В., Хныкин А.В., Шевела И.А., Юсупов Р.Г. Моделирование радиолокатора с синтезированной апертурой при решении задач его внутреннего и внешнего проектирования // Труды МАИ. 2011. № 43. URL: http: //trudymai .ru/published.php?ID=24734
17. Сучков А.В. Имитационное моделирование эффекта нормали в антенных решетках с частотным сканированием // Труды МАИ. 2017. № 92. URL: http://trudymai.ru/published.php?ID=76809
18. Gebert N., Krieger G., Moreira. A. Digital Beamforming on Receive: Techniques and Optimization Strategies for High-Resolution Wide-Swath SAR Imaging // IEEE Transactions on Aerospace and Electronic Systems, 2009, vol. 45, no. 2, pp. 564 - 592.
19. Krieger G., Younis M., Gebert N., Huber S., Bordoni F., Patyuchenko A., Moreira A. Advanced digital beamforming concepts for future SAR systems. Proceedings of 2010 // IEEE International Geoscience and Remote Sensing Symposium (IGARSS). 26 - 30 July
2010, Honolulu, Hawaii, available at: http://elib.dlr.de/64961A/igrass2010 dbf 09Dec09.pdf
20. Capece P. Active SAR Antennas: Design, Development, and Current Programs. Review Article // International Journal of Antennas and Propagation, Article ID 796064, Hindawi Publishing Corporation, 2009, vol. 2009, 11 p., available at: http://doi: 10.1155/2009/796064