Метод последовательной смены стационарных состоянийдля плоскоодномерной задачи фильтрации с предельным градиентом давления Текст научной статьи по специальности «Физика»

Вестник СамГУ — Естественнонаучная серия. 2014■ № 7(118)

МЕХАНИКА

УДК 532.546

МЕТОД ПОСЛЕДОВАТЕЛЬНОЙ СМЕНЫ СТАЦИОНАРНЫХ СОСТОЯНИЙ ДЛЯ ПЛОСКООДНОМЕРНОЙ ЗАДАЧИ ФИЛЬТРАЦИИ С ПРЕДЕЛЬНЫМ ГРАДИЕНТОМ ДАВЛЕНИЯ

© 2014 О.В. Белова,1 В.Ш. Шагапов2

С учетом нелинейных эффектов, наблюдаемых в опытах с низкопроницаемыми пластами, при малых градиентах давления (например, порядка 105 Па/м), предложено уточнение закона Дарси. На основе такой модели методом последовательной смены стационарных состояний и численно решена задача о плоскоодномерной фильтрации. Установлено, что приближенные решения, полученные методом последовательной смены стационарных состояний, для описания распределения давления в пласте и дебита скважины согласуются с численным решением уравнения фильтрации в полной постановке. Проведен анализ влияния предельного градиента давления д и показателя степени, определяющего темп выхода нелинейного закона фильтрации к линейному закону Дарси с ростом градиента давления 7, на особенности гидродинамических полей и дебит скважины.

Ключевые слова: фильтрация, сверхнизкая проницаемость, предельный градиент давления, метод последовательной смены стационарных состояний.

Введение

В последнее время в связи с сокращением запасов углеводородного сырья возросла проблема разработки пластов со сверхнизкой проницаемостью к = 10- 16 + 10- 14 м2. Результаты экспериментов, опубликованные в источниках [1; 2] показывают, что при фильтрации флюидов в низкопроницаемых коллекторах при низких градиентах давления появляются отклонения от линейного закона Дарси.

В работе [3] рассмотрено три режима течения — это фильтрация воды на водо-насыщенном керне, фильтрация нефти в керне с остаточной водонасыщенностью, фильтрация воды в керне с остаточной нефтенасыщенностью. Из анализа экспериментов следует, что нелинейное отклонение происходит не только для нефтей, но

хБелова Ольга Владимировна ([email protected]), кафедра высшей и прикладной математики, Бирский филиал Башкирского государственного университета, 452453, Российская Федерация, г. Бирск, ул. Интернациональная, 10.

2Шагапов Владислав Шайхулагзамович ([email protected]), лаборатория "Моделирование

технологических процессов", Институт механики Уфимского научного центра РАН, 450025, Рос-

сийская Федерация, г. Уфа, пр. Октября, 71.

также и при фильтрации воды. Причем величины характерных градиентов давления, при которых наблюдается нелинейность, зависят также от флюида, каким пористая среда насыщена в исходном состоянии.

Результаты измерений величин скорости фильтрации |и| от градиента давления |Ур| лабораторных исследований [4] представлены в виде зависимости коэффициента эффективной проницаемости от градиента давления. По результатам экспериментов видно, что полученное отклонение от линейного закона Дарси при фильтрации в низкопроницаемых коллекторах характеризуется наличием предельного градиента давления. Видно, что фильтрационное течение вообще отсутствует, если градиент давления становится ниже этого предельного значения. В соответствии с вышесказанным примем, что течение через пористую среду вообще отсутствует, если величина градиента давления становится ниже некоторой предельной величины д. Таким образом, фильтрация жидкости происходит только в тех зонах, где градиент давления выше этой предельной величины. При достаточно высоких значениях градиента давления закон фильтрации становится линейным.

Нарушение линейного закона при низких градиентах давления, связанное с реологией нефтей, описывалось ранее в [5]. Закон фильтрации с предельным градиентом применительно к вязкопластичным жидкостям был предложен А.Х. Мир-заджанзаде [6]. Далее теория течения в пластах неньютоновских жидкостей с такими реологическими свойствами развивалась в его школе [7-8], а также в работах Г.И. Баренблатта с его учениками [9].

1. Модель фильтрации с предельным градиентом давления

Для учета нелинейных эффектов обобщенный закон Дарси запишем в виде

к др / дх

где п — безразмерный коэффициент, учитывающий фактор нелинейности. Фактор нелинейности примем как

п = 1 - д

др у др

дх при дх

> д,

(1.1)

(1.2)

П = 0 при

др дх

< д.

Здесь д — предельный градиент, 7 — показатель степени, определяющий темп выхода закона фильтрации к линейному закону Дарси с ростом величины градиента давления. Закон фильтрации, предложенный А.Х. Мирзаджанзаде, соответствует случаю 7 = 1 . Также для 7 рассмотрены значения 7 = 1/2 и 7 = 2 , характеризующие низкий и высокий темпы выхода к линейному закону фильтрации. На рис. 1.1 представлены характерные зависимости безразмерного коэффициента п от градиента давления для различных значений 7 . Пунктирные, штрихпунктирные и точечные линии соответствуют значениям 7 =1/2, 7 =1 и 7 = 2 соответственно. Для этих значений можно построить некоторые приближенные аналитические решения.

Л

•' *** *

: / .' /' , ' ' - - ~ ~

■ /' •'/' / ;; / у /

• / :.'/ :"

V \

1 ■ ■ ■ ■ ■ ■

V , МПа/м

Рис. 1.1. Зависимости безразмерного коэффициента п от градиента давления

для различных значений 7

Сжимаемость флюида будем учитывать в линейном приближении и тогда для зависимости плотности от давления будем иметь:

Р = Ро(1 + а(Р - Ро))• (О)

Из уравнения неразрывности, записанного в пористой среде, принимая закон фильтрации в виде (1.1) с учетом (1.2) и уравнения состояния (1.3), получим основные уравнения для упругого режима фильтрации

др д дЬ X дх

др| 1 -

дх

др дЬ

— = 0 при

V др дЬ

др дх

<Я\Х

при

к

ш^а

др дЬ

> Я,

(1.4)

Здесь ц и а — коэффициенты динамической вязкости и сжимаемости жидкости, ш — пористость скелета, х — коэффициент пьезопроводности.

Отметим, что (1.4) получено без учета сжимаемости скелета пористой среды. В правой части (1.4) также пренебрегались слагаемые, связанные со сжимаемостью жидкости.

При теоретическом описании фильтрационных течений уравнением вида (1.4) в области фильтрации возникает неизвестная граница х = I (Ь), разделяющая зоны течения и покоя (где и = 0 ). Причем на этой границе х = I (Ь), помимо непрерывности давления, должно выполняться второе условие

др

дх

= Я.

(1.5)

Отметим, что при (1.5) автоматически также выполняется условие, следующее из закона сохранения масс.

1

2. Установившаяся фильтрация

В этом случае (dp/dt = 0) из (1.4) следует

4 =с, (2.1)

|l i

dp dx

где c — произвольная постоянная (причем она может быть функцией времени, если распределение давления постоянное, dp/dx = const ).

Пусть фильтрация происходит через пористую среду толщиной l при постоянном перепаде давления. Тогда для давления можем записать следующие граничные условия:

p = pe при x = 0,p = po при x = l. (2.2)

Решение уравнения (2.1) с учетом (2.2) имеет вид:

p = (p0 — pe) J + pe. (23)

Видно, что для установившегося режима распределение давления точно такое же, как и в случае линейного закона Дарси. Что касается величины скорости фильтрации (или расхода жидкости через единичную площадь), то для нее, в соответствии с (1.1), будем иметь:

- k AI ( 1 -

Л l

ql

Ap

(Ap = po - pe) при |Ap| > ql, 0 |Ap| < ql. (2.4)

Y

u

u

3. Решение методом последовательной смены стационарных состояний

Решение задачи для неустановившегося режима фильтрации будем строить методом последовательной смены стационарных состояний (ПССС). Отметим, что метод ПССС является частным случаем метода интегральных соотношений [9; 10], в котором при построении приближенного решения ограничиваются лишь моментом нулевого порядка.

Пусть в исходном состоянии Ь < 0 в полубесконечной пористой среде (х > 0), насыщенной жидкостью, имеет место однородное распределение давления р = ро, и в момент времени Ь = 0 начинается отбор жидкости через границу х = 0 . Полагаем, что распределение давления в зоне, где происходит фильтрационное течение, описывается решением вида (2.3), причем координата границы возмущенной области х = 1(Ь) является функцией времени. Полагая, что скорость снижения упругого запаса жидкости в возмущенной области пласта равна массовому расходу через границу х = 0, можем записать следующее уравнение баланса массы:

i

d

ma— (po — p)dx = ue, (3.1)

o

где ие — скорость фильтрации на границе.

Рассмотрим режим фильтрации, когда на границе х = 0 поддерживается постоянное значение давления ре . Принимая выражение (2.3) для распределения

давления в пористой среде на основе (3.1) с учетом (1.1) при х = 0, можем по-

лучить:

Ш2 АН

л" = 4х( 1"

ч1

Ар

(3.2)

Отсюда получим закон движения границы возмущенной области в неявном виде:

2xt =

М1

1 - №/|Ар|)

7 '

Закон движения границы Щ) имеет следующий неявный вид:

1 3 2 / -

х=-2 АО2<2--2(Ар)211 -ы- г

(3.3)

л =1/2,

2х^ = - АР1 - (АР я \ я

ъч -< £

Ар / ( 4хдН

1 = Т I - ехЧ- ~Ар2~

1/2

,7 = 2.

(3.4)

(3.5)

(3.6)

Из (3.3) видно, что существует верхняя граница возмущенной области, определяемая как IФ = Ар/я , при достижении которой процесс фильтрации прекращается. В дальнейшем эту величину будем называть предельным расстоянием зоны влияния скважины. Характерное время фильтрации до остановки границы возмущенной области составляет величину порядка £ = Ар2/4хя2 .

Для скорости фильтрации ие на границе х = 0 при 7 = 2 согласно (2.4) следует:

пе = -^ А ехр (-4хя^/Ар2) . (3.7)

4. Сравнение аналитических решений с численными

Для обоснования применимости метода ПССС был проведен сравнительный анализ приближенных решений (3.4)—(3.6) и (3.7) с численным решением уравнений (1.4). Чтобы получить численное решение по модели с предельным градиентом, проведена следующая регуляризация постановки задачи. Принято, что при низких градиентах давления закон фильтрации имеет линейный вид

ек др ¡л дх'

(4.1)

где е — малый безразмерный параметр (е << 1). При некотором значении градиента давления (др/дх)ф величина скорости фильтрации согласно нелинейному закону (1.1) (с учетом (1.2)) и линейному (4.1) равны. Из условия такого "сращивания" при 7 = 2 получим уравнение

к(ддр\ и -

¡ дх

я,

др дх

ек др ¡ дх

Отсюда будем иметь

др дх

v/Г—

(4.2)

(4.3)

7

0

1

7

*

*

я

*

В области фильтрации, где \др/дх\ < д/\/1 — £, уравнение фильтрации примем как

др д 2р

дЬ X дх2

(4.4)

При высоких градиентах давления (\др/дх\ ^ д/л/Г—£) уравнение фильтрации примем в виде (1.4). Решение такой регуляризированной задачи при £ ^ 0 соответствует решению по модели с предельным градиентом.

Рассмотрим нестационарную задачу о фильтрации жидкости через пористый образец длины хо. Пусть в исходном состоянии давление в пористой среде р = = ро (Ь ^ 0, 0 < х ^ хо). Далее, в момент времени Ь = 0, давление на левой границе (х = 0) мгновенно изменяется до величины ре и поддерживается постоянным, а на правой границе (х = хо) поддерживается начальное значение р = ро.

Для удобства сравнения результатов расчетов введем безразмерные переменные и параметры:

X = х/хо, т = хЬхо, Р = (ро — р)/Ар, ие = пе^хо/(кАр), Q = дхо/Др, где ие — величина скорости фильтрации на границе х = 0.

Для численных значений параметров, отвечающих за предельный градиент и регуляризацию, приняты следующие величины: £ = 10~2, Q = 1.

На рис. 4.1, а представлена эволюция по времени распределений давления для 7 = 2. Сплошные линии соответствуют решениям полученным методом ПССС. Пунктирные линии соответствуют численным решениям регуляризированной задачи. Числа на линиях соответствуют безразмерным моментам времени т. На рис. 4.1, б представлены безразмерные расходы от времени т для 7 = 2. Видно, что в плане определения расходов жидкости расхождение между кривыми, полученными методом ПССС и численным решением по общему уравнению фильтрации с предельным градиентом, не превосходит десяти процентов. Причем это различие с приближением характерной протяженности I возмущенной области к предельному значению I* несколько повышается. Аналогичное сходство для распределений давлений и расходов имеет место также при 7 =1/2 и 7 =1.

и

0,0 0,2 0,4 0,6 0,8 1,0 X

а б

Рис. 4.1. Сравнение приближенных решений, полученных методом ПССС с численным решением нелинейного уравнения фильтрации

Таким образом, приближенные решения, полученные методом ПССС в плане описания поля течения, а также расхода жидкости через границу пористой фа-

зы, неплохо согласуются с численным решением уравнения фильтрации. Все это, в свою очередь, дает основание на возможность использования метода ПССС для построения приближенных решений.

5. Влияние нелинейности

В случае линейного закона фильтрации Дарси (я ^ 0) формулы (3.6) и (3.7) переходят в известные [10] законы движения границы возмущенной области и расхода жидкости через единичную площадь:

I = и ие = -кЛр. (5.1)

М I

Введем коэффициент п(1) = ие/ие, определяющий безразмерный расход жидкости для принятого нелинейного закона фильтрации, как отношение к расходу жидкости, когда справедлив закон Дарси. На основе (3.7) и (5.1) с учетом (3.6) для п(1) в случае 7 = 2 можем получить

(1) _ 2я ехр (-4Хд2уЛр2) ^Е

п — "л ТТт". (5.2)

ЛР (1 - ехр {-4хч2Е/ЛР2)у'

На рис. 5.1 приведена зависимость безразмерного расхода п(1) от времени при Я = 0,1 МПа/м. Для параметров, определяющих физические свойства пористой среды и жидкости, приняты следующие величины: т = 0,1; к = 10_ 15 м2; р = = 860 кг/м3; м = 1, 6• 10~3 Пас. Числа на кривых соответствуют величине перепада давления в МПа. Из численных расчетов следует, что с увеличением перепада давления увеличивается расход жидкости. Из графиков также видно, что для любого значения перепада давления Лр существует некоторое значение времени Е(т), при котором расход жидкости для разных значений 7 имеет различие не более чем на четыре процента. Причем при Е < Е(т) с увеличением значения 7 расход увеличивается, а при Е < Е(т) с увеличением значения 7 расход уменьшается.

Можно рассмотреть эту задачу, когда фильтрация происходит в режиме постоянного объемного расхода жидкости через границу пористой среды х = 0. Тогда скорость фильтрации на этой границе постоянна и пусть равна ие. Тогда, согласно (2.4), полагая, что для распределения давления имеет место (2.3) при 7 = 2, можем записать уравнение:

к Лр ( я212 '

V Лр2

С другой стороны, из уравнения баланса массы (3.1) с учетом (2.3) можем записать:

"^(1 -^2) = -«е. (5.3)

I

та— Лр (1--) —х = ие. (5.4)

—Ь ] V и

о

Отсюда будем иметь

Лр1 = -(5.5) ат

Исключая Лр из (5.3) и (5.5), получим уравнение для определения закона движения границы возмущенной зоны. Разрешая такое уравнение, можем получить:

I = л/хУ^Ь (/1 + 4/Ув - 1)1/2 ,Ув = (мие/кд)2 . (5.6)

Г 1,0

0,6

0,4-

0,2

\\ ■ • У4'. ■ ■

и \ V • п \ V. - ' Ч'ч '■ ч •

\ V " \ \ ч ч " N • " ч • .15

4 ч\. "

10 1 1

100

200

300

400

500

600 t, Ч

Рис. 5.1. Динамика безразмерного коэффициента расхода п(1) от времени. Числа на кривых соответстуют величине перепада давления в МПа

Подставляя это выражение в (5.4), получим решение для закона роста перепада давления:

\ 1/2

(5.7)

/-

(лД+4/Уз - 1) Для случая д ^ 0 из (5.6) и (5.7) следует

I =у/2&, Ар = -

(5.8)

Введем безразмерные протяженность возмущенной области и перепад давления как

ь = 1, АР = Ар.

I Ар

(5.9)

Используя решения (5.5)—(5.8), для этих введенных безразмерных величин получим

Ь =((^1+4/Уз - Уз/2^ 1/2 , АР =1/1. (5.10)

Заключение

Предложена теоретическая модель, описывающая процессы фильтрации в низкопроницаемых пластах с проявлением предельного градиента давления и степенного закона выхода к закону Дарси при повышении величины градиента давления, в плоскоодномерной постановке. Методом последовательной смены стационарных состояний и численно решена задача о притоке жидкости к скважине. Выявлено, что приближенные решения, полученные методом ПССС, в плане описания полей давления и дебита хорошо согласуются с численным решением уравнения фильтрации в полной постановке. Установлено, что картина гидродинамических

полей, а также дебит скважины характеризуются не только величиной предельного градиента давления, но и показателем степени, определяющим темп выхода нелинейного закона фильтрации к линейному закону Дарси с ростом градиента давления.

Литература

[1] Behavior of Flow trough Low-Permeability Reservoirs / Qun Lei [et al.] // SPE 113144. 2008. P. 1-7.

[2] Отклонения от закона Дарси при фильтрации в низкопроницаемых пористых средах / В.А. Байков [и др.] // Нефтепромысловая химия: тез. докл. VI Всероссийской научно-практ. конф. (г. Москва, 23-24 июня. 2011 г.). М.: Изд-во РГУ нефти и газа им. И.М. Губкина, 2011. С. 37-40.

[3] Нелинейная фильтрация в низкопроницаемых коллекторах. Анализ и интерпретация результатов лабораторных исследований керна Приобского месторождения / В.А. Байков [и др.] // Вестник ОАО НК "Роснефть". 2013. Т. 31. № 2. С. 8-12.

[4] Нелинейная фильтрация в низкопроницаемых коллекторах. Лабораторные фильтрационные исследования керна Приобского месторождения / В.А. Байков [и др.] // Вестник ОАО НК "Роснефть". 2013. Т. 31. № 2. С. 4-7.

[5] Коробов К.Я., Антилин Ю.В. О нарушениях линейного закона фильтрации при низких градиентах давления // Нефтяное хозяйство. 1968. № 8. С. 26-28.

[6] Мирзаджанзаде А.Х. Вопросы гидродинамики вязкопластичных и вязких жидкостей в применении к нефтедобыче. Баку: Азернефтнешр, 1959. 409 с.

[7] Алишаев М.Г., Розенберг М.Д., Теслюк Е.В. Неизотермическая фильтрация при разработке нефтяных месторождений. М.: Недра, 1985. 271 с.

[8] Мирзаджанзаде А.Х., Аметов И.М. Технология и техника добычи нефти. М.: Недра, 1986. 409 с.

[9] Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в пористых пластах. М.: Недра, 1984. 211 с.

[10] Подземная гидравлика: учебник для вузов / К.С. Басниев [и др.]. М.: Недра, 1986. 303 с.

References

[1] Behavior of Flow trough Low-Permeability Reservoirs. Q. Lei [et al.]. SPE 113144, 2008, pp. 1-7.

[2] Deviation from Darcy law at filtering in low permeability porous media. V.A. Bikov [et al.]. Neftepromyslovaia khimiia: tez. dokl. VI Vserossiiskoi nauchno-prakt. konf. g. Moskva, 23-24 iiunia. M., Izd-vo RGU nefti i gaza im. I.M. Gubkina, 2011, pp. 37-40 (in Russian)

[3] Nonlinear filtering in low permeability reservoirs. Analysis and interpretation of results of laboratory researches of a core of Priobsky field. V.A. Bikov [et al.]. Vestnik OAO NK "Rosneft" [Vestnik of OJSC Oil Company "Rosneft"], 2013, no. 2, pp. 8-12 (in Russian)

[4] Nonlinear filtering in low permeability reservoirs. Laboratory filtrational researches of a core of Priobsky field. V.A. Bikov[et al.]. Vestnik OAO NK "Rosneft" [Vestnik of OJSC Oil Company "Rosneft"], 2013, no. 2, pp. 4-7 (in Russian)

[5] Korobov K.J., Antilin Yu.V. On violations of linear filtration law at low pressure gradients. Neftianoe hoziaistvo [Oil Industry], 1968, no. 8, pp. 26-28 (in Russian)

[6] Mirzadzhanzade A.Kh. Viscoplastic and viscous fluids problems in hydrodynamics applied to oil production. Baku: Azerneftnesh, 1959, 409 p. (in Russian)

[7] Alishaev M.G., Rosenberg M.D., Teslyuk E.V. Nonisothermal Filtration in the development of oil fields. M., Nedra, 1985, 271 p. (in Russian)

[8] Mirzadzhanzade A.Kh., Ametov I.M. Technics and technology of oil production. M., Nedra, 1986, 386 p. (in Russian)

[9] Barenblatt G.I., Entov V.M., Ryzhik V.M. Movement of liquids and gases in porous formations. M., Nedra, 1984, 211 p. (in Russian)

[10] Underground hydraulics: college textbook. K.S. Basniev [et al.]. M., Nedra, 1986, 303 p. (in Russian)

Поступила в редакцию 11//V/2014;

в окончательном варианте — 11/IV/2014.

METHOD OF SEQUENTIAL CHANGING OF STATIONARY STATES FOR ONE-DIMENSIONAL FILTRATION PROBLEM WITH LIMITING GRADIENT

OF PRESSURE

© 2014 O.V. Belova? V.Sh. Shagapov4

Taking into account nonlinear effects observed in experiments with low-permeability layers, at low pressure gradients (e.g., about 105 Pa/m), refinement of Darcy law is proposed. On the basis of this model, by means of method of sequential change of stationary states and the problem of one-dimensional filtering is numerically solved. It is established that approximate solutions received by the method of sequential change of stationary states, for the description of distribution of pressure in layer and a well production, will be agreed with the numerical solution of the equation of a filtration in full statement. The analysis of influence of pressure gradient q and limiting exponent defining the rate of yield of the nonlinear filtration law to the linear Darcy's law with increasing pressure gradient y, on the features of hydrodynamic fields and well production is carried out.

Key words: filtration, ultra-low permeability, limiting pressure gradient, method

of sequential change of stationary states.

Paper received 11//V/2014. Paper accepted 11/IV/2014.

3Belova Olga Vladimirovna ([email protected]), the Dept. of Higher and Applied Mathematics, Birsk branch of Bashkir State University, Birsk, 452453, Russian Federation.

4Shagapov Vladislav Shayhulagzamovich ([email protected]), laboratory "modeling of technological processes", Institute of Mechanics, Ufa Scientific Center, Russian Academy of Sciences,

Ufa, 450025, Russian Federation.