Файл: Закиров, С. Н. Проектирование и разработка газовых месторождений учебное пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 164

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ГЛАВА VI

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

§ 1. О методах расчета продвижения контурных или подошвенных вод в газовые залежи

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

При проектировании разработки газового месторождения в усло­ виях водонапорного режима необходимо бывает знать закон измене­ ния во времени границы раздела газ— вода (газонасыщенного объема порового пространства). Тогда согласно, например, следующему уравнению материального баланса (см. § 4 глава И)

( 1 >

можно определить изменение во времени средневзвешенного по газо­ насыщенному объему порового пространства газовой залежи пласто­ вого давления.

Присоединяя к уравнению (1) уравнение притока газа к средней скважине

Р г (0 — P c(t) — A q (t) -f- В ф (t)

(2)

или, в случае учета реальных свойств газа, — уравнение

 

Р 2 (0 — Pi (0 = А *z)cyq (t) - f B zcyq2 (t),

(2a)

а также уравнение технологического режима эксплуатации средней скважины (например, уравнение технологического режима эксплуа­ тации скважины при заданной депрессии на пласт)

P ( t ) — p c (t) = b { t )

(3)

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

n ( t )

<?(0

(4)

9 ( 0 ’

получаем систему из четырех уравнений с пятью неизвестными. Совместное решение системы уравнений (1)— (4) и задачи о продвиже­ нии воды в залежь позволяет определить изменение во времени сле­ дующих основных показателей разработки месторождения в усло­ виях водонапорного режима (в предположении равномерного раз­ мещения скважин на площади газоносности):

P = P ( t ) , Pc = P c ( t ) , a Q = f ( t ) , q — q (t), n = n ( t ) ,

(5)

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

Расчеты продвижения контурных или подошвенных вод в газо­ вые залежи относятся к наиболее сложным задачам подземной газогидродинамики — к задачам с подвижной границей раздела (типа Стефана — в математической физике [61]). А так как при этом важно учитывать неоднородность пласта вдоль осей х, у, z, то эти задачи

еще более усложняются.

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

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

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

Простой приближенный метод расчета продвижения контурной воды в га­ зовую залежь с учетом вязкости воды предложен Б. Б. Лапуком [83]. Этот метод был развит и для расчетов продвижения границы раздела газ—вода в неоднород­ ных пластах (для одномерных фильтрационных потоков). Метод, изложенный в работах Б. Б. Лапука, основан на последовательных приближениях при нахо­ ждении истинного количества поступающей в газовую залежь воды и соответ­ ствующего изменения во времени среднего пластового давления в газовой за­ лежи. В следующем параграфе на идее метода Б. Б. Лапука остановимся под­ робнее.

172


Применение метода последовательной смены стационарных состояний по­ зволило И. А . Парному приближенно учесть упругие свойства водоносного пласта и противодавление поступающей в газовую залежь воды.

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

ния воды в галерею, расположенную в полубесконечном газоносном

пласте,

и закачки воды в скважину исчезающе малого радиуса, находящуюся

в беско­

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

Задачами с подвижной границей раздела при учете неоднородности пласта по проницаемости занимались М. Т. Абасов, С. И. Алекперов и М. В. Филинов. Полученные решения могут использоваться как эталонные для сопоставления с ними приближенных решений.

Некоторые исследователи (М. Т. Абасов, Ф. Г. Гасанов, Т. Гольдфрахт, П. Ионеску, П. К. Страдымов, Е. Н. Храменков, А. И. Ширковский, П. Т. Шмыгля, К. Мюллер п др.) решение задачи о движении границы раздела газ—вода находят путем совместного рассмотрения уравнения материального баланса и соответствующей формулы теории упругого режима для падения давле­ ния на стенке укрупненной скважины при пуске ее в работу с постоянным де­ битом в бесконечном (или ограниченном) однородном по коллекторским свой­ ствам пласте. В ряде других исследований или упрощаются существующие реше­ ния, или показывается возможность использования какого-либо метода или под­ хода к решению задачи о движении границы раздела газ—вода (М. Т. Абасов, К. Н. Джалилов, А. М. Кулиев, А. С. Малых, Н. Мухидинов, Ф. Г. Темпель, М. В. Филинов, И. А. Парный).

Так как в последние годы лабораторными и промысловыми исследованиями установлено, что газ не полностью вытесняется водой, в работе [81 ] известный метод Баклея— Леверетта использован для подсчета остаточной газонасыщенности за фронтом вытесненпя. Наиболее совершенная методика расчета посту­ пления воды в залежь, учитывающая особенности проявления водонапорного режима, изложена в работе [71].

Задачи о разработке газонефтяных месторождений при использовании энер­ гии сжатого газа газовой шапки аналогичны задачам о вторжении воды при раз­ работке газовых месторождений. Разница состоит в том, что отбираемое и зака­ чиваемое количества газа входят в уравнение материального баланса с различ­ ными знаками. В связи с этим полученные решения относительно расширения газовой шапки могут быть использованы при расчетах разработки газовых месторождений в условиях водонапорного режима и наоборот. В этом смысле представляют интерес исследования М. Д. Розенберга и М. М. Глоговского про­ цесса разработки газонефтяных месторождений при заданных давлениях на забоях системы инжекционных (газовых) и эксплуатационных (нефтяных) скважин, а также при заданных дебитах эксплуатационных скважин и расходах нагнетаемого в пласт газа. Результаты этих работ могут быть использованы при рассмотрении возможности поддержания пластового давления в газоконденсат­ ных месторождениях путем законтурного заводнения и разработки газовых месторождений в условиях водонапорного режима, а также при создании и эксплуатацпп подземных газохранилищ в водоносных пластах. Методики прогноз­ ных расчетов применительно к подземным газохранилищам в водоносных пла­ стах приводятся в работах [9, 73, 75].

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

173


в залежь воды и показателей разработки залежи учитываются наи­ более интересные подходы. В § 2, 3, 4 и 5 изложены методы расчета продвижения воды в газовую залежь при замене ее укрупненной скважиной. Результаты расчетов могут быть использованы для определения показателей разработки в условиях равномерного (близкого к равномерному) размещения скважин на площади газо­ носности.

Методам расчета продвижения воды в газовые залежи при разме­ щении скважин в центральной зоне или при неравномерном размеще­ нии скважин на площади газоносности посвящается последний параграф. Задача с п о д в и ж н о й границей раздела газ — вода при не­

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

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

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

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

на основе теории упругого режима фильтрации

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

Принимаем следующую схематизацию. Газовая залежь предста­

вляется укрупненной скважиной радиусом R

3 Г Радиус

укрупнен­

ной скважины определяется из равенства S =

nRI (здесь

S — пло­

щадь газоносности). Если возмущение, вызванное разработкой газо­ вой залежи, за рассматриваемый период времени не достигает внеш-1

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

174


ней границы, то водоносный пласт принимается бесконечным по протяженности. В противном случае водоносный пласт предста­ вляется круговым с радиусом R K.

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

Необходимо определить показатели разработки газовой залежи при водонапорном режиме, при которых обеспечивается получение заданного отбора газа Q = Q (t). Согласно изложенному в преды­

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

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

тока воды

к укрупненной скважине.

 

 

 

Пусть

укрупненная

скважина

радиусом

R 3

эксплуатируется

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

Ар

— р а Рс (Рн

начальное

давление в водоносном пласте, р с

забойное давление

на стенке

укрупненной

скважины).

Требуется

найти зависимость

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

уравнения теории упругого

режима

фильтрации

 

1

др

.

д12р

 

1

др

(1)

г

дг

1

dr2

 

к

dt

при следующих условиях:

 

 

 

 

 

 

 

 

t =

0,

Р =

Рп,

 

(2)

r = R 3,

Ар =

р н — р с =

const;

(3)

или

г - + о о ,

р = рн

(4а)

 

 

 

 

 

 

 

или

r =

R K,

Р =

Рн

 

(46)

 

 

 

 

 

 

 

 

г =

К

,

4 г

“ °-

(4в)

 

 

В уравнении (1) и — средний для водоносного пласта коэффи­ циент пьезопроводности.

Уравнение (1) описывает изменение давления в любой точке водоносного пласта в любой момент времени. Начальное условие (2) характеризует невозмущениость водоносного пласта до начала разработки залежи. Граничное усло­ вие на стенке укрупненной скважины (3) показывает, что она эксплуатируется при постоянном во времени перепаде давления Др. На внешней границе водо­ носного пласта задается одно из условий (4а)—(4в). Если водоносный пласт мо­ жет рассматриваться как бесконечный, то используется условие (4а). При нали­ чии области питания на расстоянии Дк решается задача (1)— (3), (46) Н При

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

175


непроницаемости внешней границы водоносного пласта задача решается при со­ блюдении условия (4в).

В результате решения задачи (1)— (4) определяется зависимость изменения во времени дебита воды, поступающей в газовую залежь (укрупненную скважину):

9в =

(0 -

(5)

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

 

 

t

 

<?в (0

=

(6)

 

 

О

 

Интегрирование (1)— (4),

осуществленное Ван

Эвердингеном и

Херстом [96], приводит к следующему выражению для QB (t):

& (*) =

-^

f - A p § (fo>.

(7)

Здесь к — средний для водоносного пласта коэффициент прони­ цаемости; h — средняя для водоносного пласта эффективная мощ­

ность; [хв — коэффициент динамической вязкости

воды при пласто­

вой температуре; fo =

— параметр Фурье.

 

Функция Q (fo) записывается по-разному в зависимости от того,

при каком из условий —

(4а), (46) или (4в) —

решается задача.

В случае, например, пласта, бесконечного по протяженности,

функция Q (fo) будет

 

 

 

( l - e “ 2f0) Ли

(8)

 

из i q (M) _ y 2 {tt)]

Здесь / 0 и У 0 — функции Бесселя соответственно первого и вто­ рого рода нулевого порядка.

Расчеты по формуле (7) с учетом выражений типа (8) (для различ­ ных условий на внешней границе водоносного пласта) очень гро­ моздки. Поэтому в работе [96] для функций Q = Q (fo) составлены

таблицы (см. также [35]). Помимо таблиц для функции Q (fo), можно пользоваться и приближенными, но достаточно точными формулами Ю. П. Борисова, В. Е. Влюшина и В. Н . Щелкачева.

Естественно, что укрупненная скважина (газовая залежь) экс­ плуатируется с переменными во времени перепадами давления на забое. Поэтому воспользуемся принципом суперпозиции. Однако прежде необходимо задаться темпом изменения забойного давления на стенке укрупненной скважины (пластового давления в газовой залежи).

Так как изменение давления в газовой залежи зависит от темпа поступления воды, а поступление воды определяется изменением

176