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

ТРАНСПОРТ И ХРАНЕНИЕ НЕФТИ И ГАЗА

УДК 574+502.1:622.692.48

В.В. Пшенин1, e-mail: [email protected]; Н.А. Зарипова1, e-mail: [email protected]; А.К. Николаев1, e-maii: [email protected]

1 Федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский горный университет» (Санкт-Петербург, Россия).

Моделирование процесса вытеснения нефтяных паров при наливе железнодорожных цистерн

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

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

Ключевые слова: испарение нефти, налив нефти, железнодорожные цистерны, потери нефти, установки рекуперации паров.

V.V. Pshenin1, e-mail: [email protected]; N.A. Zaripova1, e-mail: [email protected]; A.K. Nikolaev1, e-mail: [email protected]

1 State Federal-Funded Educational Institution of Higher Professional Training «Saint-Petersburg Mining University» (Saint Petersburg, Russia).

Modeling of the Process of Oil Vapors Displacement during Railroad Tanks Loading

The article proposes the mathematical model of the oil vapor displacement process in loading of tank cars which allows approximate estimation of hydrocarbon emission parameters with minimum controlled variables and under operational production project. The Model is not hard binding to loading operations of tank cars and characterized by sufficient commonness which makes its application in loading of any tank cars possible taking into account a number of limitations in estimating emissions. The research carried out being of current importance is conditioned by considerable hydrocarbon emissions resulting from oil loading operations of tank cars. Intensive oil (product) evaporation in loading causes air pollution and losses of a valuable product. Vapor regeneration units are designed to solve the problem but their introduction into loading facilities is limited by their high price, operation complexity and, generally, by the absence of a positive economic effect. The situation being aggravated can be explained by incomplete or no control over the hydrocarbon emission quantitative variables even at the facilities equipped with the above units.

The central objective of the Research was to find the ways of higher environmental safety and performance optimization of vapor regeneration units at operational oil tank loading facilities in view of the current level of automation and

78

№ 6 июнь 2019 ТЕРРИТОРИЯ НЕФТЕГАЗ

OIL AND GAS TRANSPORTATION AND STORAGE

instrumentation. It has been established that application of the Model proposed with regard to installation of the oxygen analyzer at the line of controlled vapor removal from loading tank cars allows dynamic monitoring of hydrocarbon emission variables. The Model efficiency has been proved by numerical calculations of changes in volume concentration of air-gas components and the excess factor during a loading operation on the railway route.

Keywords: evaporation of oil, oil loading, railway tanks, oil losses, vapor recovery units.

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

1 - подводящий трубопровод; 2 - площадка обслуживания; 3 - стояк налива; 4 - рукав отвода паров; 5 - цистерна; 0н - объемный расход налива нефти, м3/ч; 0ГВС - неучтенный объемный расход газовоздушной смеси, м3/ч; QrBC - объемный расход газовоздушной смеси, м3/ч; P - давление, Па; J - плотность массы потока газовоздушной смеси, кг/м3

Fig. 1. Basic diagram of the system to remove the resultant vapor phase from the oil loading tank car: 1 - supply line; 2 - operating platform; 3 - loading pipe riser; 4 - vapor removal tube; 5 - tank; QB - volume oil loading flow rate, m3/h; 0ГВС - unaccounted volume air-gas flow rate, m3/h; QrBC - volume air-gas flow rate, m3/h; P - pressure, Ра; J - air-gas flow mass density, kg/m3

1

Обеспечению промышленной и экологической безопасности объектов налива железнодорожных цистерн, морских нефтеналивных терминалов,станций автоналива как ключевых элементов нефтетранспортной системы уделяется большое внимание [1-3]. Экологические риски при эксплуатации подобных объектов сопряжены прежде всего с загрязнением атмосферного воздуха за счет эмиссии углеводородов при заполнении транспортных емкостей. Точное определение количественных параметров эмиссии затруднено ввиду большого количества влияющих на эти параметры факторов (тип налива, условия окружающей среды, физико-химические свойства продукта, конструктивные особенности и т. д.), часть из которых носит случайный или сложнопрогно-зируемый характер. В российской практике при расчете потерь от испарения в пунктах налива железнодорожных цистерн принято пользоваться формулами П.Р. Ривкина [5-7], к числу достоинств которых можно отнести относительную простоту. В то же время данные формулы недостаточно точны и опираются на эмпирические коэффициенты, не имеющие достаточного физического обоснования и полученные на основе весьма ограниченной выборки. В зарубежной практике в основном применяется нормативный подход [8]. Отсутствие должного внимания к данной проблематике в последние десятилетия связано с отсутствием запроса на исследования - сложившаяся практика инструментальных замеров и расчетов по упрощенным методикам оказалась крайне удобной для всех участников производственного процесса.

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

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

Ссылка для цитирования (for citation):

Пшенин В.В., Зарипова Н.А., Николаев А.К. Моделирование процесса вытеснения нефтяных паров при наливе железнодорожных цистерн // Территория «НЕФТЕГАЗ». 2019. № 6. C. 79-82.

Pshenin V.V., Zaripova N.A., Nikolaev A.K. Modeling of the Process of Oil Vapors Displacement during Railroad Tanks Loading. Territorija «NEFTEGAS» [Oil and Gas Territory]. 2019;6:79-82. (In Russ.)

TERRITORIJA NEFTEGAS - OIL AND GAS TERRITORY No. 6 June 2019

79

ТРАНСПОРТ И ХРАНЕНИЕ НЕФТИ И ГАЗА

Кислород 02 Oxygen 02

- Углеводороды Ct-C1( Carbohydrates Ц-С^

— —-------—

О ОД 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0

Относительное время налива т Nondimensional time of loading т

Рис. 2. Изменение объемных долей компонентов газовоздушной смеси за время операции налива железнодорожного маршрута

Fig. 2. Volume concentration changes of air-gas components during a loading operation on the railway route

МОДЕЛИРОВАНИЕ ПРОЦЕССА ВЫТЕСНЕНИЯ ПАРОВ НЕФТИ ПРИ НАЛИВЕ ЖЕЛЕЗНОДОРОЖНОЙ ЦИСТЕРНЫ

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

Изменение во времени массы воздуха в газовом пространстве (ГП) можно рассчитать с помощью дифференциального уравнения:

с//77

где V - начальный объем газового

" нач

пространства, м3 (поскольку железнодорожные цистерны заполняются практически полностью, можно принять допущение Кач = 0н^н); 0н - объемный расход налива нефти, м3/ч; V - объем газового пространства, м3; Ь - время, прошедшее от начала налива, ч; тн -продолжительность операции налива, ч; т - относительное время операции налива, безразмерное. Величины С и Смв в уравнении (1) связаны соотношением:

Ргвс'

где рГВС - плотность ГВС в цистерне, кг/м3, при давлении Р, Па. Плотность ГВС при давлении Р и при атмосферном давлении Ратм определяется соответственно по формулам:

Ргвс = -

Р-Мг

ГВС

RT

Р ■М

РгВС(а) - ^ j ,

(4)

(5)

где рГВС, рГВС(а) - плотность ГВС при давлении Р и при атмосферном давлении Ратм соответственно, кг/м3; МГВС - молярная масса газовоздушной смеси, кг/кмоль; R - универсальная газовая постоянная 8314 Дж/(кмоль.К); T - температура ГВС, К.

Поскольку продолжительность налива железнодорожного маршрута относительно невелика, изменением температуры ГВС допустимо пренебречь T~ const, или dT/dx « 0. Если ввести параметр относительного давления, y « Р/Ратм, то из уравнений (4) и (5) следует РгВ/РГВОД = Р/Ратм. Коэффициент превышения, или коэффициент относительного расхода [4], K по определению равен отношению объемного расхода ГВС при давлении в конечной точке отвода паров (в рассматриваемом случае - при атмосферном давлении Ратм) к объемному расходу налива:

(3) к-

Q-BC(a)+QrBC(a)

о, ,

(6)

где тв - масса воздуха в ГП, кг; Смв -массовая концентрация воздуха в газовоздушной смеси (ГВС), д. ед., равная отношению массы воздуха к общей массе ГВС в ГП; СГВС - массовый расход ГВС, поступающий на УРП, кг/ч; бгвс - неучтенный массовый расход ГВС, обусловленный негерметичностью системы отвода паров (иными словами, расход утечки), кг/ч. Массовое содержание воздуха в единице объема газовоздушной смеси Св, кг/м3, по определению равно:

1,4 i V £ 1,35 -x i ¡5 1-3 S 3 1,25 ■ si с С 12- X <D £ ™ 1,15 t § 1,1 --в-"-1 о 1,05 1

0

0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 Относительное время налива т Nondimensional time of loading т

т.

т.

т.

V (2)

Рис. 3. Изменение коэффициента превышения за время операции налива железнодорожного маршрута

Fig. 3. Excess factor changes during a loading operation on the railway route

80

№ 6 июнь 2019 ТЕРРИТОРИЯ НЕФТЕГАЗ

OIL AND GAS TRANSPORTATION AND STORAGE

где 0ГВС(а) - объемный расход ГВС при атмосферном давлении, м3/ч; ОГВС(а) - неучтенный объемный расход ГВС при атмосферном давлении, обусловленный негерметичностью системы отвода паров, м3/ч.

Подставив (2) и (3) в уравнение (1), после алгебраических преобразований получаем:

dC, di

(1 -т) = С.

1 У

(7)

Массовое содержание воздуха в единице объема ГВС определяется по формуле:

Св Рв 5Св У Рв(a) ^

(8)

где рв, рв(а) - плотность воздуха при давлении Р и при атмосферном давлении Р соответственно, кг/м3; у - объем-

атм ' Л<в

ная доля воздуха в ГВС, д. ед. Подставив (7) в (8), получим уравнение для мгновенно значения коэффициента превышения в общем виде:

¿Мхв-У))

К = у

1-(1-т )•

dt

Плотность воздуха при атмосферном давлении Ратм = 101325 Па и фиксированной температуре постоянна. Например, при Ратм = 101325 Па и Т = 293 К она составляет р « 1,2 кг/м3. Объемная доля кислорода в воздухе составляет 20,9476 %. С учетом этих допущений можно записать (8) в виде:

Св = Рв^в = 4,7738 У * РвЫ* ,,

(10)

где ^ - объемная доля кислорода в ГВС, д. ед.2

С учетом (10), сократив константы в левой и правой части уравнения (7), получим:

с/т

а-тИ^)!1—) (11)

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

(9) К = у

(1-т)

1-

'd%o, . dy

dx

Хо-У

. (12)

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

■3 атм

относительного давления и dy/dт « 0 и у « 1, и уравнение (12) значительно упростится:

(1-х)

К = 1-

dx

5Со,

4пМ

di

(13)

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

XXVI специализированная выставка

Кшшж, I марка

АО «Тюменская ярмарка»

Адрес: Россия, 625013, г. Тюмень, ул. Севастопольская, 12, Выставочный зал телефакс: (3452) 48-55-56, 48-66-99, 48-53-33; e-mail: [email protected]. www.expo72.ru

2019

ТРАНСПОРТ И ХРАНЕНИЕ НЕФТИ И ГАЗА

ЧИСЛЕННЫЕ ПРИМЕРЫ РАСЧЕТА

Приведем численный пример расчета коэффициента превышения по разработанной методике. В качестве исходных данных приняты экспериментальные значения, полученные при наливе железнодорожного маршрута на действующей нефтеналивной эстакаде. На рис. 2 представлен характер изменения объемных долей компонентов ГВС за время операции налива, принятый по показаниям поточного кислородо-мера на трубопроводе организованного отвода газовой фазы. На рис. 3 представлен график изменения коэффициента превышения за время операции налива железнодорожного маршрута, рассчитанный по формуле (13) с использованием приемов дифференцирования функций,заданных таблично. Из рис. 3 видно, что характер кривой для коэффициента превышения во многом совпадает с теоретическими

ожиданиями. Функция принимает значение 1 в начале и конце операции налива, при этом она имеет явно выраженный максимум (Ктах = 1,31). Максимум функции смещен вправо, поскольку для рассматриваемых условий интенсивное испарение углеводородов начинается приблизительно с середины операции налива. В рассматриваемом примере в точке максимума объемный расход ГВС составлял 2620 м3/ч при расходе налива 2000 м3/ч, среднее интегральное значение коэффициента превышения составило К = 1,12. По графику ясно видно, сколь существенное влияние оказывает явление превышения на общий объем вытесняемой из железнодорожных цистерн ГВС.

ВЫВОДЫ

1. Объемный расход ГВС, вытесняемый из железнодорожных цистерн, может существенно (до 30 % в рассматри-

ваемом случае) превышать объемный расход наливаемой нефти (нефтепродукта).

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

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

4. Уравнения (12) и (13) имеют при учете некоторых допущений ряд аналитических решений, однако поскольку они представляют научный интерес, их целесообразно рассмотреть в следующих работах.

References:

1. DavLetjarov R.R., Lapina O.A., Polovkov S.A. Ensuring Environmental Safety around the Oil Loading Port Primorsk. Nauka i tekhnologii truboprovodnogo transporta nefti i nefteproduktov [Science & Technologies: Oil and Oil Products Pipeline Transportation]. 2015;2(18):84-87. (In Russ.)

2. Vykhodtseva N.A., Kostina E.A., Ukolova E.S. Biomonitoring of the Offshore Area of Koz'mino Oil Port Bay in the Area of the Oil Loading Terminal of Spetsmornefteport Kozmino LLC. Nauka i tekhnologii truboprovodnogo transporta nefti i nefteproduktov [Science & Technologies: Oil and Oil Products Pipeline Transportation]. 2014;1(13):86-90. (In Russ.)

3. Vakaev A.Y., Roslyakov D.A., Kinyashov M.A. Analysis of the Current State and Prospects of Reconstruction of Moscow Region Transneft JSC Loading Depots. Nauka i tekhnologii truboprovodnogo transporta nefti i nefteproduktov [Science & Technologies: Oil and Oil Products Pipeline Transportation]. 2014;3(15):80-85. (In Russ.)

4. Pshenin V.V., Gaysin M.T. Gas Dynamics Modeling during Tanker Loading Operations. Gornyi informatsionno-analiticheskii byulleten' (nauchno-tekhnicheskii zhurnal) [Mining Informational and Analytical Bulletin (Scientific and Technical Journal)]. 2017;S28:3-12. (In Russ.)

5. Rivkin P.R., Kurganskiy M.I., Bronshtein I.S. Vapor Losses when Loading Tank Cars in Different Ways. Neftyanoye khozyaistvo [Oil Industry]. 1974;2. (In Russ.)

6. Rivkin P.R. Investigation on Evaporation Losses of Oil and Its Products in Loading of Tank Cars. Thesis for a Candidate Degree in Engineering Sciences. Moscow, 1975. (In Russ.)

7. Rivkin P.R., Battalov A.Z., Bronshtein I.S., et al. On the Excess Factors of Vapor Volume Leaving the Tank over the Oil Volume| Pumped. In collected research papers "Technology of oil transfer, prevention of losses, automated pipeline transportation management system". Ufa: All-Russian Research Institute for Gathering, Treatment and Transportation of Oil and its Products. 1985. P. 75-79. (In Russ.)

8. Mulder T. VOC recovery systems. Hydrocarbon Engineering. 2007;12(6):37-40.

Литература:

1. Давлетяров Р.Р., Лапина О.А., Половков С.А. Обеспечение экологической безопасности в районе нефтеналивного терминала Порта Приморск // Наука и технологии трубопроводного транспорта нефти и нефтепродуктов. 2015. № 2 (18). С. 84-87.

2. Выходцева Н.А., Костина Е.А., Уколова Е.С. Биомониторинг акватории бухты Козьмина в районе нефтеналивного терминала «Спецморнефтепорт Козьмино» // Наука и технологии трубопроводного транспорта нефти и нефтепродуктов. 2014. № 1 (13). С. 86-91.

3. Вакаев А.Ю., Росляков Д.А., Киняшов М.А. Анализ современного состояния и перспектив реконструкции наливных пунктов московского региона ОАО «АК «Транснефть» // Наука и технологии трубопроводного транспорта нефти и нефтепродуктов. 2014. № 3 (15). С. 80-85.

4. Пшенин В.В., Гайсин М.Т. Моделирование газодинамических процессов при наливе танкеров // Горный информационно-аналитический бюллетень (научно-технический журнал). 2017. № S28. С. 3-12.

5. Ривкин П.Р., Курганский М.И., Бронштейн И.С. Потери нефтепродуктов от испарения при различных способах налива железнодорожных цистерн // Нефтяное хозяйство. 1974. № 2.

6. Ривкин П.Р. Исследование потерь нефти и нефтепродуктов от испарения при наливе в транспортные емкости: дисс. ... канд. техн. наук. М., 1975. 252 с.

7. Ривкин П.Р., Батталов А.З., Бронштейн И.С. и др. К определению коэффициентов превышения объема выходящей из резервуара паровоздушной смеси над объемом закачанной нефти // Технология перекачки нефти, предотвращение потерь, автоматизированная система управления неф-тепроводным транспортом: сб. научных трудов Всероссийского научно-исследовательского института по сбору, подготовке и транспортировке нефти и нефтепродуктов. Уфа, 1985. С. 75-79.

8. Mulder T. VOС recovery systems // Hydrocarbon Engineering. 2007. Vol. 12. Iss. 6. P. 37-40.

82

№ 6 июнь 2019 ТЕРРИТОРИЯ НЕФТЕГАЗ