СПОСОБ АВТОМАТИЧЕСКОГО ОПРЕДЕЛЕНИЯ В УСЛОВИЯХ ОКЕАНА ПАРАМЕТРОВ СОСТОЯНИЯ МОРСКОЙ ВОДЫ
METHOD FOR AUTOMATIC DETERMINATION OF PARAMETERS OF THE STATE OF SEAWATER IN OCEAN CONDITIONS

Изобретение относится к области океанологических исследований и может быть использовано при построении автономных и зондирующих устройств и для определения комплексных гидрофизических параметров морской воды океана, а также может быть использовано для контроля различных примесей, пузырьков газа, взвешенных твердых частиц не только в морской, но и в речных и сточных водах и для исследования физических свойств воды, смесей жидкостей и уточнения эмпирических выражений при определении термодинамических величин и молекулярных свойств жидкостей.
Известны способы определения взаимосвязанных гидрологических параметров морской воды путем измерения текущих значений электропроводности и температуры, вычисления других параметров, например, солености по различным эмпирическим формулам [Маклаков А.Ф. и др. Океанографические приборы. Л.: Гидрометеоиздат, 1975. С. 261-270].
Известен также способ автоматического определения в условиях океана параметров состояния морской воды: температуры, солености и плотности, путем синхронного измерения текущих значений двух из совокупности взаимосвязанных параметров и последующего вычисления значений других параметров [Маклаков А.Ф. и др. Океанографические приборы. Л.: Гидрометеоиздат, 1975. С. 54-59. С. 281-283].
Известен способ автоматического определения в условиях океана параметров состояния морской воды с целью повышения быстродействия и разрешающей способности при комплексном определении в квазиреальном масштабе времени средних и пульсационных значений температуры, солености и плотности, синхронно измеряют текущие значения оптического показателя преломления и скорости распространения звука [А.С. СССР 717633, Кл. G01N 21/46. Способ автоматического определения в условиях океана параметров состояния морской воды. Брамсон М.А., Гончаров Э.Г., Красовский Э.И., Митрофанов И.А., Морозова А.А., Федоровский А.Д.; Опубл. 25.02.80. Бюл. №7.].
К недостаткам известного способа относится то, что в указанном способе не определяются коэффициенты акустического поглощения, затухания, рассеяния и суммарный диссипативный коэффициент.
Целью предлагаемого изобретения является повышение разрешающей способности при комплексном определении в квазиреальном масштабе времени параметров состояния морской воды.
Задача решена тем, что в известном способе, где измеряется температура и скорость звука, дополнительно синхронно измеряют текущие значения акустических коэффициентов затухания и рассеяния на фиксированной частоте и стабильной измерительной базе, а также измеряют гидростатическое давление, после чего последовательно определяют плотность, соленость, коэффициент поглощения и суммарный диссипативный коэффициент
где b - суммарный диссипативный коэффициент;
п
п
з
p
k - коэффициент поглощения, k=k-k;
з
k - коэффициент затухания;
p
k - коэффициент рассеяния;
ρ - плотность морской воды;
с - скорость звука в воде;
ƒ - частота ультразвукового сигнала.
Для реализации предлагаемого способа должны быть использованы измерители температуры, скорости звука, давления, акустических коэффициентов затухания и рассеяния и частоты, которые расположены на хорошо промываемой стабильной измерительной базе, выполненной, например, из ситала, благодаря чему обеспечивается необходимое пространственно-временное разрешение. Структурно-функциональная схема измерителя представлена на фигуре, где изображены: 1 преобразователь временных и амплитудных сигналов; 2 - преобразователь ультразвука; 3 - половинчатый отражатель; 4 - полный отражатель; 5 - дополнительный акустический приемник; 6 - стабильная измерительная база; 7 - усилитель сигналов; 8 - блок управления; 9 - датчик температуры; 10 - датчик давления; 11 - ПЭВМ (персональная электронно-вычислительная вычислительная машина).
Представленная на фигуре схема работает следующим образом. Электрический сигнал с фиксированной частотой с преобразователя 1 поступает на преобразователь ультразвука 2, где преобразуется в акустический сигнал, который проходит через водную среду, отражается от первого половинчатого отражателя 3, затем отражается от полного отражателя 4 и оба отраженных акустических сигналов поступает на ультразвуковой преобразователь 2, где преобразуется в электрический сигнал, который поступает на преобразователь 1. По разности времен между излученным и отраженным сигналами определяется скорость звука, амплитуда приемного сигнала характеризует коэффициент затухания. Рассеянный акустический сигнал от неоднородностей поступает на дополнительный акустический приемник 5. Все отражатели 3, 4, акустический преобразователь 2, акустический приемник 5 расположены на стабильной измерительной базе (ситал) 6. Электрический сигнал с акустического приемника 5, амплитуда которого пропорциональна коэффициенту рассеяния, через усилитель 7 поступает в блок управления 8. В зоне прохождения акустических сигналов расположен малоинерционный датчик температуры 9, сигнал с которого поступает на преобразователь 1. Одновременно датчиком давления 10 измеряется гидростатическое давление, электрический сигнал с которого поступает в блок управления 8. Блок управления 8 формирует сигналы управления и одновременно осуществляет прием и преобразование сигналов в последовательный код, который через линию связи поступает в ПЭВМ 11, где происходит преобразование последовательного кода в значение температуры, гидростатического давления, скорости звука, затухания и рассеяния ультразвука, используя измеренные величины с учетом эталонной базы и фиксированной частоты, определяют плотность, соленость, коэффициент поглощения и суммарный диссипативный коэффициент.
Измерение скорости звука основано на времяпролетном методе измерений. Каждое измерение скорости звука выполняется с помощью одного акустического импульса, проходящего стабильную измерительную базу. В зависимости от длины измерительной базы, которая формируется за счет высокостабильных стержней из ситала или карбонового композитного материала, можно получить хорошее разрешение и точность по скорости звука.
H
0
з
При измерении коэффициента затухания необходимо знать начальную амплитуду U и амплитуду отраженного сигнала U при известном расстоянии, и далее, используя формулу  рассчитывают коэффициент затухания k.
0
Датчик скорости звука имеет два отражателя, поэтому определить амплитуду U проблем не составляет. Разность амплитуды сигналов между приемными сигналами от двух отражателей и будет характеризовать величину коэффициента затухания ультразвука в среде на длине базы между двумя отражателями.
Но учитывая то, что половинчатый отражатель имеет полупрозрачную конструкцию, коэффициент его отражения учитывается при градуировке прибора.
Схема измерения амплитуды сигнала осуществляется с помощью преобразователя TDC-GP30. Алгоритм измерения амплитуды работает следующим образом. После принятия отраженного сигнала пиковым детектором фиксируются пиковая амплитуда. Затем после окончания действия пакета импульсов приемного сигнала с помощью аналого-цифрового преобразователя с учетом наклона импульса определяется время его разряда в границах амплитуды от максимальной до установления уровня, в нашем случае 0,7.
По результатам измерения времени вычисляют амплитуду акустического сигнала по следующим формулам [Ultrasonic-Flow-Converter Data Sheet TDC-GP30 June 27th, 2019 Document-No: DB_GP30Y_Vol1_en V0.3 System-Integrated Solution for Ultrasonic Flow Meters Volume 1: General Data and Frontend Description (Available at https://ams.com/documents/20143/36005/TDC-GP30_DS000391_3-00.pdf/f96f8c8b-87e5-ac8d-a26c-65756fd240fa) (дата обращения: 03.06.2019)]:
где
Для измерения коэффициента рассеяния используют акустический приемник, который расположен под определенным углом на стабильной измерительной базе. В зависимости от количества находящихся примесей в воде уровень рассеянного акустического сигнала изменяется, и при соответствующей градуировке получаем коэффициенты апроксимирующей функции.
В качестве примера для расчета плотности приведем разработанное авторами новое уравнение
0
где γ=(ρ-ρ)/ρ*;
0
τ=(Т-Т)/Т*;
абс
0
π=(Р-Р)/Р*;
0
ω=(с-с)/с*;
0
3
ρ=990 кг/м;
3
ρ*=100 кг/м;
0
Т=-10°С;
Т*=50°С;
0
Р=0,101325 МПа;
Р*=120 МПа;
0
с=1300 м/с;
с*=500 м/с.
Для уравнения (1) использована индексная матрица, приведенная в табл. 1.
ijk
Здесь знак «+» означает, что данная комбинация индексов использована, а знак «-» означает, что данная комбинация индексов пропущена. Всего в данной индексной матрице использовано 80 комбинаций. Таким образом, уравнение (1) содержит всего 80 ненулевых коэффициентов. Коэффициенты b уравнения (1) были рассчитаны авторами по методу наименьших квадратов, и совместно с индексами приведены в табл. 2.
Для определения солености используется разработанное авторами уравнение скорости звука для соленой воды, включая пресноводную и солевую части
где для вычисления пресноводной части скорости звука используется уравнение
где Т - температура, °С;
Р - абсолютное давление, МПа.
Диапазон применимости: по давлению: от 0,101325 до 120,0 МПа; по температуре: от кривой плавления до 300°С.
ij
Коэффициенты уравнения b приведены в табл. 3.
Для вычисления солевой части скорости звука используется уравнение
km
где d
Следовательно, используя синхронно измеренные значения температуры T, давления Р, скорости звука в воде с на фиксированной частоте ƒ, определяют по представленным алгоритмам плотность ρ и соленость S.
p
з
п
з
p
Измеренные коэффициенты рассеяния k и затухания k позволяют определить коэффициент поглощения k=k-k, и далее определяется суммарный диссипативный коэффициент
п
где k - коэффициент поглощения;
ρ - плотность морской воды;
с - скорость звука в воде;
ƒ - частота ультразвукового сигнала.
Все измеренные значения определяются на стабильной эталонной базе для малого объема жидкости, что позволяет повысить пространственную, а за счет синхронного измерения и временную разрешающую способность.