СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ТУРБУЛЕНТНОЙ АТМОСФЕРЫ
METHOD OF DETERMINING PARAMETERS OF TURBULENT ATMOSPHERE

Изобретение относится к физике атмосферы и может быть использовано при определении таких ее параметров как структурная характеристика показателя преломления, параметр Штреля и радиус Фрида.
Известен способ определения внутреннего масштаба атмосферной турбулентности (радиус Фрида) [1], заключающийся в том, что регистрируют интенсивность прошедшего турбулентную атмосферу излучения, определяют дисперсию ее флуктуации, и после определения дисперсии флуктуации интенсивности определяют внутренний масштаб атмосферной турбулентности.
К недостаткам известного способа следует отнести значительную удаленность в пространстве приемной и передающей излучение аппаратуры (на расстояние порядка нескольких километров), что существенно ограничивает возможности оперативного изменения как направления, в котором проводятся измерения, так и протяженности исследуемой атмосферы.
Известен также способ [2], согласно которому пропускают импульсное лазерное излучение через турбулентную атмосферу, фильтруют излучение обратного рассеяния на мелкомасштабных неоднородностях, измеряют коэффициент усиления обратного рассеяния одиночных импульсов, с частотой следования, не превышающей частоты флуктуации турбулентной атмосферы. Затем определяют коэффициент усиления обратного рассеяния импульсного излучения усреднением коэффициентов обратного рассеяния одиночных импульсов и восстанавливают профиль структурной характеристики показателя преломления турбулентной атмосферы обратным преобразованием коэффициента усиления обратного рассеяния импульсного излучения.
Данный известный способ требует использования мощного лазерного источника излучения (в силу малого значения величины обратного рассеяния), что существенно ограничивает возможности его широкого применения, особенно на больших расстояниях (превышающих 1000 м).
Техническим результатом предполагаемого изобретения является проведение экспресс-анализа состояния турбулентной атмосферы в масштабе реального времени при одновременном включении в состав определяемых параметров турбулентной атмосферы помимо структурной характеристики показателя преломления радиуса Фрида и параметра Штреля.
Указанный технический результат достигается тем, что используют адаптивную оптическую систему, состоящую из телескопа, сопряженной с ним видеокамеры и компьютера, обрабатывающего получаемые с видеокамеры изображения точечного (для данной апертуры) объекта, регистрируют распределение интенсивности искаженного турбулентной атмосферой изображения, с его помощью определяют структурную характеристику показателя преломления, радиус Фрида и аппаратную функцию атмосферы между объектом и входной апертурой оптической системы, восстанавливают искаженное турбулентной атмосферой изображение, регистрируют распределение интенсивности восстановленного изображения и по полученным данным определяют параметр Штреля, как отношение соответствующих интенсивностей в центрах искаженного и неискаженного изображений точечного объекта.
Возможность достижения технического результата основывается на следующем. Атмосферные неоднородности показателя преломления можно разделить на две части по степени их влияния на изображение [3]. Крупномасштабные неоднородности (характерный размер которых больше, чем диаметр входной апертуры оптической системы) приводят к случайным смещениям изображения как целого (дрожанию). Мелкомасштабные неоднородности вызывают размытие мелких деталей изображения и ухудшение вследствие этого разрешающей способности. Разрешающая способность оптической системы из-за влияния атмосферных искажений не превосходит разрешающую способность оптической системы с апертурой, равной радиусу Фрида, величина которого вблизи подстилающей поверхности порядка нескольких сантиметров.
n
2
2
Как известно [4-6], структурная характеристика показателя преломления C  выражается через дисперсию дрожаний (средний квадрат углового смещения центра тяжести изображения точечного объекта) σследующим образом:
<mrow><msub><mi>C</mi><mi>n</mi></msub><msup><mrow /><mn>2</mn></msup><mo>=</mo><msup><mi>σ</mi><mn>2</mn></msup><mo>/</mo><mn>5,68</mn><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>R</mi><mo stretchy="false">)</mo></mrow><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mi>L</mi><mtext>                          (2)</mtext></mrow>
где 2R - диаметр приемной апертуры;
L - длина трассы наблюдения.
0
n
2
Радиус Фрида (r) выражается через C  согласно [5, 6] как:
<mrow><msub><mi>r</mi><mn>0</mn></msub><mo>=</mo><mn>1,68</mn><msup><mrow><mo stretchy="false">(</mo><msup><mi>k</mi><mn>2</mn></msup><msub><mi>C</mi><mi>n</mi></msub><msup><mrow /><mn>2</mn></msup><mi>L</mi><mo stretchy="false">)</mo></mrow><mrow><mo>−</mo><mn>3</mn><mo>/</mo><mn>5</mn></mrow></msup><mo>,</mo><mtext>                               (3)</mtext></mrow>
где k=2π/λ - волновой вектор;
λ - длина волны излучения;
Мелкомасштабное расплывание усредненного изображения при условии, что дрожание устранено, определяется оптической передаточной функцией, зависящей от радиуса Фрида и диаметра входной апертуры оптической системы [7, 8]:
<mrow><mi>H</mi><mo stretchy="false">(</mo><mi>Ω</mi><mo stretchy="false">)</mo><mo>=</mo><mi>exp</mi><mo>{</mo><mo>−</mo><mn>3,44</mn><msup><mrow><mo stretchy="false">(</mo><mi>λ</mi><mi>Ω</mi><mo>/</mo><msub><mi>r</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mrow><mn>5</mn><mo>/</mo><mn>3</mn></mrow></msup><mo stretchy="false">[</mo><mn>1</mn><mo>−</mo><mi>α</mi><msup><mrow><mo stretchy="false">(</mo><mi>Ω</mi><mo>/</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo stretchy="false">]</mo><mo>}</mo><mo>.</mo><mtext>          (4)</mtext></mrow>
Здесь Ω - угловая пространственная частота,
0
Ω=D/λ - частота обрезания спектра пространственных частот оптической системы, D - диаметр апертуры. Параметр α принимает значение 1 для «ближнего поля» (когда существенны только фазовые эффекты) и значение 0,5 для «дальнего поля» (применимое, когда одинаково существенны и амплитудные и фазовые искажения).
На чертеже представлена схема системы экспресс-анализа состояния атмосферы, поясняющая сущность изобретения.
n
0
2
Посредством компьютерной программы сначала регистрируют дрожание изображения, вычисляют его дисперсию, параметры C  и r, a также функцию H(Ω), затем корректируют дрожание, производят усреднение и цифровую фильтрацию с использованием выражения (4), устраняющую мелкомасштабное расплывание. В результате такой обработки получают исправленное (неискаженное) дифракционное изображение [9], и, сравнивая его с исходным, определяют параметр Штреля.
Литература
4
1. А.с. 1840633 СССР, МКИ G01W 1/00. Способ измерения внутреннего масштаба атмосферной турбулентности / П.А.Бакут, И.В.Безденежных, К.Н.Свиридов, Ю.П.Шумилов (СССР) - №3183005/28; заявлено 28.10.1987; опубл. 27.06.2007.
3
2. А.с. №1840481 СССР, МКИ G01S 17/95. Способ измерения структурной характеристики показателя преломления турбулентной атмосферы / П.А.Бакут А.Б.Александров, В.А.Логинов, В.П.Логинов, И.Н.Матвеев, Ю.П.Шумилов (СССР) - №2220623/09; заявлено 08.06.1977; опубл. 27.03.2007.
3. Кандидов В.П., Чесноков С.С., Шленов С.А. // Оптика атмосферы и океана. - 11, 522 (1998).
4. Щеглов П.В. Проблемы оптической астрономии - М.: Наука, 1980.
5. Татарский В.И. Распространение волн в турбулентной атмосфере - М., Наука, 1967.
6. Гурвич А.С., Кон А.И., Миронов В.Л., Хмелевцов С.С. Лазерное излучение в турбулентной атмосфере - М.: Наука, 1976.
7. Гудмен Дж. Статистическая оптика - М.: Мир, 1988.
8. Fried D.L. J. Opt. Soc. Am. // 56, 1372 (1966).
9. Аверин А.П., Пряничников B.C., Тяпин В.В. Квантовая электроника // №40, 5 (2010), С.418-420.