Метрологические особенности определения вертикальных градиентов температуры при CTD–зондировании в море

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Метрологические особенности определения вертикальных градиентов температуры при CTD-зондировании в море

Степанюк Иван Антонович

доктор физико-математических наук, профессор

Российский государственный

гидрометеорологический университет

Аннотация

Рассматриваются погрешности определения вертикальных градиентов температуры в море при зондировании, обусловленные рассогласованием динамических свойств каналов измерения температуры и давления используемых зондов. Анализируются возможности коррекции получаемых искажений.

Ключевые слова: зондирование морской среды; вертикальные градиенты температуры; метрологические особенности; рассогласование динамических свойств измерительных каналов; возможности коррекции

Abstract

The errors of determination of the vertical temperature gradients in the sea during sensing are considered in the article. These errors result from mismatch of the dynamic properties of the sonde's channels of pressure and temperature measurements. The possibilities of correction of the obtained errors are analyzed.

Keywords: sensing of the marine environment; vertical temperature gradients; metrological features; mismatch of the dynamic properties of the measurement channels; the possibility of correction

Введение

Вертикальные градиенты температуры на различных горизонтах в море необходимы для последующих океанологических расчетов в задачах термодинамики океана (например, [1]). При этом очень важна «привязка» каждого значения к определенному горизонту, поскольку такой горизонт определяет положение границ квазиоднородных слоев, например, при TS-анализе водных масс. Соответственно, становится необходимым оценивать погрешности определения градиентов, вызванных средствами и методиками CTD-зондирования.

При CTD-зондировании (conductivity, temperature, depth - электрическая проводимость, температура, давление) используются первичные измерительные преобразователи (датчики) с различными значениями постоянных времени. Это обусловливает рассогласование («нестыковку по глубине») измеряемых характеристик, что, в свою очередь, приводит к искажениям результатов последующего анализа.

Ранее в нашей работе [3] рассматривались особенности рассогласования каналов температуры и электрической проводимости при условии, что датчик температуры является инерционным (постоянная времени >0), а датчик проводимости безинерционен, что для большинства зондов характерно. Оценивались погрешности конечных океанологических расчетов условной плотности и условного удельного объема. Было показано, что пренебрежение таким рассогласованием постоянных времени может приводить к неприемлемым погрешностям при расчетах этих характеристик.

В данной работе анализируется ситуация рассогласования каналов измерения температуры и давления, приводящая к значительным погрешностям определения вертикальных градиентов. Рассматриваются перспективы совершенствования CTD-зондов.

Погрешности определения вертикальных градиентов температуры, обусловленные динамическими свойствами каналов

Рассмотрим случай измерений градиентов температуры при зондировании с постоянной скоростью Vз=const. Градиент вычисляется по выражению

(1)

метрологический градиент температура зондирование

где zi=0,5(z2+z3). Обозначения соответствуют рис. 1.

Будем считать, что определения глубин являются практически безинерционными, а измерения температуры проводятся по каналу с динамическими свойствами системы первого порядка. Для этого случая типичное дифференциальное уравнение, описывающее динамические свойства канала температуры преобразуется к виду

(2)

где Tпр - измеренноe (приборное) значение; Tис - истинное значение; Vз - скорость зондирования; фe - постоянная времени.

Выражение (2) получается из типичного уравнения для случая изменчивости по времени путем использования известной взаимосвязи:

(3)

Из выражений (1) и (2) легко получить

 (4)

Таким образом, интересующая наблюдателя истинная разность значений температуры на двух соседних горизонтах соответствует измеренной (приборной) с некоторой знакопеременной добавкой (±), зависящей не только от Vз и фe, но и от различий градиентов температуры на этих горизонтах.

На рис. 1 проиллюстрированы искажения, вносимые динамическими свойствами канала измерений температуры в результаты определений разности Tис.

Как следует из рис. 1 разность Tиcв окрестностях горизонта zi здесь является отрицательной величиной, в то время как разность (Tпр3 - Tпр2) положительна. Т. е. обеспечивается результат измерений "с точностью до наоборот". Если же учитывать выражение (4), то получается (см. рис. 1), что измеренная производная на горизонте z2 существенно положительна, а на горизонте z3 близка к нулю. В этой ситуации второй член в выражении (4) является отрицательным. Это обеспечивает уменьшение разности (Tпр3 - Tпр2) - для наглядности полезно сравнить наклоны прямых 1 и 2.

Рисунок 1 - Искажения за счет динамических свойств канала измерения температуры

Таким образом, использование выражения (4) дает возможность некоторой коррекции вычисляемых разностей Tис по данным прямых измерений Tпр. Если же коррекция не производится, то результат может получаться таким, как на приводимом рисунке.

Сложность же этой коррекции состоит в том, что, во-первых, приборные градиенты в выражении (4) не измеряются непосредственно, а в свою очередь, вычисляются по дискретным значениям. Это может приводить к значительным искажениям. Во-вторых, значения Vз и фe только считаются константами, на самом деле они изменчивы. Это требует некоторых пояснений.

В первую очередь следует отметить, что зондирование не может производиться со скоростью свободного падения зонда в воде - это приведет к обрыву троса потере дорогостоящего прибора. Обычно задается некоторая вынужденная скорость (примерно 1 м/с). Однако судно, с которого производится зондирование, находится в условиях качки (штиль маловероятен). Качка приводит к вертикальным перемещениям точки подвеса зонда, что вызывает изменения скорости Vз. Подробнее эта задача рассмотрена в работе [4].

Во-вторых, постоянная времени фe также зависит от ряда факторов. В идеальном случае (если датчик не облепила нефтяная пленка, либо планктон, и т.д.) значение фe зависит от скорости зондирования, т.е. от скорости обтекания датчика окружающей водой.

К сожалению, высказанные замечания, хотя и не отвергают возможности расчетных методов коррекции, но не позволяют считать эти методы идеальными.

Из приведенного рисунка также видно, что в целом измеренный профиль сохраняет свой вид, лишь смещен по глубине и сглажен. Соответственно, следует полагать, что текущие измеренные градиенты эквивалентны истинным, но сдвинуты по глубине ("заглублены) и их изменчивость также сглажена.

Действительно, при предыдущем условии безинерционности канала измерений глубины получается

(5)

где FT(ikT) - комплексная передаточная функция канала измерений температуры; kТ - пространственная "круговая частота". Эту частоту можно выразить через масштаб вертикальной изменчивости температуры lТ:

Здесь lТ - масштаб колебаний в метрах.

Вид передаточной функции определяется выражением

, (6)

Сглаживание текущих градиентов характеризуется амплитудно-частотной характеристикой

, (7)

Сдвиг по глубине измеренных градиентов может быть выражен через фазо-частотную характеристику

(8)

Кажущаяся простота выражения (8), как бы позволяющая легко определить текущий сдвиг дzт, на самом деле применима на практике весьма ограничено, поскольку сюда входит масштаб lТ. Для реальных профилей этот масштаб является существенно переменным по глубине, соответственно также изменчив сдвиг измеренных градиентов.

Из приведенных рассуждений можно сделать следующий вывод. Если решаемые при исследованиях физические задачи позволяют примириться с динамическими искажениями, то о них просто следует помнить и не делать скоропостижных выводов о характеристиках какой-либо "сверхтонкой вертикальной структуры". Если же примириться нельзя, то либо надо поменять используемые приборы, например, поменяв характеристики канала измерения глубины.

Применение инерционных измерителей глубины, несомненно изменяет рассмотренную ситуацию. Действительно, в таком случае выражение для измеренного (приборного) градиента может быть записано в виде

(9)

где zпр(z) - показания канала измерения глубины; Fz(ikz) - передаточная функция канала измерения глубины, zис3 и zи2 - истинные значения глубины.

На первый взгляд, уравнивая динамические характеристики обоих измерительных каналов, легко получить полную коррекцию искажений. К сожалению, это не так. У зависимостей Tис(z) и zис(z) различна изменчивость. Глубина изменяется практически линейно (конечно, при отсутствии качки судна, с которого производится зондирование), а Tис(z) - с весьма сложными закономерностями, которые, собственно, и подлежат изучению в выполняемом зондировании.

Тем не менее, преимущества все же появляются. Действительно, за счет инерционности возникает запаздывание измеренных значений глубины по сравнению с истинными. Т.е. измеренные (приборные) профили температуры (см. рис. 1) как бы смещаются вверх на эту величину запаздывания.

Рассмотрим уравнение, описывающее динамические свойства канала глубины

, (10)

где фez - постоянная времени канала глубины.

Зададим zис(z) изменяющейся по линейному закону

, (11)

где zи0 - некоторая начальная глубина, например, "начало" пикноклина; b- постоянный коэффициент.

Будем искать решение уравнения (10) в виде

(12)

Подставляя решение (12) в уравнение (10), получим

(13)

Таким образом, уравнивая дzz по выражению (13) и дzТ по выражению (8) при некотором среднем ожидаемом значении масштаба lT, можно получить для измеренных значений градиента коррекцию среднего сдвига по глубине. При этом остается некоторая остаточная знакопеременная погрешность, зависящая от диапазона вариаций масштаба lТ.

Естественно, что погрешности за счет сглаживания вертикального профиля температуры остаются теми же, что и ранее, т.е. в соответствии с выражениями (5) и (7). Конкретные оценки в силу изменчивости масштаба lТ могут быть выполнены лишь для функций пространственной спектральной плотности

(14)

где Sпр(kT) - измеренная (приборная) функция спектральной плотности изменчивости температуры по вертикали; Sис(kT) - истинная функция спектральной плотности.

В заключение следует отметить, что прямые измерения градиентов, в частности, рассмотренные в работе [2], при одинаковых фe,Vз и lT практически не дают преимуществ по сравнению с вышерассмотренным вариантом косвенных измерений. Действительно, в случае прямых измерений выражения (5) и (9) остаются такими же. Достоинства прямых измерений проявляются преимущественно по отношению к статическим погрешностям [2].

Рассогласование динамических свойств измерительных каналов при CTD-зондировании в настоящее время никак не учитывается при использовании получаемых конечных данных. Однако, как здесь показано, результаты могут получаться «с точностью до наоборот», что приводит к существенным искажениям информации о реальных гидрофизических процессах.

При отсутствии технических возможностей согласования динамических свойств каналов температуры и давления целесообразно воспользоваться расчетным способом коррекции данных с использованием выражения (4) с контролем скорости зондирования

Библиографический список

1. Доронин Ю.П. Физика океана.- СПб.: изд. РГГМУ, 2000.- 340 с.

2. Ковчин И.С., Степанюк И.А. Методы специальных океанологических измерений.- СПб.: Изд. РГГМУ, 2002.- 271 с.

3. Степанюк И.А., Лебедева И.К. Погрешности вертикального зондирования из-за рассогласования динамических свойств измерительных каналов // Космогеофизические и гидрофизические факторы в морских технологиях. Под ред. проф. И.А. Степанюка.- СПб.: Астерион, 2008.- С.79-83

4. Степанюк И.А. Информационно-измерительные системы в океанологии.- СПб.:Изд. РГГМУ, 1998.- 90 с.

Размещено на Allbest.ru