Точная информация об осадках и особенностях динамических процессов, которые происходят в атмосфере при их выпадении, представляет большой практический интерес в таких областях, как наземная и спутниковая радиосвязь, электроэнергетика и авиация. Она востребована при проведении научных исследований в области физики облаков и осадков, геофизики, метеорологии и климатологии, а также для мониторинга динамических процессов в атмосфере при решении задач, связанных с контролем и обеспечением заданных режимов полета различных летающих аппаратов в наиболее неблагоприятных погодных условиях, сопровождающихся осадками. Тем не менее, в настоящее время ученые и инженеры не располагают в нужном объеме точными и достоверными данными об особенностях динамических атмосферных процессов в условиях выпадения осадков и в частности дождей [53, 106, 107].
Арсенал технических средств, при помощи которых можно было бы получать информацию о динамике атмосферы в дождях, достаточно узок. Практически единственным средством, обеспечивающим оперативные измерения на больших площадях (результаты которых в первую очередь интересуют современного потребителя) является радиолокация. Однако развитие и широкое практическое применение радиофизических средств исследования атмосферы с осадками сдерживается отсутствием надежных и достоверных методик измерения, которые допускали бы сопоставление и однозначную геофизическую интерпретацию результатов, полученных разными авторами.
Первые работы, показавшие возможность использования доплеровских радиолокационных методов для измерений в осадках, были проведены в Англии [1, 108]. Вслед за этим Роджерс и Баттан [2, 89] предложили два независимых метода измерения вертикальных воздушных потоков в атмосфере при выпадении дождей. Метод, предложенный Роджерсом, основывался на использовании двух параметров, измеряемых радиолокационно: средней доплеровской частоты ^ и отражаемости Ъ.
Метод, предложенный Баттаном, основывался на определении минимального размера капель, регистрируемых радиолокатором с заданным потенциалом на заданном расстоянии. Разность между фактическим положением на оси доплеровских скоростей сигнала от таких капель и скоростью их падения в спокойной атмосфере позволяла оценить скорость вертикальных потоков. Подобный метод сильно зависел от технических параметров радиолокатора, воздействия турбулентных пульсаций на каплю и линейности приемника, что не позволяло надеяться на высокую точность измерений.
Средняя частота доплеровского спектра, используемая в методике Роджерса в меньшей степени зависела от технических параметров локатора и> турбулентности, поэтому дальнейшее развитие получил именно это метод.
Предполагалось, что средняя частота фактического доплеровского спектра под воздействием восходящего или нисходящего потока сдвигается по оси частот без существенного изменения формы спектра (это справедливо при равномерном увлечении капель потоком). Для экспоненциального распределения капель дождя по размерам, которое было предложено и экспериментально обосновано Маршалом и Пальмером [3], а также квадратичной аппроксимации зависимости Ганна-Кинцера [4] справедливо следующее выражение для «теоретической среденй скорости гравитационного падения капель дождя» в зависимости от радиолокационнои отражаемости — Разность между «теоретической средней скоростью» гравитационного падения капель Ув, и «фактической средней скоростью» падения которая получалась по данным фактических измерений средней частоты доплеровского спектра, считалась скоростью вертикального потока м>= F (FЛ-3.84Z1'/I4, где «и> - скорость вертикального потока, ^ - средняя частота доплеровского N спектра, полученного за время измерения, г = ]Г.06 — радиолокационная 1 отражаемость, N — концентрация капель в единице объема, Б — диаметр капли Дождя.
Развитие метода шло по пути увеличения числа информационных параметров. Предполагалось, что оа — ширина доплеровского спектра, может выступать одним из таких параметров, и его использование позволит повысить точность и достоверность измерений. Однако экспоненциальное распределение капель по размерам и квадратичная зависимость скорости падения от размера капли, которые использовал Роджерс, оказались слишком грубой моделью, справедливой лишь при значительном усреднении результатов измерений.
Выяснилось, что в рамках «классического подхода» к моделированию микроструктуры дождя и на периодах времени, которые соответствуют радиолокационным измерениям (минуты), возможно использовать лишь среднюю частоту доплеровского спектра и отражаемость Ъ, а известные попытки связать ^ с шириной спектра <ус1 и Ъ, не получили необходимого экспериментального подтверждения [5].
В 60-х и начале 70-х годов в Советском Союзе широкий комплекс совместных радиолокационных, наземных и самолетных контактных измерений микрофизических характеристик облачных систем и происходящих в них динамических процессов при выпадении дождей, был выполнен научными учреждениями гидрометеослужбы [64]. В ходе этих работ была доказана необходимость дальнейшего совершенствования методов обработки и анализа радиосигналов, полученных от дождей, а так же аппаратуры радиолокационных комплексов пусть даже путем их усложнения.
Было показано, что дистанционные измерения с использованием только радиолокационной отражаемости Ъ и «априорных моделей» микроструктуры дождя (неизменных на всем протяжении эксперимента), приводит к значительным ошибкам. Использование контактной информации, получаемой у поверхности земли, в большинстве случаев не приводит к уменьшению ошибок. Были также выявлены основные факторы, определяющие точность радиолокационных измерений в осадках: временная и пространственная динамика микроструктуры, особенности рассеяния на отдельных гидрометеорах и динамические атмосферные процессы (воздушные потоки и турбулентные пульсации).
Результаты работ показали первостепенное значение разработки «полностью дистанционных» методов измерений, но сложности, которые возникли на этом пути — не смотря на большие финансовые затраты и энергичные усилия ученых разных стран в течении последних 50-ти лет — не удалось разрешить окончательно. Разработка «полиостью дистанционных методов» требовала создания теоретической базы, способов проведения и достоверной геофизической интерпретации результатов радиолокационных измерений.
Основная трудность, которую следовало преодолеть на этом пути, заключалась в том, что вышеперечисленные «основные факторы, определяющие точность радиолокационных измерений в осадках», в то же время являлись факторами, которые в основном определяют характеристики спектра рассеянного поля и его интенсивности. Выделение воздействия одного из них на параметры радиоэхо и даже оценка его в отдельности от остальных, на базе существовавших в то время моделей формирования отраженного сигнала от осадков, оказалась достаточно сложной задачей, что затрудняло геофизическую интерпретацию результатов радиолокационных измерений и снижало точность получаемых оценок.
От успешного решения задачи о разделении воздействия различных геофизических и метеорологических факторов на спектры рассеянного поля и его интенсивности принципиально зависит возможность использования дополнительного «источника» информации, которым являются характерные особенности доплеровского спектра и спектра флуктуаций квадратичного детектора. В противном случае использование спектральной информации будет ограничено только средней частотой доплеровского спектра, что не позволиткак было показано выше — проводить интерпретацию радиолокационной информации с требуемой точностью и достоверностью.
Таким образом, современные задачи дистанционных измерений в дождях требуют использования математических моделей, которые позволяют корректно и точно разделять влияние динамики атмосферы и микроструктуры осадков на характеристики рассеянного поля и его интенсивности. В этом случае точность оценки микрофизических параметров осадков будет определять точность измерения динамических характеристик атмосферных процессов и наоборот: точность оценки динамических характеристик атмосферных процессов будет тесно связана с точностью измерения микрофизических параметров осадков.
В нашей стране основа для разработки таких моделей, была сформирована в 50-х годах прошлого столетия в теоретическом отделе ИРЭ АН СССР работами по рассеянию радиоволн на блуждающих неоднородностях, которые были выполнены Г. С. Гореликом и его сотрудниками A.B. Францисон и М. И. Родак [6−9]. Основной недостаток работ этих авторов состоял в том, что окончательные формулы, которые были предложены ими, не учитывали различную отражающую способность и скорость гравитационного падения отдельных рассеивателей, входящих в ансамбль облученных частиц и не позволяли учесть взаимосвязь между движением рассеивающих частиц и их отражаемостью.
Модели, позволяющие установить связь между параметрами движения отдельных частиц в рассеивающем объеме и параметрами сигнала на выходе PJTC, были разработаны в 1959;61 годах в Советском Союзе [10]. В этом же направлении работы проводились и во Франции [11]. Опираясь на работы Г. С. Горелика и его сотрудников, А. Г. Гореликом и Ю. В. Мельничуком была разработана достаточно строгая модель, описывающая рассеяние радиоволн на блуждающих неоднородностях [12]. Были получены зависимости, связывающие основные спектральные параметры рассеянного поля и его интенсивности с параметрами рассеивающих неоднородностей. При выводе этих зависимостей полагалось, что отражающая способность неоднородностей (которыми могли выступать частицы осадков) может быть функционально связана со скоростью их движеиия.
В дальнейшем удалось связать спектральные параметры рассеянного поля и его интенсивности с метеорологическими факторами, характеризующими динамическое состояние атмосферы в объеме зондирования. Было показано [10, 12], что спектр рассеянного поля в масштабе (К/2)Р повторяет по форме распределение проекций скоростей рассеивателей на направление зондирования, с учетом того вклада, который вносит каждый рассеиватель в отраженный сигнал. Спектр интенсивности рассеянного поля на выходе квадратичного детектора повторяет по форме (в масштабе (К/2)¥-) распределение разности проекции скоростей рассеивателей на направление зондирования с учетом того, какой вклад вносят различные рассеиватели с данным значением разности скоростей в отражаемость.
В ходе экспериментальной проверки полученных результатов был предложен и проверен метод разделения вклада мелкомасштабной турбулентности и гравитационных скоростей падеиия капель различного размера (имеющих распределение п (Б), где Б — диаметр капли) в спектр интенсивности радиоэха [13]. Для оценки параметров мелкомасштабной турбулентности, подчиняющейся условиям закона Обухова-Колмогорова, предлагалось использовать спектры интенсивности рассеянного поля от «вложенных объемов» пространства. Масштаб наибольшего из «вложенных объемов» должен был быть меньше внешнего масштаба инерционного интервала атмосферной турбулентности. С использованием этого метода удалось существенно повысить точность и надежность интерпретации данных, получаемых от метеорологического радиолокатора, а также оценить скорость диссипации турбулентной энергии в облаках различных форм, осадках в виде дождя и снега [12,14]. ,.
Полученные результаты позволяли вплотную приблизиться к решению задачи о разделении влияния микроструктуры осадков и динамических процессов в атмосфере на параметры рассеянного поля и его интенсивности. Однако в рамках классической модели микроструктуры дождя — как было отмечено выше — имелись определенные трудности с использованием спектральной информации. Они заключались в расхождении ширины «теоретически ожидаемых» и «фактически регистрируемых» спектров.
Поэтому в начале 70-х годов в ЦАО А. Г. Гореликом и В. Ф. Логуновым [15, 16], опираясь на ранее полученные результаты, был выполнен цикл теоретических работ, направленных на повышение точности радиолокационных измерений скорости вертикальных воздушных потоков, параметров мелкомасштабной турбулентности и распределения капель по размерам в дождях.
В ходе этих работ была предложена «усовершенствованная модель микроструктуры» дождя, которая включала в себя: трехпараметрическое гамма-распределение и уточненную аппроксимацию зависимости скорости падения капель дождя в атмосфере от их размера. «Усовершенствованная модель» предполагала большую гибкость при интерпретации радиолокационных данных. В итоге, для радиолокатора, работающего в режиме вертикального зондирования при выпадении дождей, был предложен новый метод определения скорости вертикальных воздушных потоков.
Для определения «теоретической средней скорости» гравитационного падения капель было получено параметрическое семейство графиков вида: = ?(1(^(2/Л), ст), где Я — интенсивность дождя, Ъ — радиолокационная отражаемость. В качестве параметра семейства графиков выступала ста — ширина спектра рассеянного поля.
Предполагалось, что аи фактически регистрируемых спектров будет предварительно скорректирована (уменьшена) на ширину «турбулентного размытия». Для определения величины подобной коррекции спектра мог быть использован метод «вложенных объемов», о котором уже упоминалось выше. Способ получения скорости вертикального потока м? при известной теоретической средней скорости свободного гравитационного падения Уя и средней частоте фактического доплеровского спектра ^ не отличался от предложенного Роджерсом.
К сожалению в работах [15, 16], область применимости методики не была определена с необходимой точностью. Для формирования «оперативной» и «полностью дистанционной» версии методики следовало исключить «контактный» и «долгоизмеряемый» параметр Я — интенсивность дождя, входящий в соотношение ogiZIR). Способ его исключения, предложенный в [16] предполагал использование априорной (не изменяемой в ходе эксперимента), зависимости между радиолокационной отражаемостью Ъ и интенсивностью дождя Я. Но такая зависимость справедлива только на длительных периодах усреднения и в рамках классической модели микроструктуры. Следовательно, применение «полностью дистанционного» варианта, предложенного в [16] может обеспечить требуемую достоверность и точность оценок лишь на длительных периодах усреднения [104, 105]. Достоверность и точность оценки геофизических параметов на малых периодах времени измерения с использованием методики [15, 16] требует специальной оценки, которая не была проведена и это не позволяет в полной мере использовать одно из основных преимуществ радиолокационных методов: возможность проведения измерений в больших объемах пространства на малых промежутках времени.
Экспериментальная проверка полученных результатов была выполнена на доплеровских радиолокаторах, работающих как в импульсном режиме (длина волны Л = 3 см и, А = 30 см) так и непрерывном режиме с ЛЧМ модуляцией (длина волны, А = 30 см) [17, 25]. Последний позволял исследовать микроструктуру осадков, протекающие динамические процессы и их пространственно-временную изменчивость в малых масштабах осадков.
В ходе этих экспериментальных работ было показано, что далеко не всегда «мгновенные» (полученные на коротких периодах усреднения) распределения капель по размерам можно описать гамма-распределением. Достаточно часто в дождях (особенно интенсивных) присутствуют сверхкрупные капли. Подобные же эффекты следует ожидать при зондировании облаков, снега, града.
Появление «сверхкрупных» рассеивающих частиц вызывает существенные отклонения формы распределения рассеивателей по размерам от формы теоретически ожидаемого трехпараметрического гамма-распределения. Следовательно, применение методики [15−16], априорно полагающей существование трехпараметрического гамма-распределения частиц по размерам в зондируемом объеме, требует разработки способов оперативного обнаружения «сверхкрупных» частиц и компенсации их влияния при интерпретации измерений.
Очевидно, что задача обнаружения сверхкрупных рассеивателей и компенсации их влияния увеличивает требования к производительности системы обработки радиолокационных данных. Однако, появление мощной вычислительной техники позволяет, используя ранее накопленный опыт, вновь вернуться к решению проблем, связанных с проведением точных радиолокационных измерений в облаках и осадках, в том числе с учетом влияния крупных капель.
Как результат завершения этих работ можно ожидать появления новых методик для определения в режиме вертикального радиолокационного зондирования различных микрофизических параметров дождя и динамических процессов в атмосфере, протекающих в широком интервале пространственно-временных масштабов [83].
С использованием подобных методик станет возможной проверка гипотезы Колмогорова-Обухова для свободной атмосферной турбулентности в условиях выпадения осадков различного типа и интенсивностиизучение тонкой структуры воздушных потоков (одновременно на различных высотах) в кучевой и мощно-кучевой облачности на различных стадиях её развитияизучение связи между структурой воздушных потоков, интенсивностью и типом выпадающих на землю осадков, атмосферным электричеством.
Цель предлагаемого исследования состояла в разработке и теоретическом обосновании методики радиолокационных измерений скорости и знака вертикальных воздушных потоков, а также параметров атмосферной турбулентности в дождях, которая была бы применима для доплеровского локатора, работающего в режиме вертикального зондирования в дециметровом, сантиметровом и миллиметровом диапазоне длин волн, с учетом вариации микроструктуры дождей на различных пространственно-временных масштабах. Предполагалось, что разрабатываемая методика должна быть применима без использования «контактной информации» о микроструктуре дождя у поверхности земли.
Для достижения поставленных целей потребовалось:
1. Провести изучение современных публикаций и сформировать требования к методике дистанционного измерения динамических параметров атмосферы, применимой в широком интервале пространственно-временных масштабов и длин волн зондирующего излучения без использования «контактной информации».
2. Провести анализ особенностей доплеровских спектров, полученных от осадков, и специфических особенностей распределения капель дождя по размерам, которые получены по данным контактных измерений на малых пространственно-временных масштабах.
3. Опираясь на результаты проведенных исследований предложить методику измерений динамических параметров атмосферы и микроструктуры осадков с учетом особенностей проведения измерений в широком диапазоне пространственно-временных масштабов и в различных диапазонах зондирующего излучения.
Для разработки рекомендаций п. З потребовалось рассмотреть целый ряд задач, решение которых, предлагаемое в настоящем исследовании, может иметь самостоятельную ценность. К ним можно отнести: разработку «динамической модели» стандартного дождясравнительную оценку точности различных способов измерения динамических параметров атмосферы и микроструктуры дождя на её основеуточнение понятия «эффективной радиолокационной отражаемости», что оказалось необходимым для расширения области применимости полученных результатов на коротковолновые диапазоны зондирующего излучения (характеризующиеся наличием условий рассеяния Ми на каплях дождя).
3.5. Выводы.
1. На основе анализа экспериментальных данных о спектральных измерениях в дождях показано, что одним из основных отклонений формы фактически регистрируемых спектров рассеянного поля в дождях от теоретически ожидаемых является немономодалыюсть экспериментальных спектров, и проведен качественный анализ возможных причин возникновения «немономодальных» спектров.
2. С использованием математических моделей, предлагаемых в настоящей работе, впервые установлено, что наличие локальных отклонений от закона гамма-распределения в области крупных капель при использовании только ширины и средней частоты доплеровского спектра может приводить к переоценке турбулентности на коротких периодах времени измерения в 2−3 раза.
3. На основе разработанных в настоящей работе математических моделей, показано, что в условиях интенсивной турбулентности исследование формы «релеевского» спектра в области частот, соответствующих крупным каплям, осложнено сглаживанием спектрального пика от крупных капель. Предложено в этих условиях дополнительно исследовать спектры рассеянного поля, полученные на длинах волн с рассеянием Ми на каплях дождя. Возможность достоверного определения характерных точек на спектрах рассеянного поля и локальных отклонений от закона гамма-распределения в области крупных капель в случае интенсивной турбулентности требует дальнейшего исследования.
4. С использованием результатов, полученных в настоящей работе и математических моделей из методики Горелика-Логунова [15, 16] предложен новый способ измерения скорости вертикальных воздушных потоков, в составе которого предполагается использование информации о ширине распределения капель по размерам вместо ширины доплеровского спектра или спектра интенсивности. Преимущество такого способа измерений состоит в том, что в случае достоверной оценки максимального размера капли в зондируемом объеме, предлагаемый алгоритм обработки данных бдует нечувствителен к воздействию атмосферной турбулентности.
5. На основе результатов детального исследования возможных причин формирования бимодальных доплеровских спектров с выбросом в области высоких доплеровских частот предложены длины волн, характеризующиеся рассеянием Ми (X, = 4.5 мм — 5 мм) которые наиболее предпочтительны для проведения измерений в дождях.
6. Расчетным путем установлен порядок количества «крупных капель» в зондируемом объеме, которые могут приводить к бимодальной форме доплеровских спектров наиболее часто наблюдаемой экспериментально.
Глава 4. Способы интерпретации Z-R зависимости.
4.1. Общие положения.
Исследование радиолокационной отражаемости Z в зависимости от интенсивности дождя Я, было одной из первых задач радиолокационной метеорологии. За прошедшее время накоплен массовый экспериментальный материал, наработана значительная библиография теоретических исследований. Подавляющее большинство исследований сделано в рамках классичесих подходов сформулированных на базе сильноусредненных данных. Это частично объясняет отсутствие надежных полностью дистанционных методов оценки интенсивности дождя И. на коротких периодах времени измерения. Как было показано выше, отсутствие дистанционных методов оценки интенсивности дождя Я является существенным препятствием в развитии радиолокационных методов измерения динамических параметров атмосферы в условиях выпадения дождей.
Классический подход к измерениям в дождях предполагал детерминированную связь между интенсивностью дождя Я и его радиолокационной отражаемостью Z. Результаты, полученные в ходе экспериментальных работ по проверке классической модели дождя, не подтвердили наличие детерминированной связи между К и Z: она явно носила вероятностный характер и соответствовала «классическим» представлениям только при сильном усреднении данных.
Основные идеи вероятностной модели, которая пришла на смену детерменированной, заключались в следующем: в любом дожде N0 и Г)0 (где N0 — концентрация капель в единице объема, Б0 — медианный диаметр капель) ведут себя как случайные, независимые1 переменные, следовательно, как Ъ, Отмстим, что «независимость» N0 и D0 (в строгом смысле) не доказана. Некоррелированность этих параметров была экспериментально подтверждена в ходе эксперимента TOGA COARE [60]. так и Я имеют склонность быть распределены по логонормальному закону. Когда это справедливо, то Ъ и II, соотносятся в среднем по степенному закону.
Чтобы корректно учитывать стохастический характер 7,-К. зависимости на коротких периодах времени усреднения помимо Z требуются дополнительные информационные параметры, очевидным источником которых являются спектры рассеянного поля или его интенсивности. Однако до сих пор, несмотря на значительное количество публикаций, единый подход к учету стохастического характера 2−11 зависимости в оперативной обстановке радиолокационных измерений не согласован.
Одной из причин его отсутствия является то, что в 70-е годы основное внимание исследователей было обращено к изучению иных вопросов: статистических свойств и микроструктуры осадков, выпадающих из различной облачности, разработке теории конденсации в дождевой облачности различных форм и физическим процессам, которые происходят в осадках на высоте падения. В это время иследования, касающиеся Z-R соотношения ведутся лишь силами отдельных ученых, среди которых Фудживара [61, 62] в Японии, а также Атлас и Чимела в США [63], Боровиков и Костарев [64] в СССР.
Последнее десятилетие в публикациях наблюдается возобновление интереса к исследованию микроструктуры дождя. Вновь стали обсуждаться проблемы, связанные с формированием физически обоснованных позиций для объяснения изменчивости коэффициентов, А и Ь в Ъ — II соотношении [65−71], а так же с их экспериментальной проверкой [60, 55]. Однако значительных успехов в этом направлении пока не получено.
В одной из наиболее цитируемых западных работ [71], заложена основа для разработки алгоритма определения параметров Ъ-К соотношения, а в работе [72] она доработана с учетом требований поляриметрических измерений. В [71] для расчета оптимальной концентрации капель по данным дистанционных измерений N0 предлагается использовать значение водности.
W. Использование этого показателя в практических приложениях затруднено отсутствие прямых методов его измерения. Поэтому, в идейном плане [71] мало чем отличается от работы [73], опубликованной на 20 лет раньше. Тем не менее, [71] содержит одно существенное новое положение: использование постоянного показателя степени b — 1.5 при интенсивности дождя в Z-R соотношении. Такое значение b получено в результате статистической обработки массового материала эксперимента TOGA GOARE. Физических оснований для выбора указанного значения показателя степени в [71] не представлено.
В этой главе на основе работ отечественных авторов будет представлено теоретическое обоснование использования постоянного показателя степени b = 2 при R. Будет показано, что такое значение показателя степени является физически обоснованным в случае мгновенных измерений, а также представлены результаты использования предлагаемой квадратичной регрессии с.
Рис.
4.1.
Взаимозависимость обобщенным мультипликативным коэффициентов, А и ь в соотношении г = коэффициентом" для моделирования АК" ' полУченная различными авторами для разных типов дождей [74]. соотношения на конечных периодах времени параметров R и Ъ,.
Ясное физическое содержание и.
I — интенсивные дожди, Zmax > 104 (I = измерения 15мм/час);
II — слабые радиоэхо при однородных осадках Z ~ 5мм6/мЗ;
III — ячейки малых размеров со слабым радиоэхопростота интерпретации выгодно IV — малое радиоэхо от слабых осадков при отличает предлагаемую методику в ряду аналогичных, опубликованных сухой атмосфереV — радиоэхо, полученные для верхней части облака. VIИзолированные радиоэхо в грозах. до настоящего времени методик интерпретации Z-R соотношения.
Говоря об экспериментальных работах, следует отметить, что наиболее полным инструментальным исследованием Z-R зависимости в настоящее время все еще остается работа Фудживары [61] в которой отмечалось слабое соответствие между, А и Ь, такое, что большим значениям, А соответствовали меньшие значения Ь. Причем, явной зависимости от типа осадков найдено не было.
В исследованиях Фудживары мультипликативный коэффициент, А принимал значения в диапазоне от, примерно, 60 до 1100, Ь от 1 до 2. В целом исследования Японских ученых не позволили продвинуться дальше качественной диаграммы изображенной на рис. 4.1. [74], что отражает сложность поставленной задачи использования Z-R соотношения для классификации дождей и интерпретации данных радиолокационных измерений.
4.2. Мультипликативный коэффициент, А на подинтервалах.
Работа [75] была одной из первых исследований, направленных формирование теоретических оценок коэффициентов в формуле Z = АЯЬ, связывающей интенсивность и радиолокационную отражаемость дождя. Деление интервала скоростей на три подинтервала не позволяли оценить вариацию мультипликативного коэффициента, А с достаточной точностью для дождей различной интенсивности. Недостатком цитируемой работы было также отсутствие рекомендаций для обработки экспериментальных данных, выходящих за теоретические границы.
В [75] предлагалось, что статистика распределения капель по размерам соответствует трехпараметрическому Гамма-распределению вида п (О) = С0Оте~ро, чтобы не сильно проиграть в точности использовалась кусочная аппроксимация зависимости скорости свободного гравитационного падения от диаметра. На первом отрезке (минимальные размеры капель) закон скорости падения был линейным у = а-,£>, на следующем отрезке диаметров — квадратичным у8=а2-Л5. На последнем отрезке диаметров (максимальные размеры капель), скорость принимала постоянное значение.
V = аъ.
Подставляя последовательно кусочные аппроксимации скорости в формулу водности дождя, можем получить теоретические выражения мультипликативного коэффициента для дождей различных интенсивностей (см. Таблицу 4.1).
Заключение
.
В настоящем исследовании впервые были получены следующие результаты:
1. Представлена и теоретически обоснована полностью дистанционная методика измерения динамических параметров атмосферы для доплеровского радиолокатора, работающего в режиме вертикального зондирования, которая применима в различных пространственно-временных масштабах и на различных длинах волн зондирующего излучения. В методике используется концепция «обобщенного мультипликативного коэффициента», который требует учета параметров формы доплеровского спектра. Форма доплеровского спектра используется также для повышения достоверности измерений.
2. Теоретически обоснована концепция «обобщенного мультипликативного коэффициента» для измерений Z — Я на малых («мгновенных») периодах времени. Получены окончательные формулы для «обобщенного мультипликативного коэффициента», использующие параметры формы доплеровского спектра.
3. Сформулированы основные положения «динамической модели дождя», и с её использованием предложен механизм применения «обобщенного мультипликативного коэффициента» на периодах времени измерения конечной длительности. Показано удовлетворительное совпадение теоретических и экспериментальных результатов. Полученные зависимости позволяют проводить более точную и достоверную интерпретацию дистанционных измерений интенсивности дождя в широком интервале временных периодов с использованием радиолокатора в режиме вертикального зондирования.
4. Предложен механизм применения концепции «обобщенного мультипликативного коэффициента» в условиях рассеяния Ми, что позволяет расширить диапазон длин волн зондирующего излучения при проведении измерений динамических параметров атмосферы в режиме вертикального зондирования. Получены окончательные расчетные формулы. С использованием «динамической модели дождя» показано удовлетворительное совпадение теоретических и экспериментальных результатов для условий рассеяния Ми на конечных периодах времени измерения.
5. Предложена концепция методики измерения динамических параметров атмосферы с использованием информации о ширине распределения капель по размерам. Получены основные формулы. Особенностью предлагаемой методики является простой способ исключения воздействия атмосферной турбулентности — как наиболее трудноизмеримого параметра — на результаты радарных оценок динамических параметров атмосферы. Методика требует дальнейших исследований в части способов достоверного определения верхней границы распределения капель по размерам в зондируемом объеме.
6. На основе анализа характера нижней и верхней границы диапазона контактах измерений у земли, предложен метод улучшения оценок среднего диаметра капель с использованием классических контактных инструментов, основывающийся на исключении статистического смещения.
Приведенные выше результаты основывались на следующих основных предпосылках полученных в ходе настоящего исследования и учитывают их:
1. С использованием массовых экспериментальных данных установлено, что наиболее частым отклонением распределения капель по размерам от формы гамма-распределения на коротких периодах измерения, является локальные отклонения формы распределения в области крупных капель. Наибольшее количество спектров с «крупнокапельными отклонениями» следует ожидать на периодах времени измерения, равном времени пересечения наиболее крупной каплей высоты зондируемого объема.
2. Расчетным путем установлено, что наличие крупных капель в зондируемом объеме может приводить к переоценке интенсивности турбулентности в 2−3 раза. Этот факт говорит о необходимости использования критериев применимости методик, априорно допускающих гамма-распределение капель по размерам в дождях. Для случая крупных капель такие критерии могут быть разработаны с использованием параметров формы доплеровских спектров.
3. Расчетным путем получены оценки количества «крупных капель» в зондируемом объеме, которые могут приводить к бимодальной форме доплеровских спектров наиболее часто наблюдаемой экспериментально. Установлено, что для крупнокапельных дождей и дождей с широким распределением — это отдельные капли во всем зондируемом объемедля мелкокапельных дождей с узким распределением — это отдельные капли в кубическом метре зондируемого объема.
4. На основе результатов анализа возможных причин формирования бимодальных доплеровских спектров с выбросом в области высоких доплеровских частот предложены оптимальные длины волн, характеризующиеся рассеянием Ми (к = 4.5 мм — 5 мм), которые наиболее предпочтительны для проведения измерений в дождях.
Замечания по реализации.
Замечания по точности. Использование параметра Я, который достоверно измеряется контактными средствами на периодах времени измерения порядка 10 мин. не позволяет перейти к полностью дистанционным исследованиям в рамках методики [10, 12, 13, 15, 16].
Оценка точности для [15, 16] без непосредственного измерения параметра Я, (но с использованием «классической зависимости», которая соблюдается лишь для сильно усреднённых данных) составляет — 0.2 м/с.
Принимая подход [16] к оценке точности на больших периодах времени усреднения (порядка 60 мин), можно ожидать, что при переходе к шестимннутным измерениям точность ухудшится (измениться доверительный интервал) не менее чем в л/Го раз. Для негаусового процесса, которым является микроструктура дождя, оценки будут хуже.
Переход к одноминутным измерениям расширит доверительный интервал, ухудшив теретнческую оценку точности в примерно 8 раз — для гаусового процесса. Можно считать, что ассимтернчный пегауссовский процесс будет иметь погрешность в 10 раз больше. Таким образом, на одноминутных интервалах измерения замена Я статическим выражением, полученным на основе «классических Z-R зависимостей» приводит к теоретической оценке погрешности порядка 2 м/с.
Использование подходов, рассмотренных в настоящей работе, позволяет полностью исключить методическую погрешность измерений, так как предполагает — с использваиием информации о форме доплеровского спектра — вести оперативный контроль и учет фактического Z-R соотношения, сложившегося на времени измерения.
Замечания по достоверности. Спектральные имзмерения, необходимые для реализации предлагаемых методов обработки данных, и в частности информация о форме спектра, позволяют дополнительно контролировать в оперативном режиме соблюдение гипотезы о соответствии формы распределения капель по размерам форме трехпараметрического гамма-распредения. Для оценки степени соответствия может использоваться любой из критериев соответствия, вычисление которого с учетом производительности подсистемы обработки данных будет возможно на требуемом времени измерения [114]. Вопрос о допустимом значении критерия соответствия в зависимости от требуемой точности измерений является предметом дальнейших экспериментальных исследований.
Выход значения критерия соответствия за границы допустимых значений будет означать, что фактически измеряемый процесс находится вне области определения базовой методики, использующей базовую гипотезу о наличии гамма-распределения капель по размерам в зондируемом объеме пространства осадков. В таком случае для получения желаемых оценок динамических параметров атмосферы потребуется использовании другой методики обработки данных или коррекция (предобработка) данных, полученных от радиолокатора с целью их использования в рамках базовой методики.
Простейшим и наиболее частым отклонением формы распределения капель дождя по размерам от формы гамма-распределения является локальное отклонение формы распределения в области крупных капель, рассмотренное в этой работе. Такое оклонене формы па пятисекундных измерениях может наблюдаться в 20%.30% измерений [17, 19, 25], но оно допускает коррекцию и использование базовой методики.
В целом применение статистических критериев соответствия позволит ввести в практику и контролировать в оперативном режиме не только оценки точности, но и достоверности измерений.
Возможная схема. Изложенные выше соображения определяют следующую перспективную схему построения системы обработки данных, получаемых от радиолокатора, работающего в режиме вертикального зондирования.
2: ш ш о х о м.
— 1 >. О т о. <ц о О. ш ш 2.
Рис. 5.1. Схема реализации методики измерений с оперативной оценкой точности и достоверности.
Предлагаемая схема является, по сути, схемой системы распознавания образов, которые определяются соответствующими критериями, применяемыми к сигналу, поступающему с радиолокатора. Совокупность критериев определяет выбор наиболее оптимальной методики обработки радиолокационных данных в рамках выбранной базовой модели.