Моделирование процессов статического конусообразования при разработке нефтяных, газовых и нефтегазовых залежей
Изучение существующих работ, связанных с разработкой указанных залежей, показывает, что эта проблема находится в стадии интенсивного теоретического и промыслового исследования.
2.2 Моделирование процесса статического конусообразования
2.2.1 Общие представления
Во многих случаях при разработке нефтегазовых залежей (НГЗ) вскрываются газоводонефтяные зоны или нефтяные оторочки при разработке газоконденсатнонефтяных залежей (ГКНЗ) с подошвенной водой. Разработка таких залежей обусловливается следующими характерными особенностями, полной гидродинамической связью нефтяной залежи с газовой шапкой и водоносным пластом и вероятной подвижностью газонефтяного и водонефтяного контактов в окрестности скважин в процессе разработки залежи; практически неподвижностью контуров газоносности и нефтеносности в пласте; равномерным распределением пластовой энергии по площади нефтеносности; равенством начального пластового давления и давления насыщения; относительной близостью расположения к забоям скважин водонефтяного и газонефтяного контактов при дренировании нефтяной оторочки; неустойчивостью процесса вытеснения нефти газом, приводящей к быстрому прорыву газа к забоям добывающих скважин и их загазованности и в конечном счете к значительной потере пластовой энергии и снижению нефтеотдачи; возможностью проявления ретроградной конденсации из-за снижения давления в газонасыщенной зоне пласта, предопределяющей пластовые потери конденсата; трудностью регулирования перемещением ГНК и ВНК и др.
Как видим, указанные особенности требуют создания технологии разработки НГЗ и нефтяных оторочек, совершенно отличных от технологии разработки как нефтяных залежей обычного типа, так и нефтегазовых залежей с краевой водой.
При разработке НГЗ и ГКНЗ с подошвенной водой темп отбора нефти обусловливается деформацией контактов и прорывом газа и воды к забоям скважин. При этом весьма важным параметром при установлении режима работы скважин и прогнозировании технологических показателей разработки является анизотропия пласта [ 1 ], обоснование которой необходимо для каждой конкретной залежи. М. Маскет также указывает [ 1 ], что анизотропность коллектора существенно влияет на эффективность размещения скважин. Низкая проницаемость по вертикали препятствует быстрому поднятию вершины конуса и способствует выполаживанию поверхности раз дела вода-нефть. Высокая проницаемость по вертикали (малая анизотропия пласта) способствует быстрому продвижению вершины конуса к забою скважины, что обусловливает концентрированную деформацию поверхности раздела вблизи скважины с низким коэффициентом охвата вытеснения нефти подошвенной водой. Поэтому М. Маскет утверждает, что критерием размещения скважин с напором подошвенной воды должен быть параметр размещения, представляющий собой отношение половины расстояния между скважинами Ro к произведению толщины продуктивного пласта h0 и анизотропии æ*, т.е.
p0=Ro/æ*h0.
Эффективность вытеснения нефти, очевидно, может быть улучшена в условиях непрерывной эксплуатации скважин при малых дебетах, когда снижение забойного давления не на много превосходит напор столба жидкости P=hH(ρB-ρH)g, а также при периодической эксплуатации скважин с продолжительными интервалами ее закрытия, приводящей к опусканию и выполаживанию возникшей конусообразной поверхности раздела вода-нефть или газ-нефть под действием сил тяжести.
2.2.2 Математическая модель статического конусообразования Маскета-Чарного. Методы расчета предельных депрессий и дебитов несовершенных скважин
Модель предполагает установившийся приток нефти или газа к открытому забою скважины, частично вскрывшей однородный или однородно-анизотропный по проницаемости ограниченный горизонтальный пласт постоянной толщины, подстилающийся подошвенной водой. На контуре пласта и на забое скважины поддерживается постоянное давление, фильтрация происходит по закону Дарси, капиллярными силами пренебрегается, вытеснение нефти или газа водой предполагается поршневое. Решение для распределения потенциала в пласте, вызванного работой несовершенной скважины, принимается для условий невозмущенной границы раздела двух жидкостей, т.е. первоначальный ВНК и ГНК предполагаются непроницаемыми.
При эксплуатации нефтяных или газовых скважин с подошвенной водой, а также при дренировании нефтяной оторочки в определенных условиях проявляется тенденция к деформированию поверхности раздела двух фаз, которая принимает холмообразный вид, образуя конусы воды, конусы нефти или конусы воды и газа. При некоторых установившихся условиях отбора соответствующие деформированные поверхности раздела находятся в равновесии и не оказывают существенного влияния на приток добываемого флюида к скважине. Если же превысить депрессию и, соответственно, отбор нефти или газа сверх некоторой предельной величины, то вода прорвется в скважину, что может привести к ее прогрессирующему обводнению, а при дренировании нефтяной оторочки возможен прорыв подошвенной воды и верхнего газа. Таким образом, существует предельная высота вершины конуса, которой соответствуют предельная депрессия и безводный или безгазовый дебит.
Точной теории конусообразования ввиду сложности процессов не имеется. Приближенная теория этого явления, основанная М.Маскетом и И.А.Чарным [1,2] и позволяющая рассчитывать предельные дебиты и депрессии, исходит из допущения, что отклонение поверхности раздела двух фаз от первоначальной плоской формы не влияет на распределение потенциала скорости фильтрации в нефтяной (газовой) части пласта.
Дальнейшее развитие приближенной теории устойчивых конусов Маскета-Чарного и ее практическое использование нашли отражение в работах как отечественных, так и зарубежных исследователей (Б.Б.Лапук, Б.Е Сомов, А.Л.Брудно, Д.А.Эфрос, Р.Г.Аллахвердиева, А.К.Курбанов, П.Б.Садчиков, А.П.Телков, Ю.И.Стклянин, З.С.Алиев, А.П.Власенко, Е.С.Абрамов, С.Н.Закиров, Р.Чаней, Д.Сирси и др.) Здесь изложен более универсальный метод решения задач конусообразования, основанный на двухзонной схеме притока.
Предельный безводный дебит нефтяной скважины определяется по формуле
Q1 = Q0(ρ0,Ђ); Q0 = 2πKrh02 Δρ; Δρ=ρB-ρr, (2.1)
где q(ρ0,Ђ) - безразмерный безводный дебит, определяемый по соответствующим графикам или таблицам (Прил.1) для параметра ρ0≤1 [79,82] и по графикам рис.2.3 для ρ0>1. Безразмерная ордината конуса ξ0=z0/h0 для ρ0>1 определяется по графикам рис.2.4, для ρ0<1 - по таблицам (см. Прил.1).
Для безразмерной депрессии при р0> 1 имеется формула:
ΔPпр = = (ε0+Δε)q(ρ0,Ђ), (2.2)
где
ε0=[ln - - ψ(ρ0,Ђ)] ; (2.3)
ψ(ρ0,Ђ) - некоторая функция, определяемая по таблице (Прил.2). Для Де построены графики (рис.2.5). Возможно другое, наиболее полное представление для функции фильтрационных сопротивлений