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

книги / Разработка и эксплуатация нефтяных месторождений

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

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

Из подземной гидрогазодинамики известен принцип электрогидродинамической аналогии (ЭГДА), согласно которому сила тока I соответствует расходу жидкости (дебиту Q), разность на­ пряжений ДU разности давлений (депрессии Ар), электрическое сопротивление проводника /?эл фильтрационному сопротивлению пласта /?ф. Принцип ЭГДА легко доказывается из анализа фор­ мул закона Дарси или Дюпюи и закона Ома:

<7=

пл Рз)

Др

дР

 

(2.38)

IxL

 

 

 

 

kSn

 

 

 

 

 

 

Др

 

2nkh (рПЛ— рз) _

Др

Лк

(2.39)

 

р In Лк

In

 

 

 

 

 

 

2nkh

 

 

 

 

I = MJ/RM,

 

 

(2.40)

где k — проницаемость пласта; S n, L — площадь поперечного се­ чения и длина полосообразного пласта; р — динамическая вяз­ кость жидкости; h, RK— толщина и радиус контура кругового пласта; гс— радиус скважины; Ар = р Плрз — депрессия давле­

ния (разность пластового p njl и забойного р 3 давлений); Яф=

— —^— L — фильтрационное сопротивление в полосообразном kSn

пласте; Лфг 2лkh In- Я* — фильтрационное сопротивление

в круговом пласте.

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

2nkh (рпл Рз)

 

 

Рпл — Рз

Др

Ц

\ о

+ In — —)

 

kh (2о)

In а/я

Q -f со

 

 

ягс /

 

2nkh

(2.41)

где

 

 

 

 

 

 

 

Й =

 

 

 

In а/я

 

 

 

kh (2а)

U

СО=

 

 

 

 

2яkh

 

Первое слагаемое й в знаменателе, как нетрудно заметить из сопоставления с формулой закона Дарси, равно фильтраци­ онному сопротивлению в полосообразном пласте на участке дли­ ной L от контура пласта до галереи, расположенной на линии ряда («галеризация» ряда). Площадь поперечного сечения пла­ ста, приходящегося на данную скважину из ряда, равна произ­ ведению толщины пласта h на ширину 2а, равную расстоянию между скважинами.

Второе слагаемое со равно фильтрационному сопротивлению в круговом пласте с радиусом контура а/п. Таким образом, сложный фильтрационный поток можно разбить на два про­ стейших: прямолинейнопараллельный поток от контура пласта до галереи, расположенной на линии ряда скважин; плоскора­ диальный поток внутри галереи в круговом пласте с длиной контура 2я/?к = 2а, т. е. Як=а/я. Величину й принято называть

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

Аналогично для кругового пласта дебит одной скважины в концентричном круговом ряду

<7= -

 

kh (Рпл — Рз)

 

 

ц (л In

Як

+

1п

Ri

 

 

 

 

 

 

 

 

Ri

 

 

пгс

 

 

---------------------- р з------------------- = —

(2.42)

ц

J

ЯК

,

Ц

1п

p/я

Q +

со

(2alRi) kh

"

tfi

kh

"

гс

 

 

где n = 2nR\l(2a)=nR\la — число скважин

в ряду;

R i — радиус

линии размещения кругового ряда скважин.

Первое слагаемое й в знаменателе представляет собой внеш­ нее фильтрационное сопротивление части кругового пласта (сек­ тора с углом 2a/R] радиан) от контура до круговой галереи дли­ ной 2а и радиусом /?1( а второе слагаемое ш — внутреннее филь­ трационное сопротивление притоку к скважине внутри галереи в круговом пласте с длиной контура 2nRK=2a, т. е. RK=a/n. В данном случае сложный поток к одной скважине в круговом ряду можно разбить на плоскорадиальный поток от контура до круговой галереи и плоскорадиальный поток к скважине внутри галереи.

Q = Y,qi = qn = n —^

АР

i=l

£2 +

со

(Q + со)

 

 

 

 

п

 

____________ Ар_

 

 

 

(2.43)

 

 

 

 

 

 

------ К------ L + —

■In ■

ягс

kh (2а) п

п

 

2nkh

 

 

Q = Z <7г = <7«

Ар

 

 

 

 

Ар

Q +

со

 

 

 

 

1=1

 

(Q + со)

 

 

 

 

п

 

 

 

'Ар

 

 

 

(2.44)

 

1

 

 

 

—К— 1п-^Г

 

■In

 

ЯГс

л

2я£Л

 

2JT£/I

 

 

Отсюда следует, что приток ко всем скважинам можно рас­ сматривать как параллельное соединение проводников с одина­ ковыми сопротивлениями (£2 + со). Таким образом, фильтраци­ онный поток к скважинам можно представлять эквивалентной схемой электрических сопротивлений и для расчета использо­ вать законы Ома и Кирхгофа (первый или второй закон), под­ разумевая в соответствии с принципом ЭГДА под силой тока, разностью напряжений и электрическими сопротивлениями их аналоги — расход жидкости, перепад давлении, фильтрационные

сопротивления.

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

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

фориои

Р

 

ной или круговой. Реал^

иХ

между

«проницаемыми»

гале-

рядов заменяется фильтрацией

между

 

 

 

 

леями с RHVTDPHHHMH Фильтрационными сопротивлениями сква-

реями с внутренними филыр^д

 

внеШние фильтрационные

жин внутри галереи, дополняют«ми

в

п

стФавляяР фильтра-

сопротивления между

галер

•_

 

^ электрической

схемой

ционную схему пласта

эквивалентно

 

законы

0ма и

_

сопротивлении и применяя

* " ° ^ ед

 

ии рядов скважин для

гофа, составляют уравнения

п]плений. Составим эти уравне-

расчета дебитов илизабойных Д

 

 

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

ния для кольцевого (круговог )

 

еНТричНыми

рядами сква-

и толщине пласта с круговыми ко ц

 

в

ой закон Кирхгофа,

жин (рис. 2.1). Для этого

испч г д^д

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

между

согласно которому на

основе, ^ в произведений дебита жидко

двумя точками схемы равен сумм в

 

 

сопротивление ЭТОГОэтого

сти в пределах участка на фильтрационное

I

 

 

63

g

Рис. 2.1. Схема кругового пласта (а)

и эквивалентная схема

сопротивле­

ний (б)

 

 

 

 

 

 

 

участка. Получим

систему

уравнений-

интерференции

(взаимо-

действия)

рядов скважин:

 

 

 

 

 

 

 

з

 

 

 

 

 

Рк

p3! =

X Qi 4" %Qi

 

 

 

 

 

 

1=1

 

 

 

 

 

 

 

3

3

 

 

 

(2.45)

Рк —Рзг = & 1 X Qi ”Ъ ^2 2

 

Qi ~Ь W2Q2

 

 

I=1

i=2

 

 

 

 

 

3

3

 

 

 

 

Рк— Рзз =

2 Qi +

й 2 X

 

Qi +

Q3Q3 + W3Q3

 

 

 

i=l

i=2

 

 

 

где рк — давление на контуре питания пласта; р3; — забойные давления скважин i-ro ряда; Q ,— дебит всех скважин i-горяда;

(о£ = —|— 2^kh^n ~

-----внутреннее фильтрационное сопротивле­

ние, одинаковое

для круговых

и

полосообразных пластов;

Qi = —^ I n ----внешнее

фильтрационное

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

/-го ряда; гс/, /г,-,

а — радиус

и

число скважин,

половина рас­

стояния между скважинами i-го

ряда

(oi = nRi/rii) ; Ri — радиус

/-го ряда скважин.

 

 

 

 

 

Залежь можно схематизировать сектором, тогда у выраже­ ния со/ и Q/ вместо 2л необходимо подставить величину угла сек­ тора в радианах. Число уравнений в системе должно равняться числу рядов скважин. Аналогичные уравнения интерференции составляется для полосообразной залежи, только внешнее

64

фильтрационное сопротивление и расстояние между скважи­

нами соответственно равны Q,= —^—

L,; 2а* = — , где sn — ши-

khsn

т

рина залежи; L;— расстояние между

соседними рядами или

между контуром питания и первым рядом.

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

основе системы (2.45) можно записать

 

Рк p3i — Л (&j 2

Q/1 + ®lQh

(2-46)

/—i V

/=i

/

 

где /, N — соответственно число

рядов, предшествующее i-му

ряду и общее число одновременно работающих рядов.

 

Уравнения интерференции еще составляют путем обхода

схемы сопротивлений от р3 j-{ до р3j, тогда

 

N

 

 

 

Рз/-1 —Рз/ = О/ £ Q. + Ю/Q/ —<0/_iQ/_x-

(2.47)

<=/

 

 

 

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

N

бит залежи) QCyM=EQi> дебиты скважин в рядах qi = Qi/ni.

1 = 1

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

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

3 В. С. БоЯко

65

а третьего— 15—20 % от дебита первого ряда или дебиты ря­ дов составляют соответственно 60—70, 30—20 и 5—10 % от сум­ марного отбора. Если в скважинах внешних рядов поддержи­ вать более высокие забойные давления, чем во внутренних ря­ дах, то дебиты внешних и внутренних рядов в значительной сте­ пени выравниваются, однако общий отбор из залежи умень­ шается.

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

аналогично системе (2.45) .для левой и правой частей схемы;

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

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

В заключение отметим, что рассмотрение дано примени­ тельно к фильтрации однородной (одинаковой плотности и вяз­ кости) несжимаемой жидкости в однородном по проницаемости и анизотропном плоском пласте к совершенным скважинам при одинаковых условиях (дебитах и забойных давлениях) работы скважин в пределах каждого ряда. Метод эквивалентных филь­ трационных сопротивлений можно применять при граничных по­ ложениях ВНК, т. е. когда в пласте движется только нефть или только вода (после полного обводнения). На использовании ме­ тода эквивалентных фильтрационных сопротивлений основаны методики ВНИИ-1, ТатНИПИнефть и др.

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

Рассмотрим процесс поршневого вытеснения нефти-водой из прямолинейного однородного пласта проницаемостью k, порис­ тостью т , толщиной h} шириной В и длиной LK. Начальное по­ ложение ВНК определяется координатой L0, а текущее в мо­ мент времени t — координатой x(t), где соответствующие дав­ ления составляют р0 и р. На пласт создан постоянный перепад давления Ар = Рк—Рг, где рк, рг — постоянные давления соответ­ ственно на контуре пласта и на галерее (остальные поверхности непроницаемые). Жидкости считаем несжимаемыми, взаимно нерастворимыми и химически не реагирующими одна с другой

ис пористой средой. Полагаем, что плоскость контакта нефти

иводы вертикальная. Это справедливо для случая либо пре­ дельно анизотропного пласта (проницаемость в вертикальном направлении равна нулю), либо равной плотности нефти и воды. Различны только вязкости нефти р,н и воды рв. В пласте выде­ ляются водяная, заводненная и нефтяная зоны. В первых двух движется вода, а в третьей — нефть. До начала вытеснения на­ сыщенность неподвижной связанной водой в нефтяной зоне со­ ставляла SCB- В заводненной зоне остаточная нефтенасыщенность остается постоянной и равной sOH, а связанная вода не­ подвижна и смешивается с закачиваемой водой. Тогда в силу несжимаемости и неразрывности потока скорости фильтрации

во всех трех зонах будут одинаковыми, т. е.

k (рк — р0) _

kkB(Рор) __

kkH(р — Рг)

(2.48)

Ив^о

M-ВLo)

Цн (^к х)

 

где kBi kH— относительные проницаемости для воды и нефти. Применяя правило производных пропорций, исключаем не­ известные давления ро, р и окончательно имеем выражение ско­

рости фильтрации

 

 

(РкРг)

(2.49)

а также расхода жидкости

k (Рк — Рг) Bh

. (2.50)

q = vBh

-------——

 

 

Отсюда следует, что скорость фильтрации и расход изменя­ ются с перемещением ВНК, т. е. во времени. Следовательно, не­ смотря на постоянство перепада давления Ар движение жидко­ сти неустановившееся. При (хн>Цв точнее Ци/кн>рв/кв, скорость v и расход q увеличиваются во времени. Это объясняется умень­ шением знаменателя (в общем фильтрационного сопротив­ ления).

Допустим, что положение ВНК не параллельно галерее (ис­ кривлено). Из формул (2.49) и (2.50) следует, что чем больше Lo, тем больше о и q. Значит, в тех сечениях, где LQ больше или граница раздела ближе к галерее, будет происходить опере­ жающее перемещение ВНК и дальнейшее искривление линии раздела. Отсюда приходим к выводу, что если на границе раз­ дела образовался «язык обводнения», то в дальнейшем он не только не исчезает, но еще больше вытягивается, продвигаясь с большей скоростью. Искривленное, вернее горизонтальное по­ ложение ВНК по отношению к галерее, отмечается в наклонных пластах, что приводит к более быстрому обводнению галереи по подошве пласта. В реальных условиях неизбежны возмуще­ ния на границе раздела (например, изменение проницаемости) и образование «языков обводнения», т. е. проявляется вязкост­ ная неустойчивость вытеснения. Если движение образовавшихся «языков обводнения» замедляется, то такое перемещение гра­ ницы раздела называют устойчивым. Как известно из подзем­ ной гидрогазодинамики, оно возможно в первом приближении при условии, что скорость фильтрации нефти на границе раз­ дела

— — — (рв— pH)g s in а,,,

(2.51)

Ри —Ив

 

 

где рв, рн — плотность соответственно

нефти и воды;

g — уско­

рение свободного падения; ан — угол

наклона пласта к гори­

зонту.

 

 

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

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

Определим закон движения границы раздела x ^ x (t) . Не на­ рушая общности рассуждений, с методических позиций в даль­ нейшем примем L0 = 0. Это соответствует случаю, например, бло­

кового заводнения. Тогда формулы скорости фильтрации и рас­ хода упростятся и примут вид:

v

 

к (рк — Рг)

 

;

 

(2.52)

Рн

т

(

Рн

 

Рв \

 

 

 

 

 

 

 

К

V kH

 

kB J X

 

 

 

п

 

к (рк

Рг) Bh

 

 

(2.53)

Ч

Рн

J

(

Рн

 

Рв Ч

 

 

 

 

 

 

 

 

кн

к

V

 

кв ) Х

 

 

 

соотношения скорости фильтрации v и средней скорости

движения w находим

 

 

 

 

 

 

 

 

 

V tTl {\

SCBSQII) ^ — tTl{ \ SCB SQH) -----»

(2.54)

 

 

 

 

 

 

 

dt

 

 

откуда

 

 

 

 

 

 

 

 

 

dt-- m (1 — 5CB

SOH)

dx

tTt {\ SCB

$он)

 

 

 

 

 

 

 

 

k (Рк— Pr)

 

 

 

 

 

 

 

 

 

 

 

( 2 ' 5 5 )

Интегрируя уравнение

(2.55)

в пределах от 0 до / и от 0 до

х, имеем

 

 

 

 

 

 

 

 

 

t = _ m ( l — sCB- s 0H)

Г_Рн_

1 лХ _ Г

_MIH---------/ 2 .5 б )

г 2

 

k (Рк-- Рг)

L

 

 

\

/

2

J

откуда

 

 

 

 

 

 

 

 

 

X= -

Рн^к

 

 

[ 1 —V 1—

]>

 

(2.57)

 

\ «н

 

«в /

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

а —

 

 

О

 

 

 

 

 

 

 

 

 

S0 H)

 

 

 

 

 

 

 

- S CB -

 

 

 

Перед квадратным корнем принят знак минут из физических соображений, что легко установить при t = 0, когда х = 0. Из уравнения (2.56) при х = LKполучаем формулу для определения времени полного обводнения пласта (мгновенного обводнения продукции пласта)

/обв:

’ ( ^

SCB S0H )

(2.58)

2k

(Ркр — Ррг)

( - C - + - S - ) -

 

Изложенное можно распространить на слоисто-неоднород­ ный пласт, полагая, что рассмотрен один из пропластков такого пласта. Допустим, что слоистый пласт состоит из множества изолированных пропластков (не сообщающихся между собой вдоль пласта). Мысленно сложим их в «штабель», начиная с пропластка с наибольшей проницаемостью. Тогда в соответ­ ствии с формулой любого закона распределения проницаемости суммарную толщину hKпропластков, проницаемость каждого из которых не меньше значения k, можно записать

h jh = F (k),

(2.59)

где h — общая толщина всех пропластков р «штабеле»; F(k) — закон (или функция) вероятностно-статистического распределе­ ния проницаемости.

Дифференцируя уравнение (2.59), имеем

- ^ - = F'{k)dk = f(k)dk,

(2.60)

h

где f(k)= F '(k) — плотность вероятностно-статистического рас­ пределения проницаемости.

Считаем, что расход жидкости dq через слой с проницае­ мостью k и толщиной dhK в соответствии с формулой (2.53) можно записать

dq =

ft (Рк — Рг) BdhK

(2.61)

С учетом уравнений

(2.57) и (2.60) окончательно

находим

dq _

k (рк — Рг) BkHhf (k) dk

(2.62)

 

рн£к Vl — akt

 

Принимаем, что к моменту времени t — t0бв обводнились слои с проницаемостью k>-k0tB- Согласно теории поршневого вытес­ нения из них поступает только вода. Из слоев с проницаемостью k< ko6B пока еще добываем нефть. Тогда для расхода нефти из

слоистого пласта на основе формулы

(2.62) можно записать

 

*обв

kf(k)dk

 

Я»(0 =

(Рк — рг) Bhkн Г

(2.63)

о

V l — akt

 

 

 

 

 

Для определения расхода воды через обводнявшийся слой необходимо в формулу (2.62) вместо t подставить t0бВ. Тогда

Соседние файлы в папке книги