Способ измерения массы жидкости в резервуаре

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

$$M_X=b_1{q}_{\mathrm{1,}X}+\dots +b_k{q}_{k,X}+\dots +b_n{q}_{n,X}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k,X}}$$

где M_X - искомое значение массы контролируемой жидкости, q_k,X - совокупность измеренных параметров контролируемой жидкости. При этом коэффициенты bk определяются по формуле:

$$b_k=\frac{\Delta b_k}{\Delta }$$ ,

где $$\Delta =\left|\begin{array}{c}{q}_{\mathrm{1,1}}\cdots q_{i\mathrm{,1}}\cdots q_{n\mathrm{,1}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}k}\cdots q_{i,k}\cdots q_{n,k}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}n}\cdots q_{i,n}\cdots q_{n,n}\end{array}\right|$$ , $$\Delta b_k=\left|\begin{array}{c}{q}_{\mathrm{1,1}}\cdots M_1\cdots q_{n\mathrm{,1}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}k}\cdots M_k\cdots q_{n,k}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}n}\cdots M_n\cdots q_{n,n}\end{array}\right|$$ .

M_k - совокупность эталонных значений массы контролируемой жидкости, используемой при калибровке, q_i,k - совокупность физико-химических параметров жидкости, измеренных при калибровке, i - соответствует номеру измеряемого параметра (i=1…n), k - номеру калибровочного замера (k=1…n). Технический результат, достигаемый от реализации заявленного способа, заключается в значительном повышении точности измерений, в том числе может быть эффективно использован для оперативного определения массы нефтепродуктов в резервуаре при изменяющихся внешних условиях, по различным параметрам, косвенно связанным с искомой величиной. 2 ил.

 

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

Аналогами изобретения являются способы и устройства, выполненные по книге «Контроль количества и качества нефтепродуктов» / под ред. Новоселова В.Ф. - М: Недра, 1994, с 127-134; а также по патентам № RU 2188396 «Способ измерения уровня, плотности, границы раздела и температуры жидкости в резервуаре и устройство для его реализации», № RU 2343426 «Устройство для измерения массы жидкого продукта в резервуаре», № RU 99152 «Автоматизированная система контроля и учета массы нефтепродуктов «Гамма» при их хранении в вертикальных и/или горизонтальных резервуарах на топливных базах» и заявкам на изобретение № 2006126629 «Система и способ мониторинга уровня топлива в топливном резервуаре», № RU 2007110252 «Устройство для контроля массы нефтепродуктов в цистерне», № RU 95107394 «Установка для измерения количества нефти и нефтепродуктов», № RU 95104019 «Устройство для обнаружения утечек в резервуарах», № RU 2006135822 «Способ автоматизированного контроля и учета массы нефтепродуктов на топливных базах при их хранении в вертикальных и горизонтальных резервуарах», а также реализованные в системе измерения массы светлых нефтепродуктов УИП 9602, руководство по эксплуатации АТУШ. 400000.01.РЭ, Королев, 2003; в системе измерительной «СТРУНА», ТУ 4210-001-2434704-2004, заключающиеся в том, что с помощью датчиков определяют уровень, температуру, плотность жидкости, уровень подтоварной жидкости количество, закаченной жидкости и количество израсходованной жидкости, а масса или объем оставшегося продукта определяется устройством обработки как разность между текущими показаниями уровнемеров или входного и расходного расходомеров с учетом перечисленных измеренных параметров.

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

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

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

Поставлена задача: повысить точность определения остаточной массы жидкого товарного продукта в резервуаре при изменяющихся внешних и внутренних условиях.

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

$$M_X=b_1{q}_{\mathrm{1,}X}+\dots +b_k{q}_{k,X}+\dots +b_m{q}_{m,X}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k,X},\left(1\right)}$$

где M_X - искомое значение массы контролируемой жидкости, q_k,X - совокупность измеренных параметров контролируемой жидкости. При этом коэффициенты b_k определяются по формуле:

$$b_k=\frac{\Delta b_k}{\Delta },\left(2\right)$$

где $$\Delta =\left|\begin{array}{c}{q}_{\mathrm{1,1}}\cdots q_{k\mathrm{,1}}\cdots q_{n\mathrm{,1}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}i}\cdots q_{k,i}\cdots q_{n,i}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}n}\cdots q_{k,n}\cdots q_{n,n}\end{array}\right|$$ , $$\Delta b_k=\left|\begin{array}{c}{q}_{\mathrm{1,1}}\cdots M_1\cdots q_{n\mathrm{,1}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}i}\cdots M_i\cdots q_{n,i}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}n}\cdots M_n\cdots q_{n,n}\end{array}\right|$$ .

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

q_k,i - совокупность физико-химических параметров жидкости, измеренных при калибровке, k - соответствует номеру измеряемого параметра (i=1…n), i - номеру калибровочного замера (k=1…n).

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

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

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

$$M=F\left(q_1,\cdots q_k,\cdots q_n\right),\left(3\right)$$

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

В общем случае функция F неизвестна. Ее поиск является важнейшей задачей математического описания измерительного процесса. Учитывая, что любой измерительный процесс по определению подразумевает в себе, кроме сбора и обработки информации, также операцию калибровки по эталонам, для решения поставленной задачи предложен алгоритм измерений, показанный на фигуре 1. Сущность алгоритма состоит в том, что для нахождения функции преобразования при неизвестной функции F составляется калибровочная модель процесса измерения. Для этого в резервуар i - раз поэтапно наливают количество жидкости известной массы M_i, измеряют соответствующую ей совокупность контролируемых физико-химических параметров q_k,i, k - соответствует номеру измеряемого параметра (k=1…m), i - номеру калибровочного замера (i=1…n), значения которых запоминаются в центральном устройстве обработке. Количество эталонных замеров n должно быть не меньше числа контролируемых параметров m, то есть n≥m.

Пусть имеется n калибровочных наливов жидкости с известными значениями массы M₁,…M_i,…M_n. При этом следует помнить, что M_i=M_i-1+ΔM, где ΔM - фиксированное приращение массы на каждом этапе калибровки. При каждом наливе контролируются совокупность параметров мониторинга резервуара q_k,i. Отметим, что при каждом очередном наливе физико-химические параметры жидкости будут меняться, так как процесс налива при больших объемах резервуара длительный, то будет меняться температура Т, соответственно плотность ρ и вязкость η основного и подтоварного продукта.

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

$$\begin{array}{c}{M}_1=b_1{q}_{\mathrm{1,1}}+\cdots +b_k{q}_{k\mathrm{,1}}+\cdots +b_m{q}_{m\mathrm{,1}}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k\mathrm{,1}}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ M_i=b_1{q}_{\mathrm{1,}i}+\cdots +b_k{q}_{k,i}+\cdots +b_m{q}_{m,i}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k,i}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ M_m=b_1{q}_{\mathrm{1,}m}+\cdots +b_k{q}_{k,m}+\cdots +b_m{q}_{m,m}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k,m}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ M_n=b_1{q}_{\mathrm{1,}n}+\cdots +b_k{q}_{k,n}+\cdots +b_m{q}_{m,n}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k,n}.}\end{array}\left(4\right)$$

Система имеет n уравнений для m неизвестных b_k. Так как число эталонных замеров всегда можно сделать больше числа контролируемых параметров n≥m, то из системы (4) всегда можно выбрать достаточное количество уравнений, из которых можно найти коэффициенты b_k. Желательно, чтобы выбранные для совместного решения уравнения охватывали возможно больший диапазон контролируемых параметров. Решение этой системы дает совокупность коэффициентов b_k, которые в средневзвешенной форме определяют искомый показатель качества. Система (4) имеет единственное решение, если главный определитель отличен от нуля. Решением этого уравнения являются соотношения (2). Учитывая независимость получения измерительной информации и возможность варьирования сигналов датчиков всегда можно исключить равенство нулю главного определителя в заданном диапазоне измерения.

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

Численный эксперимент показывает, что при любых значениях min<q_k,i<max искомое значение находится в интервале M_min<M_x<M_max, что полностью подпадает под определение измерительного процесса, как фактора уменьшения неопределенности. Методом численного эксперимента доказано, что при любом количестве калибровочных отсчетов, равных числу используемых для мониторинга резервуара параметров, искомый параметр будет находиться в интервале M_min÷M_max. Таким образом, предложенный способ измерений не требует точного знания функции, связывающей искомую массу жидкости с измеряемыми косвенными параметрами. Он предполагает создание математической модели в процессе калибровки. Чем больше параметров мы измеряем, и чем больше число калибровочных отсчетов, тем выше точность измерений. Одним из измеряемых параметров может и должна быть температура.

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

$$M_X=S{H}_X{\rho }_X,\left(6\right)$$

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

При вычислении массы по патентуемой методике для рассматриваемого случая калибровку достаточно проводить по двум эталонным массам, полностью перекрывающим диапазон ожидаемых значений M_min÷M_max. В формуле (2) будем иметь M₁=M_min, M₂=M_max, q₁₁=H₁, q₁₂=H₂, q₂₁=ρ₁, q₂₂=ρ₂. Тогда выражение (1) примет вид:

$$M_X=b_1{H}_X+b_2{\rho }_X,\left(7\right)$$

где $$b_1=\frac{{\rho }_2{M}_1-{\rho }_1{M}_2}{H_1{\rho }_2-H_2{\rho }_1},$$ $$b_2=\frac{H_1{M}_2-H_2{M}_1}{H_1{\rho }_2-H_2{\rho }_1}.\left(8\right)$$

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

$$b_1=\frac{M_1-M_2}{H_1-H_2},$$ $$b_2=\frac{H_1{M}_2-H_2{M}_1}{\rho \left(H_1-H_2\right)}.\left(9\right)$$

$$M_X=b_1{H}_X+b_2{\rho }_X=\frac{M_1-M_2}{H_1-H_2}{H}_X+\frac{H_1{M}_2-H_2{M}_1}{\rho \left(H_1-H_2\right)}{\rho }_X\left(10\right)$$

Для подтверждения правильности предложенного способа, покажем, что при ρ_X=ρ выражения (6) и (10) идентичны. Учитывая, что M₁=ρSH₁, M₂=ρSH₂ и подставляя эти выражения в (10), получим (6).

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

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

$$M_X=b_1{H}_X+b_2{\rho }_X+b_3,T_X\left(11\right)$$

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

$$\Delta =\left|\begin{array}{ccc}{H}_1& {\rho }_1& T_1\\ H_2& {\rho }_2& T_2\\ H_3& {\rho }_3& T_3\end{array}\right|$$ , $$\Delta b_1=\left|\begin{array}{ccc}{M}_1& {\rho }_1& T_1\\ M_2& {\rho }_2& T_2\\ M_3& {\rho }_3& T_3\end{array}\right|$$

$$\Delta b_2=\left|\begin{array}{ccc}{H}_1& M_1& T_1\\ H_2& M_2& T_2\\ H_3& M_3& T_3\end{array}\right|$$ $$\Delta b_3=\left|\begin{array}{ccc}{H}_1& {\rho }_1& M_1\\ H_2& {\rho }_2& M_2\\ H_3& {\rho }_3& M_3\end{array}\right|$$

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

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

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

Способ измерения массы жидкости в резервуаре, основанный на калибровке и заключающийся в том, что с помощью датчиков измеряются параметры, характеризующие физико-химические характеристики хранимой жидкости, величины которых на каждой емкости передаются в базовый блок, связанный информационным каналом передачи данных с центральным устройством обработки информации, отличающийся тем, что в процессе калибровки в резервуар поэтапно наливают количество жидкости известной массы, измеряют соответствующую ей совокупность контролируемых физико-химических параметров, значения которых запоминаются в центральном устройстве обработки, а в процессе измерения искомую массу определяют по формуле:
$$M_X=b_1{q}_{\mathrm{1,}X}+\dots +b_k{q}_{k,X}+\dots +b_n{q}_{n,X}={\displaystyle \sum _{k=1}^m{b}_k{q}_{k,X}}$$ ,
где М_х - искомое значение массы контролируемой жидкости, q_k,X - совокупность измеренных параметров контролируемой жидкости, при этом коэффициенты b_k определяются по формуле:
$$b_k=\frac{\Delta b_k}{\Delta }$$ ,
где $$\Delta =\left|\begin{array}{c}{q}_{\mathrm{1,1}}\cdots q_{i\mathrm{,1}}\cdots q_{n\mathrm{,1}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}k}\cdots q_{i,k}\cdots q_{n,k}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}n}\cdots q_{i,n}\cdots q_{n,n}\end{array}\right|$$ , $$\Delta b_k=\left|\begin{array}{c}{q}_{\mathrm{1,1}}\cdots M_1\cdots q_{n\mathrm{,1}}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}k}\cdots M_k\cdots q_{n,k}\\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ q_{\mathrm{1,}n}\cdots M_n\cdots q_{n,n}\end{array}\right|$$ .
M_k - совокупность эталонных значений массы контролируемой жидкости, используемой при калибровке,
q_i,k - совокупность физико-химических параметров жидкости, измеренных при калибровке, i соответствует номеру измеряемого параметра (i=l…n), k - номеру калибровочного замера (k=l…n).