Материал: Разработка нефтяных месторождений.

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Продолжение таблицы 1

 

ег f(x)

 

er f(x)

 

er f(x)

 

er f w

0 ,1 8

0 ,2 0 0 9

0 ,4 3

0 ,4 5 6 9

0 , 6 8

0 ,6 6 3 8

0 , 9 3

0 ,8 1 1 6

0 ,1 9

0 ,2 1 1 8

0 ,4 4

0 ,4 6 6 2

0 , 6 9

0 ,6 7 0 8

0 , 9 4

0 ,8 1 6 3

0 , 2 0

0 ,2 2 2 7

0 , 4 5

0 ,4 7 5 5

0 , 7 0

0 ,6 7 7 8

0 , 9 5

0 ,8 2 0 9

0 ,2 1

0 ,2 3 3 5

0 , 4 6

0 ,4 8 4 7

0 ,7 1

0 ,6 8 4 7

0 , 9 6

0 ,8 2 5 4

0 , 2 2

0 ,2 4 4 3

0 , 4 7

0 ,4 9 3 7

0 ,7 2

0 ,6 9 1 4

0 ,9 7

0 ,8 2 9 9

0 , 2 3

0 ,2 5 5 0

0 , 4 8

0 ,5 0 2 7

0 , 7 3

0 ,6 9 8 1

0 , 9 8

0 ,8 3 4 2

0 , 2 4

0 ,2 6 5 7

0 , 4 9

0 ,5 1 1 7

0 , 7 4

0 ,7 0 4 7

0 ,9 9

0 ,8 3 8 5

1 ,0 0

0 ,8 4 2 7

1,25

0 ,9 2 2 9

1 ,5 0

0,9661

1 ,7 5

0 ,9 8 6 7

1,01

0 ,8 4 6 4

1 ,2 6

0 ,9 2 5 2

1,51

0 ,9 6 7 3

1 ,7 6

0 .9 8 7 2

1 ,0 2

0 ,8 5 0 8

1 ,2 7

0 ,9 2 3 5

1 ,5 2

0 ,9 6 8 4

1 ,7 7

0 , 9 8 7 7

1 ,0 3

0 ,8 5 4 8

1 ,2 8

0 ,9 2 9 7

1,53

0 ,9 6 9 5

1 ,7 8

0 ,9 8 7 2

1 ,0 4

0 ,8 5 8 6

1 ,2 9

0 ,9 3 1 9

1 ,5 4

0 ,9 7 0 6

1 ,7 9

0 ,9 8 8 6

1 ,0 5

0 ,8 6 2 4

1 ,3 0

0 ,9 3 4 0

1 ,5 5

0 ,9 7 1 6

1 ,8 0

0 ,9 8 9 1

1 ,0 6

0 ,8 6 6 1

1,31

0 ,9 3 6 1

1 ,5 6

0 ,9 7 2 6

1,81

0 ,9 8 9 5

1 ,0 7

0 ,8 6 9 8

1 ,3 2

0 ,9 3 8 1

1 ,5 7

0 ,9 7 3 6

1,82

0 ,9 8 9 9

1 ,0 8

0 ,8 7 3 3

1,3 3

0 ,9 4 0 0

1 ,5 8

0 ,9 7 4 5

1 ,8 3

0 ,9 9 0 3

1 ,0 9

0 ,8 7 6 8

1 ,3 4

0 ,9 4 1 9

1 ,5 9

0 ,9 7 5 5

1 ,8 4

0 ,9 9 0 7

1,10

0 ,8 8 0 2

1 ,3 5

0 ,9 4 3 8

1 ,6 0

0 ,9 7 6 3

1 ,8 5

0 ,9 9 1 1

1,11

0 ,8 8 3 5

1 ,3 6

0 ,9 4 5 6

1,61

0 ,9 7 7 2

1 ,8 6

0 ,9 9 1 5

1 ,1 2

0 ,8 8 6 8

1,37

0 ,9 4 7 3

1 ,6 2

0 ,9 7 8 0

1 ,8 7

0 ,9 9 1 8

1 ,1 3

0 ,8 9 0 0

1 ,3 8

0 ,9 4 9 0

1 ,6 3

0 ,9 7 8 8

1 ,8 8

0 ,9 9 2 2

1 ,1 4

0 ,8 9 3 1

1 ,3 9

0 ,9 5 0 7

1 ,6 4

0 ,9 7 9 6

1 ,8 9

0 ,9 9 2 5

1 ,1 5

0 ,8 9 6 1

1 ,4 0

0 ,9 5 2 3

1 ,6 5

0 ,9 8 0 4

1 ,9 0

0 ,9 9 2 8

1 ,1 6

0,8 9 9 1

1,41

0 ,9 5 3 9

1 ,6 6

0,9 8 1 1

1,91

0,9 9 3 1

1 ,1 7

0 ,9 0 2 0

1,42

0 ,9 5 3 8

1 ,6 7

0 ,9 8 1 8

1,92

0 ,9 9 3 4

1 ,1 8

0 ,9 0 4 8

1 ,4 3

0 ,9 5 6 9

1 ,6 8

0 ,9 8 2 5

1 ,9 3

0 ,9 9 3 7

1,1 9

0 ,9 0 7 6

1 ,4 4

0 ,9 5 8 3

1,69

0 ,9 0 8 2

1 ,9 4

0 ,9 9 4 0

1 ,2 0

0 ,9 1 0 3

1 ,4 5

0 ,9 5 9 7

1 ,7 0

0 ,9 8 3 8

1 ,9 5

0 ,9 9 4 2

1,21

0 ,9 1 3 0

1,46

0 ,9 6 1 1

1,71

0 ,9 8 4 4

1 ,9 6

0 ,9 9 4 4

1 ,2 2

0 ,9 1 5 5

1 ,4 7

0 ,9 6 2 4

1 ,7 2

0 ,9 8 5 0

1 ,9 7

0 ,9 9 5 0

1 ,2 3

0 ,9 1 8 5

1 ,4 8

0 ,9 6 3 7

1 ,7 3

0 ,9 8 5 6

1 ,9 8

0 ,9 9 5 0

1 ,2 4

0 ,9 2 0 5

1 ,4 9

0 ,9 6 4 9

1 ,7 4

0,9861

1 ,9 9

0 ,9 9 5 0

ПРИЛОЖЕНИЕ 2

Интегральная показательная функция

Расчет распределения давления в бесконечном пласте в случае притока жидкости к точечному стоку при упругом режиме про­ изводится по формуле (11.132), в которую входит интегральная показательная функция —Ei(—х). Эта функция, как и erf(x) не выражается через конечное число алгебраических функций. Значения функции —Ei(—л:) приведены в табл. 2.

31G

Т а б л и ц а

2

 

 

 

 

 

~ E i (-*>

 

~ E i (-* )

 

- Е i ( - а)

0,00

ОО

0,26

1.0139

0,52

0,5362

0,01

4,0379

0,27

0,9849

0,53

0,5250

0,02

3,3547

0,28

0,9573

0,54

0,5140

0,03

2,9591

0,29

0,9309

0,55

0,5034

0,04

2,6813

0,30

0,9057

0,56

0,4930

0,05

2,4679

0,31

0,8815

0,57

0,4830

0,06

2,2953

0,32

0,8583

0,58

0,4732

0,07

2,1508

0,33

0,8361

0,59

0,4636

0,08

2,0269

0,34

0,8147

0,60

0,4544

0,09

1,9187

0,35

0,7942

0,61

0,4454

0,10

1,8229

0,36

0,7745

0,62

0,4366

0,11

1,7371

0,37

0,7554

0,63-

0,4380

0,12

1,6595

0,38

0,7371

0,64

0,4197

0,13

1,5889

0,39

0,7194

0,65

0,4115

0,14

1,5241

0,40

0,7024

0,66

0,4036

0,15

1,4645

0,41

0,6859

0,67

0,3959

0,16

1,4092

0,42

0,6700

0,68

0,3883

0,17

1,3578

0,43

0,6546

0,69

0,3810

0,18

1,3098

0,44

0,6397

0,70

0,3738

0,19

1,2649

0,45

0,6253

0,71

0,3668

0,20

1,2227

0,46

0,6114

0,72

0,3599

0,21

1,1829

0,47

0,5979

0,73

0,3532

0,22

1,1424

0,48

0,5848

0,74

0,3467

0,23

1,1099

0,49

0,5721

0,75

0,3403

0,24

1,0762

0,50

0,5598

0,76

0,3341

0,25

1,0443

0,51

0,5478

0,77

0,3280

0,78

0,3221

1,00

0,2194

3,2

0,01013

0,79

0,3163

1, Ю

0,1860

3,3

0,008932

0,80

0,3106

1,20

0,1584

3,4

0,00789

0,81

0,3050

1,3

0,1355

3,5

0,00697

0,82

0,2996

1,4

0,1162

3,6

0,00616

0,8d

0,2943

1,5

0,1000

3,7

0,003448

0,84

0,2891

1,6

0,08631

3,8

0,004820

0,85

0,2840

1,7

0,07465

3,9

0,004267

0,86

0,2790

1,8

0,06471

4,0

0,003779

0,87

0,2742

1,9

0,05620

4,1

0,003349

0 ,8й

0,2694

2,0

0,04890

4.2

0,002969

0,89

0,2647

2,1

0,04261

4,3

0,002633

0,90

0,2602

2,2

0,03719

4,4

0,002336

0,91

0,2557

2,3

0,03250

4,5

0,002073

0,92

0,2513

2,4

0,02844

4,6

0,001841

0,93

0,2470

2,5

0,02491

4,7

0,001635

0,94

0,2429

2,6

0,02185

4,8

0,001453

0,95

0,2387

2,7

0,01918

4,9

0,001291

0,96

0,2347

2,8

0,01686

5,0

0 j 001148

0,97

0,2308

2,9

0,01482

6,0

0,00036

0,98

0,2269

3,0

0,01304

7,0

0^000116

0,99

0,2231

3,1

0,01149

8,0

0,00004

317

ПРИЛОЖЕНИЕ 3

Решение задачи о притоке жидкости из бесконечного пласта к скважине конечного радиуса

Распределение давления в случае притока жидкости к скважи­ не конечного радиуса из бесконечного пласта при упругом ре­ жиме с постоянным дебитом q рассчитывается по формуле (3.15), получаемой на основании формулы (3.14). Исходным для получения этой формулы является дифференциальное урав­ нение упругого режима (3.11) с граничным и начальным усло­ виями (3.12). Для решения указанного дифференциального уравнения используется преобразование давления по Лапласу, определяемое формулой (3.13).

 

Функция f(p, т), определяемая формулой (3.14), также удов­

летворяет

уравнению

теплопроводности. При

этом

получаем

следующую постановку задачи для определения /(р, т):

 

а*/

 

_L 1

L

I L

 

 

 

 

 

 

 

 

(3.1>

dp*

 

 

 

 

 

 

 

 

 

 

Р

Ф

 

дх

 

 

 

 

 

 

 

 

/ = 0 при х =

О,

 

 

 

 

-> оо,

( * L -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем преобразование

Лапласа

функции / (р, т) в виде

/(Р. s) = J/(p, т)е

"’dr,

 

 

 

 

 

 

(3.2)

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

где s — некоторый параметр.

 

 

 

 

 

 

Умножим левую и правую часть дифференциального урав­

нения

(3.1)

на e-ST

и проинтегрируем

его от 0 до оо. Получим

Г

d2/

 

—sx

1

 

(*

д /

—sx

 

Г д /

—ST

dx-

 

(3.3>

J a ^ e

 

dx+ T

J-aTe

dx = ) 1 Г е

 

b

Для

первого

 

o

 

 

o

членов,

входящих

в правую часть

 

 

и

второго

(3.3), имеем следующие выражения:

 

 

 

 

Г

д/

—st

д

 

г*

—sx

 

df

 

 

 

 

 

j -

^

e

 

<'г = " Ф

. ) /е

dT— & •

 

 

 

 

о

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

О

Член,

 

 

 

О

 

 

части (3.3), преобразуется

следую­

 

стоящий

в правой

щим образом

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

оо

 

 

 

+ s j /е

sxdx=sJ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.5)

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

318

так

Первый член в правой части выражения (3.5) равен нулю,

как

при т = 0,

в соответствии

с начальным условием (3.1),

/ = 0.

граничного

условия (3.1)

при р= 1

получаем, применяя

 

Из

к нему преобразование Лапласа

 

 

df____ _1_

 

 

 

 

 

(3.6)

dp

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая (3.4) и (3-5), из (3.1), получаем обыкновенное

дифференциальное уравнение

 

 

d2/

 

1

df

-

 

 

 

 

dp2

+

р

dp

— s^-

 

 

 

(3,7)

Решение этого уравнения имеет вид

 

Т = А К 0 (1/ф ),

 

 

 

 

(3.8)

где

/Со — функция

Бесселя второго рода;

А — постоянная ин­

тегрирования дифференциального уравнения.

 

Выполним граничное условие (3.6).

 

 

В результате получим

 

 

 

df

 

dK0(x)

dx

AKi (У Г ) y T ,

 

dp |P = I

 

 

dx

dp

 

 

 

 

 

 

 

 

 

1

 

 

 

 

(3.9)

‘4 = s3/J /c„(yr)'

 

 

 

 

 

 

 

где K\ — также функция Бесселя.

 

 

 

Отсюда

 

 

 

 

 

_

Ко (/sp)

 

 

 

 

(3.10)

 

 

Кх(УГ)'

 

 

 

 

 

 

 

 

 

Переход от изображения по Лапласу функции /(р,т) к ори­ гиналу осуществляется с использованием формулы обращения

Меллина, т. е. Y~J-foO

1

■Уо (У»р)

 

/(Р. %) = 2JU !•

S3/2 Кг (1/Г) ds.

(3.11)

V—ioо

 

осуществляется на основе тео­

Вычисление интеграла (3.11)

рии функций комплексного переменного. В результате получа­ ем, что

 

 

оо

 

 

 

 

 

 

 

-

е “ Т) (и) Уо (up) - Yx (и) J0 M l du

(3.12)

Л(Р. х) - ■ И '- 1

« W

(«)+lV («)]

 

где Jo Ju

Yo,

Yi — соответствующие обозначения функций Бес­

селя.

р=1, т. е. на контуре скважины

 

При

 

,

2

( *

( 1

— е

 

(3.13)

 

Я

J

 

« W

( “) + 1V(«)1

 

 

 

319

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

А =

0

при т =

0,

р ----- »- оо ,

 

А =

1

при т ^

0,

р =

1,

(3.14)

f

Роо — Р

 

 

 

 

fl

Роо — Рс

’

 

 

 

где

рс— давление

в скважине, /?«, — давление при р— *оо.

В результате использования преобразования Лапласа функ­

ции

/i(p, т),

выполнения начального

и граничного условия

(3.14)

получаем

 

 

 

=

sKo (Vs)

 

 

 

(3.15)

 

 

 

 

 

После применения к f\, определяемой формулой (3.15), обращения Меллина и вычисления соответствующего интеграла с помощью функций комплексного переменного, имеем для /1 (р, т) следующее выражение:

-ц2т Jо (up) Vо (и) — Y0 (up) J0 (и)

du

(3.16)

/1( p , t ) = i + 4 'J е

l.’ W + V W

и

 

о

f(p, т) и /i(p, т) определяются численным

Значения функций

путем на ЭВМ.

 

 

 

ПРИЛОЖЕНИЕ 4

Течение жидкости в трещиновато-пористом пласте при упругом режиме

В главе II (см. пример 4) дано решение задачи о распределе­ нии давления при упругом режиме в прямолинейном нефтяном пласте с обычной пористостью конечной длины I. Рассмотрим ту же в принципе задачу для трещиновато-пористого пласта.

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

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

320

Источник: https://tut-files.ru/previewfile/82937