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

Проектирование разработки нефтяных и газовых залежей. Практикум учеб

.pdf
Скачиваний:
18
Добавлен:
15.11.2022
Размер:
1.4 Mб
Скачать

2πkh(Pн Рд ) =

µвq

=12z

 

+ 25z2

+1,7ln

rф

+

πLµ

0

 

ln

σ2

,

(120)

 

r

2σ

 

 

π2r r

 

ф

ф

 

 

 

 

0

 

 

 

 

 

 

 

сн

 

 

 

 

 

 

cд cн

 

 

где σ – половина расстояния между скважинами в ряду. Радиус текущего положения фронта водонефтяного кон-

такта

r =

3L

(2σ+0,05L),

(121)

ф 2π

где L – расстояние между рядами.

Значение времени подхода фронта воды к добывающим скважинам определяется по формуле

 

 

 

 

 

 

 

 

µ

в

mδr2

 

 

 

 

 

 

 

 

 

 

 

t =

 

 

ф

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

 

 

 

 

 

 

 

 

 

 

2k (Pн Рд )

 

 

 

 

 

 

 

 

 

 

πLµ0

 

 

 

 

 

σ

+(1,7

 

 

)ln

r

 

 

× 12z

 

+ 25z2

+

+ 2µ

 

ln

 

 

−µ

 

ф

 

. (122)

 

 

 

 

 

 

 

rсн

 

 

ф

ф

 

2σ

 

0

 

π

 

rr

 

0

 

е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание № 2.1 (исходные данные – в табл. П.3.1)

Система расположения скважин пятиточечная, залежь находится на начальной стадии разработки; давления на забоях добывающих и нагнетательных скважин Рд и Рн; динамическая вязкость воды и нефти µв и µн; проницаемость коллектора k; толщина пласта h; кратчайшее расстояние от нагнетательной до добывающей скважины d; радиусы добывающих и нагнетательных скважин rсд и rсн; остаточная нефтенасыщенность SНО; доля связанной воды Sсв и пористость m. Вытеснение нефти водой непоршневое. Определить дебит одной добывающей скважины

61

непосредственно перед ее обводнением и время подхода к ней фронта воды.

Задание № 2.2 (исходные данные – в табл. П.3.2) Определить разность забойных давлений нагнетательных

идобывающих скважин непосредственно перед началом обводнения добывающих скважин и время достижения фронтом вытеснения середины расстояния между добывающими и нагнетательными скважинами при следующих исходных данных: система размещения скважин по залежи – рядная; дебит добывающей скважины q; динамические вязкости воды и нефти µв

иµн; проницаемость коллектора k; толщина пласта h; расстояние между скважинами в ряду σ; радиусы добывающих и нагнета-

тельных скважин rсд и rсн; остаточная нефтенасыщенность SНО; доля связанной воды Sсв; пористость m; расстояние между рядами L. Вытеснение нефти водой непоршневое.

Пример выполнения работы

Приведем пример выполнения задания № 2.1 для варианта № 1.

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

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

(118):

q =

 

 

 

2πkh(Pн Рд )

 

 

 

 

 

.

 

 

 

+25z2

 

rф

 

 

 

d

2

 

µ

12z

+1,7ln

 

ln

 

 

r

 

4r r

 

 

в

ф

ф

 

 

0

 

 

 

 

 

 

 

 

сн

 

 

 

ф cд

Коэффициент zф для условий работы определяется по уравнению

62

100zф3 +87,5zф2 1,789 = 0.

Решим уравнение с помощью компьютерной программы

Derive:

zф = 0,1584,

q= 2 3,14 21 1012 11(24 15) 106 /1,9 103 ×

×12 0,1584 + 25 0,15842 +1,7ln 0,68 960 +0,12

 

 

 

 

 

+

5,1

ln

 

 

 

960

2

 

 

 

 

= 0,176 м3/с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,9

4 0,68 960 0,1

 

 

 

 

 

 

 

 

 

 

 

 

Определяется значение времени

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

mδr2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t =

 

 

 

в

 

ф

 

 

×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2k (Pн Рд )

 

 

 

 

 

 

 

 

 

 

× 12z

 

+ 25z2

+

πLµ0

 

 

 

ln

d 2

 

+(1,7 −µ

 

)ln

rф

 

,

 

 

 

 

 

 

 

 

4rr

 

 

r

 

 

 

 

 

 

ф

 

 

 

ф

 

 

2σ

 

 

0

 

 

 

 

 

0

 

 

е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сн

 

 

 

 

где

δ =1S

 

S

 

2 z

=10,03 0,095 2 0,1584 = 0,769, t =

 

 

 

 

св

 

 

 

НО

 

3

ф

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

=

1,9 103 0,235 0,769 (0,68 960)2

·

(12 0,1584 +25 0,1584

2

+

 

2 21 1012 (24 15) 106

 

 

 

 

 

 

5,1 103

 

 

 

 

 

 

9602

 

 

 

 

 

 

 

 

5,1 103

 

 

0,68

960

 

=

+

1,9 103

 

·

ln

 

 

 

 

 

 

 

 

+

 

1,7

 

 

 

·

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 0,12 0,1

 

 

 

 

 

 

1,9 103

 

 

0,1

2,72

 

 

=15,25 106 с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

63

Работа № 3

Определение распределения пластового давления для нефтяной залежи при упругом режиме

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

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

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

P = −

qµн

E

 

r2

,

(123)

 

 

 

 

 

4πkh

i

 

 

 

 

 

 

4χ t

 

 

где ∆P – изменение давления в точке пласта;

q – дебит (приемистость) возмущающей скважины;

64

r – расстояние от скважины до точки пласта, в которой определяется ∆P;

t – время с начала пуска или остановки скважины;

Ei – интегральная показательная функция, табулированная

всправочниках;

χ– пьезопроводность,

χ =

k

 

µн (βжm п ),

(124)

где βж, βп – коэффициенты сжимаемости жидкости и пористой среды.

Работа выполняется с помощью программных комплексов

Excel Microsoft Office и Surfer.

Задание № 3.1 (исходные данные в табл. П.3.3) Построить в программном комплексе Excel Microsoft Office

график распределения давления вдоль экрана, ограничивающего изотропный пласт (рис. 8), в котором работают две скважины (А – добывающая, В – нагнетательная). Скважины работают с переменным расходом жидкости q(t). Известны также параметры пласта и жидкости: толщина пласта h, проницаемость k,

Y, м

1000

800

А (доб.)

600

400

200

В (наг.)

0 200 400 600 800 1000 X, м

Рис. 8. Схема изотропного пласта

65

пористость m, коэффициенты объемного сжатия жидкости βж и среды βср, вязкость жидкости µ = 1 мПа с, начальное пластовое давление Рпл0 (экран проходит вдоль внешней границы координатной сетки сверху и справа).

Задание № 3.2 (исходные данные – в табл. П.3.4)

В изотропном пласте пробурены скважины 1, 2 и 3, работающие с переменным дебитом. Построить карты изобар на начало процесса нагнетания и на конец работы системы скважин в программном комплексе Surfer. Геолого-физические параметры пласта и его размеры принять по условию предыдущего задания.

Пример выполнения работы

Приведен пример выполнения задания № 3.1 для варианта № 1.

Разместим фиктивные скважины по площади залежи (рис. 9):

 

В′′

 

 

 

 

rВ′′

 

 

 

 

 

Y, м

 

rА′′

 

А′′

 

 

 

 

 

 

1000

 

 

rА

rА

 

ж

 

 

 

 

 

 

 

800

 

 

 

 

 

 

 

А

А

 

 

 

 

 

600

rВ

 

 

 

 

 

 

 

 

rВ

400

 

 

 

 

 

 

 

 

В

 

 

 

 

 

200

В (наг.)

 

 

 

 

 

 

 

 

 

0

200

400

600

800

X, м

Рис. 9. Схема размещения фиктивных скважин

66

При оценке распределения давления вдоль оси 0Х:

1. Определим коэффициент пьезопроводности по формуле

(124):

 

k

 

0,48 1012

2

 

χ =

 

=

 

= 0,8242 м

 

/с.

µн (mβж с )

0,0028 (0,14 7 1010 +1,1 1010 )

 

2.Определим расстояние между скважинами и точкой М

схарактеристиками (rA , rA, rA′′, rB , rB, rB′′):

 

 

r = X 2

+Y 2 .

 

 

 

 

 

i

i

 

i

 

 

Полученные результаты занесем в табл. 5.

Таблица 5

 

 

 

 

 

 

 

 

Х, м

 

т. А

 

 

 

т. В

 

rA , м

rA, м

 

rA′′, м

 

rB , м

rB, м

rB′′, м

 

 

 

0

824,62

1216,55

 

824,62

 

824,62

1969,77

824,62

100

728,01

1118,03

 

728,01

 

806,23

1878,83

806,23

200

632,46

1019,80

 

632,46

 

800,00

1788,85

800,00

300

538,52

921,95

 

538,52

 

806,23

1700,00

806,23

400

447,21

824,62

 

447,21

 

824,62

1612,45

824,62

500

360,56

728,01

 

360,56

 

854,40

1526,43

854,40

600

282,84

632,46

 

282,84

 

894,43

1442,22

894,43

700

223,61

538,52

 

223,61

 

943,40

1360,15

943,40

800

200,00

447,21

 

200,00

 

1000,00

1280,62

1000,00

900

223,61

360,56

 

223,61

 

1063,01

1204,16

1063,01

1000

282,84

282,84

 

282,84

 

1131,37

1131,37

1131,37

3. Определим изменение давления в точке М, исходя из принципа суперпозиции:

PM = Pi ,

где Рi изменение давления в точке М от реальных и отображенных скважин.

67

 

 

 

 

 

 

Q1aµ

 

 

 

 

 

2,246χt1a

 

 

 

 

(Q1a Q2a )µ

 

 

 

 

2,246χ(t1a

 

t2a )

 

 

PM =

 

 

 

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Q

 

 

Q

 

 

)µ

 

 

 

 

 

 

2,246χ(t

2a

t

 

)

 

 

Q

µ

 

 

 

 

 

 

2,246χt

 

 

 

 

 

 

+

 

 

 

 

 

 

 

2a

 

 

3a

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3a

 

 

 

 

 

 

 

 

1a

 

 

ln

 

 

 

 

 

 

 

 

 

1a

+

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Q

Q

 

 

 

)µ

 

 

 

 

 

 

2,246χ(t

t

2a

)

 

 

 

 

 

 

(Q

 

 

 

Q

 

)µ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

1a

 

 

 

 

2a

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1a

 

 

 

 

 

 

 

+

 

 

 

 

 

 

2a

 

 

 

 

 

3a

 

 

 

 

×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rA2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,246χ((t2a t3a ))

 

 

 

 

Q

 

µ

 

 

 

 

 

2,246χt

 

 

 

 

 

(Q

 

Q

 

 

)µ

 

 

× ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1a

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1a

 

+

 

 

 

 

 

 

 

1a

 

 

 

 

 

2a

 

 

 

 

×

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

r′′2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,246χ(t

 

 

t

 

 

 

)

 

 

 

 

 

 

(Q

 

 

 

Q

 

 

)µ

 

 

 

 

 

2,246χ(t

 

 

t

 

 

)

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

3a

 

 

 

 

 

 

× ln

 

 

 

 

 

 

 

 

 

1a

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

2a

 

 

 

3a

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r′′2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r′′2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

µ

 

 

 

2,246χt

 

 

 

 

 

 

 

 

(Q

 

 

Q

 

 

)

µ

 

 

 

 

 

 

2,246χ(t

 

 

t

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2b

 

 

 

 

 

 

 

 

 

1b

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1b

 

 

+

 

 

 

 

 

 

 

1b

 

 

 

2b

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

1b

 

 

 

 

 

 

+

 

 

 

4πkh

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Q

 

Q

 

 

)µ

 

 

 

 

 

 

2,246χ(t

2b

t

 

)

 

 

Q

µ

 

 

 

 

 

 

2,246χt

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

2b

 

 

3b

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3b

 

 

 

 

 

 

 

 

 

1b

 

 

ln

 

 

 

 

 

 

 

 

 

 

1b

+

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Q

 

Q

 

 

 

)µ

 

 

 

 

 

 

2,246χ(t

t

2b

)

 

 

 

 

 

 

(Q

 

 

 

Q

 

 

)µ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

1b

 

 

 

2b

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1b

 

 

 

 

 

 

 

+

 

 

 

 

 

 

2b

 

 

 

 

3b

 

 

 

×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rB2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,246χ(t

2b

 

t

3b

 

)

 

 

 

 

Q

 

 

µ

 

 

 

2,246χt

 

 

 

 

 

 

(Q

 

Q

 

)µ

 

 

 

× ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1b

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

1b

 

+

 

 

 

 

 

 

 

1b

 

 

 

 

2b

 

 

 

 

×

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r′′2

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,246χ(t

 

t

 

 

 

)

 

 

 

(Q

 

 

Q

 

 

 

 

)µ

 

 

 

 

 

2,246χ(t

 

 

t

 

)

 

 

 

 

 

 

 

 

 

 

2b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2b

 

3b

 

 

 

 

 

×

ln

 

 

 

 

 

 

 

 

 

 

1b

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

2b

 

 

 

3b

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

r′′2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πkh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r′′2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

68

4. Полученные расчетные данные заносятся в табл. 6.

 

Таблица 6

 

 

X

∆РМ, Па

0

200560,8

100

210378,6

200

220847,8

300

232064,9

400

244150,3

500

257256,8

600

271580,1

700

287369,2

800

304923,1

900

324509,2

1000

345933,0

5. Определяются значения давления в каждой точке границы экрана: Рi = P0 ∆Pi. Результаты заносятся в табл. 7, на основании чего строится график распределения давления.

 

Таблица 7

 

 

X

Р, МПа

0

17,80

100

17,79

200

17,78

300

17,77

400

17,76

500

17,74

600

17,73

700

17,71

800

17,70

900

17,68

1000

17,65

69

Распределение давления вдоль оси 0Y определяется аналогичным образом. По результатам строится совмещенный график распределения давления вдоль экрана.

Работа № 4

Подсчет запасов и оценка коэффициентов извлечения газа и конденсата

Различают начальные геологические и извлекаемые запасы газа. Извлекаемые запасы – это запасы газа, которые можно извлечь до достижения экономически рентабельного отбора газа из месторождения.

Геологические запасы газа объемным методом определяются по формуле (2).

Извлекаемые запасы газа отличаются от геологических на коэффициент извлечения газа, что есть отношение извлекаемого количества газа Qи к общему количеству газа в пласте до начала эксплуатации Q:

η =

Qи

=1

Qо

=1

Рк

 

zн

,

(125)

Q

 

 

Р

 

 

Q

z

0

 

 

 

 

 

 

 

 

 

 

н

 

 

где Qо – остаточный запас газа в пласте;

Рк – давление в пласте на момент окончания разработки. Приведенная формула справедлива для неизменного объе-

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

Задание № 4.1 (исходные данные – в табл. П.3.5) Определить изменение потенциальных запасов газа извест-

ного состава, наполняющего пласт площадью F, толщиной h, пористостью m при пластовом давлении Рпл, пластовой температуре Тпл, коэффициенте газонасыщенности порового простран-

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]