Применение метода совокупно-косвенных измерений массы нефтепродуктов при мониторинге параметров резервуарных парков

УДК 681.2.083

ПРИМЕНЕНИЕ МЕТОДА СОВОКУПНО-КОСВЕННЫХ ИЗМЕРЕНИЙ МАССЫ НЕФТЕПРОДУКТОВ ПРИ МОНИТОРИНГЕ ПАРАМЕТРОВ

РЕЗЕРВУАРНЫХ ПАРКОВ

© 2013 А.В. Солнцева, С.А. Борминский, А.Н. Малышева-Стройкова, Е.А. Силов

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

Поступила в редакцию 04.10.2013

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

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

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

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

Учет нефтепродуктов может производиться в различных единицах, так, к примеру, на нефтебазах и наливных станциях учет ведут в единицах массы, а на автозаправочных станциях - в единицах объема.

Для обеспечения достоверности и единства измерений массы нефтепродуктов, а также контроля их качества нефтебазы и АЗС должны иметь необходимое оборудование и средства измерений, допущенные к применению Госстандартом и имеющие

Солнцева Александра Валерьевна, аспирант кафедры электротехники. E-mail: als063@mail.ru Борминский Сергей Анатольевич, кандидат технических наук, доцент. E-mail: b80@mail.ru

Малышева-Стройкова Александра Николаевна, аспирант кафедры электротехники. E-mail: aps@ssau.ru Силов Евгений Альбертович, кандидат технических наук, старший научный сотрудник НИЛ "Аналитические приборы и системы". E-mail: aps@ssau.ru

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

В соответствии с межгосударственным стандартом [1] приняты следующие основные определения:

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

- масса брутто товарной нефти - масса товарной нефти, показатели качества которой соответствуют требованиям ГОСТ Р 51858;

- масса балласта - общая масса воды, солей и механических примесей в товарной нефти;

- масса нетто товарной нефти - разность массы брутто товарной нефти и массы балласта.

Существуют следующие методики выполнения измерений массы нефти и нефтепродуктов:

- прямой метод динамических измерений;

- косвенный метод динамических измерений;

- прямой метод статических измерений;

- косвенный метод статических измерений;

- косвенный метод, основанный на гидростатическом принципе.

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

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

а) плотности с помощью поточных преобразователей плотности, давления и температуры;

б) объема продукта с помощью преобразователей расхода, давления и температуры или счетчиков жидкости [1].

Результаты измерений плотности и объема продукта приводят к стандартным условиям или результат измерений плотности продукта приводят к условиям измерений его объема.

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

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

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

- уровня продукта - стационарным уровнемером или какими-либо другими средствами измерений уровня жидкости;

- плотности продукта - переносным или стационарным средством измерений плотности или ареометром по ГОСТ 3900, ГОСТ 31072 или лабораторным плотномером в объединенной пробе, составленной из точечных проб, отобранных по ГОСТ 2517;

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

- объема продукта - по калибровочной таблице меры вместимости с использованием результата измерений уровня продукта;

б) в мерах полной вместимости (под мерой полной вместимости подразумевается средство измерений объема продукта, имеющее свидетельство о поверке и оснащенное указателем уровня наполнения):

- плотности продукта - переносным средством измерений плотности или ареометром в лаборатории по ГОСТ 3900, ГОСТ 31072 или лабораторным плотномером в точечной пробе продукта, отобранной по ГОСТ 2517;

- температуры продукта - переносным преобразователем температуры или термометром в точечной пробе продукта, отобранной по ГОСТ 2517;

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

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

- гидростатического давления столба продукта - стационарным измерителем гидростатического давления;

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

В соответствии с требованиями ГОСТ 8.587-

2006 предел допускаемой относительной погрешности измерений массы брутто товарной нефти и массы нефтепродукта при косвенном методе статических измерений продукта массой до 120 т должен составлять 0,65%, а при и измерении продукта массой свыше 120 т - 0,5% [1].

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

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

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

м = ^

где qk - совокупность измеряемых параметров.

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

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

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

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

1 г ' п

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

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

Mi = bi Ч1.1 + +bA л + •••+bmqm i = £

к=1

m

Mi = bl4l,i +  + ЬкЧкi +  + bm4m,i = £ bk4t,i

к=1

m

M = b,q, +  + bkqk +  + b q = £bkqk

m 1 m к lk, m m lm,m ^^ к J-k, n

к =1

m

M = b.q. +  + b.q. +  + b q = £ b.q.

n 1 H,n к J-к,n m ±m,n / i к J-к, n

q и qt. qn ,1  M   Чп л

где А = q 1,i qt q„,, , А bk = ЧЦ   Mr   Ч n,,

q 1, n qt, n qn, n ЧЩ  m„   q n ,n

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

АЬ,,

bk =

А '

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

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

m

Мх = ьл,х + +Kqk,x + -+bmqm, x=£hqk,x,

к =1

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

При любых значениях min<qk<max искомое значение находится в интервале М . <М <М , что пол-

1 min x max7

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

min max

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

Рассмотрим вариант измерение массы товарного продукта в резервуаре при условии контроля трех изменяющихся параметров: уровня заполнения резервуара, плотности и температуры продукта. Измерение уровня может быть проведено согласно [3]. Тогда неизвестная масса продукта будет измеряться согласно формуле:

Mx = bflx + b2Px + b3Tx,

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

А =

H1 P T M1 P1 T1

H 2 Pi Ti Ab1 = M 2 Pi Ti

H 3 Рз Тз M 3 P3 T3

к=1

H M T H P Mi

Ab2 = H2 M2 T2 3 = H2 A M2

H M3 T3 H3 P3 M3

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

Пример. Были проведены экспериментальные исследования в рамках подтверждения разработанной теории. Были проведены измерения массы товарного продукта - дизельного топлива на танкере 1М07611212 класса река-море. Данные приведены в табл. 1.

Таблица 1. Пример производимых изм

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

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

Исследования выполнены при поддержке Министерства образования и науки Российской Федерации

ений согласно разработанному методу

Данные калибровочной модели Измерение

Н, м p, кг/м3 Т, град С М, кг

H1=0.10 Pj=860 Tj=25.1 Mj=1832 Hx=0.14

H2=0.13 p2=860 T2=25.2 M2=2460 px=860

H3=0.15 p3=860 T3=25.1 M3=2898 Tx=25.1

При измерении, характеризующемся парамет- СПИСОК ЛИТЕРАТУРЫ рами, приведенными в столбце "Измерение", из

памяти калибр°вочн°й модели были взяты дан- t ГОСТ Р 8.595-2004. Масса нефти и нефтепродуктов.

ные, приведенные в первых трех столбцах. В ре- Общие требования к методикам выполнения изме-

зультате применения разработанного метода из- рений. М., 2005. II, 27 с. (Государственная система

мерения массы, было получено значение Mx=2685 обеспечения единства измерений).

кг. После анализа имеющейся калибровочной таб- 2. Методы измерений количественных и качественных

лицы было установлено, что разработанный ме- характеристик жидких энергоносителей / С.А.Бор-

тод дает более точный результат в определении минский, Б.В.Скворцов, А.В.Солнцева. Самара:

массы, чем использование заранее подготовлен- СНЦ РАН, 2012. 222 с.

ной калибровочной таблицы, в силу того, что дает 3. Скворцов Б.В., Блинов Д.И., Солнцева А.В., Борминс-возможность учесть быстроменяющиеся факторы, кий С.А. Устройство для измерения уровня жидков данном примере - температуру (по калибровоч- сти в резервуарах //Патент России № 115886. 2012. ной таблице Mx=2674 кг). На практике примене- Бюл. №13.

APPLICATION OF METHOD OF COLLECTIVELY INDIRECT MEASUREMENT OF MASS PETROLEUM PRODUCTS MONITORING PARAMETERS IN TANK PARK

© 2013 A.V. Solntceva, S.A. Borminsky, A.N. Malysheva-Strojkova, E.A. Silov

Samara State Aerospace University named after S.P. Korolyov

Propose a new method measuring the mass of liquid cargo in tank farms. The method consists in determining the desired mass based on the parameters characterizing the physical and chemical characteristics of liquid energy, which values are transferred to the base unit tank, a related information transmission channel with the central data processing device.

Keywords: storage tanks, measurement of mass, liquid energy carrier, liquid cargo, accuracy, calibration model, method

Alexandra Solntceva, Graduate Student at the Electrical Engineering Department. E-mail: als063@mail.ru Sergey Borminsky, Candidate of Technics, Associate Professor. E-mail: b80@mail.ru

Alexandra Malysheva-Strojkova, Graduate Student at the Electrical Engineering Department. E-mail: aps@ssau.ru Evgeny Silov, Candidate of Technics, Senior Research Fellow at the Research Laboratory "Analytical Devices and Systems". E-mail: aps@ssau.ru