УДК: 532.54
A.B. Трусов1, М.Н. Овчинников1, Е.А. Марфин2
'Казанский федеральный университет 2Исследовательский центр проблем энергетики КазНЦ РАН
Alexander.trusov@ksu.ru
ОСОБЕННОСТИ РАСПРОСТРАНЕНИЯ И ХАРАКТЕРИСТИКИ ФИЛЬТРАЦИОННЫХ ВОЛН ДАВЛЕНИЯ ПРИ ИСПОЛЬЗОВАНИИ ЛОКАЛЬНО-НЕРАВНОВЕСНЫХ МОДЕЛЕЙ
В работе рассматриваются основные линейные математические модели фильтрации жидкостей в пористых средах, приводятся расчеты фазовых и групповых скоростей для используемых уравнений эволюции, результаты численного моделирования значений этих скоростей при различных соотношениях релаксационных параметров, сопоставляются полученные результаты с экспериментальными данными по фильтрации жидкости в однородных высокопористых средах в условиях волн давления.
Ключевые слова: гидродинамика, фильтрационные модели, волны давления, фазовая скорость.
Введение
Гидродинамические исследования скважин являются на сегодняшний день необходимым инструментом для получения информации о фильтрационных параметрах в межскважинном интервале, таких как гидропроводность и пьезопроводность. Информация о значениях этих параметров важна для создания технологии разработки всего месторождения, поэтому определению этих параметров, в частности нестационарными гидродинамическими методами уделяется большое внимание.
Существует несколько модельных подходов, используемых для описания фильтрации жидкостей в пористых средах, но все они имеют свои пределы применимости. Одним из способов идентификации моделей фильтрации может послужить метод фильтрационных волн давления (ФВД), суть которого заключается в создании периодических возмущений давления или дебита в скважине и отслеживании отклика на данное возмущение в соседних скважинах (гидропрослушивание) либо в самой скважине (самопрослушивание) (Бузинов, Умрихин, 1964, Чека-люк, 1948).
С 60-х годов прошлого века под руководством профессора Непримерова H.H. было проведено огромное количество исследований методом ФВД различных коллекторов. Казанской школой подземной гидродинамики накоплен огромный промысловый опыт исследований пластов методом фильтрационных волн давления.
При этом вопросы идентификации моделей фильтрации остаются актуальными и, метод ФВД здесь может быть также использован. Действительно, по дисперсионным соотношениям, зависимостям фазовой скорости распространения волн давления от частоты, представляется возможным диагностировать адекватность конкретному случаю тех или иных моделей.
Рассмотрим наиболее распространённые модели в случае плоскорадиальной фильтрации:
Модель 1
Классическая модель упругого режима фильтрации (Щелкачев, 1959):
к
W(r,t) = -~Vp(r,t),
где "(гД) - скорость фильтрации, к - проницаемость пористой среды, ц - вязкость фильтрующейся жидкости, У(р(гД) - градиент давления.
Уравнение для давления (пьезопроводности):
dp(r,t)
1 д dp(r,t)
= У--г- .
дь Л г дг дг
где х=к/цв - пьезопроводность среды, в = швж + вс - сжимаемость пласта, ш - пористость среды, в - сжимаемость фильтрующейся жидкости, в - сжимаемость пористого скелета.
Данное уравнение относится к диффузионному типу и для случая точечного стока (источника) решение для изменения давления в неограниченном пласте, начиная с момента 1=0, при задании произвольного дебита д(1) на длительно простаивавшей одиночной узкой скважине запишется в виде свёртки:
1 г* ехР (ш7)
где Е=кЬ/ц - гидропроводность среды, Ъ - толщина пласта.
В случае создания периодических изменений давления в одной или нескольких скважинах, в пласте будут распространяться волны давления с фазовой и групповой скоростью: ___
Уфаз = л/2ХШ. Кр = где ю - циклическая частота распространения волны.
Модель 2
Модель фильтрации в однородных пористых средах, учитывающая инерционность движения флюида (Хрис-тианович, 1981). Закон фильтрации:
, „ дУУ к
где тч - коэффициент релаксации по дебиту (скорости фильтрации), имеющий размерность времени.
Уравнение для давления:
др(г,0 д2р(г, г) 1 д др(г, О
+ т,
dt ' dt2 ХгдгГ~ дг
Выше представленное уравнение относится к телеграф-
ному типу с конечной скоростью распространения возмущений. Решение данного уравнения в виде свёртки:
2cosh
9 Г
£--z
r2j ~х
2т„
\
dt'
/
В условиях реализации метода ФВД фазовая и групповая скорости равны:
^фаз
V =
"гр
(1 + co2Tq2y/4cos (jarctg-^y )1/4 cos ^ arctg
(1 + ш2т2
1 со
2J*
+ ^ I-(1 + 0)2T42)
шт.
я -1
3/4 cos (— arctg-)
i Jf(i+(°2T"2r3/4sin^arctg^h
Модель 3
Модель фильтрации однородной капельно сжимаемой жидкости в трещиновато-пористых, кавернозных коллекторах (Баренблатт, 1984). Закон фильтрации:
где тр - коэффициент релаксации по давлению, имеющий размерность времени, определяемый вязкостью жидкости, упругоёмкостью блоков пористой среды и характеризующий обмен жидкостью между блоками и трещинами.
Уравнение для давления:
dp(r,t) dt
Эр(г,С)>
Решение данного уравнения в виде свёртки:
1 f'flft-t'H" /
t'xz2
rZ + ^Z2Tp
dzdt'
Фазовая и групповая скорость будут рассчитываться исходя из следующих уравнений:
У^(1 + <о2тр2)4
Kp =
1 + 2
'фаз д П
C042arCt5^J
1 /1 1 \ —== (1 + CO2T„2)~1/4 cos -arctg-
— (1 + <u2t„2)"5/4 cos (^ arcto —) от 2
X v V2
(1 + oj2tp2)~5/4 sin (^ arctg —) тр \2 й)тру p
Модель 4
Двухрелаксационная модель в трещиновато-пористых пластах учитывающая сжимаемость трещин (Молокович и др., 1980). Закон фильтрации:
дш к_ / др\
, „ aw к / ор\
Уравнение для давления:
dp(r,t) d2p(r,t) 1д д
dt
- + т,
dt2
1 d О f , ■ч
+ т.
dpOvty, ' dt J
Решение данного уравнения в виде свёртки:
хехр{-£)
1 + z2 4х
\dzdt'
Фазовая и групповая скорость для этой модели будут равны:
т/ __«1 + цу _
'фаз / / 2
(1 + ^2^ Г
7 =
'гр
1 fl + 0)2Tр2\1/4 /1 1 + Ы2Тр
2Vxco \1 + ш2т,2/ C0S \2 arCt9 фр-тч) ,
1 (1 + id2T2)-3/4<UTp2 1
+ \2~ 1 + w2T2
(1+&)2t42)4 J cos [ — arctg
--(1 + й)2Тр2)1/4й)Т 2
й)(тр_т/ 2 hr^l + w2!,^
/1 1 + (i)2TpTq\ (Тр-Тч)(1-й)2ТрТч)
sin - arctg —-- —--—--5-
\2 «(тр_тч) / ш2(тр_тч) + (1 + o)2TpTqy
Модель 5
Модель фильтрации, учитывающая сложность структуры горных пород, релаксацию переупаковки частиц коллектора при изменении давления (Молокович, 1980). В этом случае возможно запаздывание в установлении равновесного состояния напряжения пористой среды а и давления р. При этом предположении основная гипотеза о постоянстве горного давления рг запишется в виде:
дс др
Рг = ^ + Гп- + Р + Тй-,
где тп - время релаксации напряжения пористой среды, та - время релаксации по давлению.
Закон фильтрации:
ь-т
W(r,t) = --TjL
V-tq
Vp + (1 - f \р - ро) VpdT V я/ h Ъ
Уравнение для давления:
•— научно-технический журнал я
ЯДЩД 451
др(г, t) / mßxтп + ßcTd\ d2p(r, t)
dt +Vq+ ß ) dt2 *t
dt ■ VT4+ У dt2
mßxzn + /?сТсЛэ3р(ло _ ß ) dt3
dp(r,t) d2p(r,t)
+
* г дгГ dr
p(r,t) +
Параметр Образец №1 Образец№2
ßc 15,5 10*Па1 4,5 10 -*Па'
т 0,35 0,35
к 170 мкм2 27 мкм2
X \,\м2/с 0,6м2/с
+ (тр + т.
dt
■ + TpTn
dt2
Таблица.
Ниже представлены графики зависимости фазовых и групповых скоростей от частоты для первых четырёх моделей фильтрации при различных временах релаксации (Рис. 1, 2).
Исследуя различия в зависимостях фазовых скоростей от частоты при распространении волн давления, можно провести верификацию фильтрационных моделей. Так на приведённых выше рисунках отчётливо видно, что с ростом частоты задаваемых колебаний давления, а также с увеличением времён релаксации тр и т^ модельные кривые имеют всё более значимые расхождения в фазовых скоростях, что может послужить хорошим диагностическим механизмом при дальнейшем выборе математической модели исследуемого процесса. К примеру, классическая модель упругого режима фильтрации (модель 1) и модель фильтрации однородной капельно сжимаемой жидкости в трещиновато-пористых, кавернозных коллекторах (модель 3) с ростом частоты приводит нас к бесконечной скорости распространения возмущений в пласте, что противоречит экспериментальным данным (из про-
т = 10сек, 50сек< т <100сек
т = 100сек, 50сек< т <100сек
т = 10сек,50сек< т <100сек
т = 100сек,50сек<т <100сек
Рис. 1. Зависимость фазовой скорости от частоты для различных моделей фильтрации. Поверхностижелтого цвета - модель1, поверхности красного цвета - модель 2, поверхности зеленого цвета - модель 3, поверхности синего цвета - модель 4.
мысловых исследований известно, что скорость распространения периодических возмущений в продуктивных коллекторах конечна и её значение, к примеру, для девона Ромашкинского месторождения равно 1-10 см/се-к).Итак, различия в дисперсионных соотношениях фильтрационных гармонических волн давления лежат в основе инструмента по исследованию границ применимости фильтрационных моделей.
Заметим, что при гидродинамических исследованиях пористых сред методом ФВД следует выделять частотные интервалы гидродинамических воздействий, в пределах которых можно говорить, собственно, о фильтрации жидко-сти(Овчинников, 2003). Действительно, при частотах изменения давления превышающих УсЫ1М где У - характерная скорость распространения упругих колебаний, d -размер пор и межпоровых каналов, мы сталкиваемся с нарушением условий континуального перехода к описанию фильтрации в рамках сплошной среды и имеем дело с задачей распространения ультразвуковых волн в пористой среде. Поэтому мы будем проводить исследования пористых сред методом фильтрационных волн давления с периодами единиц секунд, минут, то есть задавать колебания давлений на частотах порядка 10-2-100 Гц, имея в виду, прежде всего фильтрационно-гид-родинамические, а не упруго-акустические эффекты.
Экспериментальные исследования фильтрации в однородных пористых средах по методу ФВД проводились на сконструированной в лаборатории физической динамики гетерогенных сред Казанского федерального университета оригинальной установке, схематический вид которой приведён на рис. 3, устройство и принцип работы описан в (Трусов, 2010).
На разработанной установке было проведено более 40 экспериментов по реализации метода фильтрационных волн давления в лабораторных условиях, были рассчитаны фазовые скорости для различных частот воздействия волн давления, проведено исследование применимости диффузионного (модель 1) и телеграфного (модель 2) типов уравнений для описания нестационарной фильтрации.
Ниже представлены графики зависимости фазовой скорости распространения волн давления от квадратного корня из частоты для исследуемых пористых сред (Рис. 4, 5): кривые красного цвета - экспериментальные данные, кривые
шШШЕ-
зелёного цвета - расчётные кривые по модели 1, кривые синего цвета - расчётные кривые по модели 2. Для модели 2 подбиралось время релаксации для первого образца = 0,8 сек для второго = 3 сек, также показано, что при больших временах релаксации экспериментальные данные и модельный тренд имеют значительное расхождение. Параметры исследуемыхпористых сред представлены в таблице.
Проведённые исследования пределов применимости диффузионного и телеграфного типов уравнений для описания нестационарных гидродинамических процессов в насыщенных пористых средах показали, что при данных значениях пористости, сжимаемости и проницаемости пористых сред, а также вязкости вытесняющей фазы, уравнение диффузионного типа (модель 1) и телеграфного типа (модель 2) при времени релаксации для первого образца тч = 0,8 сек для второго тч = 3 сек с высокой точностью описывает фильтрационный процесс, но уже при значении времени релаксации вдвое больших приводит к значительному несовпадению расчётных и экспериментальных данных. Другими словами можно сказать, что релаксационными эффектами при исследовании нестационарной фильтрации жидкостей с малой вязкостью в однородных высокопроницаемых пористых средах можно пренебречь. Следует отметить, что в низкопроницаемых пористых материалах времена релаксации имеют порядок < 103
т = 10сек,10сек< т <30сек
т = 100сек,10сек< т <30сек
т
р
10сек,10сек< т <30сек
т
р
Рис. 2. Зависимость групповой скорости от частоты для различных моделей фильтрации. Поверхности жёлтого цвета - модель 1, поверхности красного цвета - модель 2, поверхности зеленого цвета - модель 3, поверхности синего цвета - модель 4.
Рис. 3. Схематический вид установки.
сек (Динариев, Николаев, 1990).
Из вышеизложенного можно сделать вывод что, метод фильтрационных волн давления может являться диагностирующим инструментом в области верификации фильтрационных моделей. В приведённом исследовании показано, что применимость локально-равновесной модели фильтрации при описании нестационарного гидродинамического процесса в высокопористой среде себя оправдывает, то есть, модель классического упругого режима фильтрации адекватно описывает данный процесс фильтрации.
Отметим, что в последнее время был получен ряд новых результатов в исследованиях волн давления. Так:
1. Предложен способ определения нелокальности фильтрационных моделей, основанный на анализе дисперсионных соотношений фильтрационных гармонических волн давления (Овчинников, 2002),
2. Рассмотрены волны давления, их дисперсионные соотношения, в насыщенных пористых пластах, описываемые уравнениями в дробных производных (Овчинников и др. 2004),
3. Предложены расчетные соотношения для оптимизации перетоков в системе блоки-трещины при периодическом изменении давления в условиях неоднородных пластов (ОусЫпшкоу & 2ау1ёопоу, 2008),
4. Систематизированы проблемы проведения экспериментов по организации волн давления и их автоматизации (Одиванов и др., 2009; Оди-ванов и др., 2010),
5. Рассчитаны фазовые скорости волн давления для различных моделей, исследованы вопросы деформации горных пород и возможности регистрации смещений и деформаций в приповерхностных слоях в виде периодических сигналов, вызванных нестационарными изменениями давления в тонких пористых и трещиновато-пористых пластах (ОусЫпшкоу е! а1., 2001),
6. Рассмотрены особенности распространения волн давления в нелинейных по проницаемости насыщенных пористых средах (Куштанова, Овчинников, 2010).
100сек, 10сек< т <30сек
научно-технический иотнал
1(46)2012 I еоресурсы
Рис. 4. Зависимость фазовой скорости от корня от часто-тыдля образца №1. Красная кривая - экспериментальные данные, кривые зелёного и синего цвета - расчётные кривые по модели 1 и 2 соответственно.
1,10
0,90
0,70
0,50
0,30
Рис. 5. Зависимость фазовой скорости от корня от часто-тыдля образца №2. Красная кривая - экспериментальные данные, кривые зелёного и синего цвета - расчётные кривые по модели 1 и 2 соответственно.
Заключение
Фильтрационные волны давления могут использоваться не только как способ определения фильтрационно-ем-костных характеристик продуктивных коллекторов, неодно-родностей в призабойной и межскважинной зонах, но и как инструмент для верификации фильтрационных моделей, определения их границ применимости при описании того или иного нефтеводонасыщенного пласта.
Работа выполнена в рамках Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2012 (гос. контракт №П1005, №П511).
Литереатура
Баренблатт Г.И., Ентов В.Н., Рыжик В.М.. Движение жидкостей и газов в природных пластах. М: Недра. 1984. 211.
Бузинов С.Н., Умрихин И.Д.. Исследование пластов и скважин при упругом режиме фильтрации. М.: Недра, 1964. 273.
Динариев О.Ю., Николаев О.В. О релаксационных процессах в низкопроницаемых пористых материалах. ИФЖ. 1990. Т.58, №1. 78-82.
Куштанова Г.Г., Овчинников М.Н. Спектральные особенности фильтрационных волн давления в нелинейных средах. Известия вузов. Поволжский регион. Физ.-мат.науки. 2010, №1. 130-135.
Молокович Ю.М., Непримеров Н.Н., Пикуза В.И., Штанин A.B.. Релаксационная фильтрация. Казань, Изд-во Казанского университета. 1980. 136.
Непримеров Н.Н. Трехмерный анализ нефтеотдачи охлажденных пластов. Изд-во КГУ, Казань. 1978. 216.
Овчинников М.Н. Интерпретация результатов исследований
пластов методом фильтрационных волн давления. Казань: Новое знание. 2003. 84.
Овчинников М.Н. Об одном методе идентификации моделей фильтрации. Известия вузов. Нефть и газ. 2002. № 4. 22 - 25.
Овчинников М.Н., Гаврилов М.Н., Непримеров Н.Н., Штанин A.B. Разработка нефтяного месторождения как комплексная междисциплинарная технология. Наукоемкие технологии. 2004. №4. т.5. 75 - 80.
Одиванов В.Л., Гаврилов А.Г., Овчинников М.Н.. Автоматизация исследований флюидонасыщенных пластов методом волн давления. Учён. Зап. Казан.ун-та. Сер. Физ. матем. науки. 2010. Т. 152, кн.1. 73-79.
Одиванов В.Л., Овчинников М.Н., Гаврилов А.Г. Автоматизация, телемеханизация и связь в нефтяной промышленности. Специфика и средства обработки данных исследований нефтяных пластов методом фильтрационных волн давления. 2009. №1. 16-23.
Одиванов В.Л., Овчинников М.Н., Гаврилов А.Г. Специфика задач автоматизации гидродинамических исследований на нефтяных месторождениях. Автоматизация, телемеханизация и связь в нефтяной промышленности. 2009. №1. 9-23.
Трусов А.В. Экспериментальное гидродинамическое исследование пористых сред методом фильтрационных волн давления. Нефть. Газ. Новации. №2. 2010. 67-69.
Христианович С.А. Механика сплошной среды. М., Наука. 1981. 493.
Чекалюк Э.Б. К анализу методов исследования скважин. Нефтяное хозяйство. 1948. №11. 27-30.
Щелкачев В.Н. Разработка нефтеводоносных пластов при упругом режиме. М.: Гостоптехиздат. 1959. 467.
Ovchinnikov M.N., Zavidonov A.Yu. Searching for an optimum frequency. Georesources. 2008. N1. 21-22.
Ovchinnikov M.N.,.Zavidonov A.YuEvtushenko, S.P.. Deformations of Rocks in Periodic Regimes of Filtration. Journal of Engineering Physics and Thermophysics. 2001. v.44. number 5. 10671072.
A.V. Trusov, M.N. Ovchinnikov, E.A. Marfin. Filtration waves of pressure distribution peculiarities and characteristics during local unbalanced models usage.
The paper is considered with the main linear mathematical models of liquids filtration in porous medium, calculations of phase and group velocities for applied equations of evolution, results of numerical modeling of these velocities values at various ratios of relaxation parameters are given, acquired results with experimental data for liquid filtration in homogeneous highly porous medium in the waves of pressure conditions are correlated.
Key words: hydrodynamics, filtration model, wave pressure, phase velocity.
Александр Валерьевич Трусов
Аспирант Института физики Казанского федерального университета. Научные интересы: гидродинамика, разработка нефтяных месторождений.
Марат Николаевич Овчинников
Д.ф.-м.н.,профессор Казанского федерального университета. Научные интересы: динамика сложных систем, молекулярно-динамическое моделирование, нестационарная фильтрация.
Евгений Александрович Марфин
К.т.н.,в.н.с.исследовательского центра проблем энергетики КазНЦ РАН. Научные интересы: интенсификация процессов добычи нефти, акустика природных сред.
420008, Россия, Казань, ул. Кремлевская, д.16А.
Тел.: (843) 233-70-02, (927) 419-35-45.
,,н/сек