CC BY
27
2014
НАУЧНЫЙ ВЕСТНИК МГТУ ГА
№ 206
УДК 544.03
ВОПРОСЫ ПОВЫШЕНИЯ ЭФФЕКТИВНОСТИ УЧЁТА УГЛЕВОДОРОДНЫХ ТОПЛИВ НА ПОЛЕВЫХ СКЛАДАХ
Ю.Н. РЫБАКОВ, С.И. ЧИРИКОВ, Р.И. КЮННАП
В статье изложены теоретические подходы к оценке технологических потерь нефтепродуктов при хранении в эластичных резервуарах с учётом стационарного и нестационарного режимов диффузии топлив через конструкционный материал эластичного резервуара. Предложенный алгоритм и математические модели определения технологических потерь, с учётом нестационарного режима изменения топливопроницаемости конструкционного материала эластичных резервуаров, позволяют существенно (до 40%) повысить точность нормирования естественной убыли углеводородных топлив на полевых складах горючего.
Ключевые слова эластичные резервуары, полевые склады горючего, естественная убыль, технологические потери, моторные и авиационные топлива, коэффициент диффузии, алгоритм и математические модели, нормирование.
Большинство известных работ, посвященных оценке естественной убыли при хранении в эластичных резервуарах, используют стационарную модель процесса диффузии углеводородов. Для получения более точных оценок технологических потерь нефтепродукта необходимо учитывать нестационарность данного процесса, т.е. рассматривать данное явление при переменной не только по толщине, но также и во времени концентрации нефтепродукта в эластичной оболочке. Это позволит осуществить более точный прогноз потерь углеводородов при их хранении в эластичных резервуарах.
Решение задач по определению нормативов технологических потерь (естественной убыли) нефтепродуктов за счет проницаемости нефтепродуктов через оболочки резервуаров может быть эффективно осуществлено с использованием дифференциального уравнения диффузии без набухания
дС (Х,тТ = с11у( - Б ■ gradC (х, т)) . (1)
дт
При постоянном коэффициенте диффузии В=еоп81 имеем уравнение
— = Б ■АС
дт . (2)
При решении практических задач чаще всего ограничиваются случаем одномерных моделей, используя известное дифференциальное уравнение параболического типа [1; 2]. Считая внешние размеры резервуара существенно большими толщины его стенок, выберем модель диффузионного процесса в виде неограниченной пластины толщиной ё. С учетом отсутствия источников массы внутри эластичного материала можно применить дифференциальное уравнение, обычно называемое уравнением теплопроводности [3]
дС = Б д^С
дт дх2
где т - время протекания диффузионных процессов; х - координата, направление которой совпадает с внешней нормалью к эластичной поверхности.
= Б ■—, (3)
Для количественного решения задачи определения естественной убыли нефтепродукта вследствие диффузии из эластичного резервуара уравнение (3) необходимо дополнить начальными и граничными условиями. Очевидно, что начальные условия в момент времени т = 0 следует принять в виде
С(0,х)=0 при 0 < X < /, что соответствует отсутствию молекул нефтепродукта в эластичном материале.
Для задания двух граничных условий в точках х = 0 и х = / примем следующие допущения:
1. В точке х = 0 в месте контакта границы жидкого объема с эластичной поверхностью во все моменты времени т >0 считаем концентрацию С(х, т ) постоянной величиной, равной максимально возможной при данной температуре
С (*,т)| х=0 = Со • (4)
2. В точке X = /, соответствующей внешней оболочке эластичной поверхности, принимаем граничные условия в виде
- ^ = ^Со, (5)
^X=/
считая плотность диффузионного потока на внешней границе оболочки постоянной величиной, пропорциональной максимально возможной концентрации нефтепродукта в эластичном материале - С0 (А = 0 - случай непроницаемой эластичной оболочки).
Для этого определения численных значений параметров С 0 и А в формулах (4) и (5) необходимо использовать экспериментальные данные по накоплению диффузанта (нефтепродукта) в эластичной оболочке. Т.е. необходимо проведение двух экспериментов по измерению концентрации при двух различных толщинах эластичных оболочек. При этом толщина этих двух оболочек должна отличаться в 2 раза.
Учитывая граничные условия (4) и (5) решение уравнения (3) можно записать в виде
8/ 2п +1
С(X, т) = С0 01 - Ах) - — X [ А + (-1)п -— ^0 ] •
П 1 2/
2п +1 Бт С08^П+1 П(/ - х) (6)
.е - (^ п)2 • Т--21-
2/ /2 2п +1
Задавая различные значения элементов времени т, можно получить значения концентрации нефтепродукта С в различных точках эластичной оболочки.
Наиболее важным является определение функции изменения во времени значения концентрации С (т) в точке х = /, т.е. на внешней границе оболочки. Тогда решение (6) примет вид
8/ х 2п +1
С(/т) = 0,(1 - А/) - -2 X [А + (-1)п —— пП
П 1 2/
1 2п+1 2 От (7)
--^ в"^ Т
(2 п +1)2
Данное решение С (/,х) может быть использовано при определении массы вышедшего нефтепродукта из объема, ограниченного эластичной поверхностью при нестационарном режиме диффузии.
Для определения последней величины примем, что стационарное значение концентрации
С(I,т) при т-соответствует известной величине I - массе вышедшего нефтепродукта в
единицу времени с единицы поверхности, измеряемой экспериментально в течение значительного промежутка времени (более 5 сут.).
Промежуточные во времени т (0<т<да) значения величины ¡(т) будут пропорциональны концентрации С (I, т) в тот же момент времени т. Таким образом, текущее значение плотности диффузионного потока ¡(т) будет определяться формулой
Кг)=сит±. (8)
С(1,») ' '
Значение массы нефтепродукта, вышедшей из единицы поверхности эластичного резервуара за время т, может быть определено в виде
т
| i (т)й?т, (9)
где т - топливопроницаемость через конструкционный материал, кг/м2.
Для определения зависимости во времени концентрации и плотности диффузионного потока на внешней границе эластичной оболочки воспользуемся формулой (7). Задавая различные значения т, получим соответствующие значения концентрации С (I, т) (таблица). р = 750 кг/м3, I = 3 мм, г = 20 0С, 1о = 145 г/(м3сут.),
Б = 6, 71- 10-12 м2/с = 5,8 • 10-7 м2/сут.
Таблица
Исходные данные и результаты расчета технологических потерь в нестационарном режиме диффузии
Бт
Т~ I / 1с
0,025 0,092
0,05 0,710
0,075 0,831
0,1 0,907
0,15 0,58
1
На рисунке приведен график зависимости безразмерной величины 1/10 (10 - значение диффузионного потока в стационарном режиме) от безразмерного времени.
Рисуиок иллюстрирует зависимость От при рассмотрении диффузии бензина АИ-92 через конструкционный материал эластичного резервуара типа МР [4; 5; 6]. Безразмерное время протекания неустановившегося процесса диффузии, равное 0,2, соответствует реальному времени около 3 сут.
Оценку массы вышедшего из эластичного резервуара нефтепродукта рассчитывают по формуле (9).
О 0,05 ОД 0Д5 0,2 0,25 1
Безразмерное время
Рисунок. Динамика технологических потерь при нестационарном режиме диффузии углеводородов
Промежуток времени определяют из расчётов по предыдущему пункту. Очевидно, что это время т - протяженность нестационарного процесса диффузии или время стабилизации диффузии. Интегрируя по (9) с учётом формул (8) и (7), получаем значение для рассматриваемого случая (материал - МР, бензин - АИ-92).
После интегрирования получаем значение убывшей массы в нестационарном режиме за безразмерное время Вт /I2 = 0,2, равное 0,5326.
Это значение убывшей массы вещества удобно выразить в долях от массы, убывающей за аналогичный период времени, считая процесс установившимся, с постоянной величиной диффузионного потока.
Резюмируя изложенное, можно констатировать, что без учёта нестационарности процесса диффузии оценка количества вышедшего нефтепродукта оказывается завышенной. Реальная убыль нефтепродукта при протекании нестационарного режима оказывается немногим более половины (56,5%) от убыли, посчитанной по известным методикам для стационарных режимов. Таким образом, неучёт этого фактора может очень сильно отразиться на результате для относительно коротких сроков хранения (до 15 сут.) нефтепродуктов в эластичных резервуарах.
Реальное время, соответствующее безразмерному времени Бт/I2 = 0,2, равно 3,1 сут. За 6,2 сут. выйдет около 80% от рассчитанной величины по стационарному, за 9,3 сут. - 86%.
Предложенный алгоритм и математические модели определения технологических потерь, с учётом нестационарного режима измерения топливопроницаемости конструкционного материала эластичных резервуаров позволяют существенно (до 40%) повысить точность нормирования естественной убыли углеводородных топлив на полевых складах горючего.
ЛИТЕРАТУРА
1. Рудобашта С.П., Карташов Э.М. Диффузия в химико-технологических процессах. - 2-е изд., перераб. и доп. - М.: КолосС, 2010.
2. Карслоу Г., Егер Д. Теплопроводность твердых тел. - М.: Наука, 1964.
3. Васин С.И., Филиппов А.Н. Проницаемость сложнопористых сред // Коллоидный журнал. - 2009. - Т.71. - № 1. - С. 32-46.
4. Отчет о НИР Исследование по разработке норм естественной убыли нефтепродуктов при хранении и транспортировке. - ФАУ «25 ГосНИИ химмотологии Минобороны РФ». Шифр «Свирск-4», 2009.
5. Пат. Российская Федерация. Автоматизированная установка определения массы нефтепродукта в эластичном резервуаре при хранении / Рыбаков Ю.Н., Ванчугов Н.А., Харламова О.Д., Пирогов Ю.Н. - № 111639; заявл. 31.08.2011; опубл. 20.12.2011; Б.И. № 35.
6. Пат. Российская Федерация. Способ определения массы нефтепродукта при хранении в эластичном резервуаре / Рыбаков Ю.Н., Ванчугов Н.А., Харламова О.Д., Пирогов Ю.Н. - № 2470264; заявл. 20.07.2011; опубл. 20.12.2012; Б.И. № 35.
IMPROVING OF HYDROCARBON FUELS ACCOUNTING IN FUEL FARMS
Ribakov Y.N., Chirikov S.I., Kyunnap R.I.
The article outlines the theoretical approaches to the assessment of technological losses of petroleum products during storage in flexible tanks with flexible account of steady and unsteady regimes fuels diffusion through structural material of tanks. The proposed algorithm and mathematical model for determining technological losses, taking into account the unsteady measurement mode structural material elastic reservoirs, can significantly (up to 40 %) improve the accuracy of the valuation of natural loss of hydrocarbon fuels at fuel depots field.
Key words: flexible tanks, fuel dumps, wastage, technological losses, motor and aviation fuel, the diffusion coefficient, the algorithm and mathematical models, valuation.
Сведения об авторах
Рыбаков Юрий Николаевич, 1961 г.р., окончил МИНХ и ГП им. И.М. Губкина (1983), кандидат технических наук, старший научный сотрудник, заслуженный изобретатель РФ, начальник 23 отдела «25 ГосНИИ химмотологии Минобороны России», автор более 200 научных работ, область научных интересов - нефтепродуктообеспечение, полимерные материалы.
Чириков Сергей Игоревич, 1990 г.р., окончил МАТИ (2013), младший научный сотрудник 23 отдела «25 ГосНИИ химмотологии Минобороны России», автор 5 научных работ, область научных интересов - полимерные материалы.
Кюннап Роман Игоревич, 1989 г.р., окончил УВВТУ (ВИ) (2011), младший научный сотрудник 23 отдела «25 ГосНИИ химмотологии Минобороны России», автор 7 научных работ, область научных интересов - технические средства нефтепродуктообеспечения.
