ТЕРМОДИНАМИЧЕСКОЕ ИССЛЕДОВАНИЕ МНОГОФАЗНЫХ ТЕЧЕНИЙ В ПРОЦЕССАХ РАЗРАБОТКИ МЕСТОРОЖДЕНИЙ УГЛЕВОДОРОДОВ
В.М. Максимов Институт проблем нефти и газа РАН, г. Москва, e-mail: [email protected]
1. Введение
Проектирование разработки месторождений углеводородов опирается в значительной степени на теорию многофазной фильтрации. Последняя включает основные законы
сохранения и постулат о справедливости закона Дарси для скорости wf фильтрации каждой фазы в качестве уравнения движения [1]:
wa =- kj^ V jPa, Pi - p2 = pc (s), (1)
где a = 1,2 (в двухфазном варианте); kj - тензор абсолютной проницаемости; pa - давление в фазах; по латинским индексам производится суммирование.
В соотношения (1) заложены две основополагающие гипотезы: для данной породы и пары фильтрующихся флюидов капиллярное давление pc и относительные фазовые
проницаемости (ОФП) ka являются универсальными, равновесными, стационарными,
скалярными функциями локальной насыщенности s одной из фаз.
Исследования последних десятилетий показывают, что во многих случаях эти гипотезы не подтверждаются. Нередко необходим учет эффектов неравновесности и нестационарности при измерениях материальных функций. В анизотропных средах, каковыми и являются реальные коллекторы, ОФП и капиллярное давление зависят не только от насыщенности, но и от направления течения и имеют тензорную природу. Большое значение имеет учет смачиваемости породы, которая меняется в процессе добычи нефти и газа. Индексом смачиваемости служит интеграл от капиллярного давления. Вместе с тем в уравнения Раппопорта - Лисса входит производная от капиллярного давления, что позволяет контролировать лишь изменение градиента насыщенности, при слишком больших значениях которого эти уравнения асимптотически переходят в модель Бакли - Леверетта.
Для обобщения и развития теории многофазной фильтрации авторами рассмотрены некоторые результаты применения методов неравновесной термодинамики к процессам, происходящим в гетерогенных средах.
2. Об обобщенном уравнении движения многофазных смесей
в пористой среде
В работах [2-4] с использованием методов термодинамики неравновесных процессов развита достаточно общая теория многофазного многокомпонентного химически активного континуума, рассматривающая ряд термодинамических, физико-химических, тепловых и других эффектов. Мы ограничимся анализом уравнений движения для многофазных течений несмешивающихся жидкостей, обобщающих закон Дарси. Для случая изотермических процессов в изотропном пласте эти уравнения имеют вид:
—в ф кав ( — па ф-1 / \ Л , ч
* =-*£ -т Р7-р") ™ (ф = 1,2,...,Ф -1), (2)
а=1 М \ Р /=1 )
где индекс ф соответствует скелету пласта, кав - многофазные коэффициенты ОФП, и ра - скорость фильтрации и истинная плотность а -фазы соответственно. Уравнения (2) служат для определения скоростей фильтрации фаз и выражают их зависимость от градиентов давления во всех фазах и гравитационных сил в фазе. Последнее слагаемое в правой части (2) - сила Рахматуллина - связано с силовым воздействием системы на выделенную фазу из-за несовпадения давления в отдельных фазах (капиллярные эффекты, эффекты прочности, инерции фаз в их мелкомасштабном движении). В соотношениях (2) пренебрегли термодиффузионными отношениями и вкладом термофоретических сил, являющихся эффектом второго порядка малости.
Альтернативная уравнениям (2) форма представления уравнений движения в виде обобщенных соотношений Стефана - Максвелла для гетерогенной среды [3] в рассматриваемом случае имеет вид:
( —в — а Л
(ра-ра^+ т 2 (рв-ра)в = ^ яав ^ - ^ (в* а), (3)
Р в=1 в=1 ^ ^
где Яав - матрица бинарных коэффициентов сопротивления, Яав = Ява (принцип взаимности Онзагера).
Показано [3, 4], что эти коэффициенты связаны с многофазными коэффициентами ОФП кар системой алгебраических уравнений:
ф ( ф Л кав (па
к II ^ +ЪRaRaS\Jh-J = т\£—«а\ (в**) (а,Г = 1,2, .....ф). (4)
р=1 V s=l Р ) № V Р
Представление уравнений движения в виде (3) имеет двойное преимущество.
Во-первых, определение бинарных коэффициентов Яав сопротивления с помощью стандартных экспериментов по совместному движению двух фаз в пористой среде дает принципиальную возможность вычисления многофазных ОФП кар путем решения системы уравнений (4). Отметим, что методы прямого экспериментального определения ОФП отсутствуют даже для трехфазных систем.
Во-вторых, уравнения фильтрации в форме (3), будучи подставленными в уравнения неразрывности, приводят к уравнениям, не разрешенным относительно старшей производной. Это создает определенные трудности в решении задач при численной реализации.
В частном случае двухфазного течения, описывающем процесс вытеснения нефти водой, задача определения материальных функций в уравнениях (2) и (3) упрощается. Матрицы коэффициентов сопротивлений Яав и ОФП кав задаются тремя функциями на-сыщенностей. Более того, для их определения возможен альтернативный подход, основанный на серии известных экспериментов.
3. Обобщенная модель двухфазной фильтрации
Для случая изотермического двухфазного течения обобщенный закон (2) многофазной фильтрации представляется в виде:
= -Л1(УА - Р§) - Л2 - р2§) » $2 = -^21 - Р§) - А22 (Ур2 - р§) ,
к^И0) , ^ =А , ^ =^1. (5)
Ир №2
Здесь пренебрегли силой Рахматуллина, учитывая капиллярные эффекты традиционной зависимостью капиллярного давления в форме Леверетта [1].
Как видно из системы (5), коэффициенты Хар имеют смысл обобщенных подвиж-
ностей фаз. При этом недиагональные коэффициенты матрицы Хар (а*в) обусловлены
вязким взаимодействием между флюидами и капиллярными силами, а диагональные коэффициенты представляют собой вклад обеих фаз в суммарный поток, как если бы каждая из них двигалась независимо в поровом пространстве, видоизмененном присутствием другой фазы.
Покажем, в каких случаях классический и обобщенный формализм двухфазной фильтрации эквивалентны друг другу.
Предполагая несжимаемость обеих фаз и скелета пласта, для одномерного вытеснения нефти водой в изотропном наклонном пласте после стандартных преобразований из системы (5) с учетом уравнения неразрывности получим определяющее уравнение для насыщенности ^ смачивающей фазы, аналогичное по структуре уравнению Раппопорта -Лиса классической модели [1]:
д 5 д _ , ч д 5 . д 5
- = -В ( 5 ) -+ А ( 5 ) -
дт £ ' д £ д£
(6)
Здесь введены безразмерные переменные т = и / тЬ, £ = х/Ь (и - постоянная суммарная
скорость фильтрации фаз, Ь - характерная длина пласта). Для обобщенной модели:
В(5) —ЛОт-^т-К , ар = Р -р.
иЬ А1 + А 2
А(5) = -((-^Ч ^^
А1 + А 2
и
А1 + А 2
Аа =Лаа+Лав, Т* А = АЛ2 - ^ ^ав = К
М^ Мо)
Мв
Мо = М
М2
(7)
Для уравнения Раппопорта - Лиса:
В(5) = р ; / (5 ) =
иЬ
ф)
К )+Мо к2 У
А(5) = -{/ ' + (Дт= (а = 1,2).
и
Ма
В уравнениях (7)-(8) символ «штрих» (') означает дифференцирование по насыщенности (например, р'к = ёрк / ёя, /' = ё/(я)/ ёя и т.д.); Ла и /(я)- соответственно
подвижности фаз и функция Бакли - Леверетта в классической теории [1].
Сопоставление обобщенной и классической моделей позволяет сделать вывод об идентичности соответствующих уравнений для двух независимых процессов: для проти-воточной капиллярно-гравитационной пропитки = 0) и однонаправленного вытеснения без учета капиллярных сил. Это приводит к следующим соотношениям:
Л2 = (Л1 -ЛП)(Л2 -Л22)
2' = А/ Л = А/ пропитка Л ~ /а , 2 ~ /д
вытеснение
Лa =Аa =Лaa+Лaв, (9)
а, в = 1,2; a Ф в .
В формулах (9) значок «штрих» соответствует процессу пропитки. Из соотношений (9) находятся выражения коэффициентов матрицы Лав через стандартные подвижности
Л'а и л для двух указанных процессов вытеснения в виде:
1 ЛЛ+Л
1 = Л1Л
2 = Л Л + Л
= л-
^ Л + 11
Л2 =ЛЛ (1-1)1 -Л1 (10)
(Л+ Л)2
Для определения стандартных подвижностей Ла и Л достаточно проведения двух
независимых экспериментов: по однонаправленному вытеснению и по противоточной пропитке.
Известно большое количество лабораторных экспериментов по определению ОФП для процессов пропитки и дренирования в однотипных образцах. Рассмотрим один из таких экспериментов [5].
Опыты проводились на образце, представленном песчаником Вегеа. Характерные параметры системы пористая среда - флюиды:
т = 0.15, к = 0.03-10 м1, и = 0.95-10-3Па• с, р1 = 1071
кг
м
кг
¡и2 = 9.5-10 Па • с , р2 = 833—- (1 - вода, 2 - нефть).
м
В результате были получены следующие зависимости для ОФП: для однонаправленного вытеснения:
к2(5)=
Об^0837^)15169,0 < 5 < 0.837, 0.837
0,0.837 <5 < 1.
к2(5) =
1.3612(°^7-5)15169,0 < 5 < 0.837,
0.837 0,0.837 <5 < 1.
(11)
для противоточной пропитки: '0,0 < 5 < 0.226,
ку( 5) =
0.2508( 5-0.22б)2,0.22б < 5 < 1. 1-0.226
к2(5):
0 729 - V 1.2034( ), 0 < 5 < 0.729,
0.729
0.729 < 5 < 1.
(12)
На основе формул (11) и (12) находятся подвижности, соответствующие процессам вытеснения Л,а(5) = кка(5)/ ¡а и пропитки А'а (5) = кк'а (5)/ ¡а (а = 1, 2), по которым по формулам (10) определяются диагональные Лаа и «перекрестный» коэффициенты
матрицы обобщенных подвижностей.
Эти результаты и зависимости матрицы обобщенных подвижностей представлены соответственно на рис. 1 и 2.
4. Количественная оценка вклада «перекрестных» фазовых проницаемостей
в показатели вытеснения
4.1. Постановка задачи. В рамках допущений, сформулированных в разд. 3, задача о вытеснении нефти водой в одномерном приближении приводит к определяющему уравнению (6) для водонасыщенности 5:
д5 = д
В(5)
д5 д?
+ А(5)
дя
(13)
записанному в безразмерных переменных т = / тЬ, ? = х/Ь .
Здесь коэффициенты В(^) и А(^) представлены соответственно формулами (7) и (8) для обобщенной модели и уравнения Раппопорта - Лиса.
Начально-краевые условия. Эти условия для обобщенной модели фильтрации могут быть сформулированы аналогично условиям для уравнения Раппопорта-Лиса [1].
В качестве начального условия задаются значения неизвестной функции я в зависимости от пространственных координат. Можно считать, что в начальный момент т = 0 насыщенность всюду постоянна, например:
я(д,0) = я0, где 0 < я0 < я, (14)
(я, — насыщенность связанной воды).
На входе в пласт задаются расход закачиваемой воды и равенство нулю скорости фильтрации нефти. На выходе из пласта (д = 1) можно пренебречь градиентом капиллярного давления по сравнению с градиентом давления в фазах. Если на входе в пласт не учитывать влияния капиллярных и гравитационных сил, а в выходном сечении пренебречь «концевым эффектом», то граничные условия формулируются следующим образом:
я(0,т) = я*: Л22(/) = 0, Л2(я') = 0; к2(я*) = 0;
% = 0 = 0 при д = 1. (15)
дд дд
Начально-краевая задача (13)-(15) решалась методом конечных разностей с соблюдением условий абсолютной устойчивости и сходимости в среднем разностной схемы.
4.2. Результаты расчетов и выводы. При расчетах использовались численные значения параметров и коэффициентов обобщенной матрицы подвижностей, соответствующие эксперименту [5] (см. разд. 4). При решении краевой задачи для уравнения Раппопорта - Лиса использовались функции относительных фазовых проницаемостей и капиллярного давления, соответствующие циклу пропитки (увеличению насыщенности смачивающей фазы). При этом капиллярное давление, выраженное в единицах Па-105, принималось в виде (в соответствии с экспериментом):
Рк = —0,276я + 0,068 / я + 0,051.
Ниже приведены результаты расчетов. На рис. 3 и 4 показаны профили водонасы-щенности на разные моменты времени для классической и обобщенной моделей, а также динамика коэффициента вытеснения п (рис. 5). Последний определяется как отношение вытесненного водой объема нефти к общему объему пор, занятых нефтью до начала вы-
теснения. После прорыва воды через добывающую галерею на графиках показаны безразмерное время прорыва воды т и коэффициент безводной нефтеотдачи п.
Рассмотрены три характерные ситуации:
1) до момента прорыва воды через добывающую галерею для обеих моделей;
2) до момента прорыва воды для обобщенной модели и при наступлении обводненности в классическом случае;
3) при наступлении обводненности в обоих случаях.
Отдельно дается зависимость коэффициента вытеснения от времени. Вертикальные пунктирные линии на нем соответствуют времени прорыва воды.
Графики показывают, что в случае обобщенной модели процесс вытеснения протекает медленнее и момент прорыва фронта воды наступает позднее. Это объясняется учетом поверхностного межвязкостного взаимодействия между фазами в обобщенной модели. На границе раздела сред более вязкая жидкость «тормозит» менее вязкую и, соответственно, создает дополнительное сопротивление в продвижении последней (в сравнении с классической моделью). На рис. 6 и 7 показана зависимость коэффициента вытеснения от времени. Вертикальные пунктирные линии соответствуют времени начала обводнения для каждой из моделей.
По расчетам, для классической модели обводнение начинается раньше (т = 241), чем для обобщенной модели (т = 342). После прорыва воды характер вытеснения нефти остается различным. Но асимптотически, при больших временах и объемах закачки воды, количество вытесненной нефти, рассчитанное по обеим моделям, будет практически одинаковым.
Таким образом, сравнение обобщенного закона фильтрации, учитывающего «перекрестные» фазовые проницаемости, с классической теорией показало различную картину изменения профилей насыщенности. Характер вытеснения в обобщенной модели более медленный, время начала обводнения наступает позже и коэффициент вытеснения ниже.
5. О новом эффекте капиллярных сил на межфазной границе в анизотропных средах
В предыдущих разделах мы ограничились случаем изотропии для выявления эффекта вязкостного взаимодействия между флюидами на межфазной границе. Капиллярное
давление по-прежнему оставалось скалярной функцией насыщенности. Эксперименты, выполненные на анизотропном керне, показали, что это неверно.
Первоначально была проведена серия комплексных экспериментов [6] на отобранном цилиндрическом керне сцементированного слоистого песчаника высотой и диаметром по 100 мм, который был подготовлен для проведения исследований. В соответствии с методикой [6] были определены тип анизотропии (ортотропная симметрия) и главные направления тензора абсолютной проницаемости. Затем из исходного керна были выпилены 4 образца меньших размеров, керны с диаметром 25 мм и длиной 30 мм. Три керна (образцы 1-3) были выпилены вдоль главных направлений; образцы 2 и 3 - вдоль осей X и У в плоскости напластования; образец 1 - в вертикальном направлении 2 (рис. 5). Четвертый образец - контрольный - по биссектрисе угла между экстремальными направлениями в плоскости напластования (рис. 5). Последний образец использовался для проверки предположения о том, что ось симметрии керна совпадает с главным направлением тензора абсолютной проницаемости и что полученные в результате эксперимента значения действительно являются компонентами тензора. Далее на всех ориентированных образцах были определены пористость и главные значения абсолютной проницаемости при фильтрации гелия: т = 18,6%, к1 = 668 мД, к2 = 689 мД, к3 = 579 мД. По полученным значениям к2 и к3
можно вычислить значение проницаемости в любом направлении к (п) = к^пп^, найти теоретическое значение к4 = 638 мД и сравнить его с полученным экспериментальным значением к4Э = 644 мД. Отличие составляет менее 2%.
Последующие исследования ориентированных образцов керна на компьютерном томографе БкуБеап 1172 привели к установлению тензорной природы функций плотности распределения пор по радиусам и введению тензора Щ характерных линейных размеров (эффективных радиусов пор).
После завершения этой трудоемкой работы была проведена дополнительная серия исследований, состоящая в определении остаточной водонасыщенности для тех же ориентированных образцов. Для пересчета лабораторных параметров была использована методика Хасслера - Брунера [7], позволяющая с высокой степенью точности определить величину капиллярного давления на внешнем торце образца, оценить соответствующую ему насыщенность и построить зависимости капиллярного давления от насыщенности. Результаты исследований приведены на рис. 6 для образцов 2 и 3 вдоль направлений X и У в плоскости напластования и для контрольного образца 4 и на рис. 7 для образца 1 вдоль
вертикального направления Z. Результаты сравнения теоретических и экспериментальных значений капиллярного давления на контрольном образце позволяют предположить, что капиллярное давление в анизотропных средах имеет тензорную природу.
Важным результатом этого исследования является установление эффекта зависимости капиллярного давления от направления измерений. Отсюда следует, что капиллярное давление не является универсальной функцией насыщенности для анизотропного образца породы, но зависит от направления течения, от направления воздействия на пласт. Подтверждение этого факта требует более глубокого экспериментального и теоретического исследования, выявления физических механизмов этого эффекта, особенностей физического и физико-химического взаимодействий системы скелет пласта - флюид, характера смачиваемости пласта и других факторов. Необходима повторяемость экспериментов для различных образцов пород (терригенных, карбонатных) с разным масштабом неоднородности.
Строго говоря, тензором являются поверхностные напряжения капиллярных сил на межфазной границе рС, связь которых с тензором г], обратным тензору характерных линейных размеров Ri¡■, может быть представлена в виде:
] пгП] )= а^уЫгк1, (16)
Р^As, Пп}■) = а^фЛи
где ас определяется коэффициентом межфазного натяжения и углом смачивания; пр; — диада; Ji]■kl — тензор 4-го ранга, симметричный по паре индексов и их перестановке; тензор г^ находится экспериментально [6].
Структура тензора Jykl в (16) требует дополнительных исследований для различного состава пород и свойств флюидов. При обработке данного эксперимента использовалась следующая аппроксимация этого тензора вдоль главных направлений тензора к:
J: =
а +
¡Р—,](—<)
( \ (г 3Г
\ -3 ^ У
1 — я 1 — я *
где ¡1(5)— первый инвариант тензора г], параметры ai, si определяются экспериментально, при этом а1=J:(s:).
Заметим, что интеграл от функции распределения капиллярного давления является наиболее информативным индикатором смачиваемости, которая изменяется в процессе фильтрации флюидов и движения межфазной границы. Поэтому зависимость капиллярного давления от направления течения имеет определяющую физическую основу.
Заключение
С использованием методов неравновесной термодинамики получены и исследованы уравнения движения многофазной фильтрации (обобщенный закон Дарси) и законы межфазного взаимодействия. Показана возможность прямого определения коэффициентов «многофазных» ОФП. Для двухфазной фильтрации исследован эффект вязкого поверхностного взаимодействия между флюидами, предложен метод определения «перекрестных» фазовых проницаемостей и их количественный вклад в показатели разработки месторождений.
На основе современных методов исследования анизотропного кернового материала установлен новый эффект зависимости капиллярного давления от направления измерений. Это означает, что капиллярное давление, так же как и ОФП, не является универсальной скалярной функцией насыщенности. Физическое объяснение этого факта может быть связано с изменением смачиваемости породы при извлечении нефти и газа. Известно также возникновение «гетерогенной» смачиваемости, когда на одних участках залежи поверхность породы гидрофильна, а на других - гидрофобна. А наиболее информативным интегральным показателем смачиваемости является интеграл от функции распределения капиллярного давления.
ЛИТЕРАТУРА
1. Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидромеханика. - М.: Недра, 1993. - 414 с.
2. Максимов В. М. Основы гидротермодинамики пластовых систем. - М.: Недра, 1994. -202 с.
3. Колесниченко А.В., Максимов В.М. Обобщенный закон фильтрации Дарси как следствие соотношений Стефана - Максвелла для гетерогенной среды: Препринт № 45 ИПМ РАН им. М.В. Келдыша. - М., 1991. - 31 с.
4. Колесниченко А.В., Максимов В.М. Термодинамика многофазной термически активной смеси. Законы фильтрации Дарси и диффузии: Препринт № 52 ИПМ РАН им. М.В. Келдыша. - М., 1997. - 32 с.
5. Jerauld G.R., Salter S.J. The effect of pore - structure on hysteresis in relative permeability and capillary pressure: Pore-level modeling // Trans. Porous Media. - 1990. - № 5. -P. 103-151.
6. Дмитриев Н.М., Дмитриев М.Н., Максимов В.М., Кузьмичев А.Н. Двухфазная фильтрация в ортотропной пористой среде: Эксперимент и теория // Изв. РАН. Механика жидкости и газа. - 2014. - № 6. - С. 94-100.
7. Kuznetsov A.M., Baishev A.V., Kuznetsov V.V. Determination of the initial water saturation and capillary curve by centrifugal method // Oil Industry. - 2010. - № 1. - P. 49-51.
ПРИЛОЖЕНИЕ
с.е \ ---однонаправленное вытеснение, для воды а ----- *--однонаправленное вытеснение, для нефти % - противоточная пропитка, для воды \ —*--й— противоточная пропитка, для нефти
£>
¡у?
и 0.2 0-3 0.1 0.5 О Б и.' О.'в 0.9 Водонасыщенность
Рис. 1. Зависимость фазовых проницаемостей от насыщенности [5]
Водонасыщенность
Рис. 2. Зависимость матрицы подвижностей от насыщенности в единицах (Па-с)"1
s Насыщенность
- с lime: 200 Обобщенный закон Q: 0.120528
0.4586 . ч _ ч ч \ \ \ ч \ \ \ \ \ \ _ Классический закон Q: 0.148506
50 = 0.3520 _
1 g
Рис. 3. Распределение насыщенности (т = 2GG)
N ос шжц тм ноеть
S
.SV=0.5652- tima: Э50 ÖfefeäHpvsotM ганан Q: 8,178694
П.4536 . ~ ~ - - .__ t: 342
Qt: 0.176111
ч \ X ч _ Классический »кон
- Q: С.211697
\ ч t: 241
Ql: 0.169022
50=0.3520.
i £
Рис. 4. Распределение насыщенности (т = 35G)
Рис. 5. Расположение вырезанных образцов керна
3 5 <3 и к Е
В 2 3 к ^ 1 и с й ы. 1
0 ■
0 0.02 ШМ Ш 0.03 Водонасыщенность: З^у -—» Кривая капиллярного давления для образца 2 - Кривая капиллярного давления для образца 3 - Кривая капиллярного давления для контрольного образца • • Экспериментальные точки для образца 2 А А Экспериментальные точки для образца 3 ■ ■ Экспериментальные точки для контрольного образца 4
Рис. 6. Зависимости капиллярного давления от водонасыщенности для образцов 2, 3 вдоль направлений X, Y и контрольного образца 4
Е
£ £
— и
—I
И
и
о
к
я:
с ®
И
3 < >
с
\ ----
005 01 015
Водонасыщенностъ: Л'и.
Кривая капиллярного давления для образца 1 • • Экспериментальные точки для образца 1
Рис. 7. Зависимость капиллярного давления от водонасыщенности
для направления 2