CC BY
47
Научный журнал КубГАУ, №99(05), 2014 года
1
УДК 620
ОБОБЩЕНИЕ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ ПО ИССЛЕДОВАНИЮ БЕНЗИНОВЫХ НЕФТЯНЫХ ФРАКЦИЙ
Харченко Павел Михайлович к.т.н., доцент, доцент кафедры [email protected]
Т имофеев Виталий Павлович студент
Кубанский государственный аграрный университет, Краснодар, Россия
В статье рассматриваются результаты экспериментальных исследований бензиновых нефтяных фракций, их основные физикохимические свойства, плотность и ДНП, жидкая фаза, двухфазная область, критическая область и пограничные кривые
Ключевые слова: БЕНЗИНОВЫЕ ФРАКЦИИ, НЕФТЯНЫЕ ФРАКЦИИ, ПОГРАНИЧНЫЕ КРИВЫЕ, КРИТИЧЕСКАЯ ОБЛАСТЬ, МЕСТОРОЖДЕНИЯ, ПАРАФИН, ПРЕДЕЛ ВЫКИПАНИЯ, ПЛОТНОСТЬ, ИЗОТЕРМЫ
UDC 620
SUMMARY OF EXPERIMENTAL RESEARCH OF GASOLINE AND OIL FRACTIONS
Harchenko Pavel Mikhailovich
Candidate of engineering science, associate professor
Timofeev Vitaliy Pavlovich student
Kuban State Agrarian University, Krasnodar, Russia
In the article we present the results of pilot studies of petrol and oil fractions; their main physical and chemical properties, density and DNP, a liquid phase, two-phase area, critical area and boundary curves are considered
Keywords: PETROL FRACTIONS, OIL FRACTIONS, BOUNDARY CURVES, CRITICAL AREA, FIELDS, PARAFFIN, BOILING LIMIT, DENSITY, ISOTHERMS
1. Плотность (удельный объём) на линиях насыщения.
При обобщении зависимости плотности от температуры и давления предполагалось, что все исследованные нами вещества являются членами одного псевдогомологического ряда и для обобщения результатов исследования можно воспользоваться методами теории подобия. В качестве параметров приведения были выбраны критические значения плотности рКР, давления ?КР и температуры ТКР.
Предварительная обработка данных проводилась с использованием крупномасштабных графиков р-т, п-т. В координатах наблюдалось расслоение экспериментальных точек, соответствующих линиям насыщения (рис.1).
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
2
Рис. 1. - Обобщённая зависимость плотности нефтяных фракций на левой пограничной кривой.
Изучение характера расслоения привело к выводу, что основной причиной расслоения является ошибка в определении критических
параметров. Анализ зависимостей р(т) с учётом изменения ^ показал, что
критические параметры могут быть скорректированы в пределах погрешности их выделения, корректировка проводилась из условий минимизации среднеквадратического отклонения расчётных и экспериментальных значений плотности.
После уточнения критических параметров все экспериментальные точки удовлетворительно размещаются на одном графике (рис.2).
http ://ej. kubagro. ru/2014/05/pdf/62. pdf
Научный журнал КубГАУ, №99(05), 2014 года
3
Рис. 2. - Обобщённая зависимость плотности нефтяных фракций на левой пограничной кривой, после
уточнения критических параметров.
Правильность выбранного нами метода подтверждается совпадением в координатах р-т прямолинейных диаметров всех фракций.
Аналитическое описание результатов осуществлялось по методу наименьших квадратов уравнением вида:
= Ол +
+ а2Р + а3/>
Коэффициенты уравнения найдены путём обработки всего массива экспериментальных данных на линиях насыщения. Коэффициенты уравнения и среднеквадратическое отклонение от экспериментальных данных представлены в таблице 1.
http:// ej. kubagro .ru/2014/05/pdf/ 62. pdf
Научный журнал КубГАУ, №99(05), 2014 года
4
Таблица 1.
Коэффициенты уравнения (1) и среднеквадратическое отклонение от
экспериментальных данных
Названия ао аі а2 аз 8, %
Левая погр. -0,03253335 0,10896239 -0,13758248 0,06489381 0,08
кривая
Правая погр. -0,21225496 1,1107828 -2,0763058 1,3189137 0,9
кривая
5.2. Плотность (удельный объём) жидкой фазе.
Как указывалось в аналитическом обзоре, для расчёта плотности нефтепродуктов под давлением наиболее подходящим с точки зрения точности и простоты расчёта является уравнения Тэйта [2]. Значительную трудность представляет корректное выделение температурных коэффициентов K и L, что приводит к большим ошибкам расчёта плотности (1-1,5%). Поэтому было признано целесообразным для расчёта плотности нефтяных фракций в жидкой фазы вблизи линии насыщения использовать уравнение Тэйта вида:
VTJ> = VT,P0
Теоретическое обоснование уравнения Тэйта дано в работах [1,3], а рекомендации по его применению для нефтепродуктов приведены в [2]. Коэффициент A считался постоянным для всей группы веществ, независящим от температуры, B - температурная функция, независящая от давления. В качестве опорных значений uT, P0 были приняты значения
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
5
плотности на линии насыщения, соответствующие каждой изотерме. В ходе обработки всего массива экспериментальных данных из условий минимизации среднеквадратической ошибки расчёта были получены оптимальные значения коэффициентов A и B. Общая среднеквадратическая ошибка расчёта плотности по всем фракциям составила 0,3%.
Однако анализ зависимостей погрешности описания экспериментальных данных показал, что при приведённой температуре т<0,95 ошибка описания лежит в пределах погрешности экспериментальных данных, а в околокритической области возрастает до 0,7-0,8%. Вероятно при температурах, близких критической, коэффициент A тоже является температурной функцией, а коэффициент B начинает зависеть от давления и индивидуальных свойств веществ. Аналитическое описание коэффициента B от показало возрастание нелинейного характера этой зависимости при т>0,95 (рис.3).
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
6
50
40
30
°Д Г
Д □ -Ф-
6° о □ Д
t Д ° dP # 8
20
10
1,5
Рис.З. - Температурная зависимость коэффициента В.
Таким образом проводить экстраполяцию B по температуре (как указывалось в [3]) признано нецелесообразным. По нашим данным было получено новое значение т при B=0, равное 0,925. Путём обработки всего массива экспериментальных данных получено уравнение:
BCD
В =А+ — + — + —
A=41,4232; B=-153,817; C=137,208; D=-27,9074.
http://ej .kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
7
5.3. Давление в двухфазной области и на пограничной кривой.
Давление двухфазной области зависит от температуры и соотношения объёмов ф. Для нахождения аналитической зависимости давления от этих параметров на крупномасштабных графиках в P-V координатах были выделены линии равных значений ф (рис.4).
При температурах ниже 230oC значения ф рассчитывались по массе жидкости, слитой из пьезометра, при этом предполагалось, что рЖ>>рП. Равные значения ф наносились на соответствующие изотермы и соединялись плавной кривой. Причём при температурах выше 230OC проводилась экстраполяция линий ф-const вплоть до критической точки.
Были сделаны попытки обобщения давления на различных изотермах при соответствующих значениях ф с использованием в качестве параметров приведения различных комбинаций из критической плотности (критического объёма), критического давления и критической
http://ej .kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
8
температуры. Наиболее удачной представляется описание разности
давлений в виде Ал(ф)=^ф; л'-л''), где
а п
л" =
р' -р,г
КР
В результате обработки массива из 350 значений получено двухпараметрическое уравнение вида:
= + b1 In + с1.
+ Ъ-2
+ Cz
г'-и'О + Ъ I
rll\\3
■+
n
— n
Таблица 2.
Коэффициенты уравнения (4)
ао bi Ci Ь2 C2 Ьз C3
71,1959 0,456101 53,345 -0,0817 12,5252 0,00493 0,98661
Среднеквадратическая ошибка описания Ал(ф) по уравнению составляет 3% при максимальных ошибках, не превышающих 5%. В пересчёте на абсолютное давление это соответствует средней ошибке расчёта не превышающей 0,2%. Необходимые для расчёта по уравнению
(4) значения л' и л" получены в результате аналитического описания зависимости приведенного давления от приведённой температуры (л-т) на пограничной кривой. На рисунке 5 представлен вид этой зависимости для двух фракций различной ширины. Для удобства иллюстрации на графике допущено искажение масштаба, вследствие которого расстояние л'-л'' увеличено в 5 раз.
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
9
0,5 0,75 т
Рис.5. - Пограничные кривые нефтяных фракций.
1 - опорная кривая, 2 - пограничные кривые узкой фракции, 3 - пограничные кривые
широкой фракции.
Анализ графической зависимости п-т для всех фракций показывает, что линии п' и л" закономерно расслаивается. Графики п'-т, и п''-т широких фракций лежат выше соответствующих линий более узких фракций. Разность л'-л" для всех фракций изменяется в функции от т, монотонно возрастая практически от нуля при т=0,5, достигая максимума при 0,85<т<0,90 и, затем уменьшается до 0 при т=1.
При нахождении параметров позволяющих описать расслоение были испытаны фактор ацентричности и различные комбинации параметров состояния и физико-химических характеристик. Наилучшие результаты были получены при использовании в качестве корелирующего параметра относительно ширины фракции Дт = 1г>~* . На рисунке 6 показана
^КР
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
10
закономерность изменения п'-п'' для фракций с различной относительной шириной, из которой следует, что для некоторой гипотетической фракции с Ах=0, п'-п'' во всех точках должно быть равно 0.
Предположительно в качестве модели такой фракции можно испытать индивидуальные углеводороды, близкие к исследованным фракциям по температуре кипения и критическим температуре и давлению графическая проверка зависимости давления от температуры некоторых индивидуальных углеводородов, отвечающих указанным выше требованиям по данным [1] в координатах п-т показала, что они существенно расслаиваются между собой.
Поэтому на основе статистической обработки наших экспериментальных данных посредством экстраполяции п'-п'' на значения соответствующие Ат=0 позволило выделить на графике 5 линию п от т модельной фракции, которая была принята в качестве опорной кривой для
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
11
обобщения зависимости л’ от т и л’’ от т исследованных бензиновых фракций.
Точки опорной кривой были описаны уравнениями, применяющимися при описании давления от температуры для чистых веществ (уравнения Антуона, Кирхгофа, Риделя [2]). Однако при этом наблюдалось систематическое отклонение кривых, описанных по этим уравнениям от опорной кривой. Причём отклонения по мере увеличения т возрастали. Поэтому значения л от т на опорной кривой были описаны уравнением вида:
П
А +
где A=-56,237871; B=116,70118; C=-93,837626; D=34,374218; E=-56,690291.
Затем в координатах Дл-т были построены отклонения левых и правых пограничных кривых от опорной. График зависимости представлен на рисунке 7.
http: //ej .kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
12
Описание Ал от т проводилось отдельно для правой и левой пограничных кривых уравнением вида:
О О д д 1
Ал=а+Ь0т+с0Ат+Ьіт +СіАт +Ь2т +с2Ат +b3lgT+C3lgAT+b4- (6) где Ат - относительная ширина фракции Дт = Д1 . Коэффициенты
тСОТК
уравнения представлены в таблице 3.
Таблица 3.
Коэффициенты уравнения (6)
Левая пограничная кривая Правая пограничная кривая
а=317,99019 b2=-71,756129 a=-300,50084 b2=54,414265
b0=-576,55377 c2=-10,578539 b0=533,62979 c2=18,791241
Gq=-1,4503019 b3=-582,41463 c0=3,7332203 bs=-583,94303
b1=298,29162 c3=0,12596067 b1=-243,85471 c3=-0,29811572
c1=6,5964931 b4=41,233186 c1=-13,533076 b4=-44,246061
Среднеквадратические отклонения Ап от экспериментальных данных по уравнению (6) составляют: для левой пограничной кривой 0,009 при максимальном 0,025, для правой пограничной кривой 0,007 при максимальном 0,02.
Список литературы
1. Харченко П. М. Результаты экспериментальных исследований бензиновых нефтяных фракций/ П. М. Харченко, В. П. Тимофеев//Научный журнал КубГАУ. -Краснодар. - 2014. - №98(04).
2. Харченко П. М. Исследование плотности и давления насыщенных пород нефтяных фракций / П. М. Харченко, В. П. Тимофеев// Труды КубГАУ. - Краснодар. -
2012. - Т1. - №39. - С. 140 - 142.
3. Харченко П. М. Экспериментальная установка и методики исследования плотности и ДНП агропромышленных сточных вод/ П. М. Харченко, В. В. Христиченко, Н. А. Блощинский// Труды КубГАУ. - Краснодар. - 2012. - Т1. - №37. -
С. 238 - 242.
4. Потапенко И. А. Способ термической обработки деталей машин/ И. А. Потапенко, П. М. Харченко; патент на изобретение RUS 2297459, 12.10.2005.
5. Андрейчук В. К. Термоадаптивный блок озонатора/ В. К. Андрейчук, П. М. Харченко; патент на изобретение RUS 2181103, 19.10.1999.
http://ej.kubagro.ru/2014/05/pdf/62.pdf
Научный журнал КубГАУ, №99(05), 2014 года
13
6. Оськин С. В. Ветроэнергетическая установка/ С. В. Оськин, Д. П. Харченко, П.
М. Харченко; патент на изобретение RUS 2299356, 22.02.2006.
7. Харченко П. М. Вентиляция производственных и коммунально-бытовых зданий/ П. М. Харченко, В. В. Христиченко, А. А. Тимофеюк// Труды КубГАУ. - Краснодар. -
2012. - Т1. - №37. - С. 271 - 275.
8. Харченко П. М. Расчёт вентиляции и отопления производственного здания/ П.
М. Харченко, В. П. Тимофеев// Труды КубГАУ. - Краснодар. - 2013. - Т1. - №42. - С. 152 - 155.
References
1. Harchenko P. M. Rezultati eksperimentalnih issledovaniy benzinovih i neftyanih frakciy/ P. M. Harchenko, V. P. Timofeev//Nauchniy zhurnal KubGAU. - Krasnodar. - 2014.
- №98(04).
2. Harchenko P. M. Issledovanie plotnosti i davleniya nasischennih porod neftyanih frakciy/ P. M. Harchenko, V. P. Timofeev// Trudi KubGAU. - Krasnodar. - 2012. - T1. -№39. - S. 140 - 142.
3. Harchenko P. M. Eksperimentalnaya ustanovka i metodiki issledovaniya plotnosti i DNP agropromishlennih stochnih vod/ P. M. Harchenko, N. A. Bloschinskiy// Trudi KubGAU. - Krasnodar. - 2012. - T1. - №37. - S. 238 - 242.
4. Potapenko I. A. Sposob termicheskoy obrabotki detaley mashin/ I. A. Potapenko, P. M. Harchenko; patent na izobretenie RUS 2297459, 12.10.2005.
5. Andreychuk V. K. Termoadaptivniy blok ozonatora/ V. K. Andreychuk, P. M. Harchenko; patent na izobretenie RUS 2181103, 19.10.1999.
6. Oskin S. V. Vetroenergeticheskaya ustanovka/ S. V. Oskin, D. P. Harchenko, P. M. Harchenko; patent na izobretenie RUS 2299356, 22.02.2006.
7. Harchenko P. M. Ventilyaciya proizvodstvennih i kommunalno-bitovih zdaniy/ P. M. Harchenko, V. V. Hristichenko, A. A. Timofeyuk// Trudi KubGAU. - Krasnodar. - 2012.
- T1. - №37. - S. 271 - 275.
8. Harchenko P. M. Raschet ventilyacii i otopleniya proizvodstvennogo zdaniya/ P. M. Harchenko, V. P. Timofeev// Trudi KubGAU. - Krasnodar. - 2013. - T1. - №42. - S. 152 -155.
http://ej.kubagro.ru/2014/05/pdf/62.pdf
