CC BY
28
ФИЗИКА ЗЕМЛИ, АТМОСФЕРЫ И ГИДРОСФЕРЫ
Моделирование формы отраженного импульса радиоальтиметра
К. В. Показеев1^, А. С. Запевалов1, В. В. Пустовойтенко2,6
1 Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра кафедра физики моря и вод суши. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
2 Черноморский филиал МГУ имени М. В. Ломоносова.
Украина, 99001, Севастополь, ул. Героев Севастополя, д. 7.
E-mail: a [email protected], 6 [email protected] Статья поступила 20.02.2013, подписана в печать 07.06.2013.
Анализируются возможности и ограничения использования распределения Грама-Шарлье возвышений взволнованной морской поверхности для описания формы отраженного радиоимпульса при вертикальном зондировании. Показано, что для характерных значений асимметрии и эксцесса распределения возвышений морской поверхности расчеты на основе распределения Грама-Шарлье приводят к искажениям формы радиоимпульса, в том числе могут приводить к появлению отрицательных значений. Предложен новый подход, основанный на использовании комбинированной модели распределения возвышения морской поверхности, свободный от недостатков, связанных с распределения Грама-Шарлье.
Ключевые слова: космический аппарат, радиоальтиметр, морская поверхность, модель Брауна, форма отраженного импульса, распределение Грама-Шарлье.
УДК: 537.877+532.5.013. PACS: 41.20.Jb, 47.35.Bb.
Введение
При вертикальном радиозондировании с космического аппарата (КА) всю информацию об уровне морской поверхности вдоль его трассы и о ее локальных параметрах несет передний фронт отраженного радиоимпульса [1]. В настоящее время наряду с измерениями уровня моря [2] по данным альтиметрических измерений определяются еще два параметра морской поверхности: значимая высота волн [3] и уровень шероховатости морской поверхности, по которому рассчитывается скорость приводного ветра [4]. Исследуется возможность дистанционного определения асимметрии распределения возвышений взволнованной морской поверхности [5].
Установленные на КА радиоальтиметры работают на длинах волн от 2 до 6 см. Они имеют узкую диаграмму направленности, ширина которой не превышает 3.4° (исключением является альтиметр, установленный на КА «Епу18а1», с шириной диаграммы направленности 5.25°). Зондирование осуществляется импульсами длительностью порядка одной наносекунды [2].
При вертикальном зондировании морской поверхности отраженный радиосигнал формируется в результате зеркального отражения (этот физический механизм является основным при углах падения до 20° [6]. Форма отраженного импульса альтиметра описывается сверткой трех функций (модель Брауна [7]). Две функции определяются характеристиками альтиметра, установленного на космическом аппарате, третья функция определяется плотностью вероятностей возвышений морской поверхности [8]. Зависимость формы отраженного радиоимпульса от распределения возвышений морской поверхности позволяет решать обратную
задачу, восстанавливать параметры этого распределения [9].
Основной моделью, описывающей распределение возвышений морской поверхности в приложениях, связанных с дистанционным зондированием, является модель, построенная на основе рядов Грама-Шарлье [10, 11]. Ее достоинством является то, что коэффициенты ряда рассчитываются по статистическим моментам изучаемого параметра. Недостаток заключается в том, что в натурных экспериментах, как правило, определяются статистические моменты не старше четвертого порядка. В результате при построении модели плотности вероятностей могут быть использованы только первые пять (включая нулевой) членов ряда Грама-Шарлье, что приводит к искажениям на «крыльях» распределения, в том числе к появлениям отрицательных значений [12, 13].
Некорректное описание распределений возвышений морской поверхности в свою очередь приводит к искажениям расчетной формы отраженного радиоимпульса. Соответственно появляются ошибки в измерении расстояния от космического аппарата до морской поверхности и в определении ее характеристик [14]. В настоящей работе предложен новый подход к описанию формы отраженного радиоимпульса, основанный на использовании комбинированной модели плотности вероятностей возвышений морской поверхности.
1. Модель Брауна
При квазивертикальном зондировании средняя форма отраженного от морской поверхности радиоимпульса описывается моделью Брауна [7]
V(t)= x(t) * s(t) * q(t), (1)
где x(t) — форма импульса, отраженного от плоской поверхности; s(t) — форма зондирующего импульса; q(t) — функция, связанная с плотностью вероятностей высот точек зеркального отражения; символ * означает свертку; t — время. Модель Брауна была построена в рамках следующих предположений:
— отражающая поверхность содержит достаточно большое число независимых отражающих элементов;
— статистики возвышений поверхности являются постоянными в пределах освещаемой радаром поверхности;
— отражение является процессом без эффектов поляризации;
— нормированное сечение обратного рассеяния является постоянным для всей освещаемой площади, т. е. в пределах диаграммы направленности не зависит от угла падения;
— эффектом Доплера можно пренебречь.
Физический процесс формирования формы отраженного от морской поверхности короткого (наносекунд-ного) импульса рассматривался в ряде работ (см., например, [1, 2, 15]). Из этих работ, в частности, следует, что влиянием боковых лепестков в диаграмме направленности антенны альтиметра на формирование переднего фронта короткого (наносекундного) отраженного импульса, можно пренебречь.
Следует также отметить, что зеркальное отражение имеет место при выполнении условия, что длина поверхностной волны много больше длины радиоволны. Мелкомасштабная рябь, не удовлетворяющая этому условию, вызывает диффузное рассеяние, которое приводит к уменьшению мощности квазизеркального отраженного радиолокационного сигнала. Для того, чтобы учесть диффузное рассеяние коэффициент Френеля заменяют эффективным коэффициентом отражения [16]. Согласно [17] эффективный коэффициент отражения меняется в пределах от 0.45 до 0.50 при изменении скорости ветра от 1.5 до 15 м/с.
Будем рассматривать случай, когда форма зондирующего радиоимпульса является гауссовой:
s(t) =
1
лДПД
exp -
( 2Д.)
(2)
где параметр Бг определяет ширину радиоимпульса. Будем также полагать, что зондирование осуществляется при нулевом угле падения, и что диаграмма направленности описывается функцией Гаусса. В этом случае форма отраженного от плоской поверхности импульса имеет вид [8]
x(t) = a exp
ln 4c
cos2(0/2)h
H (t),
(3)
где a — амплитуда; c — скорость света; 0 — ширина луча антенны, определенная по уровню половинной мощности; h — высота орбиты КА; H(t) — единичная функция Хевисайда.
Если форма зондирующего импульса s(t) и форма отраженного от плоской поверхности импульса x(t) описываются соответственно выражениями (2) и (3), их свертка
I(t)= x(t) * s(t) =
2 exp
52 Д -«
1 - erf( (5 - D)
2
(4)
где 5 =
ln 4c cos2(0/2)h
; erf(x) = -jn J exp(-y2) dy — функция 0
ошибок.
Учитывая (1) и (4), получаем, свертку двух функций, где первая функция определяется характеристиками альтиметра, вторая — характеристиками морской поверхности
V«) = I(О * ?(*).
Функция ц{1) рассчитывается путем преобразования плотности вероятностей возвышений морской поверхности Р(п) с помощью соотношения, связывающего возвышение поверхности п и время I [8]:
t = n/(c/2).
(5)
Поскольку преобразование (5) является линейным, асимметрия и эксцесс распределений ц{1) и Р(п) численно равны между собой. Отметим, что понижение уровня морской поверхности приводит к увеличению времени возврата радиоимпульса альтиметра. Следовательно, нечетные статистические моменты распределения ц{1) должны иметь противоположный знак, чем моменты возвышения морской поверхности, положительно направленного вверх.
Далее при численном моделировании формы отраженного радиоимпульса будем полагать, что параметры функций (2)-(4) соответствуют значениям параметров альтиметра, установленного на КА «БеазаЬ-! ет = 1.6°, = 1.327 нс и Н = 8 • 105 м [2].
2. Модели плотности вероятностей возвышений морской поверхности
Первоначально форма отраженного радиоимпульса рассчитывалась в рамках предположения, что распределение возвышений морской поверхности описывается распределением Гаусса. В настоящее время основной моделью плотности вероятностей возвышений морской поверхности, используемой при анализе данных альти-метрических измерений, является модель, построенная на основе рядов Грама-Шарлье [5, 9, 10]. Данная модель позволяет учитывать нелинейные эффекты в поле морских поверхностных волн, приводящие к отклонению распределения возвышений поверхности от распределения Гаусса.
Как правило, при моделировании плотности вероятностей возвышений морской поверхности используется эджвортова форма типа А рядов Грама-Шарлье [18]
Poc(n) = ЕaH [ | -¡= exp (-2L
1
n=0
V2n
(А V 2 W
(6)
где Нп — ортогональные полиномы Чебышева-Эрмита п -го порядка; ап — коэффициенты ряда (6). Предполагается, что среднее значение возвышения п = 0.
На практике расчеты плотности вероятностей возвышений морской поверхности ограничивают первыми членами ряда Грама-Шарлье. При анализе данных альтиметрических измерений используются несколько
a
t
форм модели РЪс(п), которые после введения норми- лей (7)-(9), удовлетворяет условию
\Щ <2.5.
ровки ^ = цш п2 можно записать в виде
Рпс(п) =
Р<(сС(п) =
P<G3C (п) =
exp (-¥)
2пп2 exp (- ¥ )
А
1 + g H (ч)
2щ2 exp(-|)
4 F
1 + 6 Яз ® + 24 Н4
(7)
(8)
4 F 42
1 + 4Н3 (v) + 24 Н4 + 4 Не fô)
(9)
где Л и £ — асимметрия и эксцесс распределения возвышений морской поверхности. Модель (9) была получена в рамках анализа нелинейных эффектов в поле поверхностных волн, обусловленных слабыми межволновыми взаимодействиями [18].
Известным недостатком моделей Грама-Шарлье является то, что они позволяют описать плотность вероятностей только в ограниченной области изменения исследуемого параметра [13]. Для возвышений морской поверхности область, где возможно применение моде-
Следствием искажения моделью Грама-Шарлье распределения возвышений морской поверхности за пределами области (10) является искажение рассчитанной с ее помощью формы отраженного радиоимпульса Vqç(î) (i = 1,2,3). При некоторых значениях асимметрии А и эксцесса F в моделях ), как видно из рис. 1, появляются отрицательные значения. В верхней части рис. 1 показаны общая форма отраженного импульса, рассчитанная по моделям Р<о£(п), а также сами модели Pq£(^) . В нижней части рис. 1 приведены
фрагменты функций V(C(t) и Pq£(^) в областях, где имеют место отрицательные значения. Приведенные графики наглядно показывают ограничения возможности использовать модели Грама-Шарлье для описания плотности вероятностей возвышений морской поверхности при моделировании отраженного радиоимпульса.
Проблема ограничения области изменения исследуемого параметра, в которой возможно применение моделей Грама-Шарлье, существует также при описании других характеристик поля поверхностных волн, в частности при описании распределения уклонов морской поверхности [19]. Для ее решения в работе [20] предложено использовать комбинированную модель. В об-
vgc(Л)
-0.01
Рис. 1. Формы радиоимпульсов V(C(t), отраженных от морской поверхности, рассчитанные по модели Грама-Шарлье Рос© : Расчеты проведены для поля поверхностных волн, имеющего следующие характеристики:
Hs = 5 м, А = 0.3 и F = -0.3
ласти (10) комбинированная модель близка к модели Грама-Шарлье, за ее пределами приближается к распределению Гаусса.
В общем виде комбинированную модель можно представить как
1
2пп2
Pc (n) = exp - ^ 1 + F (fi) £ anHn (fi)
n=3
где Р — фильтр. Фильтр Р должен быть близок к единице в области (10) и стремиться к нулю за пределами этой области. В качестве фильтра используется двухпараметрическая функция
F (n) = exp - (|n|/d)
(11)
где параметр й определяет область, внутри которой Р (п) « 1, параметр п определяет скорость, с которой функция Р стремится к нулю за пределами этой области. При моделировании распределений характеристик морской поверхности оптимальным является выбор п = 3.5 [21].
Как было отмечено выше, в некоторых ситуациях в модели Грама-Шарлье появляются отрицательные значения при относительно небольших величинах |г?|. Так, например, при А = 0.17 и Е = -0.18 отрицательные значения в модели Грама-Шарлье появляются в области г] < -3.1 (рис. 2). Отметим, что комбинированная модель может также иметь отрицательные значения, если неправильно выбран параметр й. В частности, при указанных значениях асимметрии и эксцесса в модели Рс (п) появляются отрицательные значения, если й = 4.5 и й = 4. При уменьшении й отрицательные значения в комбинированной модели не появляются. При й < 3 в области > 3 комбинированная модель везде положительна.
P(rj) 0.005
-0.001
Рис. 2. Плотности вероятностей возвышений взволнованной поверхности Р(г}): сплошные кривые — комбинированная модель (1 — й = 3, 2 — й = 3.5, 3 — й = 4, 4 — й = 4.5); пунктир — модель Грама-Шарлье; штриховая линия — распределение Гаусса. Расчет при А = 0.17 и Е = -0.18
Введем параметр 'пР, соответствующий частоте среза в частотных фильтрах. Аналогично принятому в ра-
диотехнике для частотных фильтров условию его можно определить из соотношения F (fi = fiF) = 0.707. При n = 3.5 связь параметров fiF и d можно аппроксимировать выражением fiF = 0.74d. Поскольку модель Грама-Шарлье справедлива только в пределах области (10), значение параметра d естественно принять равным 3, что соответствует fiF « 2.5.
Комбинированная модель должна удовлетворять условию PC (fi) >0 при всех наблюдаемых для распределения возвышений морской поверхности значениях асимметрии и эксцесса. По измерениям в натурных условиях получены следующие результаты. При значимой высоте волн Hs <4 м по данным измерений, выполненных на Черном море, оценки асимметрии A и эксцесса E меняются соответственно в диапазонах -0.2 < A <0.4, -0.4 < E <0.9 [3]. При значимой высоте волн Hs > 4.5 м, по данным измерений, выполненных в Северном море, оценки асимметрии A и эксцесса E меняются соответственно в диапазонах -0.05 < A <0.4 и -0.4 < E <0.4 [22]. При сильных штормах ( Hs > 10 м) в ситуациях, когда наблюдались анормальные волны, т. е. волны, высота которых более чем в два раза превышала значимую высоту, изменения статистических моментов наблюдались в пределах 0.20 < A <0.51, 0.23 < E < 1.53 [23].
Расчеты показывают, что условие PC (n) > 0 выполняется при любых значениях асимметрии и эксцесса, лежащих внутри наблюдаемых в натурных условиях диапазонов их изменений.
3. Расчет формы отраженного радиоимпульса на основе комбинированной модели PC (n)
Если плотность вероятностей возвышений морской поверхности P(n) и соответственно функция q(t) описаны некорректно в области, соответствующей высоким гребням волн, то возникают искажения переднего фронта отраженного радиоимпульса. В частности, если имеют место отрицательные значения P(n) при n>0, то могут иметь место отрицательные значения в модели V(t).
Некорректное описание P (n) в области впадин волн практически не влияет на модель переднего фронта отраженного радиоимпульса. Связано это с тем, что когда передний фронт зондирующего радиоимпульса достигает впадин, основной вклад в отраженный сигнал дают участки поверхности, для которых выполняется условие (10). На этой стадии изменение интенсивности отраженного радиосигнала определяется изменением площади облучаемой поверхности [24].
Рассмотрим ситуацию, когда распределение возвышений морской поверхности описывается комбинированной моделью PC (n). Функция qC (t) имеет вид
qc (t)=Tibexp (- 2d) x
1 + F T,
ш t (-1)nan4 TD)
где D = (4/с2)п2. Индекс С показывает, что данная характеристика рассчитана для комбинированной модели.
Формы отраженных радиоимпульсов, построенных с помощью комбинированной модели, представлены на
n
X
0.20
0.15 -
0.10
0.05
Рис. 3. Формы радиоимпульсов V, отраженных от морской поверхности, рассчитанные по комбинированной модели РС^(г}): Расчеты проведены для поля поверхностных волн, имеющего следующие характеристики:
Н = 5 м, Л = 0.3 и Е = -0.3
рис. 3. Расчеты выполнены при тех же значениях асимметрии А и эксцесса F, при которых рассчитывались формы отраженных радиоимпульсов по моделям Грама-Шарлье, представленные на рис. 1. Сравнение рис. 1 и 3 показывает, что в области, где при расчетах по моделям Грама-Шарлье имели место отрицательные значения интенсивности отраженного радиоимпульса, при расчетах по комбинированным моделям они отсутствуют.
Заключение
Проведен анализ возможностей использования распределения Грама-Шарлье возвышений взволнованной морской поверхности для описания формы отраженного радиоимпульса при вертикальном зондировании. Показано, что ограничение области возвышений морской поверхности, в пределах которой справедливы модели
Грама-Шарлье Р<зС)(п) (п = 1,2,3), приводит к искажениям формы отраженного радиоимпульса. Наглядным проявлением этих искажений является появление отрицательных значений в расчетной форме отраженного радиоимпульса.
Предложен новый подход, основанный на использовании комбинированной модели распределения возвышения морской поверхности. В области (10) комбинированная модель близка к модели Грама-Шарлье, за ее
пределами приближается к распределению Гаусса. Расчеты с помощью комбинированной модели позволяют корректно описать отраженный радиоимпульс во всем диапазоне величин асимметрии и эксцесса распределения возвышений морской поверхности, наблюдающихся в натурных условиях.
Список литературы
1. Callahan P.S., Rodriguez E. // Marine Geodesy. 2004. 27. P. 391.
2. Пустовойтенко В.В., Запевалов А.С. Оперативная океанография. Спутниковая альтиметрия: современное состояние, перспективы и проблемы. Севастополь, 2012.
3. Queffeulou P. // Marine Geodesy. 2004. 27. P. 495.
4. Cheng Y., Xu Q., Liu Y. et al. // Intern. J. Remote Sensing. 2008. 29, N 19. P. 5699.
5. Запевалов А.С., Пустовойтенко В.В. // Исследования Земли из космоса. 2012. № 5. С. 12.
6. Valenzuela G. // Boundary Layer Meteorology. 1978. 13, N 1-4. P. 61.
7. Brown G.S. // IEEE Trans. Antennas Propagat. 1977. AP-25. P. 67.
8. Hayne G.S. // IEEE Trans Antennas Propagat. 1980. AP-28. P. 687.
9. Gomez-Enri J., Gommenginger C.P., Srokosz M.A., Chal-lenor P.G. // J. Atmos. Ocean. Technol. 2007. 24. P. 1102.
10. Rodriguez E., Chapman B. // J. Geoph. Res. 1989 94, N C7, P. 9761.
11. Gomez-Enri J., Srokosz M.A., Gommenginger C.P. et al. // Marine Geodesy. 2007 30. P. 217.
12. Кендалл М.Дж., Стьюарт А. Теория распределений. М., 1966.
13. Запевалов А.С., Большаков А.Н., Смолов В.Е. // Океанология. 2011. 51, № 3. С. 432.
14. Запевалов А.С. // Изв. РАН. Сер. Физика атмосферы и океана. 2012. 48, № 2. С. 224.
15. Mailhes C., Tourneret J.-Y., Severini J., Thibaut P. // Eur. Signal and Image Processing Conference (EUSIPCO 2008). Lausanne (Switzerland), 2008.
16. Wu S.T., Fung A.K. // J. Geophys. Res. 1972. 77, N 30. P. 5917.
17. Freilich M.H., Vanhoff B.A. // J. Atmos. Ocean. Technol. 2003. 20, N 4. P. 549.
18. Huang N.E., Long S.R. // J. Fluid Mech. 1980. 101, N 1. P. 179.
19. Cox C, Munk W. // J. Mar. Res. 1954. 13, N 2. P. 198.
20. Запевалов А.С., Пустовойтенко В.В. // Изв. вузов. Радиофизика. 2010. 53, № 2. С. 110.
21. Запевалов А.С. Статистические модели взволнованной морской поверхности для задач дистанционного зондирования. SaarbrUcken, 2012.
22. Jha A.K. Winterstein S.R. // Proc. 19th Intern. Offshore Mech. Arctic Eng. Symp. ASME. Paper No. OMAE 00-6125.
23. Guedes Soares C., Cherneva Z., Antao E.M. // Applied Ocean Res. 2003. 25. P. 337.
24. Barrick, D.E., Lipa B.J. // Satellite Oceanic Remote Sensing. Adv. in Geophys. 1985. 27. P. 61.
Simulation radar altimeter return waveform
K. V. Pokazeev1,a, A.S. Zapevalov1, V. V. Pustovoytenko2 b
1 Department of Marine and Inland Water Physics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
2 Black Sea Branch, M. V. Lomonosov Moscow State University. Sevastopol 99001, Ukraine. E-mail: a [email protected], b [email protected].
The possibilities and limits of the Gram-Charlier distribution rough sea surface elevations to describe the shape of the return waveform at vertical sounding. It is shown that for typical values of skewness and kurtosis distribution of sea surface elevation calculations based on the distribution of the Gram-Charlier distort the return waveform, including the appearance negative values. A new approach based on the use of a combined model of distribution of elevation of the sea surface that is free of the disadvantages associated with the Gram-Charlier distribution.
Keywords: satellite, radioaltimetr, sea surface, model of Brown, the Gram-Charlier distribution, return waveform. PACS: 41.20.Jb, 47.35.Bb. Received 7 June 2013.
English version: Moscow University Physics Bulletin 5(2013).
Сведения об авторах
1. Показеев Константин Васильевич — доктор физ.-мат. наук, профессор, профессор; тел.: (495) 939-16-77; e-mail: [email protected].
2. Запевалов Александр Сергеевич - доктор физ.-мат. наук, ст. науч. сотрудник, вед. науч. сотрудник; тел.: (0692) 41-05-73, e-mail: [email protected].
3. Пустовойтенко Владимир Владимирович - канд. физ.-мат. наук, ст. науч. сотрудник, вед. науч. сотрудник; тел.: (0692) 24-01-83, e-mail: [email protected].
