Результат интеллектуальной деятельности: СПОСОБ ИЗМЕРЕНИЯ ВЛАГОСОДЕРЖАНИЯ ДИЭЛЕКТРИЧЕСКОЙ ЖИДКОСТИ

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

Известны различные способы и устройства для измерения влагосодержания жидкостей, основанные на измерении электрофизических параметров (диэлектрической проницаемости или(и) тангенса угла диэлектрических потерь) жидкостей с применением радиоволновых ВЧ и СВЧ резонаторов, содержащих контролируемую жидкость (монографии: Брандт А.А. Исследование диэлектриков на сверхвысоких частотах. М.: Физматгиз. 1963. 403 с. С. 37-144; Викторов В.А., Лункин Б.В., Совлуков А.С. Радиоволновые измерения параметров технологических процессов. М.: Наука. 1989. 208 с. С. 168-177). Недостатком таких способов и реализующих эти способы измерительных устройств является их ограниченная область применения, обусловленная невозможностью контроля малых изменений влагосодержания жидкостей ввиду невысокой точности измерения соответствующих малых изменений информативных параметров (резонансной частоты, добротности резонатора и др.). Для обеспечения возможности проведения таких измерений применяют двухканальные измерительные схемы с независимыми измерительным и эталонным каналами. В эталонном канале чувствительный элемент содержит жидкость с известными физическими свойствами (монография: Брандт А.А. Исследование диэлектриков на сверхвысоких частотах. М.: Физматгиз. 1963. 403 с. С. 258-268).

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

Техническим результатом настоящего изобретения является расширение функциональных возможностей способа измерения за счет повышения точности измерения.

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

Предлагаемый способ поясняется чертежом. На фиг. 1 приведена схема устройства, поясняющая принцип измерения с применением способа.

Здесь показаны волновод 1, генераторы 2 и 3, коммутатор 4, элементы связи 5 и 6, детектор 7, функциональный преобразователь 8, регистратор 9.

Способ реализуется следующим образом.

В данном способе реализуют структурный подход к достижению инвариантности к диэлектрической проницаемости ε_н контролируемой жидкости, в частности, к ее сортности, изменения которой имеют место, в частности, при контроле нефти и нефтепродуктов в какой-либо емкости или в процессе их транспортирования по трубопроводу. Этот подход связан с организацией двух измерительных каналов (двух последовательных тактов измерений на двух разных фиксированных частотах ƒ₁ и ƒ₂) и совместным функциональным преобразованием их выходных величин с целью получения результата этого преобразования, который не зависит от возмущающего фактора, в данном случае - от величины диэлектрической проницаемости ε_н контролируемой жидкости и изменений ε_н.

Предлагаемый способ заключается в возбуждении электромагнитных волн в волноводе, используемом в качестве измерительной ячейки (при измерениях в стационарных условиях) или в качестве отрезка измерительного волновода, встроенного в трубопровод с перекачиваемой жидкостью (при измерениях в трубопроводах). Электромагнитные волны в волноводе возбуждают последовательно на двух разных фиксированных частотах ƒ₁ и ƒ₂, каждая из которых ниже критической частоты ƒ_кр для волны низшего типа. При этом на каждой из этих частот вдоль волновода существует только ослабевающее реактивное поле, убывающее при удалении от возбуждающего элемента у одного из торцов емкости. Электромагнитные волны принимают после их распространения вдоль данного волновода на другом его торце и измеряют на каждой из частот ƒ₁ и ƒ₂ амплитуду напряженности электрического поля. По результату совместного преобразования измеренных амплитуд судят о влагосодержании жидкости.

Условием распространения электромагнитных волн по любому волноводу является выполнение неравенства: ƒ>ƒ_кр, которому должны удовлетворять рабочая частота ƒ и критическая частота ƒ_кр для волны низшего типа, например, для волны Н₁₁ в круглом волноводе. При ƒ<ƒ_кр имеет место режим, при котором распространения волн по волноводу не происходит, а существует только ослабевающее реактивное поле, убывающее при удалении от возбуждающего элемента (Лебедев И.В. Техника и приборы СВЧ. Т.1. М.: Высшая школа. 1970. 440 с. С. 132-136).

При этом электрическое поле (как и магнитное поле) изменяется вдоль координаты z (оси волновода) по закону:

а постоянная ослабления α есть

В этих формулах Е₀ - амплитуда напряженности электрического поля при z=0; ε - диэлектрическая проницаемость диэлектрической жидкости в волноводе, с - скорость света. Выбирая соотношение между ƒ и ƒ_кр, можно управлять величиной ослабления α.

Если частота ƒ генератора меньше критической частоты ƒ_кр данного волновода, то амплитуда напряженности Е электрического поля, являющаяся информативным параметром, в точке приема есть

где , Е₀ - амплитуда напряженности поля в области возбуждения электромагнитных волн в рассматриваемом волноводе (т.е. в области расположения связи 3), λ_кр - критическая длина волны для данного волновода, - длина измерительного участка, т.е. расстояние вдоль волновода между элементами возбуждения и съема электромагнитных колебаний. Напряженность электрического поля Е при удалении от элемента связи, служащего для возбуждения и приема электромагнитных колебаний, спадает в соответствии с соотношением (3). При этом значение Е зависит от влагосодержания контролируемой жидкости в волноводе.

Для волн низшего типа Н₁₁ имеем λ_кр=3,41d/2 и, соответственно, , где d - внутренний диаметр волновода (Лебедев И.В. Техника и приборы СВЧ. Т.1. М.: Высшая школа. 1970. 440 с. С. 132-136). Например, при d=30 мм для волн типа Н₁₁ будем иметь λ_кр=3,41d/2=5,115 см; тогда ƒ_кр=5,865/√ε ГГц. Если, например, ε=2,0, то должно быть ƒ<ƒ_кр=5,865/√ε ГГц=4,148 ГГц.

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

Согласно данному способу, определение влагосодержания производят при проведении двух последовательных тактов измерений на фиксированных частотах ƒ₁ и ƒ₂, каждая из которых меньше критической частоты ƒ_кр волновода.

В первом такте измерений, при возбуждении в волноводе на одном из его торцов электромагнитных волн на фиксированной частоте ƒ₁, меньшей критической частоты ƒ_кр волновода, принимают электромагнитные волны после их распространения вдоль волновода на другом его торце и измеряют амплитуду Е₁ напряженности электрического поля:

где , E₁₀ - амплитуда напряженности поля в области возбуждения электромагнитных волн в рассматриваемом волноводе, - длина измерительного участка, т.е. расстояние вдоль волновода между элементами возбуждения и съема электромагнитных колебаний, λ_кр - критическая длина волны в волноводе, ε(ƒ₁) - диэлектрическая проницаемость жидкости на частоте ƒ₁.

Во втором такте измерений, при возбуждении в волноводе на одном из его торцов электромагнитных волн на фиксированной частоте ƒ₂, меньшей той же критической частоты ƒ_кр волновода (или иной критической частоты для какого-либо другого рабочего типа волн, если его выбирают для измерений), принимают электромагнитные волны после их распространения вдоль волновода на другом его торце и измеряют амплитуду Е₂ напряженности электрического поля:

где , Е₂₀ - амплитуда напряженности поля в области возбуждения электромагнитных волн в рассматриваемом волноводе, - длина измерительного участка, т.е. расстояние вдоль волновода между элементами возбуждения и съема электромагнитных колебаний, λ_кр - критическая длина волны в волноводе, ε(ƒ₂) - диэлектрическая проницаемость жидкости на частоте ƒ₂.

Величина диэлектрической проницаемости ε влагосодержащего диэлектрика, имеющего диэлектрическую проницаемость ε_н и влагосодержание W, описывается при малых W формулой Винера (Теория и практика экспрессного контроля влажности твердых и жидких материалов / Кричевский Е.С., Бензарь В.К., Венедиктов М.В. Под общ. ред. Кричевского Е.С. М.: Энергия. 1980. 240 с. С. 55-66):

где , ε_в(ƒ) - диэлектрическую проницаемость воды, являющаяся функцией частоты ƒ в СВЧ-диапазоне частот.

Например, если производить измерения на частоте ƒ₁=10 ГГц, то D(ƒ₁)=1,095, а на частоте ƒ₂=37,5 ГГц, то D(ƒ₂)=1,383, если считать ε_н=2.

Выражения для D(ƒ₁), D(ƒ₂) в формулах (4) и (5) можно упростить, если положить D(ƒ₁)-W≈D(ƒ₁), D(ƒ₂)-W≈D(ƒ₂) и не зависят от ε_н, что допустимо при малых значениях влагосодержания (до ~5%) и реальных пределах изменения ε_н.

Постоянство величин D(ƒ₁) и D(ƒ₂) для соответствующих частот ƒ₁ и ƒ₂ вытекает из постоянства величин ε_н и ε_в, входящих в формулу для D(ƒ₁) и D(ƒ₂). Величина ε_н постоянна в широком диапазоне изменения частоты ƒ, величина ε_в постоянна на недисперсионном участке кривой ε_в(ƒ) и принимается постоянной на дисперсионном участке этой кривой. Это справедливо при проведении изменений с помощью измерительных устройств, работающих на фиксированных частотах.

Покажем теперь на реальном примере, что D(ƒ₁)-W и D(ƒ₂)-W, а, точнее, величины 3/(D(ƒ₁)-W) и 3/(D(ƒ₂)-W), входящие общем виде (при произвольном значении частоты ƒ) в формулу (1), не зависят (с некоторой допустимой погрешностью) от ε_н. Так, при реальном изменении ε_в на 10% по сравнению с первоначальным значением ε_н=2, т.е. до значения 2,2, при значениях частот ƒ₁=10 ГГц и ƒ₂=37,5 ГГц будем иметь:

при ε_н=2: D(ƒ₁)=1,095, D(ƒ₂)=1,383;

при ε_н=2,2: D(ƒ₁)=1,105, D(ƒ₂)=1,337.

Отсюда следует, что относительное изменение D(ƒ₁) есть ~0,9%, а относительное изменение D(ƒ₂) есть ~2,6%.

Оценим, как влияют эти изменения D(ƒ₁) и D(ƒ₂) на коэффициенты при W в формуле (1) при ƒ=ƒ₁ и ƒ=ƒ₂, а именно на коэффициент k₁=3/(D(ƒ₁)-W) и коэффициент k₂=3/(D(ƒ₂)-W):

при ε_н=2 имеем: k₁≈2,752, k₂≈2,31;

при ε_н=2,2 имеем: k₁≈2,727, k₂≈2,33.

Отсюда следует, что относительное изменение как k₁, так и k₂ есть ~0,9%, что в ~100 раз меньше относительного изменения ε_н, т.е. реальное изменение ε_н не влияет практически на k₁ и k₂. В формулах для D(ƒ₁) и D(ƒ₂) можно использовать для выражения ε_н значение - номинальное значение диэлектрической проницаемости обезвоженной диэлектрической жидкости.

Отметим, что для реализации данного способа измерений достаточно, если рассматриваемой частотной дисперсией обладает вода только на одном из двух рабочих частот (ƒ₁ или ƒ₂).

При проведении измерений в первом такте на частоте ƒ₁ будем иметь:

где , ε_в(ƒ₁) - диэлектрическая проницаемость воды, являющаяся функцией частоты ƒ₁ в СВЧ-диапазоне частот, - номинальное значение диэлектрической проницаемости обезвоженной диэлектрической жидкости.

При проведении измерений во втором такте на частоте ƒ₂ будем иметь:

где , ε_в(ƒ₂) - диэлектрическая проницаемость воды, являющаяся функцией частоты ƒ₂ в СВЧ-диапазоне частот, - номинальное значение диэлектрической проницаемости обезвоженной диэлектрической жидкости.

При достаточно больших значениях W следует использовать другие известные выражения для ε (Теория и практика экспрессного контроля влажности твердых и жидких материалов / Кричевский Е.С., Бензарь В.К., Венедиктов М.В. Под общ. ред. Кричевского Е.С. М.: Энергия. 1980. 240 с.).

При рассмотрении (7) и (8) как системы уравнений и ее решении относительно искомого влагосодержания W получим

С учетом формул (4) и (5) выражение (9) принимает следующий вид:

Формулу (10) запишем в следующем виде:

где k₁, k₂, k₃, k₄, k₅, k₆ - постоянные величины, причем

; ; ; ; ; .

Таким образом, осуществляя совместное преобразование измеряемых значений амплитуд Е₁ и Е₂ согласно соотношению (11), получаем значение влагосодержания W, которое не зависит от диэлектрической проницаемости ε_н контролируемой жидкости и ее возможных изменений.

На фиг. 1 приведена схема устройства для реализации данного способа.

Возбуждение электромагнитных волн в волноводе 1 осуществляется последовательно, в первом и втором тактах измерений, на частотах ƒ₁ и ƒ₂, меньших критической частоты ƒ_кр для этого волновода с помощью, соответственно, генератора 2 и генератора 3 через коммутатор 4 и элемент связи 5. Другой элемент связи (приема) 6 электромагнитных волн расположен на расстоянии вдоль волновода 1. Принимаемые элементом связи 6 волны поступают на детектор 7 и затем продетектированные сигналы подаются в функциональный преобразователь 8, подсоединенный выходом к регистратору 9 для определения искомого влагосодержания, которое не зависит от диэлектрической проницаемости ε_н контролируемой жидкости.

Для волноводов конкретных размеров выбором частот ƒ₁ и ƒ₂ генераторов 2 и 3 можно оптимизировать чувствительность такого устройства для измерения влагосодержания жидкости в рабочем диапазоне его изменения. При этом имеет место монотонность зависимости информативных параметров - значений амплитуд Е₁ и Е₁ напряженности электрического поля - от значений соответственно ε₁ и ε₂, функционально связанных с измеряемым влагосодержанием жидкости.

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