CC BY
119
ЭНЕРГЕТИКА, ЭЛЕКТРИФИКАЦИЯ И ЭНЕРГЕТИЧЕСКОЕ МАШИНОСТРОЕНИЕ • ЭЛЕКТРОТЕХНИЧЕСКИЕ КОМПЛЕКСЫ...
УДК 537.868:544.2
Г. Т. ЗАКИРЬЯНОВА, Л. А. КОВАЛЕВА, Н. М. НАСЫРОВ
ДВУМЕРНОЕ МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВОЗДЕЙСТВИЯ ВЫСОКОЧАСТОТНОГО ЭЛЕКТРИЧЕСКОГО ПОЛЯ
НА ЭМУЛЬСИЮ
Теоретически исследуется воздействие высокочастотного электрического поля на эмульсию и на разделение эмульсии на нефть и воду. При моделировании учитываются пондеромоторные силы, действующие на среду со стороны электромагнитного поля, перекрестные эффекты тепломассопереноса, возникающие при неизотермическом течении жидкости, и сила тяжести в поле Земли. Проведено сопоставление результатов расчетов с известными экспериментальными данными. Электрическое поле; пондеромоторные силы; перекрестные эффекты тепломассопереноса; расслоение эмульсии
ВВЕДЕНИЕ
Проблема разрушения водонефтяных эмульсий обусловлена необходимостью обезвоживания и обессоливания нефтей, ликвидации нефтешламовых амбаров и утилизации их продукции. Первая задача связана с прогрессирующим ростом числа обводненных месторождений: на многих крупнейших месторождениях России доля воды в продукции составляет 90%, а в целом по России превышает 83,5% [1]. Вторая задача - экологическая - вызвана огромным количеством загрязненных водно-почвенных сред на территории объектов нефтегазодобывающего комплекса и их транспортировки. Для разрушения высокоустойчивых эмульсий широко используется способ разделения эмульсии на нефть и воду внешним электрическим полем низкой (промышленной) частоты. Однако при повышении напряженности электрического поля больше некоторого критического значения возникает побочный процесс, ведущий к электрическому диспергированию капель на более мелкие капли. Явление электрического диспергирования весьма нежелательно, так как при этом затрудняется слияние капель, и эмульсия становится еще более устойчивой. При повышенном содержании воды в эмульсии может произойти электрический пробой. Поэтому в промышленные электродегидраторы для обезвоживания принимается эмульсия с содержанием воды не более 30%.
Контактная информация: (347) 252-99-77 Работа выполнена при поддержке РФФИ, грант
№ 08-01-97032
Для высокочастотного (ВЧ) электрического поля критической напряженности практически не существует [2]. В ВЧ поле может обрабатываться эмульсия с содержанием воды до 8090%. Кроме того, в высоком диапазоне частот эмульсия поглощает электромагнитную (ЭМ) энергию, в среде появляются распределенные источники тепла. Под действием неоднородного электрического поля происходит движение капель в сторону увеличения напряженности электрического поля (пондеромоторные силы) [3], а под действием неоднородного температурного поля происходит движение капель в сторону меньшей температуры (перекрестные явления - термо- и электротермодиффузия) [4]. Все это приводит к столкновению, слиянию и укрупнению капель, в результате чего капли воды, как более тяжелые, чем нефть, под действием гравитации опускаются вниз, происходит расслоение эмульсии на нефть и воду.
ВЧ способ разрушения стойких нефтяных эмульсий успешно опробован в лабораторных и промысловых условиях [5, 6]. Для этого были созданы два макетных образца электродегидра-тора трубного типа. Для питания электродегид-ратора использовался ламповый генератор ЛД-2-60М с водяным охлаждением и рабочей частотой 13,56 МГц. В результате установлено: • применение ВЧ электрического поля обеспечивает эффективное разрушение стойких эмульсий с объемным содержанием воды 2050% и массовым содержанием механических примесей 0,3-10%. Остаточное объемное содержание воды при этом на выходе составило 8-10%, а из отстойника, установленного после электродегидратора, - 1,5-4%;
• использование электродегидратора позволило снизить расход деэмульгатора и температуру термохимического обезвоживания нефти.
Разработанная установка для разрушения водонефтяной эмульсии воздействием ВЧ электрического поля защищена авторским свидетельством [7].
Для понимания механизма разложения водонефтяных эмульсий, обоснования режимов его проведения, прогнозирования результатов в зависимости от свойств среды и параметров внешнего воздействия в предлагаемой работе сформулирована математическая модель в диффузионном приближении к рассматриваемому процессу. Рассматривается воздействие ВЧ электрического поля на водонефтяную эмульсию в вертикальном цилиндрическом конденсаторе. Численно решается двумерная задача в цилиндрической системе координат (г, ф, z), ось которой направлена вдоль оси конденсатора снизу вверх. Начало системы координат расположено в центре основания конденсатора. Принимается аксиальная симметрия задачи, т. е. рассматриваемые процессы от координаты ф не зависят. В расчетах принимается, что вода является диэлектриком, так как на высоких частотах даже при наличии солей потери за счет проводимости среды значительно меньше диэлектрических потерь [2].
1. ПОСТАНОВКА ЗАДАЧИ И ОСНОВНЫЕ УРАВНЕНИЯ ПРИ ВОЗДЕЙСТВИИ ВЧ ПОЛЯ
На капли воды, движущиеся в покоящейся жидкости со скоростью и, действует сила сопротивления Стокса [3]:
Р = 6лтл уг0 и,
где г0 - радиус капли; п - вязкость дисперсионной среды (нефти).
Система уравнений, описывающая движение глобул воды в эмульсии под действием сил электромагнитного происхождения и тяжести, а также пространственно-концентрационное
распределение диспергированной фазы, имеет вид [3]:
до.
do.
do 3koor с
-+-
Pv 1 +Pv°r dr +Pv°z dz 4лг03 ' 2 1
doz do do 3koo
Pv^£ +Pv°r_r +Pv°z-f =—f -Pvg, dt dr dz 4лг3
e =i(e 0°:)
°с дс ,
(1)
(2)
дс D д і де ^ д2с
— =---------------1 r— I + D—-
dt r dr У дг J Iz
дс дс 1 д і j
-Ur IT-о z aZ+711r (a
T + aE
)d IL
’ dr
+ (3)
+5 К'
aT + aE
D £
q =
we0e'.tgSs
2
, E = U ln ■r2, rr
Pv = pvс, pv = const,
P/ =P/(1 - с), p/ = const,
P:=P /(1+Pv), * * (l ) *
= Cf (1 - с)+ с • сv, = 1 f (1 - с)+ сі
(4)
Л/ =Ло ехр(-уАТ).
Здесь V, / 5 - индексы, относящиеся к воде, нефти и их смеси; АТ = Т - Т0; Т - температура эмульсии; Т0 - первоначальная температура эмульсии; q - плотность распределенных источников тепла; аь 0 - параметры термодиффузии и электротермодиффузии; р/ и р° - истинные плотности несущей и диспергированной фаз; Е - напряженность электрического поля; е0 - электрическая постоянная; г'к, tg 55 -относительная диэлектрическая проницаемость и тангенс угла диэлектрических потерь смеси нефти и воды; п0 - первоначальная вязкость нефти; у - коэффициент, учитывающий зависимость вязкости нефти от температуры; р5, с5, - плотность, удельная теплоемкость и теплопроводность смеси нефти и воды; и - разность потенциалов между обкладками конденсатора; г1, г2 - радиусы внутреннего и внешнего обкладок конденсатора; с - объемная концентрация глобул воды; Б - коэффициент диффузионного перемешивания.
Расчет диэлектрической проницаемости смеси нефти и воды в зависимости от их объемного соотношения производится по формуле Лихтенеккера:
г _ г'с, /(1-с)
^ _ ч ь/
Краевые условия. Для расчетных исследований принято, что обкладки конденсатора непроницаемы для капелек воды и гидрофобны, поэтому концентрация частиц на них принята равной нулю. Конденсатор со всех сторон теплоизолирован:
(Г^,0)_ 0, иг (Г^,Ґ)_ 0, (Г,И,Ґ)_ 0, (5)
с(г, z,0) = с0, с(г1, z, t )= 0.
с(г2, z, t)= 0, = 0,
dz
2
Т (г,0)_ Т0, _ 0, 0,
ЭТ(г,0,ґ)_ 0 дГдТ(г,к,ґ) 0 дг
(7)
Эг Эг
где с0 - исходная объемная концентрация капель воды в эмульсии, к - высота цилиндрического конденсатора.
Рассматривался случай открытого сверху конденсатора. В соответствии с этим в точке г = к принималось граничное условие:
Эс(г, к, ґ)
Эг
_ 0.
(8)
2. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ
Данные по воздействию ВЧ электрического поля, электрического поля промышленной частоты и их совместного действия на эмульсию приведены в [8]. В этой работе модель вертикального цилиндрического электродегидратора представляла собой коаксиальную систему двух латунных труб диаметрами 40 и 12 мм. В качестве источника высокочастотного напряжения использовали генератор ВЧД 2,5/13-СД-Л01 с выходной мощностью около 1 кВт и частотой генерации 13,56 МГц. Искусственную водонефтяную эмульсию готовили циркуляцией в шестеренчатом насосе (2000 об/мин) из нефти Ромашкинского месторождения. Время обработки эмульсии полем составляло 30 с. Для сравнения эффективности действия электрического поля промышленной частоты и высокочастотного поля на степень разрушения эмульсий напряжение питания в обоих случаях было одинаковым - около 1,3 кВ. Обрабатывали эмульсии полем при комнатной температуре Т = 20 0С. Затем в стеклянных отстойниках исследовали кинетику расслоения эмульсий с 5,
10, 15, 20%-м содержанием воды при Т = 50 °С (до такой же температуры нагревалась эмульсия при высокочастотном воздействии).
3. РЕШЕНИЕ ЗАДАЧИ И АНАЛИЗ РЕЗУЛЬТАТОВ
Система уравнений (1)-(4) с краевыми условиями (5)-(8) решалась методом конечных разностей по неявной схеме.
Принятые для расчетов параметры среды: Г1 = 0,006 м; г2 = 0,02 м; к = 0,02 м; =
= 848 кг/м3; р0=1000 кг/м3; Cf = 1885 Дж/кг • К; су = 4190 Дж/кг • К; V = 0,125 Вт/м • К; \ = = 0,582 Вт/м • К; Т0 = 20 0С; п0 = 0,01561 Па • с; у = 0,02902 К-1; е'у =81; гв 5 = 0,0025; е/ =2,5;
tg 5^- = 0,1; с0 = 0,2. В расчетах использованы
значения термодиффузионных коэффициентов, полученные в работе [9] путем сравнения экспериментального и математического моделирования: ат = 0,0025 К-1; ает =0,115 К-1. Коэффициент диффузионного перемешивания определен путем сравнения расчетных и опубликованных в [8] экспериментальных данных: В = = 4* 10-5 м2/с. Ориентирами при сравнении служили количественные характеристики кинетики расслоения в результате ВЧ воздействия.
На рис. 1-4 приведены графики распределения концентрация капель воды и температуры между внешним радиусом внутренней обкладки и внутренним радиусом внешней обкладки конденсатора в различные моменты времени. Расчеты показали, что кривые распределения концентрации капель воды в различные моменты времени мало отличаются друг от друга. Таким образом, к моменту времени 10 с процесс перераспределения капель воды в пространстве практически завершается (рис. 1), тогда как температура эмульсии в результате воздействия продолжает расти (рис. 2-4).
Рис. 1. Распределение концентрации глобул воды в эмульсии при / = 10 с
Рис. 2. Распределение температуры в конденсаторе при ґ=10 с
Рис. 3. Распределение температуры в конденсаторе при ґ = 20 с
Рис. 4. Распределение температуры в конденсаторе при / = 30 с
Такое поведение системы объясняется тем, что в начальный период действия поля порожденные им пондеромоторные силы сдвигают капли воды ближе к внутренней обкладке конденсатора, где больше напряженность электрического поля. В точках максимума концентрация воды может достигать значительной величины, но все еще это концентрированная водонефтяная эмульсия. Известно также [10], что к концентрированным эмульсиям типа «вода в масле» относятся высокодисперсные системы со сравнительно большим содержанием дисперсной фазы (до 80 объемных % и даже выше). Максимальные значения температуры в этой области также связаны с величиной напряженности поля и, соответственно, наибольшими значениями плотности тепловых источников. Однако к моменту времени ^ = 10 с после начала ЭМ воздействия кривые распределения концентрации выходят на стационарный режим, т. е. наступает некоторый баланс всех действующих на систему сил, несмотря на то, что нагрев эмульсии продолжается.
4. РАСЧЕТ ДИНАМИКИ ОТСТОЯ ЭМУЛЬСИИ И СРАВНЕНИЕ С ЭКСПЕРИМЕНТОМ
Для проведения оценочных расчетов эффективности расслоения эмульсии необходимо учитывать, что реальные эмульсии являются полидисперсными системами, в которых большие капли опускаются быстрее и попутно захватывают маленькие, благодаря наличию ди-поль-дипольного взаимодействия [11]. При наложении на эмульсию внешнего электрического поля вводится константа коагуляции, характеризующая процесс укрупнения частиц:
К(Г^)_(иу -и„)Е
&1г 0г5
181V,
1-& Ррё
м
с
где Уг - объем большой капли; иу ,ьк - скорость оседания большой и маленькой капель
соответственно, определяемые по формуле Стокса
и
(г, г, ґ):
2/ V/3
6щн (г, г, ґ) &,г ’
(9)
где у, I - объем капель; индекс I относится к маленькой капле; g - ускорение свободного падения; р^ - разность плотностей воды и нефти; g1 = ^ ; В, — - диаметры большой и ма-
ленькой капель; X - коэффициент, компенсирующий завышенное значение силы диполь-дипольного взаимодействия (X < 1). Параметр X и соотношение — / В подбирались из условия максимальной близости экспериментальных и расчетных кривых кинетики отстоя эмульсии. Наилучшее совпадение результатов соответствует следующим значениям: — / В = 0,26; X = = 0,9.
Вследствие слияния капель воды, число их в эмульсии уменьшается. Время, в течение которого число капель уменьшается в 2 раза, определяется выражением:
^(г’2’* ^= к (у>)щ 0г) ’
где Щ(г, z, 0 - число капель в единице объема. Число капель в единице объема конденсатора в предположении, что эмульсия состоит преимущественно из больших капель, так как при рассматриваемых полях диспергирование незначительно, можно определить по формуле
с(г, z, /)
N (г,г,ґ)_ ■
V,
(10)
где с(г, z, 0 - объемная концентрация капель воды в эмульсии. Далее, определив число периодов п(г, z, 0 = 1ов2 Щ(г, z, 0, за время которых число капель уменьшается вдвое, можно вычислить время, за которое число капель в эмульсии уменьшится до нуля: 4её(г, z, 0 = = 4(г, z, ^ • п(г, z, ^.
Для расчета динамики расслоения эмульсии вычисляются интегральные значения количества капель воды в эмульсии в начальный Щ0 и текущий Щк моменты времени. Результаты подставляются в выражение
( ш \
V(ґ) _
1 - N.
N
• 100%.
0 У
Временная зависимость ^ґ) определяет динамику расслоения эмульсии.
Расчеты показали, в случае открытого сверху конденсатора расслоение эмульсии проис-
ходит приблизительно одновременно по высоте. Начинается расслоение на стороне внутренней обкладки конденсатора ближе к центру, затем распространяется к обкладкам. На рис. 5 приведена расчетная кривая динамики разделения эмульсии на нефть и воду, на рис. 6 приведены для сравнения экспериментальные кривые динамики разделения эмульсии под действием на нее различных полей, в том числе ВЧ электрического поля [8]. Сравнение показывает хорошее качественное совпадение результатов расчетов с экспериментальными данными (см. кривая 3 на рис. 6).
V, %
о
10
20
30
Рис. 5. Расчетная динамика расслоения эмульсии -зависимость объемной доли эмульсии У (%), расслоившейся на нефть и воду к моменту времени /
Рис. 6. Экспериментальная динамика расслоения эмульсии: 1 - не обработанная полем; 2, 3, 4 - обработанная соответственно электрическим полем промышленной частоты, высокочастотным полем и совместно этими полями
ЗАКЛЮЧЕНИЕ
Сформулированная математическая модель процесса воздействия ВЧ электрического поля на эмульсию и динамики расслоения ее на нефть и воду позволяет дать количественную оценку распределения концентрации капель воды в эмульсии, температуры и скорости разделения эмульсии на нефть и воду, что необходимо при прогнозировании использования данного метода для обезвоживания эмульсии. Рассмотрен случай открытого сверху вертикального цилиндрического конденсатора.
Распределения концентрации капель воды и температуры приблизительно одинаковы по высоте конденсатора. Расслоение эмульсии начинается с верхней части конденсатора. Сравнение теоретических расчетов с известными экспериментальными исследованиями показывает хорошее качественное совпадение расчетных и экспериментальных данных.
СПИСОК ЛИТЕРАТУРЫ
1. Хавкин А. Я. Нанотехнологии в добыче нефти и газа. М.: Нефть и газ. ПЦ «НТИС», 2008. 171 с.
2. Саяхов Ф. Л. Исследование термо- и гидродинамических процессов в многофазных средах в высокочастотном электромагнитном поле применительно к нефтедобыче: дисс. д-ра физ.-мат. наук. М., 1984. 449 с.
3. Диффузионная и многоскоростная модели двухфазных сред в электрическом поле / В. В. Гого-сов [и др.] // ПММ. 1980. Вып. 2. С. 290-300.
4. Перекрестные явления переноса в дисперсных системах, взаимодействующих с высокочастотным электромагнитным полем / Р. И. Нигматуллин [и др.] // Докл. РАН. 2001. Т. 377, № 3. С. 340-343.
5. Разрушение стойких нефтяных эмульсий на промыслах электромагнитными полями резонансных частот / В. С. Хакимов [и др.] // Научнопроизводственные достижения нефтяной промышленности в новых условиях хозяйствования: науч.-техн. инф. сб. М., 1989. Вып. 6. С. 12-16.
6. Interaction of Highly Stable Water-Oil Emulsion With Radio Frequency Electromagnetic Field / F. L. Sayakhov [et al.] // Heavy Crude and Tar Sands -Hydrocarbons For the 21-st Century. 5-th UNITAR International Conference on Heavy Crude and T ar Sands. 1991. V. 3. Recovery processes. P. 507-513.
7. А.с. СССР № 700163, кл. В 01 D 17/06, С 10 G 33/02. Установка для разрушения водонефтяной эмульсии / Ф. Л. Саяхов [и др.]. Опубл. 05.12.1979, Бюл. № 44.
8. Диэлектрические свойства и агрегативная устойчивость водонефтяных эмульсий / Ф. Л. Саяхов [и др.] // Нефтяное хозяйство. 1979. № 1. C. 36-39.
9. Влияние высокочастотного электромагнитного поля на перекрестные эффекты переноса в многокомпонентных системах / Ф. Л. Саяхов [и др.] //
Магнитная гидродинамика. 1998. Т. 34, N° 2. С. 148-157.
10. Эмульсии / под ред. А. А. Абрамзона. М.: Наука, 1972. 321 с.
11. Панченков Г. М., Цабек Л. К. Поведение эмульсии во внешнем электрическом поле. М.: Химия, 1969. 190 с.
ОБ АВТОРАХ
Закирьянова Галия Тимерга-зиевна, преп. Уфимск. ин-та (филиала) Российск. гос. торг.-эконом. ун-та. Дипл. физик (БашГУ, 1995). Иссл. в обл. электрогидродинамики многофазных систем.
Ковалева Лиала Ароновна,
проф., зав. каф. прикл. физики БашГУ. Дипл. физик (БашГУ, 1974). Д-р техн. наук по механике жидкости, газа и плазмы (М., 1998). Иссл. в обл. физики и гидродинамики насыщен. пористых сред, термодинамики необратимых процессов.
Насыров Нур Мутагарович,
доц. той же каф. Дипл. физик (БашГУ, 1974). Канд. физ.-мат. наук по теплофизике и молеку-лярн. физике (Уфа, 1992). Иссл. в обл. физ.-техн. проблем нефти и газа, электродинамики.
