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

книги / Основы разработки нефтяных и газовых месторождений

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

аг

х а

аз

при данном объеме закачанной воды и при 5и > 5^, можно выразить уравнение (10.23) в виде

а?

(1-5 )—- + Г $ а ( —

ог а$

§ =

(10.24)

 

а?

 

___ \у_

 

аз

Применяя интегрирование по частям, записываем

 

/ иау = иу - 1 уаи,

 

 

и получаем

 

 

аг 1 и*

 

 

$

а [— 1 =

5

_

{

азиг ^1-5

"

уаз\ X V /I

^ "

 

1 - $

 

 

 

 

ог

XV

 

 

 

 

ОГ

Подстановка этого выражения в уравнение (10.24) дает после упро­ щения

:$*г +

а

(10.25)

аз

Здесь и функция и ее производная определяются для насыщен­ ности на фронте 8^ И, наконец, приравниваем уравнения (10.22) и

(10.25) и получаем

 

 

 

<к,

( 1 - П )

1

(10.26)

 

____ ш1УГ- = ___

азиг

Значимость этого результата продемонстрирована на рис. 10.13. Чтобы выполнялось равенство (10.26), прямая, проведенная из

точки с координатами (5^ = 5^с; = 0), должна касаться кривой доли

воды в потоке в точке с координатами (3^ = 5 ^

), а экстра­

полированная касательная должна пересекать линию

= 1 в точке с

координатами (5^ - 5^; Г = 1).

 

Рис* 10*13* Касательная к кривой доли воды в потоке, проведенная из

точки 5 = 5

УУ УУС

Этот метод определения 5^, ^ |5 г и 5штребует построения кривой доли воды в потоке с помощью уравнения (10.12) или (10.21) во всем интервале значений водонасыщенности

Как отмечалось выше, при использовании и того и другого уравне­ ния игнорируется влияние градиента капиллярного давления ЭРс / Эх. Однако пренебрегать градиентом капиллярного давления допустимо лишь в области за фронтом вытеснения, то есть в интервале

Поэтому участок кривой доли воды в потоке для значений насы­ щенности менее 8 ^ является мнимым. Первая действительная точка на кривой имеет координаты 8^г, {^\$ , что соответствует насыщен­ ности на фронте. Этот простой графический метод Уэлджа нашел очень широкое применение в расчетах добычи, которые будут рас­ смотрены в следующем разделе.

10.5. РАСЧЕТ ДОБЫЧИ НЕФТИ

Для определения положения плоскостей с постоянной водонасыщенностью и, следовательно, профиля водонасыщенности в ин-

тервале < 5^ < 1 -5ог до прорыва воды в добывающую скважину можно применить уравнение (10.20). В момент прорыва и далее это уравнение используется несколько иначе, для контроля возрастания водонасыщенности около добывающей скважины. В этом случае х = 1 (Ь - длина элемента пласта). Таким образом, здесь х - постоянная величина, и уравнение (10.20) можно записать так:

1

 

__ 1

= \У.и*

(10.27)

Ь Аср

 

агЦТ

 

 

а$м

 

Здесь 8ше - текущее значение водонасыщенности около добывающей скважины (рис. 10.14), а \\7а - безразмерная накопленная закачка, выраженная в долях порового объема (один поровый объем равен произведению ЬАср).

До момента прорыва расчет добычи является тривиальным. На этом этапе вода не отбирается, и поэтому, при участии в процессе вытеснения несжимаемых флюидов, объем добытой нефти равен объему закачанной воды. В момент прорыва плоскость с насыщен­ ностью на фронте 8 , = 8 подходит к добывающей скважине, и доля воды в потоке резко возрастает от нуля до л с* |I $^. Это явление ча сто наблюдается на промыслах, и оно подтверждает существование скачка насыщенности. Для этого момента уравнение (10.22) можно записать с теми же параметрами, что и уравнение (10.27), получив

ЧаЬ( = ™и,к, = Я,А, = (5„к, - $«) =

(10.28)

1аы

й{

 

____ XV

 

аз

Здесь все объемы выражены, для удобства, в безразмерной форме, через объем порового пространства. В частности, безразмерный рас­ ход закачки равен / (ЬАф) (объем порового пространства, делен­ ный на единицу времени). Такой подход упрощает расчет момента,

когда происходит прорыв

 

 

Ч*

ы

(10.29)

 

Ям

Рис* 10*14* Распределение водонасыщенносги в момент прорыва и

далее при линейном заводнении

После прорыва параметр Ь в уравнении (10.27) остается посто­ янным, а водонасыщенность 5^е и доля воды в потоке около добы­ вающей скважины постепенно возрастают по мере продвижения воды по пласту (см. рис. 10.14). На этом этапе расчет добычи несколь­ ко усложняется, поскольку требуется применять уравнение Уэлджа (10.25) в виде

5 = 5

у/е

+ ( ! - ( )

(10.30)

4 «и»'

 

__^

а$V*$

Используя уравнение (10.27), его можно записать также в виде

5 = 5

+(1- { )Ж,.

(10.31)

'

^

\уе у ше' к!

4

И, наконец, вычитая 5^с из правой и левой частей уравнения (10.31), получаем зависимость для расчета добычи

N = § - 5

= (5 -5

) + (1 - Г )\У.,.

(10.32)

рс1 ш

^ с у ш е

уме*

4

ш еу М

'

^

Ниже описано возможное практическое применение уравнений (10.28) и (10.32).

а) Построить кривую доли воды в потоке по уравнению (10.12) или (10.21). При необходимости учесть действие силы тяжести, но пренебречь градиентом капиллярного давления ЭРс / Эх.

B) Провести касательную к этой кривой из точки (8^ = 8^с, ^ = 0). Как указывалось в предыдущем разделе, точка касания имеет коорди­ наты 8 = 8 , = 8 \ ={ , а экстраполяция этой прямой до {„= 1 дает значение средней насыщенности за фронтом в момент прорыва 5„ = 8 . Теперь можно рассчитать добычу и момент про­ рыва по уравнениям (10.28) и (10.29).

c)Выбрать 5^е в качестве независимой переменной и увеличивать ее значение приращениями, допустим, по 5 % относительно на­

сыщенности при прорыве. Каждая точка на кривой доли воды в потоке имеет в области 8 >8 координаты (8 =8 , { = { ). Построить по уравнению (10.30) касательную к кривой доли воды в потоке, точка пересечения которой с линией { = 1 даст текущее значение средней водонасыщенности 8^ в элементе пласта (см. рис. 10.15).

Значения 8^, соответствующие каждому новому значению 8^е, определяются графически. Добычу рассчитывают по формуле

]Мр<1= 5^ - 8ш.

(поровый объем)

Величина, обратная тангенсу угла наклона касательной к кривой доли воды в потоке, для каждого значения 8^е даст накопленную за­ качку, выраженную в поровых объемах (\У ), уравнение (10.27). Это обстоятельство позволяет проследить динамику добычи, поскольку

Применяя другой подход, можно использовать уравнение 10.32 непосредственно для расчета добычи, определив { и \Уа по кривой доли воды в потоке для каждого выбранного значения 8^е. Этот под­ ход рассмотрен в упражнении (10.2), где рассчитываются Ыра и \Уа.

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

Рис* 10*15* Применение практического метода Уэлджа для расчета

добычи после прорыва воды

ности по толщине в вертикальном направлении. В качестве примера можно привести фильтрацию с учетом гравитационной сегрегации и вытеснение в слоисто-неоднородных пластах. Тем не менее путем усреднения насыщенностей и зависящих от насыщенности относи­ тельных фазовых проницаемостей в направлении, нормальном к по­ току, можно свести большинство двумерных задач к одномерным. Можно построить кривую доли воды в потоке, используя вместо лабораторных кривых проницаемости породы кривые усредненных относительных фазовых проницаемостей, и рассчитать добычу по методу Бакли-Леверетта / Уэлджа. Эти замечания помогут понять материал, приведенный в остальной части главы.

УПРАЖНКНИВ10.1. РАСЧЕТ ДОЛИ ВОДЫ В ПОТОКЕ

Вытеснение нефти водой происходит из горизонтального прямо­ линейного пласта при равномерном распределении насыщенности. Значения относительных фазовых проницаемостей породы для воды и нефти приведены в табл. 10.1.

Давление поддерживается на первоначальном уровне. При таком давлении В = 1,3 пл. м3 / ст. м3 и В = 1,0 пл. м3 / ст. м3.

Требуется сравнить обводненность продукции (в поверхностных условиях) и накопленную добычу нефти на момент прорыва при сле­ дующих значениях параметров флюидов:

Случай

Вязкость нефти

Вязкость воды

1

50 мПа с

0,5 мПа с

2

5 мПа с

0,5 мПа с

3

0,4 мПа с

1,0 мПас

Принимается, что значения относительных фазовых проницаемо­ стей и данные РУТ действительны для всех трех случаев.

УПРАЖНЕНИЕ 10.1. РЕШЕНИЕ

1) При горизонтальном течении доля воды в потоке составляет

( 10.12)

{ =

1

 

И

к

1 +

 

к

т

и

 

 

«о

При этом обводненность продукции на поверхности ^ равна

^ _

а

_____*ИГ

АУ_____

 

а*иг

+ Тао / Во '

Здесь расходы выражены в пл. м3 / сут. Выполнив подстановку, полу­ чаем выражение для обводненности продукции на поверхности

{ =

1

(10.33)

.

В

/1

 

1 + — м| - н

В

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

8иг

кГУГ

кго

8иг

кги

кго

0,20

0

0,800

0,50

0,075

0,163

0,25

0,02

0,610

0,55

0,100

0,120

0,30

0,09

0,470

0,60

0,132

0,081

0,35

0,020

0,370

0,65

0,170

0,050

0,40

0,033

0,285

0,70

0,208

0,027

0,45

0,051

0,220

0,75

0,251

0,010

 

 

 

0,80

0,300

0

Таблица 10.1

 

 

 

 

 

Доля воды в потоке ( Е )

 

 

 

 

 

случай 1

случай 2

случай 3

8V

к

к

го

к0/к„

ци/цо =°,°1

МИ„ = 0>1

И* / И0 ” 2>25

/

 

 

 

 

 

0,2

0

0,800

сю

0

0

0

0,25

0,002

0,610

305,000

0,247

0,032

0,001

0,30

0,009

0,470

52,222

0,657

0,161

0,008

0,35

0,020

0,370

18,500

0,844

0,351

0,021

0,40

0,033

0,285

8,636

0,921

0,537

0,044

0,45

0,051

0,220

4,314

0,959

0,699

0,085

0,50

0,075

0,163

2,173

0,979

0,821

0,155

0,55

0,100

0,120

1,200

0,988

0,893

0,250

0,60

0,132

0,081

0,614

0,994

0,942

0,394

0,65

0,170

0,050

0,294

0,997

0,971

0,576

0,70

0,208

0,027

0,130

0,999

0,987

0,755

0,75

0,251

0,010

0,040

0,999

0,996

0,909

0,80

0,300

0

0

1,000

1,000

1,000

Таблица 10*2

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

 

8

Е

Е

Е

N .

Случай

поверхностные

 

Р 6 Ы

 

пластовые

 

 

 

 

условия

условия

 

 

1

0,28

0,55

0,61

0,34

0,14

2

0,45

0,70

0,75

0,55

0,35

3

0,80

1,00

1,00

0,80

0,60

Таблица 10.3

Важным параметром для оценки эффективности заводнения яв­ ляется отношение подвижностей в концевых точках кривых ОФП, определенное в главе 4 (раздел 9) как

М =

При горизонтальном потоке устойчивое поршневое вытеснение происходит при М < 1. Еще более важным параметром для характе­ ристики устойчивости вытеснения в соответствии с моделью БаклиЛеверетта является отношение подвижностей на фронте М$, опреде­ ляемое как

к ( 8 - ) / ц + к ( 8 3 / и

М =

(10.34)

к

/ и

ГО

Го

Относительные фазовые проницаемости в числителе этой форму­ лы определяются при водонасыщенности на фронте 8^. Используя теоретические и экспериментальные материалы, Хагурт (На§оог1) показал9, что вытеснение, соответствующее модели Бакли-Леверетта, можно считать устойчивым при менее строгом условии М < 1. Если это условие не выполняется, то происходит интенсивное образова­ ние языков воды, и прорыв происходит даже раньше, чем предска­ зывает метод Уэлджа10. В табл. 10.3 (а) приведены значения М и М§ для трех случаев, рассмотренных в упражнении 10.1. Используя эти данные, можно выполнить анализ результатов упражнения 10.1 сле­ дующим образом:

a)Случай 1. Отношение вязкостей нефти и воды очень велико, по­ этому вытеснение неустойчиво. В результате вода обходит нефть, и происходит преждевременный прорыв воды. Добыча на момент прорыва очень низка. Для извлечения всей подвижной нефти нуж­ но закачать количество воды, равное нескольким поровым объе­ мам. При таких обстоятельствах увеличение добычи путем закачки воды вряд ли осуществимо, поэтому целесообразно рассмотреть возможность применения тепловых методов повышения нефтеот­

дачи, позволяющих уменьшить отношение вязкостей.

B) Случай 2. Отношение вязкостей нефти и воды на порядок меньше, чем в первом случае. Поэтому вытеснение устойчиво и протекает в намного более благоприятных условиях (М§< 1). Этот случай под­ робно рассмотрен в упражнении 10.2, где добыча после прорыва определяется как функция накопленной закачки воды и времени.

c)Случай 3. Поскольку вязкость нефти очень мала (ро = 0,4 мПа с), как отношение подвижностей в концевых точках кривых ОФП, так и отношение подвижностей на фронте меньше единицы, и проис­ ходит поршневое вытеснение. Прямая, проведенная из точки с ко­ ординатами (8^ = 5^с; - 0), касается кривой доли воды в потоке в

точке с координатами (5 = 1 - 5 ; Г =1). Поэтому 5 =5

= 1

- 8ог. Накопленная добыча нефти на момент прорыва равна 8

- 8^8

= 1 - 8ог -8^с, то есть составляет полный объем подвижной не^ти.

Рис. 10.16. Кривые доли воды в потоке при различных отношениях

вязкостей нефти и воды (табл. 10.2)

Случай №

 

8*

к™(5*)

к „ (* У

М$

М

(упражнение 10.1)

к

 

 

 

 

 

 

1

100

0,28

0,006

0,520

1,40

37,50

2

10

0,45

0,051

0,220

0,91

3,75

3

0,4

0,80

0,300

0

0,15

0,15

Таблица 10.3(а) Значения относительных фазовых проницаемостей, соответствующие насыщенности на фронте, и относительных фазовых проницаемостей в концевых точках кривых О Ф П. Расчет выполнен по данным упражнения 10.1