CC BY
34
Исследование глубины нагрева высоковязкой диэлектрической жидкости в высокочастотном электромагнитном поле методом теории размерности
М. А. Фатыхов
Федеральное государственное образовательное учреждение высшего профессионального образования «Башкирский государственный педагогический университет им. М. Акмуллы», ул. Октябрьской революции, 3, г. Уфа, 450000, Башкортостан, Россия, e-mail: [email protected]
Получены критерии подобия, описывающие процесс нагрева нефтяного пласта в высокочастотном электромагнитном поле. Показана возможность применения их для определения глубины нагрева нефтяного пласта.
Ключевые слова: высоковязкая диэлектрическая жидкость, высокочастотное электромагнитное поле, метод, теория размерности, нефтяной пласт, температура.
УДК 539.1: 622.276
В настоящее время изучаются и находят применение новые методы повышения нефтеотдачи и интенсификации нефтедобычи. К ним можно отнести и высокочастотное электромагнитное воздействие. Объектом данного метода является нефтеносный пласт. Цель метода - интенсификация извлечения высоковязкой, неньютоновской, битумной, высокопарафинистой нефти из залежей с плохими коллекторскими свойствами, высокообводненных и глубокозалегающих пластов.
При высокочастотном электромагнитном воздействии происходят интенсивный нагрев нефтяного пласта вокруг скважины, уменьшение вязкости нефти. Радиус области разжижения нефти со временем увеличивается. С повышением температуры происходят физико-химические процессы. Эти процессы и явления исследованы теоретически и экспериментально, однако отсутствует их подробный анализ. Одним из методов решения этой проблемы является метод теории подобия и размерности [1], который применительно к рассматриваемым задачам не разработан.
Методы теории подобия и размерности играют особенно большую роль при моделировании физических процессов [2]. Применение этих методов к нашей задаче значительно облегчает обработку результатов, полученных при моделировании процесса нагрева пласта, а также позволяет получить достаточно точные выводы при сравнении промысловых и лабораторных данных.
Высокочастотное электромагнитное воздействие представляет практический интерес для глубокого нагрева нефтяного пласта [3]. Впервые это положение проверялось на Мордово-Кармальском битумном участке Республики
Татарстан, где осуществлялся высокочастотный разогрев призабойной зоны пласта по радиусу от скважины на глубину 5 м в течение 12 суток при мощности 15-20 кВт [4, 5]. На этой глубине температура термопарами не была зарегистрирована. С целью установления достоверности полученных данных в настоящей работе исследуются особенности пространственно-временного распределения температуры в пласте методом теории размерности и подобия [1].
В настоящее время при изучении температурных полей в нефтяных пластах используются два метода, названных «точной схемой» и «схемой сосредоточенной ёмкости» [6]. В «точной схеме» пласт и окружающие его породы считаются термически изотропными, имеющими теплофизические характеристики, совпадающие с характеристиками реального пласта, его кровли и подошвы. «Схема сосредоточенной ёмкости» предполагает, что пласт имеет бесконечно большую теплопроводность в вертикальном направлении, а теплопроводность пласта в направлении его простирания является конечной, совпадающей с теплопроводностью реального пласта. Породы считаются термически изотропными с реальным значением коэффициента теплопроводности. В «схеме сосредоточенной ёмкости» поле температуры в интервале пласта не зависит от вертикальной координаты.
В предпринятом исследовании воспользуемся «точной схемой». Следуя работе [7], рассмотрим температурную задачу в цилиндрической системе координат, где среда представлена тремя областями с плоскими границами раздела, перпендикулярными оси г (рис. 1).
В описании температурной задачи примем следующие допущения: пласт считается однородным и изотропным по гидродинамическим и
© Фатыхов М.А., Электронная обработка материалов, 2016, 52(1), 78-83.
теплофизическим свойствам; породы, окружающие пласт, предполагаются непроницаемыми и изотропными по теплофизическим свойствам. Пренебрегаем зависимостью коэффициента теплопроводности от радиальной координаты. Источник, связанный со скважиной, излучает высокочастотные (ВЧ) электромагнитные (ЭМ) волны в радиальном направлении. Вследствие объемного поглощения электромагнитной энергии вокруг скважины происходит разогрев пласта. Фильтрационным движением нефти и связанным с ним конвективным переносом тепла, а также теплообменом на поверхности скважины пренебрега-ется.
Т1 = Т2 = Т = Т0 при г = 0,
Подошва
Рис. 1. Геометрия модели. В пласте расположен излучатель электромагнитных волн, высота которого равна толщине пласта.
В рамках указанной модели процесс нагрева пласта описывается уравнением теплопроводности:
дТ I д2Т 1 дТ д2Т
-= а I —-+--+—-
дг \ дг г дг дг
аЫ
пгкср
0 е-2а(г-г0)
дТ1
"дГ
= а,
I д2Т1 1 дТ д2Т Л —- +---L +--1
ч дг2 г дг дг2
дТ2
"дГ
(
■ = ап
\
д 2Т2 1 дТ2 д 2Т2 —- +---2 +--2
у дг2 г дг дг2 у
(1)
(2)
(3)
где а, а1, а2 - коэффициенты температуропроводности пласта 1- и 2-го слоев соответственно,
м2/с;
- удельная теплоемкость, Дж/кг-К;
р - плотность, кг/м ; к - толщина пласта, м; г и г - цилиндрические координаты, м; Т, Т1, Т2 -температуры пласта 1- и 2-го слоев соответственно, К; г - время воздействия высокочастотного электромагнитного поля на диэлектрическую среду, с; Ы0 - мощность ВЧ генератора, Вт; г0 - радиус скважины, м; а - коэффициент поглощения среды, м-1.
На границах пласта с покрывающей и подстилающей породами заданы условия равенства температур и тепловых потоков, то есть краевые условия:
X
X.
дТ1
1
дТ2
к
г=+— 2
= X дТ
дг
к
г=— 2
= X
к
г=— 2
дТ
"дг
Т1 = Т2 = Т = Т0 при г ^ да,
дТ! дТ2 дТ п —^ = —1 = — = 0 при г = г0
дг дг дг
(4)
(5)
(6)
(7)
(8)
где Т0 - начальная температура до включения электромагнитного поля; X, Хь Х2 - коэффициенты теплопроводности пласта 1- и 2-го слоев, Вт/Км.
Из системы уравнений теплопроводности видно, что основными величинами, от которых зависит температура в пласте, являются коэффициент поглощения а, мощность электромагнитного поля Ы0, коэффициент теплопроводности X = сра окружающих пласт пород с удельной теплоемкостью с и плотностью р; коэффициент теплопроводности ХП, сП и рП - соответственно удельная теплоемкость и плотность, толщина пласта к, радиус пласта г0, время воздействия г.
При этих условиях, согласно системам (1)-(8), распространение тепла в пласте определяется следующими величинами:
ДТ, N 0, Хп, сп, рп, X, с, р, к, г, г, а, г,
(9)
где ДТ = Т - Т0 - изменение температуры пласта; Т0, Т - начальная температура пласта и температура в точке с координатами г, г, г.
Запишем размерности фигурирующих в задаче величин:
[Ы0] = 1}МТРЪ; [ДТ] = 0; [X] = МЬТ-'0-1;
[с] = Ь%-20-1; [р] = М173; (10)
[Xп ]= МЬТР 3©-1; [сп ] = Ь2Тр-20-'; [рп ] = МЬ-3;
[к]= Ь; [г] = Ь; [I] = Тр; [г] = Ь; [а] = I-'.
Размерности 13 величин (9) выражаются через пять основных единиц измерения: температуру - 6, длину - Ь, массу - М, количество теплоты - Qp и время - ТР. Число величин с независимыми размерностями также равно пяти:
[ДТ] = 0, [к] = Ь, [Г] = Тр, [р] = МЬЬ3.
(11)
Поэтому из 13 размерных величин (9) можно образовать 13-5=8 независимых безразмерных критериев подобия. Для их нахождения воспользуемся уравнением теплопроводности и граничными условиями (4)-(8). Используя параметры, которые входят в это уравнение, и величины (9),
с
можно составить следующие безразмерные комбинации, или критерии:
т* _ АЛпи ро _ у
N
сп Р п и
ы*_ ^ я*_ Г т01' И
I*_^, л_А, -, -Р-.
И ^П сп Рп
(12)
Исходя из некоторых физических соображений, число полученных безразмерных параметров можно несколько уменьшить. Величины сп и рп пласта входят в параметр Фурье Ро только в виде произведения спрп. То же самое можно сказать относительно с и р, так как в параметр Ро вместо спрп можно подставить ср. поскольку указанные величины в остальные параметры в другой комбинации не входят, то их можно объединить, образовав один безразмерный параметр:
С _
сР
спРп
(13)
С учетом (9), (12) и (13) функциональную зависимость для безразмерной температуры пласта можем записать в виде:
Т* _/(Ро,N,Л,С,Я*,I*).
(14)
Соответственно изменение температуры пласта АТ определяется решением:
АТ _ / (о, N *, Л, С, Я*, I *). (15) И у '
Из полученного выражения (15) следует, что изменение температуры АТ пласта прямо пропорционально мощности N источника высокочастотного электромагнитного излучения. Выбрав мощность N из соображений удобства проведения опытов, полученные результаты можно пересчитать для любого интересующего значения мощности N1. Соответствующее значение избыточной температуры АТ] определяется при этом из выражения:
N
АТ _ Т*
^п И
Условия моделирования процесса вытекают из необходимости выполнения равенства в модели и натуре безразмерных параметров, входящих в правую часть (15). Если модель пласта изготовлена из натуральных материалов, то этим обеспечивается выполнение условий:
(Л)мод _(Л)нат , (С)мод _(С)
(16)
Известно, что длина ВЧ ЭМ волны в пласте составляет десятки метров. Для соблюдения условий моделирования размеры эксперимен-
тальных моделей должны быть сравнимы с длиной ЭМ волны и больше. В лабораторных условиях модели таких размеров создать невозможно. поэтому условия моделирования процесса вытекают из необходимости равенства в модели и натуре безразмерных параметров - критериев подобия, входящих в правую часть (15). Для выполнения этого условия желательно иметь минимально допустимое число критериев, которые определяются количеством независимых размерностей и выбором основных единиц измерения.
применительно к нашей задаче условие равенства безразмерных критериев запишется в виде:
(Т *) _(Т *)
v /мод v /ш (ро)мод _(ро)и
( N *) _(N *)
v /мод v /а
(I *) _(I *)
v /мод \ /на
(17)
(18)
(19)
(20)
К перечисленным условиям добавляются требования, вытекающие из геометрического подобия внешних размеров моделируемой системы (пласта и окружающих пород), которые приняты неограниченными.
Из условия
(Я*) _(Я*) (21)
V /МОД \ /НАТ (21)
следует, что при уменьшении в модели толщины пласта в п раз во столько же раз в модели необходимо уменьшить радиус скважины. Равенства (18)-(20) служат для определения соответственных моментов времени и сходственных точек. Таким образом, появляется возможность решить следующую задачу: определить температуру пласта на заданной глубине высокочастотного электромагнитного прогрева. Используя полученные соотношения для параметров подобия с помощью анализа размерностей, мы можем рассчитать значение температуры на глубине 5 м.
Возникает необходимость проверки выполнения вышеперечисленных условий для натурных и модельных параметров. С этой целью проводились экспериментальные исследования температуры на модели пласта при сверхвысокочастотном электромагнитном воздействии.
Следует заметить, что высокочастотное и сверхвысокочастотное электромагнитные воздействия отличаются только способами их генерирования и ввода в рабочую среду, а нет принципиального различия в механизме взаимодействия с ней. Если мы пользуемся термином «сверхвысокочастотное» электромагнитное воздействие, то имеем в виду и высокочастотный диапазон.
Для экспериментальных исследований нагрева пористой среды, насыщенной высоковязкой нефтью, в сверхвысокочастотном электромагнитном поле разработана установка, схема кото -рой приведена на рис. 2.
Рис. 2. Схема экспериментальной установки. 1 - модель пласта; 2 - кровля; 3 - подошва; 4 - металлический кожух; 5 - термокарманы; 6 - «обсадная колонна»; 7 - «насосно-компрессорная труба»; 8 - диэлектрические шайбы; 9 - цанговое устройство; 10 - термопары; 11 - многоконтактный переключатель; 12 - микровольтметр; 13 - сосуд Дьюара; 14 - перфорированная диэлектрическая трубка; 15 - электромагнитный излучатель; 16 - металлическая крышка; 17 - сверхвысокочастотный генератор; 18 - фидер.
Установка включает металлический корпус, трехслойную модель пласта, «скважину» с помещенным в нее электромагнитным излучателем, а также источник электромагнитной энергии, контрольно-измерительные приборы.
Моделью пласта служил кварцевый песок, насыщенный битумной нефтью Мордово -Кармальского месторождения Республики Татарстан. В качестве кровли и подошвы использовалась бентонитовая глина, замешанная на минерализованной воде (25% поваренной соли) до тестообразной консистенции. Радиус модели опре-
делялся из условия r > r = — = 05
^ min '
Ь О
Ws'tgS
, где
^о - длина электромагнитной волны в вакууме, м; га - коэффициент поглощения электромагнитных волн в пласте, м-1; в', tgS - диэлектрическая проницаемость и тангенс угла диэлектрических потерь модели пласта.
В качестве натурных (промысловых) и модельных (лабораторных) данных использовались параметры, представленные в таблице.
Эти параметры применительны к Мордово-Кармальской залежи битумов. Значения натурных величин р, Т0 взяты из работы [2], все остальные значения, в том числе модельных величин, представлены в работе [5]. Используя
результаты лабораторных исследований распределения температуры в модели пласта, представленные в работе [6], получили, что коэффициент подобия равен 100.
В условиях эксперимента В = 7,5, tgS = 0,05, которые определялись экспериментально. Следовательно, радиус модели составляет 0,147 м.
Для того чтобы исключить влияние стенки модели на распределение электромагнитного поля из-за отражения от нее волны, радиус модели выбирался равным 0,25 м. Толщина подошвы и кровли составляла 0,28 м, что также исключает влияние дна и крышки модели пласта при отражении электромагнитных волн. Толщина пласта выбрана из условия согласования сверхвысокочастотного генератора с излучателем и составляла 0,08 м.
Одним из основных элементов установки является устройство ввода энергии электромагнитного поля в пласт через скважину. Устройство содержит тройник с короткозамыкающим поршнем, коаксиальную линию передачи, излучатель. Коаксиальная линия передачи состоит из внешней (латунной) и внутренней (медной) трубок. Все узлы устройства, по которым передавали в пласт энергию электромагнитного поля, имеют волновое сопротивление 50 Ом и сечение С = 0,008 м, Б = 0,0215 м (где С - внешний диаметр внутреннего проводника; Б - внутренний диаметр внешнего проводника). Узлы коаксиальной линии передачи центрированы с помощью фторопластовых шайб.
Электромагнитным излучателем служит система, состоящая из концов внешней трубки и выступающей за ней внутренней трубки.
Источником электромагнитной энергии является сверхвысокочастотный генератор «Парус» с вырабатываемой частотой 2375 МГц. Подаваемая на линию передачи мощность измерялась калориметрическим методом и составила около 450 Вт.
Вначале пласт полностью насыщался битумной нефтью.
Методика эксперимента заключалась в регистрации распределения температуры при сверхвысокочастотном электромагнитном воздействии. Температура в пласте и на границах пласта с кровлей и подошвой измерена медно-константановыми термопарами. Результаты измерений обрабатывались в безразмерных величинах (12).
На рис. 3 показаны температурные профили в пласте в различные моменты времени, на рис. 4 -изменение безразмерной температуры в зависимости от безразмерного времени.
Распределение температуры в пласте по его толщине неравномерное: температура в середине
а
Физические параметры пласта в натуре и модели
Параметры к, м /, Гц Яп, Вт/м-К а, м-1 СП, Дж/кг-К кг/м3 М0, Вт т0, °к
Натурные 8 13,56 1,75 0,036 2133 961 15-103 281
Модель 0,08 2400 1,1 3,6 1220 2050 95 293
г*-10
4 2 0 -2 -4 -6
Рис. 3. Распределение температуры в призабойной зоне пласта при электромагнитном воздействии на безразмерном расстоянии Я* = 0,44 при различных значениях параметра Фурье: 1 - 0; 2 - 0,009; 3 - 0,0250; 4 - 0,032.
9-103
12 10 8 6 4 2 0
Рис. 4. Изменение безразмерной температуры в зависимости ях от излучателя: 1 - 0,88; 2 - 0,44.
пласта больше, чем в кровле и подошве. С увеличением времени воздействия эта неравномерность увеличивается (рис. 3).
Из рис. 4 следует, что вначале наблюдается интенсивный рост температуры в пласте, который со временем замедляется, а в дальнейшем температура стремится к установившемуся значению.
Показанные данные (рис. 3 и 4) свидетельствуют об убывании температуры в пласте с расстоянием. Изменение температуры через равные интервалы расстояния неодинаково. Вблизи скважины это изменение незначительно, затем все более растет и к периферии прогретой зоны уменьшается. По мере удаления от скважины температура в середине пласта значительно отличается от температуры в кровле и подошве пласта.
При использовании безразмерных параметров (12), критериев подобия (17)-(21) и результатов
от безразмерного времени на разных безразмерных расстояни-
модельных исследований (рис. 3 и 4) были найдены значения температуры в реальном (натурном) пласте. На рис. 5 показано распределение температуры в пласте. Как видно из данного рисунка, на расстоянии 5 м изменение температуры составляет около 1К, что находится в пределах ошибки измерений температуры термопарой в промысловых условиях.
<
Рис. 5. Зависимость изменения температуры в пласте: 1 - с учетом уменьшения мощности ВЧ излучения по скважине и согласования генератора с ней; 2 - без их учета.
По кривой 1 на рис. 5 мы можем определить изменение температуры на произвольном расстоянии от излучателя.
ВЫВОДЫ
1. Получены безразмерные критерии подобия, характеризующие нагрев нефтяного пласта в сверхвысокочастотном электромагнитном поле.
2. Описана лабораторная установка, позволяющая экспериментально исследовать влияние сверхвысокочастотного электромагнитного поля на изменение температуры в пористой среде, насыщенной нефтью. Электромагнитное поле в модель пласта вводится через межтрубное пространство скважины. Критерии подобия позволяют перенести результаты исследований и на высокочастотный диапазон.
3. Выявлено немонотонное изменение температуры по толщине пласта: середина пласта нагревается значительно больше и быстрее, нежели вблизи кровли и подошвы пласта.
4. При использовании результатов лабораторных исследований и безразмерных критериев подобия показана возможность определения глубины теплового воздействия на пласт в высокочастотном электромагнитном поле.
5. Изложенные результаты могут быть использованы для обоснования эффективности применения энергии высокочастотного и сверхвысокочастотного электромагнитного полей с целью теплового воздействия на нефтяной пласт, интенсификации притока углеводородной жидкости в скважину, а также для теоретических исследований изучаемого процесса и обоснования достоверности их результатов.
6. Материалы данной работы можно успешно применить и для анализа результатов исследований, связанных с моделированием других тепловых методов воздействия на нефтяные пласты [8, 9].
ЛИТЕРАТУРА
1. Седов Л.И. Методы подобия и размерности в механике. М.: Наука, 1972. 440 с.
2. Малофеев Г.Е. О параметрах подобия нагрева пласта при тепловой обработке скважин.
Физическое и математическое моделирование механизмов нефтегазоотдачи. М.: Наука, 1981. 21-26.
3. Carrizales M.A., Larry W. Lake, Johns R.T. Production Improvement of Heavy-Oil Recovery by using Electromagnetic Heating. SPE Annual Technical Conference and Exhibition, 21-24 September, Denver, Colorado, USA. 2008. SPE-115723-MS. http://dx.doi.org/10.2118/115723-MS
4. Фатыхов М.А., Фатыхов Л.М. Особенности фазовых переходов в сверхвысокочастотном электромагнитном поле. Научное обозрение. 2014, (8-2), 576-580.
5. Фатыхов М.А., Худайбердина А.И. Комбинированные методы воздействия на нефтяные пласты на основе электромагнитных эффектов. Уфа: Изд-во БГПУ, 2010. 112 с.
6. Фатыхов М.А. Экспериментальное исследование фильтрации высоковязкой нефти в электромагнитном поле. Известия вузов: Нефть и газ. 2006, 6, 53-56.
7. Филиппов А.И., Михайлов П.Н., Иванов Д.В. Температурное поле радиоактивных изотопов в пористой среде. Теплофизика высоких температур. 2010, 48(1), 96-104.
8. Барышников А.В., Стрекалов А.В., Грачев С.И. Повышение нефтеотдачи пластов за счет воздействия электромагнитным полем. Научное обозрение. 2013, (4), 90-92.
9. Fatykhov M.A., Fatykhov L.M. Microwave Electromagnetic Method of Melting the Paraffin Plug in an Open Coaxial System. J Eng Phys Thermophys. 2015, 88(3), 724-729.
Поступила 13.10.14 После доработки 23.03.15 Summary
The criteria of similarity describing the process of heating of an oil layer in a high-frequency electromagnetic field are received. A possibility of their application for the determination of the depth of heating of an oil layer is shown.
Keywords: high-viscosity dielectric liquid, high-frequency electromagnetic field, method, theory of dimension, oil layer, temperature.
