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

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

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

0

10

20

30

40 ОП,

 

 

 

 

млн ст фут^/сут

Рис. 8.9. Анализ результатов исследования газовой скважины при допущении о существовании квазиустановившегося режима фильтращн в каждом периоде работы скважины График построен по данным из табл. 8 3

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

вэтой точке должны быть равны в обоих уравнениях, и

Дт (р) = 144,48 0 + 0,416 О2

При (Зст = 141,5 тыс. м3 / сут Д т (р) = 20443,92 + 8309,26 = 28773,18 МПа2 / Па с

Поскольку

А т (р) =

Ар,

взяв из табл. 8.1 данные для р. = 8,27 МПа, получаем

Др = (15,30 х 10*6 х 0,832 х 28773,18 х 106) / 2 х 8,27 х 106 = 22,15 кПа.

По уравнению материального баланса для залежи, работающей в условиях газового режима (1.35), находим

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

4) Данные для построения графика по уравнению (8.47) при приня­ том значении Р = 2,79 МПа2 / П а с / (тыс. ст. м3 / сут)2, которое было определено в первой части этого упражнения, приведены в табл. 8.4. Сам график показан на рис. 8.10.

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

Корректный график должен иметь некоторый положительный на­ клон, обусловленный тем обстоятельством, что дренируется ограни­ ченный пласт. Это означает, что следует уменьшать значение Р, действуя методом проб и ошибок, пока не будет получена прямая линия. Были по­ строены также графики при Р = 2,38,2,08 и 1,78 (см. табл. 8.4 и рис. 8.10).

 

“ о

Д т

(р) РО п %МПа2/ Пас / (тыс. ст. м3/ сут) ((фунт/

тыс. м3

часы

I

Д1.

О

дюйм2)2 / сП / (тыс.ст.фут3 / сут))

/ сут

 

 

' - о .

Р =2,79 (0,047)

Р = 1,78 (0,03)

Р =2,08 (0,035)

Р =2,37 (0,04)

283

12

12

5481,84 (3263)

5767,44 (3433)

5683,44 (3383)

5599,44 (3333)

566

24

18

5476,8 (3260)

6048 (3600)

5880(3500)

5712(3400)

849

36

24

5490,24 (3268)

6347,04 (3778)

6095,04

(3628)

5843,04 (3478)

1132

48

30

9235,77 (3271)

6637,68 (3951)

6301,68

(3751)

5965,68 (3551)

Таблица 8.4

4000

* т (р) - РО*

/

Р=0.030, А=43 АКРА

О»

7"

 

(фунт/дюйм2)2/сП/тыс стфут3/сут у

 

 

 

Р=0.035; А=60 АКРОВ

/

 

 

/

 

Р=0.040; А=100 АКРОВ

3500

Р=0.047; А ^ о о

^

а О,

Е -П Г ",

3000

10

20

30

Рис. 8.10. Анализ результатов исследования газовой скважины с использованием уравнения (8.47) при допущении о существовании квазиустановившегося режима фильтрации в каждом периоде работы скважины График построен по данным из табл 8.4

Для двух последних значений характерно небольшое искривление вверх, в то время как при Б = 2,38 график прямолинейный, с угловым коэффициентом 20,50 МПа2 / П ас / (тыс. ст. м3 / сут)/час. Отсюда следу­ ет, что дренируемый объем равен примерно 1 х 106 м3. При Ь = 12,2 миф = 0,2, принятых за средние для дренируемого объема, площадь области дренирования составит примерно 40,5 га. Судя по отрезку, отсекаемо­ му прямой на оси ординат, коэффициент В в уравнении двучленного закона фильтрации равен 5352,48 МПа2/ Пас / (тыс. ст. м3/ сут).

Нужно признать, что метод анализа при квазиустановившемся режиме фильтрации с помощью уравнения (8.47), мягко выражаясь, довольно рискованное занятие. В этом упражнении любой из графи­ ков на рис. 8.10, соответствующих различным значениям Р, может быть принят за прямолинейный, вследствие чего оценка площади области дренирования даст 24,3 га при Р = 2,08 и 17,4 га при Р = 1,78. В подавляющем большинстве случаев высокая чувствительность к значению параметра Р не позволяет сделать сколько-нибудь надеж­ ную оценку объема порового пространства или площади области дренирования.

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

к!

<ОА= 3600

<р (рс).А

Поэтому такие исследования дают хорошие результаты для не­ больших высокопроницаемых пластов, где относительно быстро возникает квазиустановившийся фильтрационный поток. Однако это не единственный критерий. Например, в упражнении 8.1 фак­ тическая проницаемость равна 0,1 мкм2, и скважина расположена в центре квадратной области дренирования площадью 40,5 га. В таких условиях значение 1ОА, при котором происходит переход к квазиуста­ новившейся фильтрации, равно 0,1 (см. графики МБХ, рис. 7.11, или таблицу коэффициентов формы Дитца, рис. 6.4).

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

1 = 0,1 х <р (рс). А / 3600 к = 0,1 х 0,2 х 0,522 х 10'6 х 10'6 х 40,5 х 10000 / (3600 х 0,1 х 1 0 12) = 11,7 ч.

Результат подтверждает, что в данном случае продолжительность отдельных периодов работы скважины 12 часов достаточна для того, чтобы можно было проводить анализ с использованием уравнений для квазиустановившейся фильтрации. Однако допустим, что сква­ жина расположена не в центре квадратной области дренирования с непроницаемой границей, а в центре одной из четвертей такой об­ ласти дренирования (график III на рис. 7.12 (а)). В таком случае зна­ чение 1оа, при котором происходит переход к квазиустановившейся фильтрации, равно 0,5. Следовательно, при тех же параметрах пла­ ста применение уравнений для квазиустановившейся фильтрации допустимо лишь тогда, когда продолжительность каждого периода работы скважины достигнет примерно 60 часов. Если же продолжи­ тельность периодов работы скважины останется прежней, то есть 12 часов, то применять уравнения для квазиустановившейся фильтра­ ции некорректно, поскольку все исследование будет проводиться в

позднем периоде неустановившейся фильтрации, и рассчитывать т 0 нужно по уравнению (8.40), а не по (8.33). При исследовании высоко­ проницаемых пластов, таких как рассмотренные в упражнении 8.1, использование в анализе неправильно выбранных математических выражений не повлечет за собой серьезных ошибок определения В и Е Однако при исследовании низкопроницаемых пластов ошибки могут весьма велики.

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

Вупражнении 8.1 эти значения равны: к = 0,1 мкм2, А = 40,5 га, СА

=30,9, 5 = 6,0. За исключением А, определенного путем довольно не­ точного анализа в четвертой части упражнения, ни один из параме­ тров не был получен непосредственно по данным исследования. Это обстоятельство может привести к нежелательным последствиям, ког­ да значение В, полученное по данным исследования, используется в уравнении притока с целью прогнозирования продуктивности сред­ нестатистической скважины в долгосрочной перспективе. Например, ожидается высокий скин-фактор в поисковой скважине, но при бу­ рении среднестатистической эксплуатационной скважины благодаря применению лучше подобранного бурового раствора для заканчивания скин-фактор будет меньше. Если удастся уменьшить скин-фактор

вскважине, рассмотренной в упражнении 8.1, до нуля, то значение коэффициента В уменьшится с 5335,68 до 2970,24 МПа2 / Па с/ тыс. м3. Далее, если пласт имеет большие размеры и для его разработки требуется больше одной скважины, то А и СА, неявно выраженные в В, определенном по данным исследования скважин, будут изменяться

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

По мере того как открывались новые месторождения и выполня­ лись исследования скважин, становилось очевидно, что метод ана­

лиза при квазиустановившемся режиме фильтрации, который рас­ сматривали до сих пор в этой главе, не подходит для определения коэффициентов В и Р в уравнении притока (8.44). Причина заклю­ чается в том, что при низкой проницаемости пласта продолжитель­ ность работы скважины до перехода к квазиустановившейся филь­ трации (1оа = 3600 к! / ср (рс). А) может стать чрезвычайно большой. В таких случаях использование для анализа уравнения квазиустановившегося притока становится неприемлемым, поскольку изменение давления в скважине вследствие изменения дебита существенно за­ висит от времени. В литературе были представлены два метода, учи­ тывающие эту зависимость. Эти методы используют допущение, что изменение давления в скважине можно выразить путем суммиро­ вания решений уравнения (5.20) при постоянном дебите в условиях неустановившейся фильтрации. Один из них - метод Оде-Джонса13 (1965 г.), в котором уравнение для анализа результатов исследования выражено через р2. Второй - метод Эссиса (Е8$1$) и Томаса (ТЪотаз)14 (1971 г.). Данный метод использует модифицированное уравнение Оде-Джонса, выраженное через псевдодавление реального газа

кЬ

( ш (р,) - ш ( ? „ » = 1 Щ, тв (Ч - ц () + С>„3;.

113 Т

Это записанное в общем виде выражение, содержащее сумму ре­ шений для постоянного дебита (8.39). И Оде, и Джонс, и Эссис, и То­ мас применяли свои методы анализа строго в условиях неустановив­ шейся фильтрации, когда т 0 в уравнении (8.39) можно выразить с помощью (8.32)

тгА ) = °>51п

ане с помощью общего выражения (8.40), предполагающего знание геометрии и площади области дренирования. Если можно использо­ вать выражения для т 0, применимые при неустановившемся режиме фильтрации, то анализ прост и позволяет легко найти значения к и 8. Эти значения в свою очередь можно применить для расчета коэф­ фициента В при любых размерах области дренирования и при любом коэффициенте формы (уравнение 8.44). Кроме того, можно непосред­ ственно определить второй коэффициент в уравнении притока (Р или О). Этот метод будет подробно изложен в упражнении 8.2.

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

т о ({'0тах) = т 0 (полное безразмерное время исследования)

по-прежнему можно было определить по уравнению, применимо­ му при неустановившемся режиме фильтрации (8.32). В следующем упражнении оценивается ошибка, которую можно допустить, при­ няв слишком продолжительные периоды работы скважины.

УПРАЖНЕНИЕ 8.2. АНАЛИЗ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЯ ГАЗОВОЙ СКВАЖИНЫ МЕТОДОМ МНОГОКРАТНОГО ИЗМЕНЕНИЯ РЕЖИМА

СДОПУЩЕНИЕМ О СУЩЕСТВОВАНИИ УСЛОВИЙ НЕУСТАНОВИВШЕЙСЯ ФИЛЬТРАЦИИ

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

Пластовая температура и характеристики газа опять такие же, как в табл. 8.1, и поэтому зависимость между давлением и псевдодавле­ нием для р > 19,31 МПа описывается уравнением (8.48). Ниже даны параметры пласта и скважины.

 

Продолжительность

 

 

Дебит (}, тыс. м3

работы скважины

МПа

ш (р^), МПа2 / Пас

 

1, часы

 

 

283

1

28,68

48,61 х 106

566

2

27,56

45,93 х 106

849

3

26,21

42,73 х 10б

1132

4

24,66

39,04 x10й

Таблица8«5

р.

-

29,58 МПа

 

(рс)1

-

0,522х 10*12 Па с / Па

А

-

81 га

 

ср 1

-

0,15

Форма -

4

1

Ь

-

15,24 м

4

 

 

 

 

г

-

91,4 мм

Требуется определить т 0(мвн) (1^) для 1од <0,01, то есть нижнего предела шкалы безразмерного времени на каждом графике МБХ, рис. 7.11-15. Значения т 0(мвн) при < 0,01 для формы 4:1 можно найти по приведенному на рис. 8.11 графику, который построен по табулиро­ ванной зависимости между т Е)(мвн и 1од, представленной Эрлафером15. (Примечание: в вышеуказанной работе, по-видимому, есть типограф­ ская опечатка, относящаяся к геометрической форме 4:1, а именно

И

в табл. 2 на стр. 203 нужно поменять местами.)

 

Прежде всего нужно выполнить анализ КВД и определить прони­ цаемость, в данном случае 50 х 10'3мкм2 (этот метод будет подробно изложен в разделе 8.11 и проиллюстрирован упражнением 8.3).

Требуется:

Выполнить анализ результатов исследования методом многократно­ го изменения режима работы скважины и определить к, 8 и Р, ис­ пользуя

1)Значение т 0, полученное для условий неустановившейся филь­ трации (уравнение 8.32)

2)Значение т 0, полученное по общему выражению (уравнение 8.40)

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

1) Используемое в методе Эссиса-Томаса уравнение (8.39) можно за­ писать в более удобном для практического применения виде

Рис. 8.11. График МБХ для формы области дренирования 4:1 при 10А< 0,01 (по Эрлаферу и др.15)

130Т

 

т

(р ,) - т

( р „ Гп) - РС>> =

кЬ

Е д о )18 а п- ^

1) +

 

 

 

 

 

)=*

 

 

 

 

 

 

 

 

 

 

к1

+ 0,35 + 0,87 8

 

 

 

 

 

+ С>„

18

 

 

 

 

 

 

<р(ИЧ

 

 

 

 

 

 

т

(р.) - т

(Р«г„) - р(3п

ДО,

, ,

 

 

 

 

И Л И

------------ —— --------- = т I

------1—1е (I: -1 ,)+

 

 

 

 

с>„

 

 

о

8 4

 

 

 

 

 

 

 

+ ш (

------------- + 0,35 + 0,87 8

^

(8.

 

 

 

\

 

 

 

 

 

 

/

 

Таким образом, график зависимости

 

 

 

 

 

 

 

ш (р .) - ш ( р и Г ) -

 

 

^

Д О ,

 

(4

- I .

 

 

 

 

 

 

 

2

------- Ч

е

,)

 

 

О,

 

 

=,

п

6

"

 

 

 

 

 

 

»-•

о,

 

 

 

 

должен быть линейным, иметь угловой коэффициент ш = 130 Т / и отсекать на оси ординат отрезок

/ к

т[ ------------+ 0,35 + 0,87 8

По этим данным можно рассчитать к и 5. В таком анализе значение Р ищется методом проб и ошибок, пока не будет получена прямая ли­ ния. Результаты анализа приведены в табл. 8.6. График, построенный по полученным значениям, показан на рис. 8.12.

Первый график построен без учета отклонения потока от за­ кона Дарси (Р = 0). Он заметно искривляется в области больших значений абсциссы, что объясняется возрастанием дебита к кон­ цу исследования. Если бы исследование проводилось с последо­ вательным уменьшением дебита, то искривление было бы более интенсивным при малых значениях абсциссы14. С целью линеари­ зации графика был выполнен расчет левой части уравнения (8.49) при различных значениях Р, результаты которого приведены в табл. 8.6. Как можно видеть на рис. 8.12, линеаризация достигается при Р = 2,97. Прямолинейный график имеет угловой коэффициент т = 825 (491,5 ) и отсекает на оси ординат отрезок 6708 (3993). По этим данным можно определить проницаемость и скин-фактор

к = 130 х 660 / (1,8 х 825 / (11,57 х 1015) х 15,24) = 44 х 10‘3 мкм2

и

отсекаемый отрезок

-18

3,35 | = 2,7

5= 1,151

 

ш

 

ф(нс),^

0, тыс.

1,

 

 

т

м3 / сух часы

 

МПа

МПа2/

 

 

 

 

Пас

283

1

0

28,68

48,61 х

 

 

 

 

10б

566

2

0,1505

27,56

45,93 х

 

 

 

 

10‘

849

3

0,2594

26,21

42,73 х

 

 

 

 

10б

1132

4

0,3450

24,66

39,04 х

 

 

 

 

106

т (р.) - т (р ) - Р(^

МПа2 / Па с/ (1000 мЧ

--------------------------- ,

сут) ((фунт / дюйм2)2

Оп

/ сП / (1000 фут3/сут))

•п

II

Р = 2,37

Р = 2,97

Р = 3,56

о

(0,04)

(0,05)

(0,06)

7549,92

6877,92

6709,92

6541,92

(4494)

(4094)

(3994)

(3894)

8509,2

7165,2

6829,2

6526,8

(5065)

(4265)

(4065)

(3885)

9446,64

7430,64

6926,64

6422,64

(5623)

(4423)

(4123)

(3823)

10352,16

7664,16

6992,16

6320,16

(6162)

(4562)

(4162)

(3762)