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

книги / Физика и гидравлика нефтяного пласта

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

Поэтому выражение (IV.48) будет:

&v dv

(IV.49)

* dxi — dt >

т. е. полученное дифференциальное уравнение имеет такой же вид, как и (IV.34). Граничные и начальные условия следующие: при х = О

 

ис= const,

 

(IV.50)

при t = 0, поскольку

давление

всюду постоянно

(др/дх = 0),

скорость равна нулю. Это дает

 

 

 

t»c =

®(0, t),

v(x,

0) = 0.

(IV.51)

Сравнивая дифференциальное уравнение (IV.49) и условия (IV.50) и (IV.51) с соответствующим дифференциальным урав­

нением (IV.34) и условиями (IV.38)

и (IV.41),

видим: p -^ v ,

Pc-^Vc, Рк-э-0.

 

 

 

дифференциального

Поэтому, заменив в решении (IV.43)

уравнения

(IV.34) р через v, рс через vc и рк через нуль, полу­

чим решение рассматриваемой задачи:

 

 

 

 

 

 

и = ис(1 — erf £)•

 

 

(IV.52)

Очевидно,

что вдали

от галереи

(|-»-оо)

и

и-»-0. Учитывая

k

др

 

 

 

 

 

 

V --- ------------------ -— , получим

 

 

 

 

 

дх

 

 

 

 

 

 

 

- ^ • ^ F = Uc( 1-

erf|)-

 

 

(IV.53)

Интегрируя (IV.53)

в пределах

от

нуля

до

х,

имеем:

 

 

 

 

X

 

 

U

 

 

 

 

S 1 — erf (

$~ЗГdx =

р (*’

0 = —

к

2 V t.t )]du.

(IV. 54)

В подынтегральном выражении х заменено через ы, чтобы подчеркнуть различие между обозначениями предела интеграла и аргументом.

Интеграл (IV.54) можно представить в виде безразмерной формы:

к[р(х,

t ) - P {0, ()_]_

(1 _ erf u)du.

(IV.55)

 

2n V ^ »c

 

о

 

Интегрируя

(IV.55)

по частям, получим окончательно:

Т

Р ^°’ 0 1

=6(1 -

erf 6) + ——l {l - e ~ f ) .

(IV.56)

2\iyxtvc

 

 

\ я

 

Вычислим максимальный перепад давления:

 

 

Р(оо,

t ) — р(0,

t) =

рК— р с.

(IV.57)

III

При JC- VOO (т. е. оо) первый член правой части выраже­ ния (IV.57) стремится к нулю, поэтому полный перепад давле­ ния будет

k(pK— P c ) _ _

2^

(IV.58)

 

V я

С

 

или

 

 

[J- V

Vc-

(IV.59)

k

 

ПРОСТЕЙШЕЕ РЕШЕНИЕ, СООТВЕТСТВУЮЩЕЕ ПРИТОКУ ЖИДКОСТИ К ТОЧЕЧНОМУ СТОКУ НА ПЛАСТЕ БЕСКОНЕЧНОЙ ПРОТЯЖЕННОСТИ

Пусть на плоском пласте бесконечной протяженности имеется точечный сток.

Для этого случая дифференциальное уравнение имеет сле­ дующий вид:

• ( - & + 7 - & - H F - <IV-60>

С помощью простой подстановки можно показать, что решением

(IV.60) является

 

 

 

 

 

 

 

 

 

 

 

£

 

 

 

C2, 5 = - ^ - .

(IV.61)

p =

p(r,t) = Cl \ ± - e - “4*du +

Примем,

что

в начальный

/ = 0

давление

всюду

равно рк,

т. е. что при

 

 

 

 

 

 

 

 

 

 

 

%=

ОО

^ =

0,

Р =

рк.

 

(IV.62)

Это дает

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рк =

с Л - 1 ~ е - “!/<с1и + С2.

 

(IV.63)

 

 

 

о

и

 

 

 

 

 

Продифференцируем (IV.61):

 

 

 

 

 

| £ . „ с , |

е -« '-< - ^ г =

с

4

 

(IV.64)

Дебит стока, очевидно, равен

 

 

 

 

Q0 = **>L(r °P)

= М

С1( г 1 е- ^

= - ^ С , .

Р

\

дг )г<= о

 

j i

1 \

г

) г ~ о

Р

 

 

 

 

 

 

 

 

 

(IV.65)

112

Подставив значение Сi из уравнения (IV.65) в (IV.63), най­ дем Сг■ Далее из (IV.61) получим распределение давления в пласте в следующем виде:

I

i e - ^ d

u .

(IV.66)

..

г

 

 

 

Введем новую переменную w =

и2/4.

 

 

 

Тогда уравнение (IV.66) преобразуется к следующему виду:

 

оо

 

 

 

 

S

 

 

<IV-67>

 

Г*

 

 

 

 

4xt

 

 

 

Интервал в (IV.67) называется экспоненциальной функцией

и обозначается так:

 

Г2

Ч

 

— E i (—JC) = —Ei ^

(IV.68)

ш

)

 

 

 

Итак, окончательно имеем

 

Г2

Ч

 

Рк~ Р =

 

(IV.69)

 

Ш ) '

 

 

 

Для малых значений аргумента x = r2/(Ant) можно пользо­ ваться приближенной формулой

 

—Ei (—х ) ------у — 1п * +

* —

я2 +

О (JC3) ,

 

(IV.70)

где у =

0,5772... — постоянная Эйлера.

 

 

 

 

 

Предыдущую формулу можно переписать в виде

 

 

/

г/2- 2 \

,

2 2,25^2 5 f-t

,. г2

г*

 

 

- n (

 

\

— (

4*г) =

п

г2

^ Ш

6 4 ( *

<

) 2

U V 6 4

( * . f )

3 ) '

 

 

 

 

 

 

 

 

 

 

(IV.71)

С очень высокой степенью точности при r2/(4xt) <

0,1 можно

пользоваться приближенной формулой

 

 

 

 

 

 

Рк — Р =

 

, Q o j n

2 , 2 5 f-t

pQo

 

jn

1,5 л/**

 

(IV.72)

 

 

4 я kh

 

2nkh

 

 

 

 

 

Эту формулу можно интерпретировать как формулу Дюпюи,

в которой числитель 1 , 5

— условный

радиус ri(f),

завися­

щий от времени:

 

 

К>о

| п

 

 

 

 

 

 

Р к~ Р

 

 

 

(IV.73)

 

 

2 я kh

г

 

 

 

 

 

 

 

 

 

 

 

 

 

Для точечного стока (или источника) не существует давле­ ния на стоке (источнике). В самом деле, примем в формуле (IV.69) г = 0, получим, что —Ei(—х) равна бесконечно боль­ шой величине.

8 Заказ № 283

ИЗ

От точечного стока мы перейдем к реальной скважине весьма малого радиуса Rc и будем считать давление на этой окружно­ сти равным рс. Итак, имеем

Рк Рс

У?0

р;

( - 4

- )

(IV.74)

4nkh

 

или

 

 

 

 

 

 

 

Рк — Рс =

^

In-

1,5 д/х/

 

(IV.75)

 

2nkh

111

Rc

 

 

при условии, что (У * 0 /Яс велико.

РЕШЕНИЕ, СООТВЕТСТВУЮЩЕЕ ПРИТОКУ ЖИДКОСТИ К ПРЯМОЛИНЕЙНОЙ ГАЛЕРЕЕ В ПОЛОСООБРАЗНОМ ПЛАСТЕ ОГРАНИЧЕННОЙ ПРОТЯЖЕННОСТИ

Поместим начало координат на галерее и рассмотрим две модели пласта ограниченной протяженности:

в пласт длины L отсутствует приток жидкости со стороны

внешней границы

=

на внешней границе (контуре питания) задано постоянное давление.

Пусть давление на галерее все время постоянно и равно р, начальное пластовое давление равно рк.

Начальные и граничные условия будут:

Р (*> 0) = рк, Р (0, t) = Рс, др% П- = 0.

(IV.76)

Решить первую задачу можно как классическим методом Фурье, так и операционным. Решение имеет вид (36)

оо

Рк — Р Рк Рс

Х ех р (

2 л — 1

cos (2 п — 1 )

я (L — х)

 

2L

(2л - ip

4

X

(IV.77)

где критерий Фурье Fo = xt/L2.

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

Рк — р

ао

» — _г -- +

2 (л — 1) + ^ ~

|

Z ( - u n+i erfc

erfc ■

-\/Fo

j

Рк рс

п= 1

2-v/Fo

2

где erfc и =

1 — erf и.

 

 

(IV.78)

 

 

 

 

114

Изучим характер изменения давления на непроницаемой гра­ нице пласта x = L при малых значениях Fo, т. е. в начале экс­ плуатации галереи.

На этой границе имеем

оо

Рк- P j t) = 2 У

(-l)"+ » e r fc

2п~ - .

(IV.79)

Рк — Рс

2 д / F o

'

При малых значениях Fo можно воспользоваться прибли­ женной формулой (IV.79), которая в данном случае при усло­ вии сохранения одного члена ряда (IV.79) имеет вид

 

PK— P{L)

я*—^ д /F o e

4Fo

(IV.80)

или

Рк -- Рс

V я

 

 

 

 

,

______ !_

 

 

Р (L) рс

(IV.81)

а

1 — -.---v/Foe

4Fo .

 

Рк Рс

у

%

 

 

При определении давления на непроницаемой границе при­

ближенная

формула

дает

вполне

приемлемые

результаты до

Fo < 0 ,1 (точность до

~ 1 0

% ). Отметим вместе с тем, что при

определении изменения давления на этой границе необходимо пользоваться точной формулой (IV.77), взяв достаточно боль­ шое число членов ряда.

Расчеты показывают, что до значения Fo<0,08 а практиче­ ски мало отличается от первоначального значения а = 1. Од­ нако из этого не следует заключение о малости изменения дав­

ления на этой границе. Действительно, примем

рк =

20 МПа,

Рс = 5 МПа, ркрс =

15 МПа;

тогда

при Fo =

0,08, а = 0,97

и p{L) = р0+ а(р„ — рс) =

5+0,97 •15 = 19,55 МПа,

т. е. дав­

ление изменилось на вполне ощутимую величину 0,45 МПа.

Для определения дебита галереи продифференцируем выра­

жение (IV.77). Объемный дебит при этом будет

 

 

Q = hb д

 

 

 

khb

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Х ( - 1 ) 2п ехр [ — (2п 4 1 — .+Fo|.

 

(IV.82)

При малых значениях F o (F o < 0 ,l)

эта

формула

неудобна,

так как ряд сходится очень медленно.

 

формулой (IV.78).

В этих случаях можно

воспользоваться

Объемный дебит, определенный по этой формуле, будет

Q = hb k_ ( дР \

_

khb

(Рк — Рс)

1

 

 

Р- V дх ) х = о

(1

 

 

L

VitFo i=E ix

 

/

(

я

- О

1

 

 

 

(IV.83)

Х ( — l)n + IVe

Fo

е

 

 

8*

115

В

случае

больших

значений F o (F o > 0 ,l)

в уравнении

(IV.82)

можно сохранить лишь первый член ряда,

и тогда оно

примет более простой вид

 

 

 

 

 

Q = 2exp ( - - ^ F o ) .

(IV.84)

При

малых

значениях

F o (F o < 0 ,l)

также можно ограни­

читься первым членом ряда (IV.83) (п =

1) и получим для этого

члена 1 — e~1/Fo« 1 .

 

 

 

Тогда окончательно имеем

 

 

 

 

 

1

 

(IV.85)

 

 

 

Q — VnFo ’

 

 

 

 

 

 

где Q — безразмерный дебит, определяемый

 

 

 

Q =

V-LQ

 

(IV.86)

 

 

khb (рк — рс)

 

Из уравнения (IV.85) видно, что в начальные периоды экс­ плуатации дебит не зависит от длины пласта. Действительно, преобразуем уравнение (IV.85) следующим образом:

Q д/Fo

(JLQ V ч-t

[J.Qл/ч-t

_

1

(IV.87)

khb (рк — рс) L

khb (Рк Рс)

-\jя

 

 

Это уравнение в точности совпадает с уравнением для объ­

емного дебита пласта бесконечной протяженности

 

[см. уравне­

ние (IV.44)].

 

 

 

 

 

Таким образом, приходим к очень важному выводу: когда время мало, процесс фильтрации не зависит от условий на про­ тивоположном конце пласта. Этот вывод согласуется с выводом, сделанным нами выше в отношении характера распределения давления.

Если построить кривые распределения давления для пластов конечной и бесконечной протяженности, то они должны прак­ тически совпадать при малых Fo (особенно в окрестности гале­ реи). При малых Fo мы имеем основание пользоваться для рас­ чета дебита формулой (IV.87) как более простой.

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

Начальное и граничное условия в этой задаче следующие:

р(х, 0) = рк, р (0, 0 = Рс, P(L, t) = pK.

(IV.88)

Решение задачи методом Фурье или операционным имеет вид

Р к -Р

L~ X

 

оо

( - Ц * +

r- ^ F n cin ПК(Ь-Х)

'

2

V

Рк — Рс

L

я

L J

п

L

'

(IV.89)

116

р к — р с
L

Для

малых

значений

F o (F o < 0 ,l) существует

другая фор­

мула:

 

X

 

 

 

 

 

 

 

 

2n+ J — erfc

 

Р к ~ Р

erfc

L

£

erfc

 

Рк --- Рс

 

2 л/FO +

 

2 + F o

2 д/Fo

 

 

 

/2 = 1

 

 

 

 

 

 

 

 

(IV.90)

Найдем объемный дебит Q по общей формуле

(IV.89)

 

 

 

 

оо

 

 

 

 

Q =

1 + 2 £

e -" W o.

(IV.91)

 

 

 

 

л=1

 

 

По приближенной формуле (IV.90) для малых Fo имеем

Q =

1

V^Fo

 

Так как Fo мало, то в послед­ ней формуле можно пренебречь вторым членом, и тогда мы полу­ чим формулу (IV.87) для слу­ чая непроницаемой границы. Приходим к аналогичному выво­ ду: при малых Fo порядка до 0,1 условия на внешней границе не оказывают заметного влияния на характер изменения дебита га­ лереи (дебит можно рассчиты­ вать по формулам для бесконеч­ ного пласта).

Из общей формулы (IV.91) следует, что с возрастанием Fo Q стремится к единице, т. е. наб­ людается беспредельное прибли­ жение к стационарному состоянию

khb

И

(IV.92)

Рис. IV.5. Кривые распределения дав* ления в полосообраэном пласте при за­ данных давлениях на контуре и галерее

(IV.93)

Пусть n2Fo = 3,

что

соответствует

F o«0 ,3 .

Тогда по фор­

муле (IV.91) имеем

оо

 

 

 

 

 

 

 

Q = 1 + 2

£

е -',г*!р0« 1 + 2

* 0 ,0 4 98 »

1,1,

Л= 1

т.е. текущий дебит всего на 10 % отличается от стационарного значения.

На рис. IV.5 приводится примерное распределение давления

по пласту для ряда значений Fo.

117

ПРИНЦИП СУПЕРПОЗИЦИИ. ТЕОРЕМА ДЮГАМЕЛЯ

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

две простые задачи.

Пусть в течение некоторого вре­ мени t\ давление на галерее посто­ янно и равно рс,- Тогда мы можем

 

воспользоваться

уравнением (IV.43)

 

для определения давления р\ в лю­

 

бой точке

пласта,

заменив в нем

 

erf через erfc:

 

 

 

 

 

Рк — Pi =

(Рк -

Pc,) erfc ~ ^ = г - .

 

 

 

 

 

 

 

(IV.94)

 

В некоторый момент времени t\

 

давление на галерее изменилось и

 

стало рс2-

 

 

 

 

 

 

Задача

решается

в два этапа.

 

Вначале подбирается такое решение,

 

чтобы

оно при сложении с

(IV.94)

 

обеспечивало бы при t>t\ на гале­

 

рее нулевой перепад давления. Та­

 

ким решением, очевидно, будет

 

— (Рк — Pc,)erfc-^-^===

tx)

 

 

 

 

 

 

 

(IV.95)

 

На рис. IV.6, а кривая соответст­

 

вует

распределению

давления по

 

формуле (IV.94) для момента вре­

 

мени t = t\ + At,

где

At — малая ве­

 

личина; кривая б построена по фор­

Рис. IV.*. К выводу формулы

муле

(IV.95) для этого же момента

(IV.97)

времени. Сложение

 

решений, оче­

 

видно,

обеспечит нулевой

перепад

давления на галерее, которой сохраняется для любого t>t\. Для того чтобы получить окончательный результат, необхо­

димо подобрать решение, соответствующее новому граничному условию, а именно:

(Р . — Ре,) eric 2 у . _ tl (IV.96)

118

Кривая рис. IV.6 , б построена по формуле для момента вре­

мени t =

fi+Д /.

 

всех трех

решений

распределение

давления

При

сложении

р2 после того, как

давление

на

галерее изменилось

с рс, до

рСг, будет:

 

 

 

 

 

 

 

 

Рк

Р

(Рк

Pci) erfc

^

-(- [(рк рС2)

 

 

 

-

( Р

. - <■«,)] «1с

3 Vx

 

(IV-97)

где индекс 2 при р опущен (см. рис. IV

в, г).

давление,

Первый член уравнения (IV.97) характеризует

обусловленное

предшествующим

перепадом давления

ркрС],

влияние которого, очевидно, сохраняется и после того, как на галерее установилось новое давление рС2. Второе слагаемое по

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

 

[(Рк — Р^) — (Рк — Рсг_0] erfc

2л/%

_ 7 J77) •

(IV.98)

 

Вычислим объемный

дебит для

t > t i

по формуле

(IV.97)

r t

ии k

( др \

__ khb

I Рк — Дс,

t

(Рк — Рс2) — (Рк — Pc,) {

4

по (i

\дх Л г-о—

и

(

 

 

Vnx ( / - # , )

Г

(IV.99)

Аналогичным способом решается задача для случая стока (источника).

Рассмотрим линейный сток. Пусть при t <. U дебит стока был Qi, тогда давление будет

 

 

 

-

w

E I ( - i ) -

<IV ' 100>

При

дебит изменился до Q2. Начиная с этого момента,

в соответствии с принципом

суперпозиции давление станет

Рк — Pi =

K h

с , /

г 2

\

Q2 — Q1 т - . /

г 2

4nkh

Bt { — Е г )

4nkh

 

 

(IV. 101)

Для малых значений r2/(4xf) будем иметь

 

 

2,25хГ

, |л (С?2 — Q i ) j n 2 , 2 5 x ( < - ^ r i )

(IV. 102)

b - p - j &

b

г2

+

4nkh

 

119

Очевидно, формулой (IV. 102) нельзя пользоваться при t7 близком к ti (тем более при f = /i), так как условие малости

аргумента ——— -----— при этом не соблюдается. 4х(< — U)

При многоступенчатом изменении Q задача решается таким же способом, как и рассмотренная выше.

Использование описанного выше расчетного приема неудобно в тех случаях, когда число ступеней изменения давления или расхода велико.

Для решения задач с многоступенчатым изменением входных параметров или при непрерывном изменении их применяют тео­

 

 

рему (формулу) Дюгамеля, к вы­

 

 

воду которой мы

сейчас

перей­

 

 

дем.

 

 

 

 

q>(f)

функ­

 

 

 

Обозначим через

 

 

цию erfc в уравнении

(IV.94) или

 

 

Ei

в уравнении

(IV. 100).

<р (f)

 

 

можно

рассматривать

как еди­

 

 

ничную

функцию,

т.

е.

как рас­

 

 

пределение

давления, соответст­

 

 

вующее

единичному

перепаду

 

 

давления

ркрс =

1,

в задаче

 

 

с

заданным

постоянным

давле­

 

 

нием на галерее или единичному

 

 

расходу

\xQol(4nkh) =

1 в задаче

 

 

с

заданным

постоянным

рас­

Рис. IV.7. Кривая

изменения перепада

ходом.

 

 

 

 

 

 

 

давления

на галерее

 

Для

определенности

рассмот­

 

 

рим случай,

когда перепад дав­

ления на галерее изменяется по некоторому закону f(t). Разо­ бьем кривую на ряд малых интервалов времени At и предста­ вим его в виде большого числа ступеней (рис. IV.7).

Высота этих ступеней будет

 

 

ПО) Л;,

V (2 At) At. . .

(IV. 103)

На первом этапе t = At+0

процесс можно

представить

в виде

 

 

/(0)ф(ДО +

/ '( ° ) д МО).

(IV. 104)

Первый член обусловлен действием постоянного перепада давления f( 0 ), он аналогичен первому слагаемому уравнения (IV.97). Второй член соответствует изменению перепада Н О ) X XAf (первая ступень) и вполне аналогичен второму слагаемому этого же уравнения.

На второй ступени (А^+ А^+ 0 ) действие постоянного пере­ пада давления сохраняется, но аргумент единичной функции уже вдвое больше, он равен cp(2A£)- Продолжительность пере­

120