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

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

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

области 2 должна быть больше, чем скорость движения грани­ цы при X =XQJ, умноженная на т вследствие перетока нефти из области 1 в 2, но меньше скорости фильтрации на границе х= = х2, примерно равно 2тйх0т/с11. Таким образом, приближенно устанавливают s2 и х2.

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

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

Я = <7пл-г 2<7тЬД*т-

(VII.26)

Здесь q — количество тепла,

вводимого в пласт в единицу вре­

мени вместе с паром; q

— изменение за единицу времени теп­

ла в нагретой области 1

(рис. 131); qT— изменение за единицу

времени тепла, отдаваемого в кровлю — подошву. В расчетной схеме Маркса — Лангенгейма использована схема теплопотерь

хт

Рис. 131. Схема распределения темпе­ ратуры в пласте согласно модели Маркса — Лангенгейма:

1 — нагретая область; 2 — область с пла­ стовой температурой

251

«Поверье. В области, содержащей насыщенный пар и остаточ­ ную нефть с насыщенностью SHOCT, температура равна темпера­ туре То нагнетаемого пара. В области 2, расположенной перед областью 1, температура равна пластовой Тпл.

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

Д ^ = = -;

АТ0 — Т„Та .

(VII.27)

~У я х тк(/ — т)

Для нагретой области 1 имеем

Яил~ cbh&T0

;

 

(VII.28)

С [СТРТ(1

(1 SH 0ст) (Свр в -{-С пр п) -|-/72CHp HSH ост].

Подставляя (VII.27) и (VII.28) в уравнение баланса тепла (VII.26) и переходя к пределу At—>-0, Алгт— Я), получим

илт

dx

^

д = ~ с Ш Т ^

(VII.29)

~[/лхтк {I— т)

Так как здесь искомая величина dxT/dt находится под зна­ ком интеграла, уравнение (VII.29) интегральное. Решение это­ го уравнения получаем с использованием преобразования «Лап­ ласа. Оно имеет следующий вид:

T (</) = e>'erfc(J/V2)+ 2 f4 - Y /2- l ;

=

erfc(;/V3) = i _ ( V I I . 3 0 )

 

6

Подставляя время t в последнюю формулу, находим соответ­ ствующее ему значение у, по у определяем ц>(у) и затем по пер­ вой формуле (VII.30) вычисляем хт.

Скорость теплового фронта wT = dxTJdt получаем дифферен­ цированием первого выражения (VII.30):

wT =

— еУ erfc 112).

(VII.31)

cbhkTо

252

Важным показателем процесса закачки в пласт теплоноси­ телей является ,Пт ^ >к~оэ ф ф и п и е н т т е п л о в снй э ф-ф-е-тц т и в н о с т и_процесса, определяемый следующим образом:

т^= IzjL = e^ * -0- T- = еу erfc (yV*)^7

(VII.32>

На рис. 132 показана зависимость т]т= г|т(у), из которой сле­ дует, что с ростом безразмерного времени у коэффициент теп­ ловой эффективности процесса воздействия на пласт путем за­ качки в него пара уменьшается, поскольку с течением времени все большее количество тепла будет уходить в кровлю и по­ дошву пласта.

Рассмотренную схему теплопереноса в пласте при закачке- в него пара можно также использовать и в случае радиальной фильтрации. Тогда вместо первого уравнения — см. выражение (VII.30) — будем иметь

5Т= 2ХаткАГо Ф (У)\ = лгт (VII.33).

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

П р и м е р VII.1. Характеристика прямолинейного пласта та же, что и

при расчете распределения температуры по

формуле

(VI 1.20). В этот пласт

закачивается

горячая вода с температурой Т\=203,2 К и расходом q=

= 150 м3/сут.

Пластовая температура Тпл =

303,2 К;

насыщенность пористой

среды пласта связанной водой sCB= 0,05. Полная толщина пласта h0= 20 м; тол­ щина, охваченная вытеснением, h= 15 м (коэффициент охвата пласта воздей­ ствием г|2=0,75). В язкость нефти в-пластовых условиях |х„=40-10-3 Па-с, вязкость воды р,в= 10 -3 Па-с.

Экспериментально определено, что остаточная нефтенасыщенность sHос-г

при закачке в пласт горячей воды зависит от перепада температур ДТ=7’—Тпл. следующим образом:

SH ост

л

—2,554- 10-ЗДГ

(VI 1.34)'

0,75е

Определим текущую нефтеотдачу при вытеснении нефти горячей водой в- момент времени, когда *от= /= 1 0 0 м; текущую нефтеотдачу при вытеснении

нефти водой с начальной пластовой температурой; затраты условного количе­ ства нефти на производство горячей воды, если общий к. п. д. системы водо­ грейная установка — водопровод — скважина составляет 60%. теплотворная: способность нефти составляет 38 -106 кДж/м3.

ост

Рис. 133. Зависимость sHост от х:

1 — при непрерывной закачке горячей воды, когда х от=1; 2 — при использо­

вании тепловой оторочки и длительно­ сти процесса 7=500 сут; 3 — при непре­ рывной закачке горячей воды, 7=500 сут

Прежде всего построим зависимость sHО С Т --ост(х) С ПОМОЩЬЮ формул

(VII.20) и (VII.34). Эта зависимость показана на рис. 133.

Начальные запасы нефти в области, охваченной вытеснением горячей во­

дой:

€ охв = тЪМ(1 — sCB) = 0,2-100-15-100 (1 — 0,05) = 28,5 -103 м3.

Остаточные запасы нефти в области, охваченной вытеснением в момент вре­ мени, когда хТ=1, т. е. при /=302,8 сут,

I

'Сост == ffibh

SHост (я) dx.

 

 

о

 

 

 

 

 

Значение

G0cT определяем по графику (см. рис. 133). Имеем G0ст=

= 18,47-103 м3.

 

 

 

 

Коэффициент вытеснения

 

 

G Q X B - G

Q C T

_

( 2 8 ,5 - 18,47)-1Q3

_

’Tll“

G0XB

 

28,5 -103

— и’Л* -

 

Коэффициент нефтеотдачи

 

т) =

)2 = 0,352-0,75 - - 0,264.

 

 

Если нефть вытесняется водой при пластовой температуре, то коэффи­

циент вытеснения составит

 

 

Gnvn

'Gn

 

0 ,9 5 — 0,75

0, 21.

«lot =

 

 

0,95

 

Коэффициент нефтеотдачи

 

% = 0 ,2 1 -0 ,7 5 =

0,1575.

 

Количество накопленной дополнительно добытой нефти за счет горячего заводнения

AQH= G0 ост— GOCT = 22,5-103— 18,47-103 » 4 -103 м3.

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

<)в= 150-320,8= 48,12-103 м3.

При Д7[=200 К на нагрев этого объема воды потребуется следующее количе­ ство тепла:

<>т = 4,19 -103-200*48,12-103= 40,27-109 кДж.

Г-

Это тепло эквивалентно'условному сжиганию QH3 нефти:

40,27-10е

•<2нэ = —зз. 2Q6— = Ю60 м3 нефти.

Под условным сжиганием нефти понимают расходование эквивалентного количества энергии на нагрев воды. С учетом к. п. д., равного 0,6, нужно сжечь

<QH3 = Q Q = 1770 м3 нефти.

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

254

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

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

§ 3. РАЗРАБОТКА МЕСТОРОЖДЕНИЙ ПУТЕМ ЗАКАЧКИ ТЕПЛОНОСИТЕЛЕЙ В ПЛАСТ МЕТОДОМ ТЕПЛОВЫХ ОТОРОЧЕК

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

вой

о т о р о ч к и .

Способ

перемещения нагретой

области в-,

глубь

пласта путем

закачки

в него холодной воды,

т. е. воды

с температурой, близкой к пластовой, был предложен в 50-х гг., но только в 60-х гг. по экспериментальным и теоретическим данным обосновали метод тепловых оторочек как способ раз­ работки нефтяных месторождений. Были разработаны методики выбора оптимальных размеров тепловых оторочек при различ­ ных геолого-физических условиях пластов, темпах нагнетания в пласт теплоносителей, их параметрах и других технологичес­ ких показателях разработки месторождений.

Использование тепловых оторочек позволяет получить не­ сколько меньшую нефтеотдачу по сравнению с этим показате­ лем при непрерывной закачке теплоносителей в пласт. Но в та­ ком случае на подготовку горячей воды или пара значительно­ меньше тратится энергии.

Если рассматривать отношение дополнительно извлеченной нефти AQH, получаемой при использовании метода тепловой ото­ рочки, к затрате тепла QT на нагрев теплоносителя, то опти­ мальные размеры оторочки и другие показатели теплового воз­ действия достигаются при условии

TIt0 = AQH/QT ---->- max. (VII.35>

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

Рассмотрим распределение температуры в прямолинейном пласте при создании в нем тепловой оторочки за счет закачки горячей воды, основываясь на решении (VII.20). Вначале за­ качивают в пласт горячую воду с начальной температурой Т = = Т1 и AT=ATi. В момент времени t = t* температура этой воды снижается скачком до Т = Тпл или становится ДГ = 0 при х = 0.

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

25&

ния (VII.20) вычесть такое же решение, но зависящее не от t, а от tt* (U — момент начала закачки в пласт воды с темпера­ турой, равной пластовой).

В результате для определения распределения перепада тем­ пературы ДT(x,t) в пласте с тепловой оторочкой получим сле­ дующую формулу:

“КткХ

AT (х, t) = A7\ erfc

 

(VII.36)

Первый член в формуле

(VI 1.36) справедлив при t>bx]a, а вто­

рой— при tt*>bx!a.

Входящие в формулу (VII.36) обозна­

чения те же, что и в предыдущих параграфах.

Как следует из (VII.36), максимальная температура в плас­ те достигается при х = хтах, причем

х шах

a( t — t,)

(VII.37)

Ъ

Рассмотрим пример вытеснения нефти водой методом тепло­

вой оторочки.

 

П р и м е р VII.2. Пусть

имеем тот же прямолинейный пласт с теми же

размерами и свойствами, что и в примере VII. 1. Закачка горячей воды в пласт

при

A7'i = 200 К ведется с

тем же расходом <7=150 м3/сут. Однако через

7*=200 сут после начала процесса вытеснения нефти переходят на закачку

холодной воды с ДГ=0. Будем считать, что процесс разработки рассматривае­

мого пласта

ведется

в течение /=500 сут. К. п. д.

системы

водогрейная уста­

новка— водопровод — скважина

составляет 0,6.

Теплота

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

сжигаемой в

котлах

нефти для

производства

горячей

воды составляет

3,8/107 кДж/т. Определим распределение температуры в пласте в различные моменты времени, а также остаточную нефтенасыщенность sHо с т , нефтеотдачу

и «чистую» дополнительно полученную нефть по сравнению с этими парамет­ рами при вытеснении нефти из пласта холодной водой к концу разработки пласта, т. е. при /=500 сут, как в случае использования тепловой оторочки, так и во время непрерывной закачки в пласт горячей воды.

Определим по формуле (VI 1.36) распределение температуры в пласте в различные моменты времени /. На рис. 134 показано распределение темпера­

туры при вытеснении

нефти горячей водой методом

тепловой оторочки при

/ = 300, 400 и 500 сут

с начала закачки воды. Видим,

что температура в теп­

ловой оторочке существенно снижается с течением времени. Кривая 4 харак­ теризует распределение температуры в пласте при непрерывной закачке в него горячей воды в течение 500 сут.

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

формуле (VII.34), будем считать, что из охлаждаемой области тепловой ото­ рочки, т. е. из области x<.vmax (xmax приближенно вычисляют по формуле (VI 1.37)) нефть уже не вытесняется. Вычисленное по формуле (VII.34), с уче­ том указанного положения, распределение насыщенности в случае применения тепловой оторочки соответствует кривой 2 (см. рис. 133), а в случае непрерыв-

256

Из уравнения (VIIIЛ) видим, что для обеспечения опреде­ ленной скорости поступления запасов из разведки в разработку необходимо заблаговременно направлять соответствующие ре­ сурсы в подготовку запасов более низких категорий. Коэффи­ циенты подтверждаемости запасов обычно сильно колеблются, но по каждому конкретному региону или стране в целом можно получить статистические закономерности изменения этих коэф­ фициентов с тем, чтобы рассчитывать скорости подготовки за­ пасов различных категорий. Поддержание величин G на опре­ деленных уровнях необходимо для обеспечения надежности за­ планированного перевода запасов. Если перевод запасов из какой-либо более низкой категории окажется недостаточным вследствие внезапного снижения коэффициента их подтверж­ даемости, можно обеспечить перевод данных запасов в более высокую категорию за счет временного снижения величины Gt..

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

В нефтегазоносном регионе или в стране в целом разведан­ ные нефтяные месторождения переводят в разработку постепен­

но, так что в момент

времени /= 0 вводят месторождение ЛИ®,

в ti — месторождение

Мi, в t2 — месторождение М2 и т. д.

Извлекаемые запасы месторождения М0 обозначим через Л®, а

темп его разработки — через Zo(t).

Следовательно, извлекаемые

запасы месторождения M-t будут

Ni, а темп его разработки

Zi(tti) (ti — время вступления i-ro месторождения в

разра­

ботку) .

 

 

 

Отсюда, помня определение темпа разработки месторожде­

ния, можно написать выражение

для добычи

нефтп

шз

месторождений рассматриваемого

региона или

страны.

Имеем

=

(VIlUg

 

О

Отметим следующее. Если рассматривать совместно развед­ ку нефтяных месторождений и их разработку как некоторую подсистему, то развитие такой подсистемы характеризуется ©б- щим оптимумом. Этот оптимум можно определять осреджемжо* считая свойства всех нефтяных месторождений одинаковыми* ш дифференцированно с учетом свойств отдельных месторожде­ ний. Будем в дальнейшем называть осредненньан оптимум об ­ щ им о п т и м у м о м п е р в о г о , а дифференцированный* — о б щ и м о п т и м у м о м в т о р о г о р о д а .

При заданных ресурсах (труб, буровых установок, машин ш т. д.), выделяемых народным хозяйством нефтяной промышлен­ ности, можно развивать подсистему разведка — разработка раз­ личными путями. Например, направить подавляющую часть ре­ сурсов в разработку, уплотняя сетку скважин на месторожде­ ниях и обеспечивая высокий темп их добычи. Однако разведка

19 ю . П. Желтое

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

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

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

сурсов

(установок

подготовки нефти, насосного оборудования

и др.)

в разведке

месторождений либо не потребляют, либа

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

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

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

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

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

290

пасы нефти ДЯ(т). Тогда, на основе формулы (VIII.2), получим

< u o = 2 A;v (т) г ((—■т) = 2 ^ г г(*— т )д т =

^

о

о

А'1

t

 

 

(VIII.3)

= CJjZ-z(i~x)dt.

 

 

о

 

 

 

Пусть M(t) — общий

пробуренный метраж в разведочном

и эксплуатационном бурении. Если

Mp(t) — метраж разведоч­

ного, а M3(t) — эксплуатационного бурения, то

M(t) = Mp(t)+M3(t).

 

 

(VIII.4)

Если заданы выделяемые нефтяной промышленности на плани­ руемый период ресурсы, то заданы и общая длина обсадных труб (в метрах) и общий объем других материалов и оборудо­ вания. Тогда

=

<vlI1-5>

где f(t) — заданная функция.

 

Далее примем, что скорость поступления запасов нефти из разведки в разработку g (0 пропорциональна произведению эф­ фективности разведочного бурения Ер на скорость dMPldt при­ ращения числа метров породы, пробуренных в разведочном бу­ рении, т. е.

g(t) = Ep- ^ .

(VIII.6)

Эффективность разведочного

бурения Ер, выражающаяся

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

Скорость ввода разведанных извлекаемых запасов в разра­ ботку обозначим через g3(t) =dN/dt. Эту скорость можно выра­ зить через параметр А. П. Крылова NKP (извлекаемые запасы, приходящиеся на скважину), скорость изменения длины экс­

плуатационных скважин

(в метрах) dM3/dt и среднюю глубину

скважины Я. Имеем

 

f t (О

Л^кр

dM3

(VIII.7)

н

dt

В некоторых случаях разведанные запасы нефти не могут быть введены немедленно в разработку. Тогда появляются за­ пасы Nop(t), ожидающие ввода в разработку, изменение кото­

19*

291

рых со временем происходит в соответствии со следующим уравпением:

- ^

= g(*)-g,W -

(VIII.8)

Если dNop/dt = 0, то g (t)= g 9(t) и из (VIII.6) и (VIII.7) по­

лучим

 

или

 

 

 

Н Е р т N кр

(VHI.9)

 

 

Можно выразить NKP через

параметр плотности сетки скважин

Sc. Имеем

 

Л'кр = /^ н о Т н ^ к (5 с)>

(VIII. 1 0)

где

h, т, 5Н0„ ун — средние

по вводимым в разработку место­

рождениям значения соответственно толщины пласта, пористо­ сти, нефтенасыщенности, удельного веса; т]к — конечная нефте­ отдача по месторождениям, зависящая от параметра плотности сетки скважин 5С. Если обозначить

CLa= htnSBOyH1

то из (VIII.10)

имеем

V„p = aBS0ri„(Sc).

(VIII.11)

Введем величину

£ rP= ^

.

(VIII. 12)

Здесь Егр— эффективность разведочного бурения, равная коли­ честву геологических запасов нефти (в тоннах), приходящихся на 1 м проходки разведочного бурения.

Из (VIII.9), (VIII.11) и (VIII.12) получим

g. ^ __

Дп^сЛк (Sc) £ г р / (0

(VIII 13)

 

НЕгр -{■onSc

 

Для

определения зависимости Т1к = г]к(5с)

можно использо­

вать формулу В. Н. Щелкачева

 

’U S c) = % e -c‘\

(VIII.14)

где гр — коэффициент вытеснения.

 

Подставляя (VIII.14) в (VIII.13), имеем

 

/^\ _

е с f (0

(VIII. 15)

 

 

НЕгр -f- anSc

292

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