НОВЫЙ ПОДХОД К РАСЧЕТУ ВЯЗКОСТИ ЖИДКИХ СМЕСЕЙ УГЛЕВОДОРОДОВ НА ОСНОВЕ МОДИФИЦИРОВАННОГО УРАВНЕНИЯ АРРЕНИУСА Текст научной статьи по специальности «Химические технологии»

ПРОБЛЕМЫ

НЕФТЕГАЗОВОЙ И УГОЛЬНОЙ ПРОМЫШЛЕННОСТИ

PROBLEMS OF OIL, GAS AND COAL INDUSTRY

Статья поступила в редакцию 31.07.12. Ред. рег. № 1380 The article has entered in publishing office 31.07.12. Ed. reg. No. 1380

УДК 544.034

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

А.Г. Окунев, Е.В. Пархомчук, А.И. Лысиков, В. С. Деревщиков

Институт катализа им. Г.К. Борескова Сибирского отделения РАН 630090 Новосибирск, ул. Лаврентьева, д. 5 Тел./факс: (383)326-94-32; e-mail: okunev@catalysis.ru

Заключение совета рецензентов: 15.08.12 Заключение совета экспертов: 20.08.12 Принято к публикации: 25.08.12

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

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

A NEW APPROACH TO THE CALCULATION OF VISCOSITY OF LIQUID HYDROCARBON MIXTURES BASED ON MODIFIED ARRHENIUS EQUATION

A.G. Okunev, E.V. Parkhomchuk, A.I. Lysikov, V.S. Derevshchikov

Boreskov Institute of Catalysis SB RAS 5 Lavrentieva str., Novosibirsk, 630090, Russia Tel./Fax: 383)326-94-32; e-mail: okunev@catalysis.ru

Referred: 15.08.12 Expertise: 20.08.12 Accepted: 25.08.12

Nowadays transportation and processing of heavy crudes are becoming of increasing importance. This raises the problem of developing appropriate ways to control and predict their viscosity. This paper presents a new approach to the calculation of viscosities of mixtures of hydrocarbons, based on Arrhenius equation. The equation is modified using the principle of linearity of free energies. The viscosities of binary mixtures of hydrocarbons are fitted using the data on the viscosity of the initial reagents and an approximation parameter that has the physical meaning of the correlation between the activation energy of motion of particles in the mixture and the energy of solvation. New approach is demonstrated to be promising to the description of mixtures of heavy oil and light solvents.

Keywords: viscosity, mixing rule, hydrocarbons.

Введение

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

Уравнение Аррениуса было предложено одним из первых для оценки вязкости смесей по известным

вязкостям исходных компонентов. В основе уравнения Аррениуса лежат представления об активированном процессе перемещения частиц в смеси: V = к0е~ЕЯТ , где V - кинематическая вязкость; к0 -предэкспоненциальный фактор; Е - энергия активации; Т - температура; Я - газовая постоянная. Уравнение Аррениуса получается в предположении, что энергия активации выхода из клетки частиц в смеси одинакова и линейно зависит от состава смеси:

EM X1E1 + x2 E2 ,

(1)

International Scientific Journal for Alternative Energy and Ecology № 09 (113) 2012

© Scientific Technical Centre «TATA», 2012

v _ k EM¡RT _ k + X2E2 )RT _

VM 0c ""0C

_{k0e-EiR )X {k0e-E2'RT f _vf v*

(2)

где х, - доли компонентов в смеси. Часто уравнение Аррениуса записывают в виде

In (VM )_ xiln (v! )+ x2¡n (v2 )

(3)

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

Поскольку уравнение (1) плохо описывает вязкость смесей нефтепродуктов, были предприняты многократные попытки его модифицировать. В недавнем обзоре [1] описаны 27 различных формул для предсказания вязкости смесей нефтепродуктов. Наибольшие трудности возникают при описании смесей, в которых исходные компоненты имеют многократно различающуюся вязкость.

В методе Грюнберга и Ниссана [2], позже развитом Ратклифом и др. [3], предлагается модифицировать уравнение Аррениуса, вводя в него дополнительный параметр, учитывающий неидеальность раствора, аналогично тому, как это делается в теории регулярных растворов:

ln (VM ) _ Xi ln (vi ) + x2 ln (V2 ) + X2 •

(4)

Теоретический анализ

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

_ x1v1 + x2v2 •

(6)

Используем уравнение Аррениуса для температурной зависимости вязкости компонентов:

_ k e-EM,l/RT

vi _ k0e

_ k e~ EM -2¡RT

V2 _ k0e

(7)

где Ем 1 - энергия активации при перемещении в

смеси 1-го компонента.

Так же, как и в формуле (1), предположим, что энергия активации выхода из клетки ЕМ1 для компонентов смеси линейно зависит от состава смеси. Фактически это утверждение эквивалентно наличию корреляционной зависимости между энергией активации выхода из клетки и энергией сольватации рассматриваемой молекулы. Подобного типа соотношения хорошо известны в физико-химической литературе как принцип линейности свободных энергий [9]. Однако энергии активации для разных компонентов теперь не обязательно одинаковы:

ЕМ ,1 _ Х1Е1 + Х2 Е2 — аХ2 (Е2 — Е1 ) ,

Подстановка этих выражений в (7) и (6) приводит к зависимости

VM _ xik0e

x1E1 + x2E2 -ax2 (E2 -E1) RT

+ x2 k0 e

x1E1 + x2E2 -axj (E2 -E1) RT

Параметр e находится из условия наилучшего приближения формулы (4) к экспериментальным данным.

Большое распространение также получили эмпирические формулы, использующие дополнительный параметр аппроксимации [4]. Например, очень хорошие результаты получаются с помощью формулы Вальтера [5, 6]

log log (vM + C) = x loglog (vj + C) + x2 loglog (v2 + C).

(5)

Значение С = 0,7 или 0,8 часто предполагают заранее известным на основании обработки большого массива экспериментальных данных [7]. В этом случае формула Вальтера позволяет предсказывать вязкость смесей только на основе вязкостей индивидуальных компонентов [8].

В данной работе предложена новая формула для аппроксимации вязкостей смесей нефтепродуктов на основе анализа взаимного влияния компонентов смеси.

_ k0e

_ vM

2 ((2 - Ei)

X (E2 - Ei) Л

(X (2/ vi )-aX2 + X2 (2/ vi )-axi),

(8)

где vм Агг - вязкость смеси, рассчитанная по формуле (3). Очевидно, что при а = 0 уравнение (8) совпадает с уравнением (3).

Методика расчетов

Аппроксимацию выражениями (4), (5), (8) проводили с использованием экспериментальных данных, опубликованных в работах [1, 10, 11]. Вязкости исходных индивидуальных компонентов принимали равными их экспериментальным значениям. Оптимизацию параметров выполняли с помощью метода взвешенных наименьших квадратов, минимизируя квадратичную форму

5 ^^Ь^Н ,

Международный научный журнал «Альтернативная энергетика и экология» № 09 (113) 2012 © Научно-технический центр «TATA», 2012

Проблемы нефтегазовой и угольной промышленности

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

Результаты и их обсуждение

Представляет интерес сравнение точности аппроксимации новой формулой (8) и ранее известными однопараметрическими выражениями (4), (5). Для проверки были использованы экспериментальные данные по вязкости смесей тяжелых нефтей с различными растворителями, приведенные в табл. 4 работы [1], табл. 5, 6, 8, 9 работы [10]. В последней статье приведены также оцифрованные данные из работы [11]. С помощью итерационной процедуры были найдены значения параметров, наилучшим образом описывающие все результаты в пределах серий данных, относящихся к однотипным смесям.

Во всех случаях наименьшее взвешенное среднеквадратичное отклонение достигалось при использовании предложенной формулы (8) и формулы Валь-

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

Качество аппроксимаций может быть также оценено из графического сопоставления некоторых экспериментальных данных из обработанного массива с расчетными кривыми (рис. 1-3). Коэффициент детерминированности Я2 для новой формулы (8) практически во всех случаях выше 0,97. Оптимальные значения параметра а, определенные для различных смесей нефтепродуктов, отличаются менее чем в два раза, в то время как значение параметра С в формуле Вальтера и значение параметра е в приближении Грюнберга-Ниссана изменяются в три раза. При этом следует отметить, что приписывание параметру в формуле Вальтера фиксированных значений (см., например, [7]) может являться источником существенных ошибок при оценке вязкости.

Результаты аппроксимации литературных данных по вязкости смесей нефтепродуктов Results of approximation of published data on the viscosity of mixtures of petroleum

Модель Смесь нефть-дизель [1] Смесь нефть-н-декан [10] Смесь битум-растворитель [11]

Оптимум S N Оптимум S N Оптимум S N

Эта работа (8) а = 0,26 0,121 125 а = 0,46 0,165 78 а = 0,45 0,301 27

Грюнберг-Ниссан (4) е = -3,4 0,276 е = -4,3 0,206 е = -9 0,405

Вальтер (5) С = 2,7 0,124 С = 0,84 0,19 С = 2,2 0,224

1000-

^ 100-

10-

Эксп R2, 298 К R2, 327 К

о Эта работа (8) 0.9992 0.998

V Грюнберг (4) 0.999 0.981

О Вальтер 0.998 0.999

O

a

Т = 298 К

327 К

0.0

—I—

0.2

—I—

0.4

—I—

0.6

1.0

х, (об.%)

Рис. 1. Экспериментальные зависимости вязкости смеси тяжелой нефти (плотность по API 15,93) и дизельной фракции от объемной доли вязкого компонента x (по данным работы [1]) и их аппроксимации с помощью уравнений (4), (5), (8)

Fig. 1. Experimental dependence of the viscosity of heavy oil mixture (API density of 15.93) with diesel fraction on the volume fraction of the viscous

component x (according to [1]) and their approximations using equations (4), (5), (8)

International Scientific Journal for Alternative Energy and Ecology № 09 (113) 2012

© Scientific Technical Centre «TATA», 2012

100-

10-

1 -

Эксп R2, 0.407 R2, 0.684 R2, 0.872

О Эта работа 0.71 0.977 0.987

V Грюнберг 0.97 0.339 0.9996

о Вальтер 0.791 0.88 0.704

280

300

320

—I—

340

—I—

360

—I—

380

—I—

400

—I—

420

V -1—

440

х=0.872

х=0.684 х=0.407

460

T (K)

1000000

100000

10000

Эксп R2, 298 K R2, 356 K

О Эта работа 0.987 0.987

V Грюнберг 0.993 0.916

О Вальтер 0.994 0.977

IZ

о

к m

1000

100т

V Т = 298 К

0.0

—I—

0.2

0.4

—I—

0.6

—I—

0.8

Рис. 2. Температурные зависимости вязкости смеси тяжелой нефти и н-декана, измеренные при разной объемной доле вязкого компонента x (по данным работы [10]) и их аппроксимации с помощью уравнений (4), (5), (8) Fig. 2. Experimental dependence of the viscosity of

heavy oil mixture with n-decane, measured at different temperatures and volume fractions of the viscous component x (according to [10]) and their approximations using equations (4), (5), (8)

—i

1.0

х, (об.%)

Рис. 3. Экспериментальные зависимости вязкости синтетических нефтей на основе битумов и растворителя от объемной доли вязкого компонента x (по данным работы [11]) и их аппроксимации с помощью уравнений (4), (5), (8)

Fig. 3. Experimental dependence of the viscosity of synthetic crude oil made of the mixture of bitumen

with light solvent on the volume fraction of the viscous component x (according to [11]) and their approximations using equations (4), (5), (8)

Заключение

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

Работа выполнена при финансовой поддержке Министерства образования и науки РФ ГК № 16.516.11.6078.

Список литературы

1. Centeno G., Sanchez-Reyna G., Ancheyta J., Munoz J.A.D., Cardona N. Testing various mixing rules for calculation of viscosity of petroleum blends // Fuel. 2011. Vol. 90, No. 12. P. 3561-3570.

2. Grunberg L., Nissan A.H. The energies of vaporization, viscosity and cohesion and the structure of liquids // Transactions of the Faraday Society. 1949. Vol. 45. P. 125-137.

3. Ratcliff G.A., Khan M.A. Prediction of the viscosities of liquid mixtures by a group solution model // Can. J. Chem. Eng. 1971. Vol. 49. P. 125-129.

4. Фукс Г.И. Вязкость и пластичность нефтепродуктов. М.; Ижевск: Институт компьютерных исследований, 2003.

5. Walther C. The evaluation of viscosity data // Erdol Teer. 1931. Vol. 7. P. 382.

6. Долматов Л.В., Кутуков Е.Г., Кутуков И.Е. Адекватность математических моделей для расчета вязкости смесей жидких нефтепродуктов // Химия и технология топлив и масел. 2001. Вып. 3. C. 43-45.

7. Аметов И.М. Добыча тяжелых и высоковязких нефтей. М.: Недра, 1985.

8. Эйгенсон А.С., Шейх-Али Д.М. Расчет плотности и вязкости пластовой нефти по данным поверхностной дегазации // Геология нефти и газа. 1989. Вып. 11. С. 45-51.

9. Эмануэль Н.М., Кнорре Д.Г. Курс химической кинетики. М.: Высшая школа, 1984.

10. Barrufet M.A., Setiadarma A. Reliable heavy oil-solvent viscosity mixing rules for viscosities up to 450 K, oil-solvent viscosity ratios up to 4x105, and any solvent proportion // Fluid Phase Equilibria. 2003. Vol. 213. P. 65-79.

11. Shu W.R. A viscosity correlation for mixtures of heavy oil, bitumen, and petroleum fractions // SPEJ (June 1984)277-282.

rxTi

— TATA — CXJ

Международный научный журнал «Альтернативная энергетика и экология» № 09 (113) 2012 © Научно-технический центр «TATA», 2012