Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Нефтегазовая гидромеханика

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
4.17 Mб
Скачать

В результате решения данного уравнения получена формула, названная основным уравнением упругого режима:

P(r,t ) = P

qμ

E

 

r2

,

(2.342)

 

 

 

 

 

 

 

 

0

 

 

i

 

 

 

 

 

 

4πkh

 

 

 

4χt

 

 

где Ei (x) – интегральная

экспоненциальная функция

(данная

функция табулирована); r – расстояние от скважины до точки, в которой определяется давление Р.

Формула (2.342) имеет широкое практическое применение, ее используют для описания неустановившейся фильтрации упругой (сжимаемой) капельной жидкости при работе скважин с постоянными дебитами. Коме того, важнейшей областью использования данного уравнения следует считать гидродинамические исследования скважин.

При малых значениях аргумента

r2

< 1 для интегральной

 

 

4χt

экспоненциальной функции можно записать

 

E

(x) = ln 1

0,5772.

 

(2.343)

i

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом (2.343) основное уравнение упругого режима записы-

вают в виде

 

 

 

 

 

 

 

 

P(r, t ) = P

qμ

ln

2,246χt .

(2.344)

 

 

0

 

 

 

r

2

 

 

 

 

 

4πkh

 

 

 

Приближенные решения

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

Метод последовательной смены стационарных состояний

(ПССС). В соответствии с разработанным И.А. Чарным методом ПССС, в каждый момент времени вся область движения условно

161

разбивается на две зоны – возмущенную и невозмущенную. В возмущенной зоне, начинающейся от стенки скважины (галереи), давление распределяется по законам стационарного движения, а внешняя граница данной области служит на этот момент условным контуром питания. В невозмущенной области давление всюду постоянно и равно Рк. Закон перемещения подвижной границы раздела (l(t) для одномерного движения и R(t) – для плоскорадиального) определяется конкретно для каждого случая.

Одномерный фильтрационный поток Случай 1. Мгновенный пуск галереи с постоянным дебитом. За-

кон перемещения зоны возмущения:

l (t) = 2χt.

(2.345)

Закон распределения давления:

P = Pк

Qμн

( 2χt x).

(2.346)

 

 

kah

 

Случай 2. Мгновенный пуск галереи с постоянным давлением. Закон движения границы возмущенной области:

l (t) = 2 χt.

 

 

(2.347)

Распределение давления в возмущенной зоне пласта

 

 

 

x

 

 

P = Pк (Pк Pг ) 1

 

.

(2.348)

 

 

 

 

 

 

 

 

2 χt

 

Дебит галереи

 

 

 

 

 

 

Q =

k

 

Pк

Pг

ah.

(2.349)

 

2

 

 

μн

χt

 

Плоскорадиальный фильтрационный поток Случай 1. Мгновенный пуск галереи с постоянным дебитом. За-

кон перемещения зоны возмущения:

R(t ) = r2

+ 4χt.

(2.350)

c

 

 

162

Закон распределения давления:

 

Qμ

н

 

r2

+ 4χt

 

 

P = P

 

ln

c

 

.

(2.351)

2πkh

 

r

к

 

 

 

 

Метод Пирвердяна. Этот метод аналогичен методу ПССС и уточняет его. Поток также условно разбивается на две зоны – возмущенную и невозмущенную. Но, в отличие от метода ПССС, распределение давления в пределах возмущенной области задается другим образом. В невозмущенной зоне давление также постоянно

и равно Рк.

Случай 1. Мгновенный пуск галереи с постоянным дебитом. Закон перемещения зоны возмущения:

 

l (t) = 6χt.

 

 

 

(2.352)

Закон распределения давления:

 

 

 

 

 

P = Pк

Qμ

н

6χt

 

 

x

2

 

 

1

 

.

(2.353)

 

 

 

 

2kah

 

 

 

 

 

 

 

 

6χt

 

Случай 2. Мгновенный пуск галереи с постоянным давлением Рг. Закон перемещения зоны возмущения:

l (t) = 12χt.

 

 

(2.354)

Распределение давления

 

 

 

 

 

 

 

 

x

2

 

P = Pк (Pк Pг ) 1

 

 

.

(2.355)

 

 

 

 

 

 

 

 

 

 

 

12χt

 

Дебит галереи

 

 

 

 

Q =

2kah

Pк Pг

.

(2.356)

 

 

μн

12χt

 

Метод интегральных соотношений. Предложен Г.И. Баренб-

латтом, основан на разделении пласта на возмущенную и невозму-

163

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

В соответствии с методом, в невозмущенной части пласта движение отсутствует, давление всюду постоянно и равно начальному. В возмущенной части пласта распределение давления представляется в виде

P = A ln

r

+ B + B

r

+ B

 

r 2

+ … + B

 

r n

,

 

 

 

 

 

 

 

 

 

R

0 1

R

2

 

n

 

 

 

 

 

 

R

 

 

R

 

где R – радиус границы возмущенной зоны; A, B0 , B1 ,…, Bn – коэф-

фициенты многочлена с коэффициентами, зависящие от времени. Задача сводится к получению коэффициентов многочлена, ко-

торые должны удовлетворять граничным условиям на стенке скважины и на границе зоны возмущения, а также особым интегральным соотношениям, полученным Г.И. Баренблаттом. Показатель степени n определяет точность полученного решения. Чем выше требуется точность, тем больше должно быть значение n.

Уравнение пьезопроводности методом интегральных соотношений, полученное для радиального движения жидкости к скважине с постоянным дебитом. При этом давление в любой точке пласта (в пределах возмущенной зоны) в любой момент времени определяется в соответствии с уравнением:

 

Qμн

 

 

2

+

12χt

 

 

 

r

 

 

 

P(r, t ) = Pк

 

 

rc

1

+

 

 

 

(2.357)

 

ln

 

 

 

 

 

 

 

,

2πkh

 

 

r

 

r2

+ 12χt

 

 

 

 

 

 

 

 

 

c

 

 

 

 

а закон перемещения зоны возмущения имеет вид

 

 

 

 

 

R(t ) =

r2

+ 12χt.

 

 

 

 

 

(2.358)

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

164

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.14.4. Неустановившаяся фильтрация жидкости при работе скважин с переменным дебитом

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

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

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

Изменение давления в любой точке пласта, вызванное работой скважины с переменным дебитом, определится по формуле

n

P = Pi ,

i=1

где n – количество фиктивных скважин; Pi – изменение давления,

вызванное работой i-й фиктивной скважины, заменяющей работу реальной скважины при изменении ее дебита с qi–1 до qi:

P =

μ

(q

q

)

 

E

 

r2

,

 

 

 

 

 

i

4πkh

i

i1

 

i

 

4χti

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ti – время работы i-й фиктивной скважины,

i=1

ti = T ti ,

1

где Т – полное время работы реальной скважины.

165

Рассмотрим данную тему на следующем примере.

Пусть в однородном изотропном пласте постоянной толщины h работает скважина – точечный сток. Работа скважины характеризуется переменным во времени дебитом (рис. 2.57). Определить давление в точке пласта на расстоянии r от скважины по истечении времени t1 + t2 + t3.

Рис. 2.57. График работы скважины с переменным дебитом

Данная скважина заменяется на три фиктивные. Первая фиктивная скважина включается в момент времени t = 0 с постоянным дебитом q1 q0 (в данном случае q0 = 0) и работает на протяжении времени t1 + t2 + t3. График работы первой фиктивной скважины изображен на рис. 2.58. Изменение давления в точке пласта, вызванное работой первой фиктивной скважины, определится в соответствии с основным уравнением упругого режима:

Pi

или

Pi

= μ (q1 ) Ei

4πkh

=

μ

(q1 ) ln

4πkh

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

4χ (t

+ t

 

+ t )

 

 

 

 

 

 

 

 

 

1

 

2

3

 

 

 

 

 

 

 

 

2,246 χ (t1 + t2 + t3 )

.

(2.359)

 

 

 

r2

 

 

 

 

 

Вмомент времени, отмеченный на рис. 2.59 точкой Тʹ, включается

вработу вторая фиктивная скважина с постоянным дебитом (q2 q1), работающая на протяжении периода времени (t2 + t3) или (T t1). В данном случае дебит второй фиктивной скважины будет величиной отрицательной.

166

Рис. 2.58. График работы первой фиктивной скважины

Изменение давления в точке пласта, вызванное работой второй фиктивной скважины, определится как

 

 

 

μ

(q2 q1 )

 

 

 

 

r

2

 

 

 

 

 

 

P2

=

 

 

Ei

 

 

 

 

 

 

4πkh

4χ (t

 

+ t

 

)

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P =

μ

 

(q

q

)ln

2,246 χ (t2 + t3 )

.

(2.360)

4πkh

 

 

r2

 

 

 

 

 

 

2

 

2

1

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.59. График работы второй фиктивной скважины

Очевидно, что Р также будет величиной отрицательной. Физический смысл этого заключается в следующем: давление определяется как Р0 Р, и снижение дебита реальной скважины (снижение отбора из пласта) приводит к повышению давления.

167

В момент времени Тʹʹ (рис. 2.60) включается в работу третья фиктивная скважина с постоянным дебитом (q3 q2), которая работает на протяжении времени t3 или (Т t1 + t2).

Рис. 2.60. График работы третьей фиктивной скважины

Изменение давления в точке пласта, вызванное работой третьей фиктивной скважины, определится как

 

 

 

μ

 

 

 

 

 

 

 

r

2

 

 

 

 

 

P3

=

 

(q3 q2 )

Ei

 

 

 

 

 

4πkh

4χ (t

 

 

 

 

 

 

 

 

 

 

 

 

3

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,246 χ (t3 )

 

 

 

P =

μ

 

(q

q

)ln

.

(2.361)

4πkh

 

 

3

 

3

2

 

 

 

 

r2

 

 

 

 

 

 

 

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

P = P1 + P2 + P3 .

И давление в точке пласта

P = P0 P.

2.14.5. Движение жидкости в экранированных пластах

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

168

непроницаемой границей (экраном). Наличие экрана существенным образом сказывается на распределении давления вокруг работающей скважины (или группы скважин), так как упругая волна давления, дойдя до плоскости экрана, отражается от нее и распространяется в обратном направлении. По этой причине все точки реального пласта находятся под действием как прямой, так и отраженной волны. На основе таких рассуждений мы можем воздействие отраженной волны заменить другой прямой волной от действия некоторой фиктивной (отраженной) скважины, т.е. заменить экранированный пласт бесконечным. В этом случае для достижения принципа эквивалентности действия, фиктивная скважина – отражение должна быть зеркальным отражением действительной скважины, т.е. иметь тот же дебит, те же размеры, находиться на одинаковом расстоянии от экрана (но с противоположной стороны от него) и работать синхронно с реальной скважиной (рис. 2.61). Тогда давление в некоторой точке А можно определить по принципу суперпозиции (наложение полей), используя основную формулу упругого режима для бесконечного изотропного пласта:

pA (r,t ) = p0

 

q μ

 

 

 

r2

 

 

q μ

 

 

 

r2

 

 

 

 

Ei

1

 

 

 

Ei

2

 

. (2.362)

 

4kt

 

4kt

 

 

4πkh

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.61. Пример экранированного пласта

Тогда давление в точке А можно определить, используя основное уравнение упругого режима:

169

P(r, t ) = P

q μ

ln

2, 246 χ t

4πkh

r2

 

 

 

 

 

 

 

 

 

 

 

1

 

q μ

ln

2,246 χ t

.

(2.363)

4πkh

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

2

 

 

 

Этот принцип можно распространить на любую группу взаимодействующих скважин. Если экран будет не прямолинейный, а иметь более сложную конфигурацию, то количество отраженных скважин будет больше. Так, в случае когда в пласте имеются два экрана, пересекающихся под углом θ = 120°, то количество отра-

женных скважин будет 2, при θ = 90° − 4. Это объясняется тем, что в таких случаях плоскости экранов отражаются друг в друге, тем самым образуя дополнительные отражения. Очевидно, что количество определяется углом между плоскостями экранов: n = 360° / θ (при этом принимается в расчет лишь целое количество отражений; при наличии в пласте параллельных экранов построение отражений соответствует правилам оптических построений).

2.14.6. Гидродинамические исследования скважин при неустановившихся режимах

(методом восстановления давления)

На использовании основного уравнения упругого режима разработан эффективный метод определения фильтрационных параметров пласта – метод восстановления давления. В скважину, работающую при постоянном режиме с дебитом q, спускается глубинный манометр. Скважина останавливается, забойное давление постепенно восстанавливается до величины пластового. Темп восстановления забойного давления регистрируется манометром. По результатам замеров строится график – кривая восстановления давления (КВД) (рис. 2.62).

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

170