Методичка для практических занятий по разработке


Министерство образования и науки Российской Федерации
Федеральное государственное бюджетное образовательное учреждение
высшего профессионального образования
«УФИМСКИЙ ГОСУДАРСТВЕННЫЙ НЕФТЯНОЙ
ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
Кафедра разработки и эксплуатации газовых и газоконденсатных месторождений
РАЗРАБОТКА И ПРОЕКТИРОВАНИЕ МЕСТОРОЖДЕНИЙ ПРИРОДНЫХ ГАЗОВ
Учебно-методическое пособие
по выполнению практических занятий по дисциплинам «Разработка и проектирование газовых, газоконденсатных и нефтегазоконденсатных месторождений» и
«Теория проектирования разработки месторождений»
Уфа
2014
Учебно-методическое пособие предназначено для бакалавров, обучающихся по направлению 131000 «Нефтегазовое дело» по профилям подготовки «Эксплуатация и обслуживание объектов добычи газа и газоконденсата и подземных хранилищ» и «Эксплуатация и обслуживание объектов нефтегазового комплекса арктического шельфа» (БГГ, БГШ), и для магистрантов, обучающихся по направлению 131000.68 «Нефтегазовое дело» по магистерской программе «Проектирование и управление разработкой и эксплуатацией газовых, газоконденсатных и нефтегазоконденсатных месторождений» (МГГ).Составители: Пономарев А.И., проф.
Калиновский Ю.В., доц.
Шаяхметов А.И., преп.
Рецензент:
© Уфимский государственный нефтяной технический университет, 2014СОДЕРЖАНИЕ
ВВЕДЕНИЕ…………………………………………………………………..……4
Практическое занятие №1 по теме «Определение дренируемых запасов газа газовой залежи по методу падения пластового давления»………………….…6
Практическое занятие №2 по теме «Расчет параметров «средней» скважины»………………………………………………………………………..17
Практическое занятие №3 по теме «Методика расчета показателей разработки газовой залежи при газовом режиме для заданной динамики темпа отбора газа и технологического режима эксплуатации «средней» скважины с постоянной депрессией на пласт»………………………………...19
Практическое занятие №4 по теме «Методика расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории укрупненной скважины»………...…28
Практическое занятие №5 по теме «Методика расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории Фетковича»……………………….…36
Практическое занятие №6 по теме «Методика расчета показателей разработки газовой залежи при упруговодонапорном режиме с учетом образования общей депрессионной воронки в зоне разбуривания»………....43
Практическое занятие №7 по теме «Методика расчета неравномерного распределения давления по площади газовой залежи, разрабатываемой кольцевыми батареями скважин»……………………………………………....50
Практическое занятие №8 по теме «Методика расчета показателей разработки слоисто-неоднородной газовой залежи при упруговодонапорном режиме»…………………………………………………………………………..58
Практическое занятие №9 по теме «Расчет производительности обводняющейся газовой скважины в слоисто-неоднородном пласте»…...….65
Практическое занятие №10 по теме «Определение показателей разработки газоконденсатной залежи при полном сайклинг-процессе»……………….…69
СПИСОК ЛИТЕРАТУРЫ…………………………………………………...…..81
Приложение А………………………………………………………………..…..82
Приложение Б………………………………………………………………..…..84
ВВЕДЕНИЕ
Настоящее методическое пособие разработано для бакалавров, обучающихся по направлению 131000 «Нефтегазовое дело» по профилям подготовки «Эксплуатация и обслуживание объектов добычи газа и газоконденсата и подземных хранилищ» и «Эксплуатация и обслуживание объектов нефтегазового комплекса арктического шельфа» (БГГ, БГШ), для магистрантов, обучающихся по направлению 131000.68 «Нефтегазовое дело» по магистерской программе «Проектирование и управление разработкой и эксплуатацией газовых, газоконденсатных и нефтегазоконденсатных месторождений» (МГГ).Цель настоящего методического пособия заключается в изучении предлагаемых математических моделей, описывающих процессы разработки газовых и газоконденсатных месторождений; методик и методов, используемых для поиска решения; в развитии навыков разработки алгоритмов и компьютерных программ для решения рассматриваемых задач, а также в закреплении отдельных тем лекционных курсов дисциплин «Разработка и проектирование газовых, газоконденсатных и нефтегазоконденсатных месторождений» и «Теория проектирования разработки месторождений».В методическом пособии приведены следующие задачи:
определение дренируемых запасов газа газовой залежи по методу падения пластового давления;
расчет параметров «средней» скважины;
расчет показателей разработки газовой залежи при газовом режиме для заданной динамики темпа отбора газа и технологического режима эксплуатации «средней» скважины с постоянной депрессией на пласт;
расчет показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории укрупненной скважины;
расчет показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории Фетковича;
расчет показателей разработки газовой залежи при упруговодонапорном режиме с учетом образования общей депрессионной воронки в зоне разбуривания;
расчет неравномерного распределения давления по площади газовой залежи, разрабатываемой кольцевыми батареями скважин;
расчет показателей разработки слоисто-неоднородной газовой залежи при упруговодонапорном режиме.
расчет производительности обводняющейся газовой скважины в слоисто-неоднородном пласте;
определение показателей разработки газоконденсатной залежи при полном сайклинг-процессе.
При решении вышеуказанных задач используются следующие методы:
метод последовательной смены стационарных состояний;
метод последовательных приближений (итераций);
метод деления отрезка пополам при решении трансцендентного уравнения;
метод наименьших квадратов.
Практическое занятие №1 по теме «Определение дренируемых запасов газа газовой залежи по методу падения пластового давления»Цель: изучить метод определения дренируемых запасов газа по падению пластового давления, определить по данным эксплуатации залежи начальные дренируемые запасы газа, начальное пластовое давление, балансовые и извлекаемые запасы, коэффициент газоотдачи и начальный газонасыщенный поровый объём.
Основные понятия и определения
В соответствии с классификацией запасов месторождений нефти и газов они подразделяются на две группы.
К I группе относятся месторождения простого геологического строения, продуктивные пласты которых характеризуются выдержанностью толщин и коллекторских свойств по площади и разрезу.
Ко II группе относятся месторождения сложного геологического строения, характеризующиеся невыдержанностью толщины и коллекторских свойств продуктивных пластов (горизонтов).
Все разведанные запасы нефти, горючих газов, конденсата, а также содержащихся в них сопутствующих компонентов, которые служат или могут служить сырьевой базой для действующих, реконструируемых и проектируемых предприятий, подлежат обязательной проверке и утверждению Государственной комиссией по запасам полезных ископаемых.
Запасы нефти, горючих газов и сопутствующих компонентов по их народнохозяйственному значению разделяются на две группы, подлежащие отдельному подсчету и учету:
Балансовые запасы – это запасы, разработка которых в настоящее время экономически целесообразна. В балансовых запасах нефти, растворенного в ней газа, а также конденсата в свободном газе выделяются и учитываются извлекаемые запасы, то есть запасы, которые можно извлечь при наиболее полном и рациональном использовании современной техники и технологии.
Забалансовые запасы – это запасы, разработка которых в настоящее время нерентабельна, но которые могут рассматриваться в качестве объекта для промышленного освоения в дальнейшем.
В балансовых запасах газа, а также конденсата в свободном газе выделяются и учитываются извлекаемые запасы, т. е. запасы, которые можно извлечь при наиболее полном и рациональном использовании современной техники и технологии добычи.
Запасы газа и газового конденсата подсчитываются на структурных планах, составленных в зависимости от размера месторождения, в масштабе, обеспечивающем необходимую точность замера площадей (1:5000 – 1:50000). Границы подсчета запасов по месторождению, отдельным залежам и тектоническим блокам принимаются по данным разведки и должны быть увязаны с геологическими особенностями месторождения.
Подсчет и учет запасов газа, газового конденсата и содержащихся в них сопутствующих компонентов должны производиться раздельно для каждой изолированной залежи. Запасы конденсата подсчитываются в тысячах тонн. Запасы природных газов подсчитываются в миллионах, а гелия – в тысячах кубических метров, приведенных к стандартным условиям (760 мм. рт. ст. или 0,101325 МПа и 20°С).
Подсчет запасов газа и конденсата месторождений (залежей) природных газов можно производить объемным методом и методами, основанными на принципе материального баланса (по методу падения давления). При этом метод падения пластового давления является дополнительным, используемым для оперативной оценки дренируемых запасов в ходе анализа разработки залежи. Комплексное применение объемного метода и метода падения давления позволяет повысить точность определения запасов залежи (месторождения).
Подсчет запасов газа газовых залежей производится объемным методом, а при наличии фактических полноценных геолого - промысловых данных - и по методу падения давления. По методу падения давления подсчет запасов производится по залежам, в которых доказано отсутствие запасов нефти промышленного значения и резко выраженного водонапорного режима, определено изменение приведенного пластового давления от суммарного отбора газа из залежи во времени, установлено снижение средневзвешенного пластового давления и оценено количество пластовой воды, поступившей в залежь за период эксплуатации.Забалансовые запасы газа, конденсата подсчитываются по степени их изученности с объяснением причин отнесения их к забалансовым [2, 6].
При подсчете запасов газа методом падения давления должны быть установлены:
а) размеры и форма залежи;
б) тектонические особенности и литологический состав продуктивного пласта, и изолированность отдельных частей (блоков) залежи;
в) начальное и текущие высотные положения газоводяного контакта;
г) характеристика газогидродинамической связи залежей месторождения;
д) начальное статическое и пластовое давление, и пластовая температура, а также изменение приведенного пластового давления по скважинам и средневзвешенного по залежи во времени;
е) отбор газа и конденсата по скважинам и по залежи в целом;
ж) графическая зависимость средневзвешенного приведенного пластового давления газа от отбора его по залежи;
з) степень дренируемости скважинами объема газовой залежи;
и) при наличии конденсата – содержание его в газе, состав и коэффициент извлечения его при текущем пластовом давлении;
к) газогидродинамические условия и режим работы залежи и отдельных ее частей;
л) дата начала вторжения пластовой воды и количество ее, поступившее в пласт за период эксплуатации, рассчитанное различными методами;
м) перетоки и потери газа.
Методика определения дренируемых запасов газа газовой залежи по методу падения пластового давления
Краткая теория вопроса
Предложенная ниже методика подсчета запасов газа по методу падения пластового давления применима для залежей, относящихся к I группе месторождений простого геологического строения, продуктивные пласты которых характеризуются выдержанностью толщины и коллекторских свойств по площади и разрезу.
При подсчете запасов газа по методу падения пластового давления для залежей, относящихся ко II группе, можно предложить некоторые модификации метода падения давления (метод удельных объёмов дренажа, метод средневесовой плоскости и т.п.) [2, 7].
Обработка промысловых данных (средние пластовые давления и соответствующие добытые количества газа на различные моменты при уверенности, что режим залежи газовый) по уравнению материального баланса с использованием метода наименьших квадратов позволяет вычислить газонасыщенный объем порового пространства , а затем и запасы газа. Однако для более обоснованного определения запасов газа по падению среднего пластового давления промысловые данные подвергают графической обработке. Это позволяет исключить из рассмотрения (после соответствующего анализа) дефектные точки или установить причины различных аномалий, отклонений. Графический метод обработки промысловых данных позволяет с большей наглядностью определить режим залежи, момент начала активного продвижения воды.
Для промышленной оценки запасов месторождений или отдельных залежей газов определяющее значение имеют: форма и площадь, а также толщина, коллекторские свойства, газонасыщенность и эксплуатационная характеристика продуктивных пластов (горизонтов) [2, 6, 7].
Толщины продуктивных горизонтов или отдельных пластов-коллекторов весьма разнообразны и колеблются от нескольких сантиметров до десятков, а иногда и сотен метров.
Различают общую толщину продуктивного пласта, включающую от кровли до подошвы все прослои проницаемых и непроницаемых пород; эффективную (полезную), состоящую из суммы толщин проницаемых пластов-коллекторов, и газонасыщенную толщину пластов-коллекторов, включающую только те прослои пород, которые содержат нефть или газ.
Пористость в зависимости от сообщаемости пор и их насыщенности нефтью или газом разделяется на общую, открытую и эффективную. При подсчете запасов принимается открытая пористость.
Пластовое давление – это давление пластового флюида, насыщающего поровое пространство горной породы.
Приведенное пластовое давление – это пластовое давление, приведенное к единой плоскости (начальное положение ГВК, середина этажа газоносности) [6].
Средневзвешенное пластовое давление – это среднее пластовое давление, взвешенное по газонасыщенному поровому объему и определяемое по формуле
.(1)
Представим уравнение материального баланса для залежи в случае газового режима в виде
,(2)
гдеРн, , – начальное, текущее средневзвешенное по объему порового пространства пластовое давление и давление стандартных условий (0,101325 МПа), МПа;
, – пластовая температура и температура стандартных условий (293,15 К), К;
, , – коэффициент сверхсжимаемости для текущего, начального давления и давления стандартных условий;
– объем добытого газа на момент времени t, приведенный к стандартным условиям, млн. м3;
– газонасыщенный объем порового пространства, млн. м3.
Будем откладывать по оси абсцисс отобранные объемы газа , по оси ординат - на разные моменты времени. Из уравнения (2) следует, что зависимость представляет собой линейную вида
,(3)
где b =, (4)
а = .(5)
Коэффициенты а и b найдем методом наименьших квадратов:
(6)
,(7)
где (8)

Q(t i) (9)
.(10)
,(10*)
, (11)
где Н – глубина залежи, км;
H = – при соответствии начального пластового давления гидростатическому закону;
– коэффициент конечной газоотдачи;
Qизвл – извлекаемые запасы пластового газа, млн. м3;
– балансовые запасы пластового газа, млн. м3;
– пластовое давление к концу разработки – давление «забрасывания», МПа;
– плотность воды (103 кг/м3);
g – ускорение свободного падения (9,81 Н/кг, или м/с2);
– пластовое давление на i-й момент времени, МПа;
– накопленная добыча газа на i-й момент времени, млн. м3;
N – количество интервалов времени;
i – номер момента времени.
Расчет коэффициента сверхсжимаемости производится либо аналитически (формула В.В. Латонова – Г.Р. Гуревича (12) [1], по уравнению состояния реального газа и т.п.), либо графически (по графикам Стендинга и Катца), по данному составу газа и определенным давлению и температуре.
(12)
(13)
(14)
(15)
,(16)
где , – критические температура и давление i- го компонента смеси газов;
, – псевдокритические температура и давление смеси газов;
, – приведенные температура и давление смеси газов;
М – количество компонентов в смеси;
i – номер компонента;
– мольная доля компонента смеси.
По полученным значениям коэффициентов с использованием уравнения (2) определяем балансовые () и извлекаемые запасы (), а затем коэффициент конечной газоотдачи [8].
Необходимо отметить, что пластовое давление не равно устьевому давлению и вышеуказанное равенство принято нами с целью упрощения расчетов. При проведении точных расчетов балансовых запасов необходимо рассчитать пластовое давление, соответствующее устьевому давлению в 1 атм. Расчеты производятся по нижеприведенным формулам [5]:
,(17)
,(18)
,(19)
,(20)
где Тср – средняя по стволу скважины температура, К;
Ту, Тпл – устьевая и пластовая температура соответственно, К;
– плотность газа в стандартных условиях, кг/м3;
– относительная плотность газа по воздуху;
zср – коэффициент сверхсжимаемости, рассчитанный для средних по стволу скважины давления и температуры;
L – глубина скважины, м.
Расчет пластового давления осуществляется с помощью последовательного приближения (итерационно).
Кроме того, по полученным коэффициентам определяют начальное пластовое давление и газонасыщенный поровый объем. При подсчете запасов газа по методу падению пластового давления необходимо учитывать, что значение пластового давления и точность его определения очень сильно
Рисунок 1 – Графики зависимости коэффициента сверхсжимаемости z от приведенных температуры и давления по данным Стендинга и Катцасказывается на достоверности запасов. Следует помнить, что применяемые на практике манометры градуированы на технические (1ат=0,0980665 МПа) атмосферы и показывают избыточное давление. Необходимо также учитывать класс точности манометра и способ определения средневзвешенного по объему порового пространства пластового давления (измерение статического устьевого или забойного давления, расчет давления в газовой залежи по напору законтурных вод, вычисление пластового давления по данным испытания скважин по методу противодавления).
Задача 1. Определить по данным эксплуатации газовой залежи начальные дренируемые запасы газа, начальное пластовое давление, балансовые и извлекаемые запасы, коэффициент газоотдачи и начальный газонасыщенный поровый объём. Исходные данные приведены в приложении Б таблицах Б.1 - Б.3.
Порядок расчета.
Рассчитываем псевдокритические параметры смеси по формулам (15) – (16). Результаты расчета заносим в таблицу 1.
Таблица 1 – Результаты расчета псевдокритических параметров смеси
Компонент Мольная доля, ηi РкрТкрηi·Ркр ηi·Ткр
СН4С2Н6 С3Н8 nС4Н10 iС4Н10 n C5 H12 Сумма Рассчитываем приведенное пластовое давление , используя рисунок 1 или формулы (12) - (14). Результаты расчета заносят в нижеприведенную таблицу. Строим графическую зависимость от . Проводим аппроксимирующую прямую, определяем приведенное пластовое давление на начало разработки , балансовые и извлекаемые запасы.Таблица 2 – Результаты расчета приведенного пластового давления
№ замера P РпрТпрz P/z Qстдоб1 N Рассчитываем коэффициенты a и b линейной зависимости (3) по формулам (4) – (9). Результаты расчета заносим в таблицу 3.
Таблица 3 – Результаты промежуточных расчетов произведений и сумм
№ замера
1 N Сумма Строим графические зависимости z от , от .
Определяем начальное пластовое давление , используя графическую зависимость от или по методике, изложенной в приложении А.
Рассчитываем газонасыщенный поровый объем из формулы (5)
Рассчитываем пластовое давление к концу разработки по формуле (11).
Определяем коэффициент конечной газоотдачи по формулам (10) .
Практическое занятие №2 по теме «Расчет параметров «средней» скважины»Цель: изучить метод определения параметров «средней» скважины. Рассчитать параметры «средней» скважины.
В некоторых методах определения показателей разработки месторождений природных газов используется понятие средней скважины, т.е. расчеты выполняются на среднюю скважину. Принимается, что «средняя» скважина имеет среднюю глубину, среднюю длину шлейфа, среднюю конструкцию, средние допустимые дебит и депрессию, средние коэффициенты фильтрационных сопротивлений А и В.
Если на месторождении имеется значительное число скважин, то параметры «средней» скважины можно определить на основе методов статистики и теории вероятностей. Однако из-за недостаточного объема информации при составлении проектов разработки часто используют другой метод, рассмотренный ниже.
Пусть на месторождении имеется N газовых скважин. По результатам исследований этих скважин определены уравнения притока газа к каждой скважине и допустимые дебиты (депрессии). Тогда параметры «средней» скважины рассчитываются следующим образом [8].
(21)
(22)
(23)
(24)
,(25)
,(25*)
где Аср и Вср – коэффициенты фильтрационных сопротивлений «средней» скважины, [Аср]=; ;аi ,bi – коэффициенты фильтрационных сопротивлений i-й скважины;
– дебит i-й скважины, тыс. м3/сут;
δi –депрессия на пласт в i-й скважине, МПа;
N – количество скважин, шт;
Рн – начальное пластовое давление, МПа;
– дебит «средней» скважины, тыс. м3/сут;
– депрессия «средней» скважины.
Задача 2. Рассчитать параметры «средней» скважины. Исходные данные приведены в приложении Б (таблица Б.4). Начальное пластовое давление определяется в задаче 1.
Порядок расчета линейный и не требует пояснений. Результаты расчета заносятся в таблицу 4.
Таблица 4 – Результаты промежуточных расчетов для определения параметров «средней» скважины
№ скважины AiBiδi qi Aiqi Biq2i (2Рн-δi)δi1 N Сумма Необходимо отметить, что дебит «средней» скважины, определенный по формуле (25), и среднеарифметический дебит скважин должны совпадать.
Практическое занятие №3по теме «Методика расчета показателей разработки газовой залежи при газовом режиме для заданной динамики темпа отбора газа и технологического режима эксплуатации «средней» скважины с постоянной депрессией на пласт»Цель занятия: изучить метод расчета показателей разработки газовой залежи при газовом режиме для заданной динамики темпа отбора газа и технологического режима эксплуатации «средней» скважины с постоянной депрессией на пласт. Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и накопленного отбора газа.
Краткая теория вопроса
Характерные периоды разработки месторождений природных газов. С точки зрения динамики отбора газа выделяют три периода: 1) период нарастающей добычи, 2) период постоянной добычи, 3) период падающей добычи.

Рисунок 2 – Динамика основных показателей разработки газового месторождения
Такая тенденция характерна для средних и крупных месторождений. При разработке мелких месторождений может оказаться, что период падающей добычи будет основным. В период нарастающей добычи осуществляется разбуривание месторождения, обустройство промысла и вывод месторождения на постоянную добычу газа. В период постоянной добычи в ряде случае отбирается около половины начальных запасов газа месторождения, продолжается дальнейшее разбуривание месторождения и наращивание мощности ДКС до тех пор, пока это экономически целесообразно. Для периода падающей добычи газа характерно неизменное число добывающих скважин, увеличивается число обводненных и выбывших из эксплуатации скважин. Этот период продолжается по достижению минимального рентабельного отбора из месторождения. При разработке месторождений различают также периоды компрессорной и бескомпрессорной эксплуатации. В настоящее время для дальнего транспорта используются трубы большого диаметра, рассчитанные на давление 7,5 и 5,5 МПа. С точки зрения последовательности выделяют период опытно-промышленной эксплуатации и период промышленной эксплуатации. В период ОПЭ газ подается потребителю и одновременно происходит доразведка месторождения. В этот период запасы категории А и В составляют 20%. В период промышленной эксплуатации основная задача – оптимальное снабжение конкретных потребителей газом и другой продукцией. С точки зрения технологии разработки газоконденсатных залежей выделяются 2 периода: 1) разработка на истощение, 2) период с поддержанием пластового давления (ППД). Рабочие агенты могут нагнетаться в пласт газоконденсатной залежи для предотвращения конденсации и снижения потерь конденсата. Для ППД ГКМ могут быть использованы: 1) вода, 2) сухой газ, 3) азот, 4) углекислый газ. При разработке газоконденсатных месторождений и газоконденсатнонефтяных месторождений может осуществляться рециркуляция сухого газа (сайклинг-процесс). Сайклинг-процесс заключается в добыче всего продукта, его отбензинивания и закачке сухого газа. Различают полный и частичный сайклинг-процесс. При полном сайклинг-процессе ведется добыча конденсата и коэффициент возврата газа в пласт равен 1. При частичном сайклинг-процессе давление не полностью поддерживается и продукцией является газ и конденсат. Это позволяет снизить пластовые потери конденсата.
К показателям разработки месторождений природных газов относится множество параметров, такие как динамика пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и годового отбора газа.
Эти показатели можно определить в результате интегрирования дифференциального уравнения неустановившейся фильтрации газа при соответствующих краевых условиях.
Определять перечисленные показатели разработки газовых месторождений (при некоторых допущениях) можно методом последовательной смены стационарных состояний. Теоретическое обоснование данного метода применительно к проектированию разработки газовых месторождений дано Б.Б. Лапуком.
Введение в расчеты понятия об удельных объемах дренирования существенно их облегчает. Нейтральные линии (поверхности) принимаются как бы непроницаемыми, и каждая скважина дренирует "свой" участок пласта. Такой участок пласта и называют удельным объемом дренирования. Введение понятия об удельных объемах дренирования позволяет проводить расчеты на одну среднюю скважину и определять изменение во времени потребного числа средних скважин.
Распределение давления в пределах удельного объема дренирования принимается таким, как при установившемся (стационарном) притоке газа к скважине. Важный момент эффективного использования метода последовательной смены стационарных состояний состоит в доказательстве положения о том, что при радиальной фильтрации газа к скважине средневзвешенное по газонасыщенному поровому пространству удельного объема дренирования пластовое давление () мало отличается от Рк ( ≈ Рк) на границе удельного объема дренирования радиусом RK (см. рисунок 3). Расчеты показывают, что при расстоянии между скважинами от 600 м до 4400 м и забойном давлении до 0,1 пластового давления (в условиях стационарной фильтрации) давление в удельном объеме дренирования отличается от контурного на 0,5%. При расстоянии между скважинами до 1000 м и почти свободном дебите газовой скважины среднее давление в удельном объеме дренирования отличается от контурного не более чем на 3%. Физически это объясняется значительной крутизной депрессионной воронки при притоке газа к скважине.
Доказательство отмеченного положения позволило в уравнении притока к скважине неизвестное контурное давление (пластовое давление в районе данной скважины) в момент t заменить средним давлением в удельном объеме дренирования, а при равномерном размещении скважин - приближенно средним давлением в залежи в тот же момент:.
Наибольшее распространение в последние годы получили расчетные формулы, использующие двучленную формулу притока газа к забою скважин, обоснованию применимости которой посвящены исследования Е.М. Минского.
РК
rP

Рисунок 3 – Распределение давления в удельном объеме дренирования
Итак, дана зависимость изменения во времени годового отбора газа из месторождения Q = Q(t). Известны запасы газа, начальные пластовые давление и температура, состав газа, допустимый технологический режим эксплуатации средней скважины, уравнение притока газа к средней скважине.
Расчеты показателей разработки будем вести для отборов газа из месторождения и дебитов скважин, приведенных к стандартным условиям.
Требуется определить изменение во времени средневзвешенного пластового и забойного давлений, дебита и потребного числа скважин. Определение этих показателей разработки газового месторождения методом последовательной смены стационарных состояний сводится к решению системы уравнений:
материального баланса для газовой залежи (2);
технологического режима эксплуатации скважины (26);
притока газа к забою скважины (27);
связи потребного числа газовых скважин, отбора газа из месторождения Q и дебита газовой скважины (28).
Изменение во времени средневзвешенного пластового давления определяется по уравнению материального баланса для газовой залежи применительно к газовому режиму. Забойное давление и дебит скважины определяются исходя из совместного решения системы уравнений, описывающих технологический режим эксплуатации скважины и приток газа к забою скважины. Потребное количество скважин определяется по формуле (28). В нашем случае расчеты будем проводить для эксплуатации скважины с постоянной депрессией на пласт, а приток газа к забою скважины описывается двучленной формулой.
.(26)
(27)
(28)
Расчеты показателей разработки месторождений природных газов, основанные на методе последовательной смены стационарных состояний, отличаются значительной простотой. Однако необходимо иметь в виду, какой ценой достигается эта простота, иметь представление об области применимости рассмотренной методики. При использовании метода последовательной смены стационарных состояний не учитывается неоднородность продуктивных отложений по коллекторским свойствам. Расчеты ведутся на среднюю скважину - на скважину со средним дебитом, с коэффициентами фильтрационных сопротивлений А и В, при средних, например, допустимых депрессиях на пласт.
Методика определения показателей разработки месторождений природных газов, изложенная в данном параграфе, предполагает, что в начальный момент скважины размещаются равномерно по площади газоносности. Несмотря на то, что действующий фонд газовых скважин увеличивается во времени, принимается допущение о равномерности сетки скважин в каждый момент. В определенной мере это оправдывается тем, что коэффициенты фильтрационных сопротивлений А и В не претерпевают больших изменений при существенных изменениях удельных объемов дренирования вследствие, например, добуривания новых скважин.
Методика не учитывает произвольность конфигурации месторождения и расположения скважин, различия продуктивных характеристик проектных газовых скважин.
Вследствие неоднородности пласта по коллекторским свойствам в результате неравномерного расположения скважин и неравномерного дренирования залежи могут возникать значительные общие депрессионные воронки.
Разнодебитность газовых скважин может существенно влиять на систему обустройства газового промысла. Возникающие в процессе разработки месторождения глубокие депрессионные воронки приводят к необходимости более раннего (по сравнению с расчетным) ввода дожимной компрессорной станции и установок искусственного холода.
Использование метода последовательной смены стационарных состояний для определения показателей разработки газовых месторождений целесообразно при ориентировочных, оценочных расчетах. Рассматриваемая методика широко применяется на начальных этапах проектирования разработки месторождений природных газов, когда из-за недостаточных количества и достоверности исходной информации не представляется возможным использовать более точные и совершенные расчетные методы. Этот метод заслуживает внимания при проведении технико-экономических расчетов для установления перспектив разработки газоносной провинции, при отыскании принципиальных решений вопросов разработки месторождения и обустройства промысла. Технико-экономические расчеты характеризуются большим числом вариантов, большим объемом вычислительной работы. Поэтому использование рассматриваемого метода целесообразно и при проведении технико-экономических расчетов на современных ЭВМ.
После нахождения принципиальных решений более точные расчеты проводятся при помощи соответствующих точных методов с учетом дополнительных геолого-промысловых факторов [8].
Задача 3. Разрабатывается круговая газовая залежь при газовом режиме с заданной динамикой темпа отбора газа на периоды нарастающей и постоянной добычи. Продуктивные отложения принимаются недеформируемыми, однородными по коллекторским и емкостным свойствам. Размещение скважин по площади залежи – равномерное. Технологический режим эксплуатации «средней» скважины – постоянная депрессия на пласт. Коэффициенты фильтрационных сопротивлений не изменяются во времени. Исходные данные: балансовые запасы, начальное пластовое давление, пластовая температура, состав и псевдокритические параметры пластового газа, параметры «средней» скважины, давление конца разработки известны по результатам решения задач 1 и 2. Темп отбора в период постоянной добычи, продолжительность периода нарастающей добычи и доля накопленной добычи газа за период нарастающей и постоянной добычи приведены в приложении Б (таблица Б.5).
Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и накопленного отбора газа с использованием метода последовательной смены стационарных состояний. Расчеты произвести по временным шагам, i – номер временного шага (i – й момент времени).
Порядок расчета
Определяем годовой отбор при стандартных условиях на i – й момент времени (задается на периоды нарастающей и постоянной добычи в % от балансовых запасов).
Определяем накопленную добычу на i – й момент времени путем суммирования годовых отборов.
Определяем отношение на i – й момент времени по уравнению материального баланса (2).
Определяем средневзвешенное пластовое давление по методике, приведенной в приложении А, или графику от.
По уравнению связи (26) определяем забойное давление для «средней» скважины на i – й момент времени.
По формуле (27) определяем дебит «средней» скважины на i – й момент времени.
Определяем требуемое число скважин на i – й момент времени по формуле (28).
Проверяем условие . Если условие выполняется, то повторяем все расчеты по пунктам 1 – 7, иначе переходим к расчету показателей разработки на стадии падающей добычи.
Определяем в j-м приближении (в первом приближении j=1) накопленную добычу газа по формуле
.(29)
Определяем отношение по уравнению материального баланса (2).
Определяем средневзвешенное пластовое давление по методике, приведенной в приложении А, или графику от.
По формуле (27) определяем дебит «средней» скважины .
Определяем в (j+1)-м приближении накопленную добычу газа по формуле (число скважин в период падающей добычи постоянно и равно числу скважин на последний год периода постоянной добычи)
.(30)
Определяем отношение по уравнению материального баланса (2).
Определяем средневзвешенное пластовое давление по методике, приведенной в приложении А, или графику от.
Проверяем условие
.(31)
Если условие выполняется, то накопленную добычу считают рассчитанной и переходят к следующему пункту, иначе расчеты повторяют по пунктам 12 – 15.
По уравнению связи (26) определяем забойное давление для «средней» скважины на i – й момент времени.
Годовой отбор газа определяем как разность накопленных отборов.
По формуле (27) определяем дебит «средней» скважины.
Проверяется условие
.(32)
Если условие (32) выполняется, то повторяем расчеты по пунктам 9 – 20, иначе расчет показателей разработки газовой залежи при газовом режиме для заданной динамики темпа отбора газа и технологического режима эксплуатации «средней» скважины с постоянной депрессией на пласт считается завершенным.
Результаты расчета заносятся в таблицу 5. На рисунках приводится динамика показателей разработки.
Таблица 5 – Результаты расчета показателей разработки газовой залежи при газовом режиме
год , % , млн м3/год , млн м3 , МПа , МПа , МПа Qср, тыс м3/сутNскв, шт1 N Практическое занятие №4 по теме «Методика расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории укрупненной скважины»Цель: Изучить методику расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории укрупненной скважины. Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и накопленного отбора газа, темпы внедрения и суммарное количество внедрившейся воды в газонасыщенный поровый объем залежи.
Краткая теория вопроса
Разработка большинства газовых и газоконденсатных месторождений протекает при определенной степени активности напорного режима. В результате происходит обводнение скважин, микро- и макрозащемление газа в обводненной зоне.
Имеется множество публикаций по теории проектирования и практике разработки месторождений газа при водонапорном режиме.
В теории водонапорного режима газовых месторождений имеющиеся публикации можно подразделить на следующие группы:
Экспериментальные исследования особенностей вытеснения газа водой из пористых, трещиновато-пористых и трещиноватых коллекторов.
Прогнозирование показателей разработки при аппроксимации газовой залежи укрупненной скважиной. В этой группе развиваются методы:
позволяющие рассчитать продвижение в залежь контурной или подошвенной воды;основанные на точной или приближенной теории укрупненной скважины;
учитывающие или пренебрегающие некоторыми характерными особенностями проявления водонапорного режима;
базирующиеся на решении краевой задачи с подвижной границей «газ-вода».
Методы второй группы позволяют рассчитывать количество поступающей в залежь воды без дифференциации ее распределения по площади газоносности и толщине пласта. Знание зависимости изменения во времени суммарного количества поступающей в залежь воды Q = QB(t) дает возможность, используя уравнение материального баланса, записанное, например, в виде
,(33)
где – количество внедрившейся воды в залежь;
определить зависимость изменения во времени среднего пластового давления при водонапорном режиме.
3Исследование двумерных и трехмерных задач теории водонапорного режима газовых месторождений. Методы третьей группы не требуют схематизации исходной задачи. Они позволяют учитывать неоднородность продуктивного пласта по коллекторским свойствам, произвольность сетки размещения скважин на площади газоносности, произвольность границ пласта и другие факторы. Поэтому удается проследить закономерности движения границы раздела «газ – вода», установить очередность и темпы обводнения скважин. Очевидно, что такие результаты можно получить, лишь прибегая к численному интегрированию на ЭВМ соответствующих двумерных и трехмерных краевых задач теории водонапорного режима газовых месторождений. Проблема нестационарного конусообразования также относится к двумерным задачам подземной газогидродинамики. Методы третьей группы отличаются друг от друга принятым для решения численным алгоритмом интегрирования многомерных задач с подвижной границей раздела «газ – вода». Однако наиболее существенное различие методов третьей группы связано с самой постановкой исходной краевой задачи.
Вероятностно-статистические подходы. Построение слоистых моделей продуктивного пласта позволяет учесть в прогнозных расчетах важнейший фактор - неоднородность коллекторских свойств в пределах газоносной толщины пласта. В результате удается приближенно учесть динамику избирательного обводнения скважин и их продукции, оценить потребное число резервных скважин.
Решение обратных задач теории водонапорного режима. Без знания параметров водоносного бассейна нельзя осуществлять достоверный прогноз особенностей проявления водонапорного режима. Фактические данные разработки месторождения несут в себе информацию о параметрах водоносного пласта. Методы решения обратных задач позволяют определить или уточнить параметры водоносного пласта на основе фактических показателей разработки.
Регулирование разработки месторождения при водонапорном режиме. Размещением, дебитами и расходами по добывающим, нагнетательным и разгрузочным скважинам, отборами газа из месторождения или отдельных горизонтов можно достичь тех или иных технико-экономических результатов. Работы рассматриваемого направления исследований посвящены созданию методик решения соответствующих задач регулирования разработки отдельных залежей газа или месторождения в целом.
Теория укрупненной скважины Ван Эвердингена и ХерстаПри исследовании проявления водонапорного режима газовая залежь часто аппроксимируется укрупненной скважиной. На теории укрупненной скважины основаны методики прогнозирования показателей разработки при водонапорном режиме. Кроме того, на ее основе предложены методики уточнения параметров водоносного пласта.
Допустим, укрупненная скважина радиусом Rз дренирует однородный по коллекторским свойствам водоносный пласт с постоянным во времени дебитом воды qB. Согласно решению Ван Эвердингена и Херста, изменение во времени давления Р(Rз) на стенке укрупненной скважины определяется по следующему уравнению:
,(34)
здесь, h, k и – толщина и коэффициенты проницаемости и пьезопроводности водоносного пласта соответственно;
μв – коэффициент динамической вязкости воды;
– табулированная функция параметра Фурье .
Пусть укрупненная скважина эксплуатируется с постоянным во времени противодавлением ΔP = Pн - P(R3) на водоносный пласт. Для вычисления суммарного количества воды QB, которое поступит в залежь к моменту t, Ван Эвердингеном и Херстом получено выражение
,(35)
где – табулированная функция параметра Фурье .
Таблицы функций и составлены для случаев бесконечного по протяженности, конечного замкнутого и открытого водоносного пласта.
В качестве бесконечного водоносный пласт может рассматриваться при условии Rк/R3 ≥20, где Rк – радиус внешней границы пласта.
Теория укрупненной скважины наиболее актуальна для газовых и для газоконденсатных месторождений, так как газовые месторождения разрабатываются в режиме истощения пластовой энергии, а большинство газоконденсатных месторождений также разрабатывается без поддержания пластового давления и рано или поздно они переходят на режим истощения пластовой энергии. В результате между водонапорным бассейном и залежью газа происходит увеличение разницы давлений, что вызывает приток воды в залежь газа (укрупненную скважину).
Решения (34) и (35), полученные для случаев qB=const и ΔP = const, и используются, благодаря принципу суперпозиции, для переменных во времени граничных условий на стенке укрупненной скважины.
Приближенная теория укрупненной скважины
С использованием метода интегральных соотношений получено решение для случая эксплуатации укрупненной скважины при переменном во времени дебите воды в виде системы уравнений, в которую входит уравнение материального баланса (33), уравнение притока воды к «укрупненной» скважине (36), уравнение допустимого технологического режима эксплуатации «средней» скважины (26), уравнение притока газа к «средней» скважине (27), уравнение связи потребного числа газовых скважин, отбора газа из месторождения Q и дебита газовой скважины (28).
(36)
(37)
(38)
(39)
(40)
(41)
С целью единообразия расчетных операций для всех периодов добычи газа остановимся на определении показателей разработки в период падающей добычи с использованием метода последовательных приближений по добытому количеству газа. Пусть на момент t-Δt периода падающей добычи газа все показатели процесса поступления в залежь контурной или подошвенной воды известны. Тогда в первом приближении добытое количество газа из залежи ко времени t оцениваем по формуле
.(42)
Это позволяет оценить показатели разработки на момент времени t периода падающей добычи (с учетом обводнения скважин). Во втором и последующих приближениях добытое количество газа из залежи ко времени t уточняется по формуле
.(43)
При этом последовательные приближения проводятся в пределах рассматриваемого интервала времени периода падающей добычи [8].
Задача 4. Разрабатывается газовая залежь при упруговодонапорном режиме с заданной динамикой темпа отбора газа на периоды нарастающей и постоянной добычи. Газовая залежь аппроксимируется укрупненной скважиной. Продуктивные отложения принимаются недеформируемыми, однородными по коллекторским и емкостным свойствам. Размещение скважин по площади залежи – равномерное. Технологический режим эксплуатации «средней» скважины - постоянная депрессия на пласт. Коэффициенты фильтрационных сопротивлений не изменяются во времени. Исходные данные: балансовые запасы, газонасыщенный поровый объем, начальное пластовое давление, пластовая температура, состав и псевдокритические параметры пластового газа, параметры «средней» скважины, давление конца разработки известны по результатам решения задач 1 и 2. Темп отбора в период постоянной добычи, продолжительность периода нарастающей добычи, доля накопленной добычи газа за период нарастающей и постоянной добычи приведены в приложении Б (таблица Б.5). Пористость m=0,17, начальная и остаточная газонасыщенность =0,82, =0,23, коэффициент пьезопроводности =1 м2/с, коэффициент проницаемости k = 0,22 мкм2, относительная фазовая проницательность =0,2, вязкость воды в пластовых условиях = 0,6 мПа∙с, толщина пласта h=12 м.
Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин, накопленного отбора газа, темпа внедрения воды и количества внедрившейся воды на периоды нарастающей и постоянной добычи с использованием метода последовательной смены стационарных состояний и метода последовательных приближений. Расчеты произвести по временным шагам i – номер временного шага (i – й момент времени).
Порядок расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории укрупненной скважины
Рассчитываем величины , на нулевой момент времени (i=0) (эксплуатация залежи еще не началась), принимаем, что темп внедрения воды и суммарный объем внедрившейся воды на нулевой момент времени равны нулю, средневзвешенное пластовое давление равно начальному пластовому давлению.
Переходим к расчету на (i+1) момент времени. Присваиваем величине j значение ноль (j-число итераций на каждый момент времени).
Рассчитываем накопленную добычу газа на (i+1) момент времени.
Переходим к расчету показателей на (j+1) итерации. При расчете показателей на 2-й и последующих итерациях на (i+1) момент времени переходим к пункту 8.
При расчете на 2-ой и последующие моменты времени и 1-й итерации переходим к пункту 7.
При расчете на 1-й итерации 1-го момента времени производим последовательно расчет отношения Р/z по формуле (33) при =0, по методике, приведенной в приложении А, или из графика Р от Р/z, R(t)=Rз, по формуле (44), по формуле (37), по формуле (36), по формуле (45)
(44)
,(45)
где – продолжительность периода времени, с.
Присваиваем давлению сравнения Рсрав значение начального пластового давления. Переходим к пункту 4.
Присваиваем величине темпа внедрения воды в залежь значение, равное темпу внедрения воды на прошлый момент времени . Рассчитываем суммарный объем внедрившейся воды по формуле (42)
. (46)
Рассчитываем отношение Р/z по формуле (33), по методике, приведенной в приложении А, или из графика Р от Р/z. Присваиваем давлению сравнения Рсрав значение . Переходим к пункту 4.
На 2-й и последующих итерациях на (i+1) момент времени производят расчет по формуле (38), по формуле (37), R(t) по формуле (40), по формуле (36), по формуле (39), Р/z по формуле (33), по методике, приведенной в приложении А, или из графика Р от Р/z
Проверяем условие (47). Если условие (47) выполняется, то переходят к пункту 10, иначе присваивают давлению сравнения величину и переходят к пункту 4.
.(47)
Рассчитываем на (i+1) момент времени забойное давление Рс по формуле (26), средний дебит скважины по формуле (27) и потребное количество скважин по формуле (28).
Проверяем условие (48). Если условие выполняется, то переходим к пункту 2 для расчета на следующий момент времени. Иначе расчет показателей разработки закончен
(48)
Результаты расчета заносятся в таблицу 6. На рисунках приводится динамика показателей разработки.
Таблица 6 – Результаты расчета показателей разработки газовой залежи при упруговодонапорном режиме
год , % , млн м3/год , млн м3 qв, тыс м3/сутQв, млн м3 , МПа , МПа , МПа Qср, тыс м3/сутNскв, штМетодики, основанные на точных решениях Ван Эвердингена и Херста теории «укрупненной» скважины по отношению к реальным процессам, являются приближенными. Это связано, в частности, с использованием упрощенного уравнения материального баланса, неучетом некоторых особенностей поведения защемленного газа в обводненной зоне пласта. Существует более строгая, но более громоздкая методика расчёта, однако расчёты и сопоставления подтвердили практическую приемлемость изложенных методик расчета до отбора около 50 % запасов газа. В прогнозных расчётах на дальнюю перспективу целесообразно пользоваться более строгой, но более громоздкой методикой расчета.
В изложенных методиках, в частности в формуле (40), используются коэффициенты и . Здесь требования к точности коэффициентов и не такие жесткие, так как они входят в подкоренное выражение. Это снижает влияние погрешности их определения на погрешность определения R(t), a R(t) стоит под знаком логарифма [8].
Практическое занятие №5
по теме «Методика расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории Фетковича»
Цель: Изучить методику расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории Фетковича. Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и накопленного отбора газа, темпы внедрения и суммарное количество внедрившейся воды в газонасыщенный поровый объем залежи.
Краткая теория вопроса
Расчет технологических показателей разработки газовых залежей при упруговодонапорном режиме достаточно сложный и трудоемкий процесс. Наиболее точные значения темпа внедрения пластовой воды в залежь можно получить, используя теорию неустановившегося притока воды Ван Эвердингена и Херста. В основе решения Ван Эвердингена и Херста лежит принцип суперпозиции. Указанный метод расчета притока пластовой воды довольно громоздкий и требует постоянного обращения к различным номограммам и таблицам. Более простой метод выполнения расчетов притока воды в залежь без применения принципа суперпозиции был предложен Фетковичем в 1971 г. Данный метод предусматривает моделирование притока воды из водоносной области в залежь точно таким же образом, как и притока нефти из залежи в скважину. При этом используется уравнение притока в виде:
qв=dQвdt=КпрPв-Pгвк (49)
где qв – темп внедрения воды;
Qв – суммарный объем внедрившейся пластовой воды в залежь;
Кпр – коэффициент продуктивности водоносной области;
Pгвк – давление на газоводяном контакте;
Pв – среднее давление в водоносной области.
Суммарная сжимаемость водоносной области β определяется из соотношения:
β=1Q0QвPн-Pв (50)
где Q0 – начальный объем воды в водоносной области;
Pн – начальное пластовое давление.
Уравнение (50) можно переписать в виде:
Qв=Qв max1-PвPн (51)
где Qв max=βQ0Pн – максимально возможный объем притока воды в залежь из водоносной области.
Продифференцировав обе части уравнения (51) получим следующее:
dQвdt=Qв maxPн∙dPвdt (52)
Подставляя уравнение (52) в (49) и интегрируя его, Феткович получил следующее выражение
Qв=Qв maxPнPн-Pгвк1-e-КпрPнtQв max (53)
Уравнение (53) получено при условии постоянства давления на внутренней границе залежи. Для использования этого уравнения в практических расчетах, когда давление на границе постоянно изменяется во времени, следует применить принцип суперпозиции. Однако Феткович показал, что можно использовать уравнение (53) в другой записи, благодаря чему отпадает необходимость в суммировании решений. Рассматриваемый период времени необходимо разбить на равные интервалы времени ∆t. Для притока на -ом интервале уравнение (53) можно записать так:
∆Qв i=Qв maxPнPв i-1-Pгвк i1-e-КпрPн∆tQв max (54)
где
Pв i-1=Pн1-j=1i-1∆Qв jQв max (55)
Pгвк i=Pгвк i-1+Pгвк i2 (56)
Феткович показал, что, последовательно применяя уравнения (54) и (55), можно рассчитать приток воды при различных геометриях водоносной области.
Формулы для расчета коэффициента продуктивности водоносной области, зависящего и от геометрии и от условий фильтрации, приведены в таблице 7.
Таблица 7 – Коэффициенты продуктивности водоносной области
Режим фильтрации Водоносная область круговой геометрии Кпр, м3/(Па·с) Водоносная область линейной геометрии Кпр, м3/(Па·с)
Квазиустановившаяся фильтрация (с выражением депрессии в виде Pв-Pгвк) 2πkhμlnRкRз-343khBμLУстановившаяся фильтрация (с выражением депрессии в виде Pн-Pгвк) 2πkhμlnRкRзkhBμLСледует отметить что уравнения Фетковича используются совместно с выражениями Кпр для квазиустановившейся фильтрации.
Квазиустановившийся режим фильтрации наступает, когда газовая залежь разрабатывается достаточно долго, для того чтобы на приток воды в газоносную область начала оказывать влияние внешняя граница водоносного пласта. Если говорить о модели плоскорадиального притока, то такой случай изображен на рисунке 4. Принимается, что внешняя граница пласта непроницаема и приток воды в пласт отсутствует.
Квазиустановившаяся фильтрация характеризуется одинаковым значением производной давления по времени ∂P∂t=const для всех координат пласта, находящихся в диапазоне от Rз до Rк, при постоянном значении темпа внедрения пластовой воды в залежь dQвdt=const и нулевым значением градиента давления на внешней границе пласта ∂P∂r=0 при r=Rк.
Несмотря на быстроту расчетов и точность получаемых методом Фетковича результатов у нее есть один большой недостаток. Метод Фетковича можно использовать только в случае ограниченной водоносной области небольшого размера (при RкRз<6...8). При бóльших размерах водоносной области начальный приток воды в залежь достаточно длительное время будет происходить при неустановившемся режиме фильтрации, пока внешняя граница залежи не начнет оказывать влияние на давление на границе раздела газ-вода. В этом случае на протяжении нескольких начальных интервалов времени следует использовать теорию неустановившегося притока Ван Эвердингена и Херста, а затем применить метод Фетковича.

Рисунок 4 – Схема квазиустановившегося плоскорадиального притока
Задача 5. Разрабатывается газовая залежь при упруговодонапорном режиме с заданной динамикой темпа отбора газа на периоды нарастающей и постоянной добычи. Газовая залежь аппроксимируется укрупненной скважиной. Продуктивные отложения принимаются однородными по коллекторским и емкостным свойствам. Размещение скважин по площади залежи – равномерное. Технологический режим эксплуатации «средней» скважины - постоянная депрессия на пласт. Коэффициенты фильтрационных сопротивлений не изменяются во времени. Исходные данные: балансовые запасы, газонасыщенный поровый объем, начальное пластовое давление, пластовая температура, состав и псевдокритические параметры пластового газа, параметры «средней» скважины известны по результатам решения задач 1 и 2. Темп отбора в период постоянной добычи, продолжительность периода нарастающей добычи, доля накопленной добычи газа за период нарастающей и постоянной добычи Кпд приведены в приложении Б (таблица Б.5). Пористость m=0,20, начальная и остаточная газонасыщенность α=0,75, αост=0,22, коэффициент проницаемости k=0,22 мкм2, вязкость воды в пластовых условиях μв= 0,6 мПа∙с, толщина пласта h=12 м, коэффициенты сжимаемости пористой среды пласта и пластовой воды βп=4·10-10 Па-1, βв=6·10-10 Па-1, радиус внешней границы пласта Rк=23000 м.
Рассчитать динамику средних пластовых давлений в водоносной и газонасыщенной зоне и забойного давления, среднего дебита скважин, потребного количества скважин, накопленного отбора газа и количества внедрившейся воды на периоды нарастающей и постоянной добычи с использованием метода последовательной смены стационарных состояний и метода последовательных приближений. Расчеты произвести по временным шагам i – номер временного шага (i – й момент времени).
Порядок расчета показателей разработки газовой залежи при упруговодонапорном режиме с использованием приближенной теории Фетковича:
Рассчитываем величины Rз по формуле (41), Кпр, β по формуле (57), Q0 по формуле (58), Qв max, КпрPн∆tQв max, 1-e-КпрPн∆tQв maxβ=1-mβп+mβв (57)
где βп, βв – коэффициенты сжимаемости пористой среды пласта и пластовой воды.
Q0=π(Rк2-Rз2)mh (58)
На нулевой момент времени (i=0) (эксплуатация залежи еще не началась) принимаем, что суммарный объем внедрившейся воды равен нулю, средневзвешенное пластовое давление в газонасыщенной и водоносной зоне равны начальному пластовому давлению.
Переходим к расчету на (i+1) момент времени. Рассчитываем накопленную добычу газа на (i+1) момент времени.
Присваиваем величине j значение ноль (j-число итераций на каждый момент времени). Рассчитываем отношение Ptz[Pt] по формуле (59), средневзвешенное пластовое давление в газонасыщенной области Pt по методике, приведенной в приложении А, или из графика P от P/z. Присваиваем давлению сравнения Pсрав значение Pt.
Ptz[Pt]=1αΩн-Qв i-1PнZнαΩн-PатTплTст∙Qдоб iст (59)
Переходим к расчету показателей на (j+1) итерации. Рассчитываем количество внедрившейся воды за (i+1) момент времени по формуле (60), суммарный объем внедрившейся воды по формуле (61), отношение Ptz[Pt] по формуле (62), Pt по методике, приведенной в приложении А, или из графика P от P/z.
∆Qв i=Qв maxPнPв i-1-0,5Pi-1+Pсрав1-e-КпрPн∆tQв max (60)
Qв i=Qв i-1+∆Qв i (61)
Ptz[Pt]=1αΩн-Qв iPнZнαΩн-PатTплTст∙Qдоб iст (62)
Проверяем условие (63). Если условие (63) выполняется, то переходят к пункту 7, иначе присваивают давлению сравнения Pсрав значение Pt и переходят к пункту 4.
Pt-Pсрав≤10-5 (63)
Рассчитываем на (i+1) момент времени средневзвешенное давление в водоносной области по формуле
Pв i=Pн1-Qв iQв max (64)
Рассчитываем на (i+1) момент времени R(t) по формуле (40), забойное давление Pс по формуле (26), средний дебит скважины по формуле (27) и потребное количество скважин по формуле (28).
Проверяем условие (65). Если условие выполняется, то переходим к пункту 2 для расчета на следующий момент времени. Иначе расчет показателей разработки закончен
Qдоб(t)<Кпд∙Qбал (65)
Результаты расчета заносятся в таблицу 8. На рисунках приводится динамика показателей разработки.
Таблица 8 – Результаты расчета показателей разработки газовой залежи при упруговодонапорном режиме
год Qст, % Qст, млн м3/год Qстдоб, млн м3 ∆Qв, тыс. м3 Qв, тыс. м3 Ptz[Pt], МПа Pt, МПа zPPв, МПа R(t), кмPс, МПа qср, тыс. м3/сутNскв, шт1 N Практическое занятие №6
по теме «Методика расчета показателей разработки газовой залежи при упруговодонапорном режиме с учетом образования общей депрессионной воронки в зоне разбуривания»
Цель: Изучить методику расчета показателей разработки газовой залежи при упруговодонапорном режиме с учетом образования общей депрессионной воронки в зоне разбуривания. Рассчитать динамику пластового давления в зоне разбуривания и во внешней зоне, забойного давления, среднего дебита скважин, потребного количества скважин и накопленного отбора газа, темпы внедрения и суммарное количество внедрившейся воды в газонасыщенный поровый объем залежи.
Краткая теория вопроса
При использовании приближенных методов расчета количества внедряемой пластовой воды в газовую залежь при упруговодонапорном режиме разработки предполагается допущение о равномерном распределении давления по площади газоносности, что возможно только в случае использования равномерной сетки скважин. На практике с целью предотвращения преждевременного обводнения добывающих скважин и снижения затрат на обустройство месторождения скважины разбуривают в пределах отдельных зон, преимущественно расположенных в центральной части залежи. Особенно это касается разработки морских месторождений с использованием стационарных платформ, когда зона дренирования (разбуривания) локализована вблизи платформы, и месторождений расположенных в районах со сложным рельефом местности и тяжелыми климатическими условиями.
Эксплуатация газовых залежей с сеткой скважин, расположенной преимущественно в центральной части залежи, приводит к образованию общих депрессионных воронок, что отрицательно сказывается на показателях разработки: преждевременное снижение дебитов скважин, необходимость бурения дополнительных скважин для поддержания заданных отборов газа, большие объемы остаточного газа за пределами зоны разбуривания и т.д.
Пусть круговая замкнутая залежь радиуса Rз дренируется расположенной в центре укрупненной скважиной радиуса Rр. Под Rр понимается радиус площади, разбуренной реальными скважинами, которые размещены равномерно в зоне разбуривания.
Если считать, что средневзвешенное давление в зоне разбуривания равно давлению на стенке укрупненной скважины Pрt, то задача сводится к определению зависимости этого давления от времени. Заменяя реальные скважины «средними», дальнейший расчет можно провести по следующей методике.
Уравнение материального баланса для зоны разбуривания запишется следующим образом
Mрt=Mрн-Mдобt+Mперt (66)
для внешней зоны
Mвt=Mвн-Mперt (67)
Сложив (66) и (67), получим уравнение материального баланса для всей залежи
Mрt+Mвt=Mн-Mдобt (68)
Здесь Mн, Mрн, Mвн – начальная масса газа во всей залежи, в зоне разбуривания и во внешней зоне; Mрt, Mвt – масса газа в зоне разбуривания и во внешней зоне на текущий момент времени t; Mдобt – накопленная массовая добыча на текущий момент времени; Mперt – масса газа, перетекшего из внешней зоны в зону разбуривания на текущий момент времени.
Если предположить, что средневзвешенное пластовое давление во внешней зоне равно давлению на контуре залежи Pк(t), то можно записать (66)…(68) через стандартные объемы
R12∙Pрtzрt=R12∙Pнzн-PатTплTстαΩн∙Qдобt-Qперt (69)
R22t-R12∙Pкtzкt=1-R12∙Pнzн-PатTплTстαΩн∙Qперt (70)
R22t-R12∙Pкtzкt+R12∙Pрtzрt=Pнzн-PатTплTстαΩн∙Qдобt (71)
Qдобt=0tQtdt (72)
Qперt=0tqперtdt (73)
где αΩн – газонасыщенный поровый объем пласта;
α – средневзвешенная газонасыщенность;
Pн, Pрt, Pкt – начальное и текущие давления в соответствующих зонах;
zн, zрt=zPрt, zкt=zPкt – начальный и текущие коэффициенты сверхсжимаемости;
Qt – темп отбора газа;
R1=Rр/Rз;
R2(t)=R(t)Rз.
Если предположить, что процесс фильтрации газа во внешней зоне квазистационарный, то тогда интенсивность перетока можно записать так
qперt=πkhTст∙Pк2t-Pр2tzкtμкtPатTпл∙lnR2(t)R1 (74)
где μкt=μPкt – вязкость газа.
Для расчета коэффициента сверхсжимаемости будем использовать зависимость Латонова-Гуревича
z=0,4∙lgTTп.кр+0,73PPп.кр+0,1∙PPп.кр (75)
где Pп.кр, Tп.кр – псевдокритические параметры газа.
Вязкость газа будем рассчитывать по зависимости Дина и Стила с уточнениями О.В. Калашникова и др.
μP,T=μатT+9,77∙10-5ξ∙е1,415ρпр-е-3,046ρпр1,684 (76)
ξ=Tп.кр16∙M-0,5∙Pп.кр0,1013-23 (77)
μатT==166,8∙10-5ξ∙0,1338∙TTп.кр-0,09325/9; TTп.кр≥1,5=34∙10-5ξ∙TTп.кр89; TTп.кр<1,5 (78)
ρпр=ρP,Tρп.кр; ρп.кр=MVп.кр; ρст=M24,01; ρP,T=ρстP∙TстzP,T∙Pат∙T (79)
где M – молярная масса газа, кг/кмоль;
Vп.кр – псевдокритический мольный объем, м3/кмоль;
μатT – вязкость газа при атмосферном давлении, мПа·с;
ξ - параметр.
В формулах (76)…(79) единица измерения давления – МПа, температуры – К, вязкости – мПа·с. Для метана погрешность расчета вязкости по сравнению с экспериментальными данными при 300 К в диапазоне давлений 10…60 МПа не превышает 2,7 %, при 400 К – 2,3 %.
Темп внедрения воды во внешнюю зону и, следовательно, положение газоводяного контакта R(t) могут быть вычислены по методу, изложенному в практической работе №4 или №5, либо по любому другому известному методу.
Данный метод расчета, также можно использовать и при газовом режиме разработки залежи. В этом случае принимается, что величина R(t) со временем не изменяется и равна
Rt=Rз или R2(t)=1Задача 6. Разрабатывается газовая залежь при упруговодонапорном режиме с заданной динамикой темпа отбора газа на периоды нарастающей и постоянной добычи. Газовая залежь аппроксимируется укрупненной скважиной. Продуктивные отложения принимаются однородными по коллекторским и емкостным свойствам. Размещение скважин в зоне разбуривания – равномерное. Технологический режим эксплуатации «средней» скважины - постоянная депрессия на пласт. Коэффициенты фильтрационных сопротивлений не изменяются во времени. Исходные данные: балансовые запасы и газонасыщенный поровый объем принимаем в 10 раз большими по сравнению с полученными значениями при решении задачи 1. Начальное пластовое давление, пластовая температура, состав и псевдокритические параметры пластового газа, параметры «средней» скважины известны по результатам решения задач 1 и 2. Темп отбора в период постоянной добычи, продолжительность периода нарастающей добычи, доля накопленной добычи газа за период нарастающей и постоянной добычи Кпд приведены в приложении Б (таблица Б.5). Пористость m=0,17, начальная и остаточная газонасыщенность α=0,78, αост=0,25, коэффициент пьезопроводности κ=1 м2/с, коэффициент проницаемости k=0,14 мкм2, относительная фазовая проницательность kв*=0,2, вязкость воды в пластовых условиях μв= 0,6 мПа∙с, толщина пласта h=16 м, величина R1=0,6.
Рассчитать динамику пластового давления в зоне разбуривания и во внешней зоне, забойного давления, среднего дебита скважин, потребного количества скважин, накопленного отбора газа, темпа внедрения воды и количества внедрившейся воды на периоды нарастающей и постоянной добычи с использованием метода последовательной смены стационарных состояний и метода половинного деления. Расчеты произвести по временным шагам i – номер временного шага (i – й момент времени).
Порядок расчета показателей разработки газовой залежи при упруговодонапорном режиме с учетом образования общей депрессионной воронки в зоне разбуривания:
Рассчитываем величины Rз по формуле (41), ξ по формуле (77), μатT по формуле (78), ρст и ρкр по формуле (79), zн по формуле (75) и μн по формуле (76) при P=Pн.
На нулевой момент времени (i=0) (эксплуатация залежи еще не началась) принимаем, что темп внедрения воды, суммарный объем внедрившейся воды и переток газа из внешней зоны в зону разбуривания равны нулю, средневзвешенное пластовое давление в газонасыщенной зоне равно начальному пластовому давлению.
Переходим к расчету на (i+1) момент времени. Присваиваем величине j значение ноль (j-число итераций на каждый момент времени).
Рассчитываем накопленную добычу газа на (i+1) момент времени.
Переходим к расчету показателей на (j+1) итерации. При расчете показателей на 2-й и последующих итерациях на (i+1) момент времени переходим к пункту 9.
При расчете на 2-ой и последующие моменты времени и 1-й итерации переходим к пункту 8.
При расчете на 1-й итерации 1-го момента времени в качестве начального условия можно принять, что Qпер i=0,1∙Qдоб i, а R2 i=1. Далее последовательно производим расчет отношения Pрizрi по формуле (69), Pрi по методике, приведенной в приложении А, или из графика P от P/z.
Рассчитываем отношение Pкizкi по формуле (70), Pкi по методике, приведенной в приложении А, или из графика P от P/z, рассчитываем zкi по формуле (75), μкi по формуле (76), fo по формуле (44), Pfo по формуле (37), qвi по формуле (80), Qвi по формуле (45), Ri по формуле (40), qперi по формуле (74), Qперi по формуле (81), производим расчет отношения Pрizрi по формуле (69), Pрi по методике, приведенной в приложении А, или из графика P от P/z.
qвi=2πkhμв∙Pн-PкiPfoi+1kв*∙lnRкRi (80)
Qперi=qперi∆t2 (81)
Присваиваем давлению сравнения Pсрав(j) значение Pрi. Переходим к пункту 5.
Присваиваем величине темпа внедрения воды в залежь qв i значение, равное темпу внедрения воды на прошлый момент времени qв i-1. Рассчитываем суммарный объем внедрившейся воды по формуле (46), Ri по формуле (40).
Присваиваем величине темпа перетока газа в зону разбуривания qперi значение, равное темпу перетока газа на прошлый момент времени qперi-1. Рассчитываем суммарное количество перетекшего газа по формуле (82). производим расчет отношения Pрizрi по формуле (69), Pрi по методике, приведенной в приложении А, или из графика P от P/z.
Присваиваем давлению сравнения Pсрав(j) значение Pрi. Переходим к пункту 5.
На 2-й и последующих итерациях на (i+1) момент времени присваиваем давлению Pрi значение давления сравнения Pсрав(j-1) и последовательно производим расчет отношения Pкizкi по формуле (71), Pкi по методике, приведенной в приложении А, или из графика P от P/z, zкi по формуле (75), μкi по формуле (76), fo по формуле (38), Pfo по формуле (37), qвi по формуле (80), Qвi по формуле (39), Ri по формуле (40), qперi по формуле (74), Qперi по формуле (82), производим расчет отношения Pрizрi по формуле (69), Pрi по методике, приведенной в приложении А, или из графика P от P/z.
Qперi=Qперi-1+qперi-1+qперi∆t2 (82)
Проверяем условие (83). Если условие (83) выполняется, то переходят к пункту 11, иначе рассчитывают новое значение давления сравнения Pсрав(j) по формуле (84) и переходят к пункту 5.
Pрi(j)-Pсрав(j-1)≤10-5 (83)
Pсрав(j)=Pсрав(j-1)+Pрi(j)2 (84)
Рассчитываем на (i+1) момент времени забойное давление Pс по формуле (26), средний дебит скважины по формуле (27) и потребное количество скважин по формуле (28).
Проверяем условие (84). Если условие выполняется, то переходим к пункту 3 для расчета на следующий момент времени. Иначе расчет показателей разработки закончен
Qдоб(t)<Кпд∙Qбал (84)
Результаты расчета заносятся в таблицу 9. На рисунках приводится динамика показателей разработки.
Таблица 9 – Результаты расчета показателей разработки газовой залежи при упруговодонапорном режиме
год Qст, % Qст, млн м3/год Qстдоб, млн м3 qв, м3/с Qв, тыс. м3 Pкzк, МПа Pк, МПа zкμк, мПа·сR, кмqпер, м3/с Qпер, млн м3 Pрzр, МПа Pр, МПа Pс, МПа qср, тыс. м3/сутNскв, шт1 N Практическое занятие №7
по теме «Методика расчета неравномерного распределения давления по площади газовой залежи, разрабатываемой кольцевыми батареями скважин»
Цель: Изучить метод эквивалентных фильтрационных сопротивлений для построения карты распределения давления по площади газовой залежи, разрабатываемой кольцевыми батареями скважин. Рассчитать значения пластовых давлений в зоне дренирования скважины и забойных давлений.
Краткая теория вопроса
В теории и практике разработки месторождений природных газов широкое распространение получили системы размещения скважин в виде кольцевых батарей или цепочек скважин. Количество эксплуатационных скважин на месторождении определяется из условия обеспечения заданного отбора углеводородного сырья. Поэтому в фильтрационных расчетах, связанных с разработкой месторождений, необходимо решать гидродинамические задачи определения давления на забоях скважин при заданных дебитах, или наоборот, дебитов при заданных давлениях.
При решении этих задач нужно учитывать, что при работе нескольких скважин наблюдается их взаимное влияние друг на друга - интерференция скважин.
Поэтому необходимо рассмотреть постановки и решения задач, когда одновременно работают не одна, а группы скважин. Наиболее простые постановки задач получаются в том случае, когда пласт предполагается плоским, а скважины считаются точечными источниками или стоками. При решении подобных задач широко используется предположение о потенциальности течения и метод суперпозиции (потенциала).
Течение называется потенциальным, если существует такая скалярная функция Ф, что градиент от нее равен вектору скорости V, т.е. выполняется равенство
V=grad Ф (85)
при этом скалярная функция Ф называется потенциалом.
Последнее равенство имеет структуру аналогичную закону Дарси
w=-kμgrad P (86)
где P – давление.
В самом деле, если k и μ – константы, то
w=-grad kPμ и w=-grad Ф (87)
где
Ф=kPμ (88)
Поэтому фильтрационные течения в недеформируемых пластах k=const жидкостей с постоянной вязкостью потенциальны.
Потенциал точечного стока на плоскости (добывающей скважины) определяется по следующей формуле:
Ф=qh2πlnr+C (89)
где qh – дебит добывающей скважины, приходящийся на единицу толщины пласта;
r – расстояние от центра скважины до рассматриваемой точки на плоскости;
C – постоянная интегрирования.
Уравнению Лапласа удовлетворяет не только давление, но и введенный потенциал Ф:
∂2Ф∂x2+∂2Ф∂y2=0 (90)
Поскольку уравнение Лапласа линейное и однородное, его решения обладают очень важным свойством: сумма частных решений уравнения и произведение частного решения на константу также является решением. Это свойство позволяет использовать при решении задач метод суперпозиции. Математический смысл метода суперпозиции сводится к тому, что если имеется N фильтрационных потоков с потенциалами
Фi=qhi2πlnr+Ci где i=1, 2,…, N, (91)
каждый из которых удовлетворяет уравнению Лапласа, то и линейная комбинация этих потенциалов Ф=i=1NciФi, где ci – произвольные постоянные, также удовлетворяют уравнению Лапласа (90).
Рассмотрим широко применяемый при проектировании разработки месторождений метод эквивалентных фильтрационных сопротивлений, предложенный Ю.П. Борисовым. Метод основан на аналогии между движением жидкости в пористой среде и течением электрического тока в проводах.
Рассмотрим без вывода решение задачи о притоке жидкости к одной бесконечной цепочке скважин, расположенных на расстоянии 2σ друг от друга и на расстоянии L от прямолинейного контура питания. Будем считать, что на контуре питания потенциал равен Фк а на стенках скважин – Фс (рисунок 5). Требуется определить дебит каждой скважины и суммарный дебит N скважин в цепочке.

Рисунок 5 – Схема прямолинейной цепочки скважин
Решение задачи находится с помощью метода суперпозиции. Цепочка скважин стоков отображается зеркально относительно контура питания в скважины-источники, и рассматривается интерференция двух цепочек скважин в неограниченном пласте. Вдоль прямой АВ, проходящей через скважину-сток и скважину-источник, частицы движутся наиболее быстро, а вдоль прямых А'В', делящих расстояние между скважинами пополам, движение будет наиболее медленным, так как линии, в силу симметрии потока, представляют собой непроницаемые границы.
Дебит каждой скважины может быть вычислен по следующей формуле:
q=2πhФк-ФсπLσ+lnσπrс=Фк-ФсL2σh+12πhlnσπrс (92)
где h - толщина пласта, м;
rс – радиус скважины, м;
или, если ввести обозначение
L2σh=ρ, 12πhlnσπrс=ρ' (93)
формулу (92) можно переписать в виде
q=Фк-Фсρ+ρ' (94)
Соотношение (94) аналогично закону Ома, поэтому Ю.П. Борисов предложил величину ρ называть внешним фильтрационным сопротивлением батареи, а ρ' - внутренним. Таким образом, приток жидкости к цепочке скважин можно представить схемой эквивалентных фильтрационных сопротивлений, показанной на рисунке 6.

Рисунок 6 – Схема фильтрационных сопротивлений при притоке к бесконечной цепочке скважин
Аналогом дебита q служит сила тока, а аналогом разности фильтрационных потенциалов – разность электрических потенциалов.
Дебит одной скважины кольцевой батареи скважин, состоящей из m скважин (рисунок 7), в круговом пласте радиуса Rз может быть вычислен по формуле
q=2πhФк-ФсmlnRзR1+lnσπrс (95)
где R1 – радиус батареи скважин.
При этом внешние и внутренние фильтрационные сопротивления определяются соотношениями
ρ=m2πhlnRзR1 и ρ'=12πhlnσπrс (96)

Рисунок 7 – Схема кольцевой батареи скважин
Суммарный дебит m скважин кольцевой батареи определяется по формуле
Q=qm=Фк-Фс12πhlnRзR1+12πhmlnσπrс (97)
В случае трех кольцевых батарей, соосных круговому контуру питания, приток к скважинам рассчитывается по схеме эквивалентных фильтрационных сопротивлений, представленной на рисунке 8.
Расчет схемы производится аналогично расчету электрических цепей по законам Ома и Кирхгофа. Составляются линейные алгебраические уравнения либо для дебитов Q1, Q2, Q3, либо для потенциалов Фк, Ф1, Ф2, Ф3, Фс1, Фс2, Фс3, в зависимости от того, что задано и что необходимо определить.

Рисунок 8 – Схема фильтрационных сопротивлений при притоке к трем кольцевым батареям скважин
Внешние и внутренние фильтрационные сопротивления определяются соотношениями
ρ1=12πhlnRзR1, ρ2=12πhlnR1R2, ρ3=12πhlnR2R3,
(98)
ρ1'=12πhm1lnσ1πrс, ρ2'=12πhm2lnσ2πrс, ρ3'=12πhm3lnσ3πrсгде R1, R2, R3 – радиусы батарей, м;
m1, m2, m3 – число скважин в батарее.
В случае установившейся фильтрации газа метод суперпозиции можно использовать для потенциалов, определенных через функцию Лейбензона P. Для вектора массовой скорости фильтрации газа выполняется равенство
ρw=-grad Ф* (99)
в котором
Ф*=kPμ (100)
где P – функция Лейбензона.
Для идеального газа функция Лейбензона определяется по формуле
P=ρатP22Pат (101)
Тогда для построения решений с помощью метода суперпозиции при установившейся фильтрации газа можно использовать ранее найденные решения для установившейся фильтрации несжимаемой жидкости. Аналогия между фильтрацией несжимаемой жидкости и газа устанавливается с помощью следующей замены переменных:
для несжимаемой жидкости: для газа:
давление, Pфункция Лейбензона, Pпотенциал, Ф=kPμпотенциал, Ф*=kPμскорость фильтрации, wмассовая скорость фильтрации, ρwдебит, qмассовый дебит, qmсуммарный дебит прямолинейной цепочки или кольцевой батареи, Qсуммарный массовый дебит прямолинейной цепочки или кольцевой батареи, QmМассовый дебит газовой скважины с объемным дебитом при стандартных условиях соотносятся как:
qm=ρcтqат (102)
где ρcт=Mсм24,01 – плотность газа при стандартных условиях, кг/м3;
Mсм – молярная масса газа, г/моль.
Задача 7. Разрабатывается круговая газовая залежь при газовом режиме пятью кольцевыми батареями скважин. Продуктивные отложения принимаются, однородными по коллекторским и емкостным свойствам. Фильтрационные потоки в пласте за рассматриваемый период времени принимаются установившимися. Исходные данные: давление на контуре питания принять равным половине начального пластового давления, а дебит «средней» скважины меньше на одну третью по сравнению с дебитом полученном при решении задачи 1 и 2. Газонасыщенный поровый объем, состав пластового газа и коэффициенты фильтрационных сопротивлений «средней» скважины известны по результатам решения задач 1 и 2. Радиусы батарей и количество скважин в каждой батареи приведены в приложении Б, таблице Б.6. Пористость m=0,18, начальная газонасыщенность α=0,75, коэффициент проницаемости k=0,22 мкм2, вязкость газа в пластовых условиях имеет одинаковое значение μг= 0,02 мПа∙с, толщина пласта h=15 м.
Рассчитать по батареям значения пластового давления в зоне дренирования скважин, забойного давления добывающих скважин и необходимых значений депрессий на пласт.
Порядок расчета.
Рассчитываем величины Rз по формуле (41) и Фк* по формулам (100) и (101).
Составляем схему фильтрационных сопротивлений для притока к -ому количеству кольцевых батарей скважин.
Для каждой батареи рассчитываем величины Ri и σi, значения внешних фильтрационных сопротивлений ρi по формулам (98), суммарные массовые дебиты кольцевых батарей Qm i, значения j=iNQm j.
Последовательно для каждой батареи рассчитываем значения потенциалов Фi* из формулы (103), значение функции Лейбензона Pi и пластового давления в зоне дренирования скважин Pij=iNQm j=Фi-1*-Фi*ρi (103)
Рассчитываем для каждой батареи значения забойного давления в скважинах Pсi из формулы (104) и депрессию на пласт δi по формуле (105)
Pi2-Pсi2=Aqср+Bqср2 (104)
δi=Pi-Pci (105)
Результаты расчета заносятся в таблицу 10. На рисунке приводится распределение пластового давления по площади залежи.
Таблица 10 – Результаты расчета показателей разработки газовой залежи при упруговодонапорном режиме
Батарея Ri, мσi, мρi, м-1 Qm i, кг/сj=iNQm j, кг/сФi*, кг/(м·с) Pi, Па·кг/м3 Pi, МПа Pсi, МПа δi, МПа
1 N Практическое занятие №8
по теме «Методика расчета показателей разработки слоисто-неоднородной газовой залежи при упруговодонапорном режиме»
Цель: Изучить методику расчета показателей разработки слоисто-неоднородной газовой залежи при упруговодонапорном режиме с использованием полосообразной модели пласта. Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и накопленного отбора газа, количества внедрившейся воды в пласт и по каждому пропластку в отдельности, количества обводнившихся скважин по каждому пропластку.
Краткая теория вопроса
Приближенные методы расчета технологических показателей разработки газовых залежей при упруговодонапорном режиме зачастую предполагают использование однородной по фильтрационно-емкостным свойствам модели пласта. Однако опыт разработки показывает, что все залежи по своей структуре неоднородны. Выделяют три типа неоднородности пластов: вертикальная слоистая неоднородность по проницаемости, изменение проницаемости по площади и по различным направлениям в связи с трещиноватостью пород.
При упруговодонапорном режиме разработки слоисто-неоднородной газовой залежи происходит неравномерное продвижение газоводяного контакта, что отрицательно сказывается на технологических показателях разработки. В результате избирательного обводнения залежи происходит микро- и макрозащемление значительных объемов газа, снижается конечная газоотдача продуктивных пластов. По мере разработки образовавшиеся “языки” пластовой воды достигают зоны расположения добывающих скважин. Обводнение интервалов перфорации в скважинах приводит к разрушению призабойной зоны пласта, увеличению выноса пластовой воды и песка на поверхность. По мере снижения дебитов газа в обводняющихся скважинах, наступает момент, когда скорость восходящего потока газа не обеспечивает вынос всей жидкости на поверхность, и часть ее, накапливаясь на забое, может привести к остановке скважины.
Рассмотрим полубесконечный полосообразный пласт толщиной h, шириной B и длиной газоносной зоны L (рисунок 9).

Рисунок 9 – Полосоообразная модель пласта
Для подсчета количества внедрившейся воды в газоносную область необходимо построить слоистую модель пласта (рисунок 10). Под слоисто-неоднородной моделью понимаем пласт, состоящий из n пропластков с различными значениями коэффициентов проницаемости ku, и толщины hu (u =1, 2, … n).

Рисунок 10 – Слоисто-неоднородная полосообразная модель пласта
Для простоты расчетов рассмотрим поршневое вытеснение газа водой. Количество внедрившейся пластовой воды в каждый пропласток Qв ut можно найти, зная расстояние, на которое продвинулся газоводяной контакт от начального положения lut.
Qв ut=αmBhulut (106)
Для приближенных расчетов линейных перемещений контакта в пропластках различной проницаемости и их сопоставления между собой, можно воспользоваться следующим уравнением движения границы раздела газ-вода
lut=akut (107)
Приближенную оценку коэффициента a для поршневого характера вытеснения нетрудно получить при ∆P=const – постоянном перепаде давления на контуре питания и перемещающейся границе раздела воды и газа:
a=2∙∆Pmμв (108)
где m – коэффициент пористости;
μв – вязкость воды.
Выделяют две модели слоистого пласта: с наличием газодинамической связи между пропластками и ее отсутствием. В данной работе будет рассмотрен первый случай, когда значение давления во всех пропластках одинаковое.
P1t=P2t=…=Pnt=Pt (109)
Следовательно, уравнение материального баланса будет иметь следующий вид
Ptz[Pt]=1αΩн-QвtPн∙αΩнzн-PатTплzстTстQдобст(t) (110)
где
Qвt=u=1nQв ut (111)
αΩн=αmhBL (112)
Слоистое строение пласта достаточно слабо сказывается на общем количестве поступившей в залежь воды, однако оно обуславливает избирательное продвижение воды по отдельным пропласткам, в результате чего происходит обводнение интервалов перфорации добывающих скважин. В случае равномерного размещения скважин на площади газоносности число обводнившихся скважин nобв u по -ому пропластку можно найти из следующего соотношения:
nобв u=Sобв uSuNскв (113)
где Sобв u – обводненная площадь -ого пропластка на момент времени t;
Su – площадь газоносности -ого пропластка;
Nскв – общий фонд скважин на момент времени t.
В случае использования полосообразной модели пласта с равномерным размещением скважин в зоне разбуривания (рисунок 10) число обводнившихся скважин nобв u по -ому пропластку можно определить по следующим формулам:
nобв u=0, при lu≤L0(114)
nобв u=lu-L0L-L0Nскв, при lu>L0где L0 - расстояние от начального положения газоводяного контакта до зоны разбуривания.
Задача 8. Разрабатывается слоисто-неоднородная газовая залежь при упруговодонапорном режиме с заданной динамикой темпа отбора газа на периоды нарастающей и постоянной добычи. При расчетах можно использовать полубесконечную полосообразную модель залежи, ограниченную с одной стороны непроницаемой перегородкой, а с другой водоносной областью. Продуктивные отложения в пределах каждого пропластка принимаются однородными по коллекторским и емкостным свойствам. Размещение скважин по площади залежи – равномерное в пределах зоны разбуривания. Технологический режим эксплуатации «средней» скважины - постоянная депрессия на пласт. Коэффициенты фильтрационных сопротивлений не изменяются во времени. Исходные данные: балансовые запасы, газонасыщенный поровый объем, начальное пластовое давление, пластовая температура, параметры «средней» скважины известны по результатам решения задач 1 и 2. Темп отбора в период постоянной добычи, продолжительность периода нарастающей добычи, доля накопленной добычи газа за период нарастающей и постоянной добычи Кпд приведены в приложении Б, таблице Б.5. Характеристики пропластков приведены в приложении Б, таблице Б.7. Пористость m=0,20, начальная газонасыщенность α=0,75, вязкость воды в пластовых условиях μв= 0,74 мПа∙с, ширина пласта B=4000 м, расстояние от начального положения газоводяного контакта до зоны разбуривания L0=1500 м.
Рассчитать динамику пластового и забойного давлений, среднего дебита скважин, потребного количества скважин и накопленного отбора газа, количества внедрившейся воды в пласт и по каждому пропластку в отдельности, количества обводнившихся скважин по каждому пропластку на периоды нарастающей и постоянной добычи с использованием метода последовательной смены стационарных состояний и метода половинного деления. Расчеты произвести по временным шагам i – номер временного шага (i – й момент времени).
Порядок расчета показателей разработки слоисто-неоднородной газовой залежи при упруговодонапорном режиме с использованием осесимметричной полосообразной модели пласта:
Рассчитываем длину газоносной зоны пласта L из формулы (112).
Принимаем, что на нулевой момент времени (i=0) (эксплуатация залежи еще не началась) газоводяной контакт находится в начальном положении, суммарный объем внедрившейся воды равен нулю, средневзвешенное пластовое давление равно начальному пластовому давлению.
Переходим к расчету на (i+1) момент времени. Присваиваем величине j значение ноль (j-число итераций на каждый момент времени).
Рассчитываем накопленную добычу газа на (i+1) момент времени. В первом приближении присваиваем количеству внедрившейся воды значение равное на предыдущий момент времени Qвt=Qвt-∆t.
Переходим к расчету показателей на (j+1) итерации. Рассчитываем отношение Ptz[Pt] по формуле (110), Pt по методике, приведенной в приложении А, или из графика P от P/z. Присваиваем давлению сравнения Pсрав(j) значение Pt.
Переходим к расчету показателей на (j+1) итерации. Определяем положение границы раздела газ-вода для каждого пропластка по формуле (115)
lut=lut-∆t+2Pн-0,5Pt-∆t+Pсрав(j-1)mμвku∆t (115)
Рассчитываем суммарный объем внедрившейся воды в каждый пропласток Qв u и во всю залежь Qв по формулам (106) и (111), отношение Ptz[Pt] по формуле (110), Pt по методике, приведенной в приложении А, или из графика P от P/z.
Проверяем условие (116). Если условие (116) выполняется, то переходят к пункту 9, иначе рассчитывают новое значение давления сравнения Pсрав(j) по формуле (117) и переходят к пункту 6.
Pt(j)-Pсрав(j-1)≤10-5 (116)
Pсрав(j)=Pсрав(j-1)+Pt(j)2 (117)
Рассчитываем на (i+1) момент времени забойное давление Pс по формуле (26), средний дебит скважины по формуле (27) и потребное количество скважин по формуле (28).
Рассчитываем
Проверяем условие (118). Если условие выполняется, то переходим к пункту 3 для расчета на следующий момент времени. Иначе расчет показателей разработки закончен
Qдоб(t)<Кпд∙Qбал (118)
Результаты расчета заносятся в таблицы 11 и 12. На рисунках приводится динамика показателей разработки.
Таблица 11 – Результаты расчета показателей разработки газовой залежи при упруговодонапорном режиме
Год Qст, % Qст, млн м3/год Qстдоб, млн м3 Qв, тыс. м3 Ptz[Pt], МПа Pt, МПа zPс, МПа qср, тыс. м3/сутNскв, шт1 N Таблица 12 – Результаты расчета показателей разработки газовой залежи при упруговодонапорном режиме
Год l1, мnобв 1, штl2, мnобв 1, штl3, мnобв 1, штl4, мnобв 1, штl5, мnобв 1, шт1 N Практическое занятие №9
по теме «Расчет производительности обводняющейся газовой скважины в слоисто-неоднородном пласте»
Цель: изучить метод определения дебита обводняющейся скважины по газу и попутно добываемой пластовой воде в слоисто-неоднородном пласте.
Представим, что скважина вскрывает газоносный пласт, состоящий из n пропластков различной проницаемости k1, k2, … kn и толщиной h1, h2, … hn соответственно, причем k1<k2<…<kn. По результатам испытаний скважины на стационарных режимах известны начальные значения коэффициентов фильтрационных сопротивлений A и B.
Для оценки долевого участия каждого пропластка в общем дебите скважины необходимо определить значения коэффициентов au и bu (u=1, 2, …, n), характеризующих фильтрационные сопротивления пропластков. С этой целью запишем соотношения, следующие из структуры коэффициентов фильтрационных сопротивлений:
au=kсрkufuA (119)
bu=Bfu2 (120)
где fu=huh - вероятность, с которой встречается пропласток проницаемостью ku;
h=u=1nhu - начальная газонасыщенная толщина пласта, м;
kср=u=1nkufu – средневзвешенное значение проницаемости, мкм2.
Дебит скважины qг есть сумма дебитов по каждому пропластку quqг=u=1nqu (121)
где
qг=A24B2+2Pt-δ(t)∙δ(t)B-A2B (122)
qu=au24bu2+2Pt-δ(t)∙δ(t)bu-au2bu (123)
Согласно схеме послойного вытеснения по мере обводнения скважины число газоотдающих пропластков уменьшается, а водопроявляющих, соответственно, увеличивается. Тогда при обводнении пропластка проницаемостью kl (1≤l≤n) текущий дебит газа равен
qt=u=1l-1qu(t) (124)
а дебит воды по этой скважине
qвt=u=lnqв u(t) (125)
Для определения текущего дебита воды воспользуемся следующим соотношением между коэффициентом фильтрационного сопротивления au и коэффициентом продуктивности этого пропластка по воде cв u:
cв u=2Pатkв*μгμв103au (126)
Соотношение (126) получено с учетом размерностей: cв=м3/сутМПа, a=МПа2тыс. м3/сут, Pат=МПа.
При известной депрессии на пласт δ дебит воды по -ому пропластку вычисляется как
qв u(t)=cв u(t)δ(t) (127)
Задача 9. Разрабатывается слоисто-неоднородная газовая залежь при упруговодонапорном режиме с заданной динамикой темпа отбора газа на периоды нарастающей и постоянной добычи. Продуктивные отложения в пределах каждого пропластка принимаются однородными по коллекторским и емкостным свойствам. Размещение скважин по площади залежи – равномерное в пределах зоны разбуривания. Технологический режим эксплуатации «средней» скважины - постоянная депрессия на пласт. Коэффициенты фильтрационных сопротивлений не изменяются во времени. Исходные данные: состав и псевдокритические параметры пластового газа, коэффициенты фильтрационных сопротивлений, депрессия на пласт известны по результатам решения задач 1 и 2. Динамика изменения во времени пластового давления, «среднего» дебита скважины, положения границы раздела газ-вода по пропласткам известны по результатам решения задачи 8 (см. таблицу 11 и 12). Темп отбора в период постоянной добычи, продолжительность периода нарастающей добычи, доля накопленной добычи газа за период нарастающей и постоянной добычи Кпд приведены в приложении Б (таблица Б.5). Характеристики пропластков приведены в приложении Б (таблица Б.5). Вязкость воды в пластовых условиях μв= 0,74 мПа∙с, относительная фазовая проницательность kв*= 0,2, расстояние от начального положения газоводяного контакта до зоны разбуривания L0= 1500 м.
Рассчитать коэффициенты фильтрационных сопротивлений и коэффициент продуктивности для всех пропластков, динамику дебита по газу и воде для скважины, наиболее близко расположенной к начальному газоводяному контакту (на расстоянии равном L0), на периоды нарастающей и постоянной добычи.
Порядок расчета
Присваиваем величине i значение ноль (i – год разработки залежи). Рассчитываем коэффициенты фильтрационных сопротивлений au и bu по формулам (119) и (120) для всех пропластков. Результаты расчета заносим в таблицу 12.
На каждый момент времени для рассматриваемой скважины определяем номера газоотдающих пропластков и водопроявляющих по результатам полученным в 8 задаче и рассчитываем коэффициент вязкости газа по формулам (76)…(79).
На каждый момент времени рассчитываем для газоотдающих пропластков дебит по газу qu по формуле (123), для водопроявляющих – коэффициент продуктивности cв u и дебит по воде qв u по формулам (126) и (127) соответственно.
Рассчитываем на каждый момент времени суммарный дебит рассматриваемой скважины по газу q и по воде qв по формулам (124) и (125) соответственно.
Результаты расчета заносим в таблицы 13 и 14.
Таблица 13 - Результаты расчета коэффициентов фильтрационных сопротивлений
Номер пропласткаКоэффициенты фильтрационных сопротивлений
au, МПа2/(тыс. м3/сут) bu, МПа2/(тыс. м3/сут)2
1 nТаблица 14 - Результаты расчета дебита скважины по газу и воде
Год Pt, МПа μг, мПа·сqu, тыс. м3/сут или
cв u, м3/сут/МПа и qв u, м3/сутq,
тыс. м3/сут qв, м3/сут1-ый пропл. 2-ой пропл. 3-ий пропл. 4-ый пропл. 5-ый пропл. 1 NПрактическое занятие №10 по теме «Определение показателей разработки газоконденсатной залежи при полном сайклинг-процессе»Цель занятия: изучить методику расчета показателей разработки газоконденсатной залежи, разрабатываемой при полном сайклинг-процессе. Рассчитать динамику накопленного отбора конденсата, коэффициента охвата вытеснения жирного газа сухим, безразмерного объема закачиваемого газа, объема закачиваемого газа, доли жирного газа в потоке добываемой продукции, среднего дебита скважин, потребного количества скважин. Рассчитать показатели разработки газоконденсатной залежи, разрабатываемой при полном сайклинг-процессе.
Краткая теория вопроса
Пластовая газоконденсатная смесь в общем случае состоит из большого числа углеводородов (метана, этана, пропана, изобутана, н - бутана, пентана, гексана, гептана, октана, нонана, декана и более тяжелых), азота, сероводорода, углекислого газа, гелия, паров воды. Эта сложная система при изменении температуры и давления ведет себя иначе, чем индивидуальные чистые углеводороды (рисунок 4).
При повышении давления и неизменной температуре или понижении температуры и постоянном давлении происходят процессы конденсации пара в жидкость. Зависимость давления от температуры для чистого углеводорода, называемая кривой испарения (линия ОК на рисунке 11), является граничной кривой, ниже которой имеется одна паровая фаза, а выше, в области повышенных давлений, – одна жидкая фаза. Линия испарения есть граница скачкообразного изменения агрегатного состояния вещества.
Конечная точка этой линии К является критической. Она характеризует максимальную температуру Тк, при которой еще существует граница раздела фаз, т. е. паровая и жидкая фазы еще находятся в равновесии. При температуре выше критической нельзя превратить паровую фазу в жидкую при сколь угодно большом повышении давления. Критическим называется давление паров вещества (Ркр) при критической температуре. Объем вещества при критических температуре и давлении (Ткр и Ркр), отнесенный к одному молю или другой единице массы вещества, называется критическим молярным или удельным объемом.
Температура, при которой средняя молекулярная кинетическая энергия становится равной потенциальной энергии притяжения молекул, называется критической, так как при более высокой температуре невозможно существование жидкой фазы.
Математическим критерием критического состояния является равенство

где Ткр — точка перегиба изотермы на плоскости р—V при критических давлении и объеме.
80
60
40
20
С
Ср
в
Б
Па
А
Скк0
И
К
б
Д
В
О
Ж+ПРккРкре
д
Ж
Скр100
Е
ТккТкрТ
Р
Рисунок 11 – Диаграмма фазовых превращений газоконденсатной системы постоянной массы и состава при изменении давления и температуры (критическая точка лежит правее крикондебары)
Фазовая диаграмма газоконденсатной смеси дана на рисунке 11. Кривая ССкр — линия кипения. Выше этой линии — одна жидкая фаза. Линия СкрБСккДИ — линия конденсации. Правее и ниже нее — одна паровая фаза. Внутри линии ССкрСккИ — двухфазная область, область одновременного сосуществования паровой и жидкой фаз. Цифры на линиях означают процентное объемное содержание жидкой фазы в смеси.
Точка Скр — критическая точка, точка Ткк — крикондентерма, максимальная температура (выше критической Ткр), при которой жидкая и паровая фазы еще могут находиться в равновесии, т. е. еще имеется граница раздела фаз пар — жидкость.Рассмотрим изотермический процесс уменьшения давления от точки А. В точке А углеводородная смесь находится в области паровой фазы. При снижении давления и увеличении объема сосуда высокого давления при неизменном составе смеси до точки Б не будет фазовых изменений. В точке Б при уменьшении давления образуется первая капля жидкости. В этом случае происходит процесс обратной конденсации, т. е. образование жидкой фазы при уменьшающемся давлении.
При дальнейшем снижении давления объем образовавшейся жидкой фазы увеличивается и в точке В достигает максимального значения. Область СкрВСккБСкр называется областью обратной конденсации, линия СкрВСкк — линией давлений максимальной конденсации. При дальнейшем снижении давления от точки В до точки Д ранее образовавшаяся жидкая фаза будет уменьшаться в объеме, испаряться и в точке Д испарится последняя капля жидкости.
При снижении давления от точки В до точки Д идет процесс испарения жидкости при уменьшении давления. Дальнейшее падение давления от точки Д до точки Е не влечет за собой фазовых превращений, смесь находится в паровой фазе. Процесс обратной конденсации наблюдается только в интервале температур Ткр - Ткк.
Рассмотрим процесс изобарического (при постоянном давлении) снижения температуры от точки а. В точке а газоконденсатная смесь находится в жидкой фазе. При ее охлаждении до точки б фазовых переходов нет. В точке б образуется первый пузырек пара. Образование паровой фазы при понижении температуры при постоянном давлении называется процессом обратного испарения. При снижении температуры от точки б до точки в объем паровой фазы увеличивается и в точке в достигает максимума. Область СкрвСрбСкр называется областью обратного испарения, а линия СрвСкр — линией температур максимального испарения.
При понижении температуры от точки в до д объем образовавшейся паровой фазы уменьшается, паровая фаза конденсируется и в точке д сконденсируется последний пузырек пара. При уменьшении температуры от точки в до точки д идет процесс нормальной конденсации. Дальнейшее снижение температуры от точки д до точки е не вызывает фазовых переходов, углеводородная смесь находится в жидкой фазе. Явление обратного испарения наблюдается только в интервале изменения давления Ркр – Ркк, Ркк – криконденбар – максимальное давление на линии точек кипения, при котором может существовать жидкость.
Образование жидкой фазы в пористой среде приводит к потерям жидкости. При разработке газоконденсатной залежи без поддержания пластового давления в условиях газового режима потери жидкого конденсата в пласте могут составлять 30 – 60 % от начального (потенциального) содержания конденсата (C5+) в пластовом газе. Эти потери могут составлять миллионы тонн.
Процесс конденсации в пористой среде с ее громадной удельной поверхностью протекает иначе, чем в сосуде PVT при большой плоской границе раздела пар — жидкость. В поровых каналах небольшого радиуса — в капиллярах – будет происходить процесс капиллярной конденсации, где граница раздела пар — жидкость будет криволинейной. В связи с проявлением капиллярных сил в пористой среде давление начала образования жидкой фазы, объем образовавшейся жидкости, объем оставшейся жидкой фазы в пористой среде при одинаковом давлении будут больше, чем в сосуде PVT.
В зависимости от характера флюидов, находящихся в продуктивном пласте, (р-V-Т) - состояния этих флюидов месторождения природных углеводородов классифицируют на нефтяные, нефтегазоконденсатные, газоконденсатные, газовые и газогидратные. На рисунке 12 представлена типичная диаграмма фазового состояния флюида для многокомпонентных углеводородных смесей. Нефтяные месторождения (I) могут существовать при температурах ниже критической (левее точки Скр). Различают нефтяные месторождения с давлением выше давления насыщения (зона I, выше точек кипения); нефтяные месторождения с давлением, равным давлению насыщения (зона I, кривая точек кипения); нефтяные месторождения с давлением ниже давления насыщения (зона I, ниже точек кипения) или так называемые двухфазные нефтяные месторождения (нефть и газовая шапка). Для всех нефтяных месторождений характерен процесс выделения газа при снижении давления и постоянной температуре.

Точка Скр — критическая точка; Ср — криконденбар — максимальное давление на линии точек кипения, при котором может существовать жидкость и пар могут существовать в равновесии, Скк — крикондентерм — максимальная температура на линии точек росы, при которой жидкость и пар могут существовать в равновесии; F — точка гидратообразования, месторождения I — нефтяные, II — нефтегазоконденсатные, III — газоконденсатные, IV — газовые, V — газогидратныеРисунок 12 – Диаграмма фазового состояния многокомпонентной углеводородной смеси (критическая точка лежит левее крикондебары).
В ретроградной области СкрСрСккСкр имеет место обратное явление. При снижении давления при постоянной температуре в ретроградной области наблюдается выделение жидкости (конденсация углеводородов). В этой области между критической точкой (Скр) и криконденбаром (Ср) лежит зона II – зона существования нефтегазоконденсатных месторождений. В этой же ретроградной области между криконденбаром (Ср) и крикондентермом (Скк) лежит зона III – зона существования газоконденсатных месторождений.
Нефтегазоконденсатные и газоконденсатные месторождения могут быть как недонасыщенными (пластовое давление выше давления точек кипения), насыщенными (пластовое давление равно давлению точек кипения), так и двухфазными (пластовое давление ниже давления точек кипения).
Зона IV правее точки Скк (крикондентерм) и ниже линии точек росы представляет зону существования газовых месторождений. В этой зоне снижение давления при постоянной температуре не приводит к возникновению углеводородной жидкой фазы. Наконец, левее точки F (точка гидратообразования) находится зона V – зона существования газогидратных залежей.
Начальное пластовое давление, как правило, соответствует гидростатическому давлению столба воды на глубине залегания продуктивного коллектора. Известны месторождения с аномально высокими пластовыми давлениями (давление в пласте выше гидростатического) и с аномально низкими пластовыми давлениями.
Сайклинг - процесс – это метод разработки газоконденсатных залежей, заключающийся в поддержании пластового давления путем обратной закачки отбензиненного (осушенного) газа в пласт. Различают полный и частичный сайклинг-процесс. При полном сайклинг-процессе поддерживается давление начала конденсации. При этом весь отбензиненный газ закачивается в пласт. Могут привлекаться дополнительные ресурсы сухого газа, поскольку отбензиненный газ по сравнению с жирным пластовым газом имеет усадку.
Приближенная методика расчета показателей разработки газоконденсатной залежи при полном сайклинг-процессе
Накопленная с начала разработки добыча конденсата рассчитывается по формуле
,(128)
где – коэффициент охвата вытеснения жирного газа сухим;
– балансовые запасы конденсата, тыс. т, рассчитываемые по формуле:
,(129)
где – начальное содержание пентанов и высших в пластовом газе, г/м3;
– балансовые запасы пластового газа, млн. м3.
Коэффициент охвата вытеснением жирного газа сухим – это отношение объема пласта, охваченного воздействием ко всему объему пласта. Он показывает эффективность процесса вытеснения. На эффективность вытеснения влияет степень неоднородности пласта по коллекторским свойствам. Коэффициент охвата вытеснением зависит от безразмерного объема закаченного газа и рассчитывается по формулам (130):
(130)
где – приведенный к пластовым условиям объем закачиваемого газа, приходящийся на 1 м3 газонасыщенного порового объема пласта;
– безразмерный объём закачиваемого газа, при котором происходит прорыв сухого газа к добывающим скважинам; величина этого объёма равна величине коэффициента охвата при прорыве.
Прорывное значение коэффициента охвата вытеснением зависит от степени неоднородности пласта по коллекторским свойствам и в некоторых случаях может быть рассчитано по эмпирической зависимости:
,(131)
где v – вариация проницаемости.
Вариация проницаемости зависит от среднеквадратического отклонения и средневзвешенного значения проницаемости и определяется по формуле
,(132)
где – среднеквадратическое отклонение проницаемости, рассчитываемое по формуле:
,(133)
где – проницаемость i-го пропластка, мкм2;
n – количество пропластков;
– средневзвешенное значение проницаемости, мкм2;
,(134)
где hi – толщина пропластка, м.
Приведенный к пластовым условиям объем закачиваемого газа, приходящийся на 1 м3 газонасыщенного порового объема пласта (безразмерный объем закачиваемого газа) на момент времени t, рассчитывается по формуле
,(135)
где – объем закачиваемого газа, приведенного к пластовым условиям, млн м3
– газонасыщенный объем порового пространства, млн. м3.
Объем закачиваемого газа на момент времени t, приведенного к пластовым условиям, рассчитывается по формуле
(136)
,(137)
где – объем закачиваемого газа на момент времени t, приведенного к стандартным условиям, млн. м3;
– темп отбора на момент времени t, %;
Pст, Tст, zст – давление (МПа), температура (К) и коэффициент сверхсжимаемости в стандартных условиях;
Pпл, Tпл, zпл – давление (МПа), температура (К) и коэффициент сверхсжимаемости в пластовых условиях.
Доля жирного газа в потоке добываемой продукции до прорыва сухого газа равна единице
при .(138)
После прорыва закачиваемого сухого газа доля жирного газа в потоке добываемой продукции будет меньше единицы и рассчитывается по формуле
при ,(139)
где , – приращение величин за предыдущий год, определяемых по формулам:
, (140)
,(141)
где – коэффициент сверхсжимаемости закачиваемого сухого газа в пластовых условиях.
Забойное давление в нагнетательных скважинах рассчитывают по формуле
.(142)
Так как газоконденсатная залежь разрабатывается с поддержанием пластового давления, то для расчета по формулам (142) и (145) соответственно используется начальное пластовое давление.
Дебиты нагнетательных скважин рассчитывают по формуле
.(143)
Количество нагнетательных скважин рассчитывают по формуле
.(144)
Забойное давление в добывающих скважинах рассчитывают по формуле
.(145)
Дебиты добывающих скважин рассчитывают по формуле
.(146)
Количество добывающих скважин рассчитывают по формуле
.(147)
Задача 5. Разрабатывается газоконденсатная залежь с поддержанием пластового давления путем обратной закачки отбензиненного газа в пласт. Разработка газоконденсатной залежи ведется при полном сайклинг-процессе. Модель залежи – слоисто-неоднородный пласт, состоящий из пропластков с различными коэффициентами проницаемости, а также толщиной пропластков. Технологический режим эксплуатации нагнетательных и добывающих скважин – постоянная депрессия на пласт. Депрессия на пласт нагнетательных скважин в 2 раза больше, чем для добывающих скважин. Коэффициенты фильтрационных сопротивлений нагнетательных и добывающих скважин равны и не изменяются во времени. Исходные данные: балансовые запасы пластового газа, начальный газонасыщенный поровый объем, пластовая температура, состав и псевдокритические параметры пластового газа, параметры «средней» добывающей скважины известны из исходных данных и по результатам решения задач 1 и 2. Темп отбора в период постоянной добычи и продолжительность периода нарастающей добычи приведены в приложении Б, таблице Б.5. Характеристики пропластков приведены в приложении Б, таблице Б.8. Начальное содержание пентанов и выше кипящих в пластовом газе составляет 300 г/м3. Разработка считается рациональной до тех пор, пока доля жирного газа больше 50%.
Рассчитать динамику пластового и забойного давлений, среднего дебита добывающих и нагнетательных скважин, потребного количества добывающих и нагнетательных скважин, доли жирного газа, добытого за год конденсата, накопленного отбора газа и конденсата. Расчеты произвести по временным шагам, i – номер временного шага (i – й момент времени).
Порядок расчета
Рассчитываем запасы конденсата по формуле (129).
Рассчитываем прорывное значение коэффициента охвата вытеснением по формулам (131) – (134).
По известному годовому отбору газа определяем накопленную добычу газа на момент времени t.
По формуле (137) определяем объем закачиваемого газа на момент времени t, приведенного к стандартным условиям.
По формуле (136) определяем объем закачиваемого газа на момент времени t, приведенного к пластовым условиям.
Рассчитываем по формуле (135) безразмерный объем закачиваемого газа на момент времени t.
Рассчитывается коэффициент охвата вытеснением на момент времени t по формуле (130).
Рассчитывается накопленная с начала разработки добыча конденсата на момент времени t по формуле (128).
Рассчитываем забойное давление, дебит и потребное количество добывающих и нагнетательных скважин на момент времени t по формулам (142) – (147).
Определяем долю жирного газа в потоке добываемой продукции на момент времени t по формулам (138) – (141).
Если полученная доля жирного газа выше предельной заданной величины доли жирного газа, то переходят к расчету показателей разработки на следующий момент времени, т.е. производят расчеты по пунктам 3 – 10, иначе прекращают расчеты.
По результатам расчета показателей разработки заполняется таблица 15 и на рисунках приводится динамика показателей разработки.
Таблица 15 – Результаты расчета динамики показателей разработки газоконденсатной залежи при полном сайклинг-процессе
год , % , млн м3/год , млн м3 Qзак, млн м3 , млн м3 Q* fохвqж, МПа , МПа qн, тыс м3/сутqд, тыс м3/сутnн штnд, шт, тыс т
СПИСОК ЛИТЕРАТУРЫ
АСУ ТП газопромысловых объектов/ А.Г. Ананенков, Г.П. Ставкин, О.П. Андреевич и др. – М.: ООО «Недра-Бизнесцентр», 2003.-343с.
Инструкция по применению классификации запасов к месторождениям нефти и горючих газов.
Козлов А. Л., Коротаев Ю. П., Фиш М. Л., Фриман Ю. М., Букреева Н. А., Тверковкин С.М.  Подсчет запасов газа по падению давления. Тематический научно – технический обзор. – М.: ВНИИЭОПТИГП, 1969.
Компьютерные технологии вычислений в математическом моделировании: Учеб. пособие. – М.: Финансы и статистика, 1999.–256 с.
Коротаев Ю.П., Ширковский А.И. Добыча, транспорт и подземное хранение газа. – М.: Недра,1984.–486 с.
Нефтегазопромысловая геология: Терминологический справочник/ Под ред. М.М. Ивановой. – М.: Недра, 1983. -262 с.
Подсчет запасов нефти, газа, конденсата и содержащихся в них компонентов: Справочник/И. Д. Амелин, В. А. Бадъянов, Б. Ю. Вендельштейн и др.; Под ред. В. В. Стасенкова, И. С. Гутмана. – М.: Недра, 1989.–270 с.:ил.
Разработка газовых, газоконденсатных и нефтегазоконденсатных месторождений/ С. Н. Закиров - М.: Струна, 1998.- 628с.
Физические величины: Справочник/ А. П. Бабичев, Н. А. Бабушкина, А. М. Братковский и др.; Под. ред. И. С. Григорьева, Е. З. Мейлихова. – М., Энергоатомиздат, 1991. – 1232 с.
Приложение АМетодика расчета пластового давления из отношения Р/zДля решения поставленной задачи можно использовать графический метод путем построения графиков Р от Р/z. Однако при вычислениях на ПЭВМ рациональнее использовать аналитический метод, который заключается в решении системы уравнений, представленной ниже:
(69)
Эту систему можно свести в общем случае к трансцендентному уравнению, которое решается любым численным методом.
,(70)
гдеС – значение , определяемое по уравнению материального баланса (2);
– известная функция зависимости коэффициента сверхсжимаемости от давления.
Полученное уравнение можно решить методом деления отрезка пополам.
Порядок решения уравнения методом деления отрезка пополам
Процесс решения нелинейного уравнения осуществляется в два этапа. На первом этапе находят такие отрезки, внутри которых находится строго один корень. Для этого используют графический (построение графика) или аналитический методы (численное или аналитическое дифференцирование). На втором этапе производят поиск корня тем или иным способом.
В нашем случае вначале необходимо определить отрезок, содержащий решение уравнения. Для этого будем определять знак функции при значениях . Величина Δ первоначально принимается небольшой (0,001 МПа). Если произведение > 0, тогда величину Δ увеличивают, т.е. расширяется отрезок, на котором происходит поиск корня уравнения. Если произведение < 0, то на отрезке содержится корень уравнения.
Сущность метода решения уравнения методом деления отрезка пополам состоит в делении отрезка, содержащего корень, на две равные части и после этого определяют, в какой из этих двух частей содержится корень. Процесс поиска корня итерационный. На первой итерации находится среднее:
.(71)
Затем определяем знаки произведений функций и . В дальнейшем ищем корень на отрезке с отрицательным произведением функций. При этом в зависимости от отрезка, содержащего корень, происходит присвоение или . Производят соответственно проверку условия:
или .(72)
Если условие (36) не выполняется, то переходят к следующей итерации, т.е. повторяют расчет для нового отрезка, содержащего корень, определяют знаки произведений функций и и т.д. Если условие выполняется, то считают, что корень найден и равен [4].
Приложение БВарианты исходных данных к задачам 1-10
Таблица Б.1 – Динамика накопленной добычи и средневзвешенного пластового давления
№ замера , млн. м3 , МПа
1 548,3 22,42
2 1118,5 21,79
3 1638,9 21,65
4 2318,1 21,01
5 2751,4 20,91
6 3310,6 20,09
7 3871,8 20,03
8 4407,7 19,27
9 5305,1 18,47
10 5772,4 18,29
Таблица Б.2 – Множители к расчету накопленной добычи и средневзвешенного пластового давления по вариантам, глубина залежи
Вариант х Пластовая температура, КВариант х Глубина залежи, мА 1,20 Т=325 1 1,20 1650
Б 1,10 2 1,10 1530
В 1,05 3 1,05 1460
Г 0,95 Т=305 4 0,95 1405
Д 0,85 5 0,9 1300
Е 0,80 6 0,85 1275
Таблица Б.3 – Состав пластового газа
Компонент Мольная доля, ηi Ркр, МПа Ткр, К
Вариант I Вариант II Вариант III CH4 0,97 0,95 0,92 4,604 190,55
C2 H6 0,01 0,015 0,015 4,880 305,43
C3 H8 0,005 0,01 0,01 4,249 369,82
n C4 H10 0,004 0,005 0,015 3,796 425,16
i C4 H10 0,001 0,004 0,007 3,647 408,13
n C5 H12 0,01 0,016 0,033 3,369 469,65
Таблица Б.4 – Коэффициенты фильтрационных сопротивлений и депрессия на пласт скважин
№ скв. δi, МПа
1 0,0135 0,41 0,8
2 0,0119 0,5 0,75
3 0,0065 0,9 0,5
4 0,011 0,4 0,6
5 0,0387 0,28 0,7
6 0,0329 0,3 0,9
7 0,0145 0,63 0,45
8 0,002 1,15 0,72
9 0,0245 0,8 0,83
10 0,0018 0,4 0,55
Таблица Б.5 – Темп отбора в период постоянной добычи и продолжительность периода нарастающей добычи
Вариант Темп отбора в период постоянной добычи, % от Кпд Вариант Продолжительность периода нарастающей добычи, год
А 4 0,6 1 7
Б 5 0,55 2 6
В 6 0,65 3 5
Г 5 0,57 4 7
Д 6 0,5 5 6
Е 7 0,52 6 4
Таблица Б.6 – Радиусы батарей и количество скважин в каждой батареи
Батарея Радиус батареи RiЧисло скважин в батарее mi, шт1 0,70∙Rз17
2 0,55∙Rз15
3 0,40∙Rз12
4 0,25∙Rз8
5 0,10∙Rз6
Таблица Б.7 – Проницаемость и толщина пропластков№ пропласткаПроницаемость ku, мкм2Толщина hu, м1 0,055 8,1
2 0,098 1,2
3 0,012 2,9
4 0,005 1,7
5 0,031 5,1
Таблица Б.8 – Проницаемость и толщина пропластков№ пропласткаПроницаемость ki, мкм2Толщина hi, м
1 0,088 2,1
2 0,167 0,4
3 0,255 0,7
4 0,533 1,5
5 0,290 2,5
6 0,180 0,9
7 0,105 0,2
8 0,308 1,7
9 0,200 0,5
10 0,067 1,0
Таблица Б.9 – Критические параметры некоторых веществ
Вещество Химическая формула υкр, см3/(г∙моль) Ркр, МПа Ткр, К zкрρкр, кг/м3
метан СН499,5 4,604 190,55 0,290 162,0
этан С2Н6 148,0 4,880 305,43 0,285 203,0
пропан С3Н8 200,0 4,249 369,82 0,277 220,0
н-бутан n-C4H10 255,0 3,796 425,16 0,274 228,0
изо-бутанi-C4H10 263,0 3,647 408,13 0,283 221,0
н-пентан n-C5H12 311,0 3,369 469,65 0,269 232,0
изо-пентан i-C5H12 310,0 3,381 460,39 0,270 236,0
гексанC6H14 368,0 3,013 507,35 0,264 234,0
гептан C7H16 426,0 2,736 540,15 0,352 235,0
октан C8H18 486,0 2,486 568,76 0,256 235,0
азот N2 90,1 3,398 126,26 0,291 311,0
водород H2 65,0 1,296 33,25 0,304 30,7
воздух 86,5 3,648 132,46 — 335,0
водяной пар Н2О 57,0 21,408 647,30 0,234 316,0
кислород O2 74,4 5,080 154,78 0,292 430,0
сероводород H2S 95,0 9,007 373,60 0,268 359,0
двуокись углерода CO2 94,0 7,381 304,20 0,274 468,0
окись углерода CO 93,1 3,499 132,93 0,294 301,0
двуокись азота NO2 82,0 10,132 431,00 0,232 561,0
окись азота NO 58,0 6,535 180,30 0,260 520,0
двуокись серы SO2 122,0 7.893 430,65 0,268 525,0
гелий He57,8 0,229 5,20 0,300 69,2
аргон Ar75,2 4,863 150,72 0,290 531,0
криптон Kr92,2 5,492 209,39 0,291 908,0
фтор F2 — 5,573 144,20 0,292 630,0
хлор Cl2 124,0 7,711 417,20 0,276 573,0
этилмеркаптанC2H5SH 207,0 5,315 499.10 0,274 __
вода H20 56,0 22,119 547,40 0,228 325,0

Приложенные файлы

  • docx 25333072
    Размер файла: 1 MB Загрузок: 1

Добавить комментарий