С увеличением значения Kife>/ возрастает миграция массы вещества к поверхности S, а следовательно, снижается темп затвердевания области Кц,о. Роль этой миграции может быть доведена до минимума при быстром затвердевании изделия, то есть с увеличением скорости продвижения поверхности S.
ЛИТЕРАТУРА
1. Акимов И.А., Акимов А.И., Тюков Н.И. Математическая модель процесса изготовления изделий из композиционных материалов в результате полимеризации // Труды Стерлитамакского филиала Академии наук Республики Башкортостан, 2001, Вып. 2. Уфа: Гилем. С. 6 - 9.
Поступила 15.01.03
Summary
Criteria of similarity have been obtained which describe the polymerization process in multi-layer structures made of composite materials in an electric field. The paper has determined the effect of individual criteria upon the process dynamics.
О.Н. Сизоненко, Г. А. Барбашова, О.В. Хвощан
К ВОПРОСУ О РАСПРОСТРАНЕНИИ ВОЛНЫ НАПРЯЖЕНИЙ В ПЛАСТЕ ПРИ ЭЛЕКТРИЧЕСКОМ РАЗРЯДЕ В СКВАЖИНЕ
Институт импульсных процессов и технологий НАН Украины, пр. Октябрьский, 43 А, г. Николаев, 54018, Украина
Теоретическими и экспериментальными исследованиями, проведенными в различных странах в последние годы, установлено, что наиболее эффективные технологии разработки трудноизвлекае-мых запасов нефти основываются на использовании различных физических полей. Накопленный нами опыт использования высоковольтного электрического разряда в скважине для воздействия на продуктивные пласты показывает, что при оптимальном выборе объектов обработки и применяемых технологических приемах можно заметно интенсифицировать фильтрационные процессы в пластах и повысить их нефтеотдачу [1 - 3].
Согласно существующим представлениям о механизме электроразрядного воздействия на пласт, превалирующее действие на изменение фильтрационных характеристик в пористой насыщенной среде нефтяного коллектора оказывают волны сжатия, создающие и развивающие радиальные и цилиндрические трещины [4]. В условиях скважины картина распространения волн сжатия из скважины в пласт усложняется ввиду наличия обсадных труб и цементного кольца между трубой и породой продуктивного коллектора. Поэтому значительный интерес представляет исследование воздействия волн сжатия при высоковольтном электрическом разряде в скважине на призабойную зону скважин, представленное в данной работе.
Для решения поставленной задачи использовались методы физического и математического моделирования. Характеристики электрического разряда получены в камере высокого давления, которая является частью экспериментального стенда для физического моделирования скважинных условий. Схема разрядной камеры электроразрядного устройства (ЭРУ) приведена на рис. 1.
При построении математической модели распространения волн сжатия в скважине приняты следующие допущения: стенки электроразрядного устройства и скважины абсолютно жесткие; канал разряда, образованный в жидкости при ее пробое [5], в начальный момент времени имеет форму прямого кругового цилиндра, высота которого равна длине межэлектродного промежутка, а ось симметрии совпадает с осью симметрии скважины и ЭРУ; стенка канала непроницаема; скважина заполнена идеальной сжимаемой жидкостью, а канал разряда - идеальной плазмой.
© Сизоненко О.Н., Барбашова Г.А., Хвощан О.В., Электронная обработка материалов, 2003, № 4, С. 51-56.
Рис. 1. Схема разрядной камеры электроразрядного устройства: 1, 2 - электроды, 3 - изолятор, 4 - "окно"разрядной камеры, 5 - стенка скважины, 6- обсадная колонна, 7 - цементное кольцо, 8 - продуктивный пласт
В соответствии с принятыми допущениями математическая постановка задачи состоит в следующем. В области, внутренней границей которой является контактный разрыв плазма - жидкость (стенка канала разряда), а внешней - неподвижная жесткая стенка (поверхности скважины и ЭРУ), необходимо решить систему двумерных нелинейных уравнений газовой динамики, записанных в цилиндрической системе координат и представляющих собой законы сохранения массы, импульса и энергии [6]:
д(г^), ) )
ы
&
дг
= ¥Л.
(1)
где
р1 = ^ р^ , Pvг, е]Т;
р2 = [РЛ^^, руг2 + p, (е + р) ]Т ;
рз = [рл ,р^г2 + р^л, (е+р К ]Т; ^ = [0,р,0,0]Т;
t - время; г, г - цилиндрические координаты; уг, у2 - радиальная и осевая компоненты вектора скоро-
+ (+ V?)/2
сти жидкости; р - давление; р - плотность; е = [6], определяемая по формуле
р - со? ( Р " Ро)
8 =
Р (к -1)
; в - удельная внутренняя энергия
(?)
где р0 - плотность покоящейся жидкости; с0 - скорость звука в покоящейся жидкости; к = 7,15.
На внутренней границе расчетной области задается давление, определяемое при решении уравнения баланса энергии в канале разряда [5]:
+=я (,),
Л у -1 Л V ;
л (3)
ж у -
где рк - давление в канале разряда; Ук - объем канала разряда; у = 1,26; N(1) - скорость ввода энергии в канал разряда.
На внешней границе расчетной области - неподвижной жесткой стенке - ставится условие непротекания [6]:
Уп= 0,
(4)
у„ - нормальная составляющая вектора скорости жидкости.
Начальные значения гидродинамических характеристик соответствуют невозмущенному состоянию жидкости, в начальный момент времени скорость расширения канала разряда равна нулю, его радиус равен 0,1 мм, а давление в канале разряда превышает давление в окружающей жидкости на 1 МПа.
Система уравнений (1) - (2) решается конечноразностным методом С.К. Годунова [6]. Используется подвижная расчетная сетка, при построении которой учитываются форма внутренней поверхности разрядной камеры и расширение канала разряда. Давление в канале разряда определяется из уравнения (3) двухшаговым методом Эйлера. Объем канала разряда и его производная по времени вычисляются по координатам поверхности канала и скорости его расширения, полученными при решении задачи о распаде произвольного разрыва [6] на контактной границе плазма-жидкость. Алгоритм решения внутренней задачи гидродинамики при электрическом разряде в воде приведен в работе [7].
Функция N(0 в уравнении (3) определяется по формуле
N(0 = ы(0 • 7(0, (5)
где н(0 - напряжение в канале разряда; ¡(() - ток в канале разряда.
Значения тока и напряжения получены экспериментально в ЭРУ, схема разрядной камеры которого приведена на рис. 1. Начальное напряжение на конденсаторной батарее, имеющей емкость 2,4 мкФ, равно 30 кВ, индуктивность электрической цепи - 3 мкГн, расстояние между электродами -0,02 м. Исследования проводились при значениях гидростатического давления 0,1; 10; 20; 30; 40 МПа. Рабочей средой в разрядной камере являлась вода и водонефтяная эмульсия (ВНЭ), обводненность которой составляла 30%.
Обводненность О и плотность р ВНЭ определялись соответственно методом Дина-Старка и пикнометрическим методом по методике [8].
Ранее [9] скорость звука с0 в ВНЭ определялась в зависимости от температуры без учета влияния гидростатического давления. В настоящей работе использованы экспериментальные данные работы [10] по зависимости скорости звука в ВНЭ от давления и температуры, поскольку значения плотности ВНЭ (р = 958 кг/м3 при 200С), используемой нами и приведенные в этой работе, совпадают. Зависимость скорости звука от гидростатического давления при температуре 200С (температура среды в условиях эксперимента) для ВНЭ с р = 958 кг/м3 приведена в табл. 1.
Таблица 1
рг, МПа 0,1 10 20 30 40
с0, м/с 1480 1580 1590 1600 1620
Для расчета продольных компонентов напряжений, возникающих в пласте при гидродинамической нагрузке, которые являются определяющими при создании трещин, использована формула [11]:
8С1С2С3Р1Р2Р3
о3, =
х/
хЯ
(С3Р3 + С2Р2 ) (( + С2Р2 ) (( + С0Р0 ) ) Г
1 (Г ~ Г3 ) ( ~ Г2 ) (Г2 -Г1 ) | (Г1 -Г0 ) С3 С2 С1 С0
(Г ~ Г3 ) ( ~ Г2 ) (Г2 -Г1 ) | (Г1 -Г0 )
Здесь с7 - скорости продольных упругих волн; р7 - плотности упругих сред; г - радиальная координата; г1, г2, г3 - соответственно внутренний и внешний радиусы обсадной трубы и внешний
радиус цементного кольца; /^) - давление в жидкости при г = г0; г0 - координата фиксированной точки в жидкости; Н (х) - единичная функция Хевисайда. Индексами 1, 2, 3 обозначены величины,
соответствующие стальной трубе, цементному кольцу и породе коллектора.
По этой же формуле определялось напряжение в стальной обсадной трубе и цементном кольце.
Продольные компоненты напряжения устанавливались в плоскости серединного сечения канала разряда по величине давления на стенку скважины, то есть г0 = гь а функция / (^) - временная зависимость давления на стенку скважины в этой же плоскости.
Характеристики упругих сред, которые использованы при выполнении расчетов по данным [12], представлены в табл. 2.
Таблица 2
Материал Скорость звука сь м/с Плотность рг-, кг/м3
Вода 1500 1000
Сталь 5000 7850
Цемент 2900 1800
Песчаник 3500 2500
Карбонат (известняк) 3600 2600
Согласно приведенным на рис. 2 и 3 временным зависимостям гидродидинамического давления на стенку скважины в плоскости серединного сечения канала разряда, амплитуда избыточного давления с ростом гидростатического давления уменьшается как в воде (рис.2), так и в ВНЭ (рис.3). Хотя необходимо отметить, что амплитудные значения давления в ВНЭ при рг свыше 20 МПа существенно (почти в 2 раза) выше, чем в воде, и зависимость амплитуды давления от рг значительно меньше. Разница между абсолютными значениями давления в воде и в ВНЭ при равных гидростатических давлениях составляет около 10%.
^р, МПа
□ 100 200 % мкс 300
Рис. 2. Давление на стенку скважины в плоскости серединного сечения канала разряда при электрическом разряде в воде. Рг, МПа: 1 -0,1; 2 - 10; 3 - 20; 4 - 30
Рис.3 Давление на стенку скважины в плоскости серединного сечения канала разряда при электрическом разряде в ВНЭ. Рг, МПа: 1 - 0,1; 2 - 10; 3 - 20; 4 - 30; 5 - 40
Та же тенденция зависимости от среды в скважине сохраняется и в величинах напряжений в упругих средах при прохождении в них упругих волн. Так, величины напряжений в обсадной трубе, цементном кольце и в породе (рис.4) при электроразряде в ВНЭ выше, чем в воде. Наибольшие потери механической энергии (в 2 раза) происходят при прохождении волной сжатия обсадной трубы (см. рис.4). В цементном кольце потери составляют около 15%, далее в породе напряжения экспоненциально затухают, и на расстоянии 0,3 м от стенки скважины разница в значениях амплитуды напряжений (в зависимости от среды в скважине) составляет не более 10% (рис.5).
Рис. 4. Максимальные значения напряжений в упругих средах при рг = 30 МПа, 1 - песчаник, разряд в ВНЭ; 2 - песчаник, разряд в воде; 3 - известняк, разряд в ВНЭ; 4 - известняк, разряд в воде
арг, МПа
3
1
□ 100 200 Ъ мкс 300
Рис. 5. Продольные напряжения на расстоянии 0,1 м от стенки скважины при рг = 30 МПа. 1 - песчаник, разряд в ВНЭ; 2 - песчаник, разряд в воде; 3 - известняк, разряд в ВНЭ;
4 - известняк, разряд в воде
(Ргг- МПа
\Ч 4 \4 \ —
0 100 200 С, мкс 300
Рис. 6. Разность между напряжениями, вызванными волной сжатия в ВНЭ и внутрипаровым давлением жидкости на расстоянии 0,1 м от стенки скважины. Рг, МПа: 1 - 10; 2 - 20; 3 - 30; 4 - 40
Напряжения, возникающие в известняках, несколько выше, чем в песчаниках, что объясняется их акустическими свойствами, однако тенденция затухания для этих пород одинакова как при электроразряде в воде, так и при электроразряде в ВНЭ (см. рис.5).
Влияние гидростатического давления в скважине на изменение напряжений в пласте сказывается более значительно, чем на давление в жидкости, при увеличении гидростатического давления с 10 до 40 МПа максимум гидродинамического давления снижается на ~ 20%, а напряжения (разность между напряжением, вызванной волной сжатия, и напряжением, вызванным в породе внутрипаровым давлением жидкости) - на ~ 30% (рис. 6).
Таким образом, выполненная работа позволяет дать оценку электроразрядного воздействия на призабойную зону скважины в зависимости от породы коллектора, пластового давления и жидкости, заполняющей скважину, и дает возможность управлять процессом воздействия в конкретных геолого-технологических условиях.
ЛИТЕРАТУРА
1. Максутов Р.А, Сизоненко О.Н., Малюшевский П.П. и др. Использование электровзрывного воздействия на призабойную зону скважины // Нефтяное хозяйство. 1985. № 1. С. 34 - 35.
2. Сизоненко О.Н., Швец И.С. Особенности электроразрядной технологии интенсификации фильтрационных процессов в нефтяных скважинах // Современные проблемы электрофизики и электродинамики жидкостей. Санкт-Петербург, 1998. С. 137 - 139.
3. Сизоненко О.Н., Швец И.С., Кучернюк А.В. Применение электроразрядного воздействия для обработки добывающих и нагнетательных скважин // Нефтяное хозяйство. 2000. № 12. С. 133 - 135.
4. Сизоненко О.Н., Любимов А.Д., Денисюк О.Н. Исследование влияния профиля волны сжатия электроразряда на изменение пористости и проницаемости осадочных пород // Физика и техника высоких давлений. 1997. № 1. С. 107 - 111.
5. Наугольных К.А., Рой Н.А. Электрические разряды в воде. М., 1971.
6. Численное решение многомерных задач газовой динамики / Под ред. С. К. Годунова. М., 1976.
7. Барбашова Г.А., Иванов А.В. Поведение тонкой цилиндрической оболочки под действием импульсной нагрузки // Электрический разряд в конденсированных средах. Киев, 1989. С. 55 60.
8. БороваяМ.С. Лаборант нефтяной и газовой лаборатории. М., 1968.
9. Сизоненко О.Н., Любимов А.Д., Денисюк О.Н., Хвощан О.Н., Петухов С.В. Исследование зависимости фильтрационных характеристик коллекторов от воздействия электрического разряда в водонефтяной эмульсии // Нефтяное хозяйство. 2002. № 1. С. 43 - 45.
10. Беляков В.Л. Автоматический контроль параметров нефтяных эмульсий. М., 1992.
11. Дворкин Я.П., Максутов Р.А., Сизоненко О.Н. О распространении волны напряжений в пласте при осесимметричном разряде в скважине // Тез. докл. III Всесоюзной научно-технической конференции «Электрический разряд в жидкости и его применение в промышленности». Киев, 1984.
12. Геофизические методы исследования скважин. Справочник геофизика / Под ред. В.М. Запорожца. М., 1983.
Поступила 23.01.03
Summary
The mathematical model of propagation and influence of stress waves on the hole bottom zone of oil wells is created and considered. Theoretical and experimental data about using the water electrolyte and water-oil emulsion as working fluid is used.