Методы повышения точности измерения массы жидкости в резервуаре Текст научной статьи по специальности «Физика»

УДК 531.7

МЕТОДЫ ПОВЫШЕНИЯ ТОЧНОСТИ ИЗМЕРЕНИЯ МАССЫ ЖИДКОСТИ

В РЕЗЕРВУАРЕ

© 2012 А. В. Солнцева, Б. В. Скворцов, С. А. Борминский

Самарский государственный аэрокосмический университет имени академика С. П. Королёва

(национальный исследовательский университет)

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

Масса жидкости, резервуар, повышение точности, калибровка, нефтепродукт, интерполяция данных.

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

при изменяющихся внешних и внутренних условиях.

На рис. 1 схематично представлен резервуар 1, заполненный контролируемой жидкостью 2 на уровне Н. В резервуар опускается измерительный зонд 3, в частном случае труба, с отверстиями для размещения датчиков 4. Размещение датчиков 4 на разных уровнях, определяемых координатой у (0<у<Н), внутри резервуара 1 позволяет получить значение градиента параметров, таких, как температура Т, давление р, плотность р, вязкость п и др.

Полученная с датчиков 4 информация о параметрах контролируемой жидкости поступает на устройство сбора информации 5, затем через блок передачи данных 6 в устройство обработки данных 7.

В общем случае масса продукта М контролируемой жидкости связана с косвенными измеряемыми параметрами обобщенной функциональной зависимостью

М = ¥(цх,..цк..цп), (1)

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

Рис. 1. Резервуар с контролируемой жидкостью: 1 - резервуар; 2 - контролируемая жидкость, параметры которой изменяются в зависимости от глубины налива; 3 - измерительный зонд; 4 - датчики; 5 - устройство сбора данных; 6 - блок передачи данных; 7 - устройство обработки

данных

В общем случае функция F неизвестна. Её поиск является важнейшей задачей математического описания измерительного процесса. Учитывая, что любой измерительный процесс по определению подразумевает в себе, кроме сбора и обработки информации, также операцию калибровки по эталонам, для решения поставленной задачи предложен алгоритм измерений, показанный на рис. 2. Сущность алгоритма состоит в том, что для нахождения функции преобразования при неизвестной функции F составляется калибровочная модель процесса измерения. Для этого в резервуар i раз поэтапно наливают количество жидкости известной массы М, измеряют соответствующую ей совокупность контролируемых физико-химических параметров qk,i, к - соответствует номеру измеряемого параметра (к=1...т), i - номеру калиб-

ровочного замера @=1...п), значения которых запоминаются в центральном устройстве обработке. Количество эталонных замеров п должно быть не меньше числа контролируемых параметров т , то есть п > т [1].

Пусть имеется п калибровочных наливов жидкости с известными значениями массы М1, ..М, ...Мп . При этом следует помнить, что М{=М-.1+АМ, где АМ - фиксированное приращение массы на каждом этапе калибровки. При каждом наливе контролируется совокупность параметров мониторинга резервуара qк,i. Отметим, что при каждом очередном наливе физико-химические параметры жидкости будут меняться, так как процесс налива при больших объемах резервуара длительный, то будет меняться температура Т, соответственно плотность р и вязкость п основного и подтоварного продукта.

41 42 ■ ■ 4к ■■

1 Г 1 ' \ ' 1 г

Калибровка

г

Запись калибровочных коэффициентов

Измерение

Запись измеренных значений

Рис. 2. Алгоритм измерений

При построении калибровочной моде- стемы (2) всегда можно выбрать достаточное

ли для произвольных значений измеряемых количество уравнений, из которых можно

параметров всегда можно подобрать такие найти коэффициенты К,. Желательно, чтобы

нормирующие коэффициенты Ьк, что будут выбранные для совместного решения урав-выполняться равенства:

М1 = Ьчхл + •••+ъкЧк,х + •••+ъшЧт,1 =^ЕькЧкл;

к=1

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

Мг = Ь^ + ••• + Ькчкг + ••• + Ьтчтг = £Ькчкг; (2) Система (2) имеет единственное решение,

к=1 '

Мт = Ь1^1, т + •••+ Ьк4к,т + ••• + Ьт4т, т = ^Ьк4к, т

к=1

Мп = Ь1^1,П + ••• + Ьк4к, п + ••• + Ьт4т, п = ^Ьк4к,п •

к=1

Система (2) имеет п уравнений для т неизвестных ЬТак как число эталонных замеров всегда можно сделать больше числа контролируемых параметров п>т, то из си-

если главный определитель отличен от нуля. Решением этого уравнения являются соотношения:

К = ^, (3)

где

4х 1 • -4к , 1 • -Чп, 1 ¡1, 1 • М1 • -Чп, 1

л = 41, ••4кг • •4п,г ; ль, = ¡1, г • М •¡п.

41, п -4к, п -•¡п. п 41 п •Мп ••¡п, п

Мк - совокупность эталонных значений массы контролируемой жидкости, используемой при калибровке; д - совокупность физико-химических параметров жидкости, измеренных при калибровке; к - соответствует номеру измеряемого параметра (¡=1...п), г - номеру калибровочного замера (к=1 ...п).

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

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

мх = Ъ&,х + ■■■+ЪкЧк,х + ■■■+Ътдт,х =

(4)

■ТЬк^к,

,Х'

к=1

где Мх - искомое значение массы контролируемой жидкости; дк,х - совокупность измеренных параметров контролируемой жидкости.

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

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

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

Для пояснения физического смысла предложенного способа проанализируем простейший случай, когда измеряется только два параметра - уровень Н и плотность жидкости р. При этом считаем, что температура жидкости не изменяется и равна температуре жидкости, наливаемой при калибровочном отсчете, а подтоварная жидкость отсутствует. Очевидно, что в этом случае масса жидкости будет вычисляться по формуле:

Мх = 8Нх рх, (5)

где £ - площадь основания емкости.

При вычислении массы по данной методике, для рассматриваемого случая калибровку достаточно проводить по двум эталонным массам, полностью перекрывающим диапазон ожидаемых значений Мт\п Мтах . В формуле (3) будем иметь: М1 = Мтп ; М2 = Мтах ; Ч11=Н1; Ц12=Н2 ; Ц21=Рь Ч22=Р2. Тогда выражение (2) примет вид:

Мх = Ъ нх+ъ2 рх, (6)

где

ъ = р 2М1 - рм 2

Н21У± 1 ~ /-У" 2 . Ъ = Н1М2 - Н2М1 (7)

Н1р2 - Н2 Р/ 2 Н1р2 - Н2 Р1 ' Если плотность жидкости в процессе калибровки не изменялась, что вполне может соответствовать действительности, то есть р1 = р2 = р, то выражения (7) и (6) примут вид:

Н М - Н М

М1 - М2

Ъ = ■ 1 2

Н1 - Н2

Ъ2 =

1 2^-4

р(Н 1 - Н 2)

(8)

Для подтверждения правильности предложенного способа покажем, что при рХ = р выражения (5) и (9) идентичны. Учи-

тывая, что М1 = р£Иг, М2 = р£И 2 и подставляя эти выражения в (9), получим (5).

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

Способ предусматривает возможности неограниченного нарастания точности измерений, связанной с увеличением числа контролируемых параметров. Например, если к контролю уровня и плотности в процессе калибровки и измерений добавить контроль температуры, то выражение (4) примет вид: Мх = ь Их + Ь2 рх + ЬЪТХ, (10)

где коэффициенты Ьк определяются по формуле (3) с определителями

а =

И1 И 2

И

Л р 2

р 3

Тх

Т2 Т

аЬ1 =

аЬ2 =

их м1 тх

И2 М2 Т2

И3 М3 Т3

М1 Р1 Т1

= М2 Р2 Т2

М3 Р3 Т

И1 Р1 М1

И 2 Р2 М2

И3 Р3 М3

резервуаре.

Приведенная в формуле (5) площадь основания емкости £ в реальных условиях эксплуатации не является постоянной величиной, а изменяется в зависимости от многих факторов, в том числе и от уровня наполнения [2]. Поскольку резервуар подвержен естественной деформации при его заполнении продуктом, при вычислении объема необходимо учитывать изменение площади основания £ по всему участку измерения у, то есть от дна резервуара у=0 до максимального уровня наполнения у=Н (рис. 1).

Отметим также, что изменение температуры одновременно приводит к изменению и плотности, и объема жидкости, поэтому масса при изменении температуры не меняется.

Значение плотности продукта р также в реальных условиях не является постоянным на всей протяженности у & Н ], поэтому

справедливо выражение для определения массы:

Мх = |Р (у )• £ (у )у.

(11)

Нетрудно показать, что если температура постоянна и при измерениях такая же, как и при калибровке, то выражение (10) превращается в (6). Если температура и плотность постоянны и такие же, как и при калибровке, то выражение (10) превращается в (5), что еще раз подтверждает правильность предложенного способа измерений и возможность учета различных изменяющихся факторов при вычислении массы жидкости в резервуаре.

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

Также предлагается еще один метод определения массы жидкости, находящейся в

Площадь основания резервуара и объем наполняемого продукта в общем случае математически связаны формулой:

где £(у) - площадь горизонтального сечения

тела с данной аппликатой у.

Согласно (12), выражение (11) перепишется следующим образом:

М.

(IV.

1

(13)

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

3

18-49 шп ЛБ 18-49шп ПБ 49-88 шп. ЛБ 49-88 шп. Пр.Б 88-127 шп. ЛБ 88-127 шп. Пр.Б 127-165 шп. Л Б 127-165 шп. Пр.Б

Уровень 1 2 3 4 5 6 7 8

0 0 0 0 0 0 0 0 0

0,01 0,19 0,19 0,71 0,71 0,74 0,74 0,7 0,7

0,02 0,38 0,38 1,43 1,43 1,48 1,49 1,41 1,42

0,03 0,58 0,58 2,15 2,15 2,23 2,24 2,12 2,14

0,04 0,78 0,78 2,88 2,89 2,98 2,99 2,84 2,86

0,05 0,99 0,99 3,62 3,63 3,74 3,76 3,57 3,6

0,06 1,21 1,2 4,36 4,37 4,5 4,52 4,29 4,33

0,07 1,43 1,43 5,1 5,11 5,27 5,29 5,03 5,07

0,08 1,66 1,65 5,85 5,86 6,05 6,08 5,76 5,82

0,09 1,89 1,88 6,62 6,62 6,83 6,86 6,51 6,57

0,1 2,13 2,11 7,38 7,39 7,62 7,64 7,27 7,32

Рис. 3. Фрагмент калибровочной таблицы на восемь емкостей

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

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

Полученные в результате калибровки значения объема и уровня необходимо представить в виде аналитической зависимости У (у). В результате применения полиномиальной регрессии, получим зависимость вида:

У(У) = . (14)

1=0

Плотность продукта может быть определена несколькими методами. Как видно из рис. 1, в измерительном зонде (трубе с открытым концом) имеется несколько отверстий, в которых установлены датчики, в том числе плотномер. Таким образом, в первом случае возможно проведение к измерений плотности на уровнях расположения отверстий в измерительном зонде. Во втором случае, по мере заполнения емкости, возможно перемещение датчика плотности на величину реперного шага, таким образом, может быть произведено к2 измерений. Способ измерения плотности определяется исходя из технического оснащения резервуарного парка.

Аппроксимируя данные, полученные о плотности жидкости, получаем зависимость (для удобства пользуемся тем же видом ре-

грессии, что и для зависимости объема от уровня):

р(у) = £ьу. (15)

3=0

С учетом (14) и (15) выражение (13)

примет вид:

>.' к-Г

,Н1 (1:У 'У/!? :"Г '

м | ' " яг т

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

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

Кроме того, выражение (16) даёт возможность создать универсальный алгоритм и программу обработки данных для ЭВМ, в которую будут вводиться только данные о коэффициентах аппроксимации и текущих значениях плотности.

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

данных, применяемых в настоящее время на всех нефтебазах РФ.

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

1 . Скворцов, Б.В. Электрофизические устройства контроля качества углеводород-

ных топлив [Текст] / Б. В. Скворцов. - Самара: Изд-во СГАУ, 2000. - 264 с.

2 . Скворцов, Б.В. Импульсные методы измерений количества и качества жидких углеводородных топлив [Текст] / Б.В. Скворцов, Борминский С. А. - Самара: Изд-во СНЦ РАН, 2010. - 220 с.

METHOD TO IMPROVE ACCURACY MASS MEASUREMENT FLUID

IN THE TANK

© 2012 A. V. Solntseva, B. V. Skvortsov, S. A. Borminsky

Samara State Aerospace University named after academician S. P. Korolyov (National Research University)

Considered two methods to improve the accuracy of measuring the mass of fluid in the tank based on the non-linear interpolation of the experimental calibration data. The first method is to find a transformation function mass and the total parameters of the fluid measured by monitoring with the help of the calibration model of the measurement process. The second method is based on the determination of an individual analytical dependence for the tank, connecting a lot of fluid, with its qualitative and quantitative parameters.

Fluid mass, reservoir, improve the accuracy, calibration, petroleum, interpolation of data.

Информация об авторах

Солнцева Александра Валерьевна, аспирант, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: аналитическое приборостроение.

Скворцов Борис Владимирович, доктор технических наук, профессор, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: аналитическое приборостроение.

Борминский Сергей Анатольевич, кандидат технических наук, доцент, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: аналитическое приборостроение.

Solntseva Alexandra Valeryevna, post-graduate student of Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Area of research: analytical device construction.

Skvortsov Boris Vladamirovitch, doctor of technical science, professor of Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Area of research: analytical device construction.

Borminsky Sergey Anatolyevich, candidate of technical science, docent of Samara State Aerospace University named after academician S. P. Korolyov (National Research University), doctoral of Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Area of research: analytical device construction.