
- •8.1. Краткий обзор существующих работ
- •8.2. Построение обобщенного дифференциального уравнения неустановившейся фильтрации однородной жидкости и газа в пористой среде при изотермическом процессе
- •(Источников) в пространстве
- •8.3. Приток к несовершенной линии стоков (скважине) в ограниченном пласте при наличии подошвенной воды
- •Прямоугольной формы за счет напора подошвенной воды
- •9. Методы расчета фильтрационных сопротивлений. Табулирование сложных функций
- •9.1. Краткий обзор существующих работ; постановка задач
- •9.2. Методы расчета фильтрационных сопротивлений при установившемся притоке жидкости и реального газа к несовершенной скважине. Табулирование функций
- •Ограниченном однородно-анизотропном пласте
- •Т абулированные значения функции
- •Экраном и относительным вскрытия пласта
- •Обусловленного нелинейным законом фильтрации
- •С1 от относительного вскрытия пласта при параметрах ρ0 и
- •9.3. Методика расчета фильтрационных сопротивлений при неустановившемся осесимметричном притоке жидкости (газа) к несовершенной скважине в неограниченном пласте.
- •При параметре
- •9.4. Методика расчета фильтрационных сопротивлений при неустановившемся притоке жидкости к несовершенной скважине в ограниченном пласте по линейному закону
- •9.5. Методика расчета фильтрационных сопротивлений, обусловленных перфорацией колонны
- •Пласта æ* при фиксированной глубине l0 пулевого канала (см)
- •Канала при фиксированном значении анизотропии пласта æ*
- •10. Интерпретация результатов исследования гидродинамически несовершенных скважин при нестационарной фильтрации
- •10.1. Общая характеристика прискважинной зоны пласта
- •10.2. Основы дифференциального и интегрального методов обработки кривых восстановления давления в пласте
- •10.3. Влияние учета несовершенства скважин на точность определения параметров пласта при интерпретации кривых восстановления давления
- •10.4. Влияние изменения проницаемости на характеристики пласта
- •Исходные данные для обработки квд
- •10.5. Определение радиуса кольцевой неоднородности по квд при дренировании однородно-анизотропного пласта несовершенной скважиной
- •Неоднородностью
- •10.6. Интерпретация кольцевой неоднородности пласта и скин-эффект в условиях плоско-радиального потока
- •Литература к гл. 8-10
- •11. Моделирование процессов статического конусообразования при разработке нефтяных, газовых и нефтегазовых залежей
- •11.1. Сущность проблемы конусообразования
- •11.2. Моделирование процесса статического конусообразования
- •Статическом равновесии границы раздела
- •11.3. Методы расчета предельных безводных и безгазовых дебитов несовершенных скважин, дренирующих нефтегазовые залежи с подошвенной водой
- •При безнапорном притоке к несовершенной скважине
- •Воды в условиях напорного притока к несовершенной скважине
- •Зависимости от расположения интервала вскрытия пласта
- •11.4. Расчет предельных безводных дебитов несовершенных сважин и депрессий в газовых залежах с подошвенной водой при линейном законе фильтрации
- •Результаты расчетов погрешности d0 по формуле (11.49)
- •11.5. Решение задач конусообразования по двухзонной схеме притока
- •Определение ординаты x0 и функции е0(x0, r, )
- •Литература к гл. 11
- •12. Моделирование процессов динамического конусообразования при разработкЕ водонефтяных и газонефтяных залежЕй
- •12.1. Краткий обзор теоретических работ по конусообразованию
- •12.2. Упрощенные и строгие методы расчета времени безводной эксплуатации скважин с подошвенной водой
- •Скважины t от относительного вскрытия пласта
- •12.3. Методика прогнозирования продвижения границы раздела и нефтеотдачи за безводный период по удельному объему дренирования
- •12.4. Уточненная методика расчета безводного периода эксплуатации несовершенной скважины при опережающей разработке нефтяной оторочки
- •12.5. Уточненная методика расчета времени прорыва нефти из оторочки к забою газовой скважины при опережающей разработке газовой шапки
- •12.6. Уточненная методика расчета времени прорыва газа из газовой шапки к забою несовершенной скважнны, дренирующей нефтяную оторочку
- •Залежи несовершенной скважиной
- •Литература к гл. 12
- •13. Установившийся и неустановившийся приток жидкости и газа к вертикальным трещинам грп и горизонтальным стволам
- •13.1. Установившийся приток к вертикальным трещинам и горизонтальным стволам скважин
- •Скважине и несовершенной щели в полосообразном пласте
- •13.2. Наиболее известные формулы дебита горизонтальных стволов нефтяных скважин при установившемся притоке
- •13.3. Определение дебита горизонтального ствола скважины по методу эквивалентных фильтрационных сопротивлений
- •Горизонтальной скважины по сравнению с дебитом вертикальной
- •13.4. Определение оптимального местоположения и дебита горизонтального ствола скважины, дренирующего нефтегазовую залежь с подошвенной водой
- •Залежи с подошвенной водой
- •Погрешность формул (13.4.1) и (13.4.2)
- •Определение безразмерного дебита 10 скважины-трещиы
- •13.5. К обоснованию оптимальной сетки горизонтальных скважин и сравнительная эффективность их работы вертикальными трещинами и скважинами
- •Расположением горизонтальной скважины
- •Результаты расчета оптимальных размеров а и b сетки размещения горизонтальных скважин и вертикальных трещин и их эффективности при исходных параметрах a, l
- •13.6. Неустановившийся приток жидкости и газа к несовершенной галерее (вертикальной трещине грп) и горизонтальному стволу скважины по двухзонной схеме
- •4.Приток к горизонтальному стволу
- •Трещины q0 от степени вскрытия пласта
- •5. Приток реального газа к вертикальной трещине грп и горизонтальному стволу по нелинейному закону фильтрации
- •13.7. Установившийся и неустановившийся приток жидкости к многозабойным горизонтальным скважинам
- •13.7.1. Некоторые типовые профили многозабойных скважин
- •Разработке нефтегазовых залежей
- •Воды горизонтальными стволами в плоскости (X, z)
- •(Y, z) при одновременно–раздельном отборе воды и нефти
- •Линиями нагнетания
- •13.8. Решение некоторых гидродинамических задач притока жидкости к горизонтальным стволам скважин на основе теории функций комплексного переменного.
- •Продуктивном блоке
- •Результаты расчета фукнкции f(ρ,
- •Литература к гл. 13
- •1.Чарный и.А. Подземная гидромеханика. Гтти, 1948.
- •Результаты расчета добавочных фильтрационных сопротивлений при
- •Табулированные значения функции фильтрационного сопротивления по формуле (9.3.4)
- •Значение безразмерных плотностей по формулам (11.25) и (11.26)
4.Приток к горизонтальному стволу
Используя известную аппроксимацию, предложенную Ю.П.Борисовым [6], Р.Супрунович и Р.Батлер [31] нашли выражение для падения двления, описываемое формулой (13.5.7). Для нашей задачи с учетом анизотропии æ* она принимает следующий вид
(19.6.31')
Авторы [31] решали плоскую задачу фильтрации жидкости к трещине (горизонтальному стволу). Поэтому, строго говоря, использование формулы (13.6.31') для трехмерного пространства будет некорректным.
Исходя из сопоставления формул для фильтрационных сопротивлений Jтр и Jгс, соответствующих притоку к вертикальной трещине и горизонтальному стволу, дренирующих полосообразный полубесконечный пласт, предлагается формула падения давления за счет конвергенции линий тока [6,31]:
(13.6.32)
где
(13.6.33)
Здесь
– безразмерная
вертикальная ордината положения
горизонтального ствола.
Таким образом, при установившейся фильтрации формула (13.6.9) для горизонтального ствола представится в виде
(13.6.34)
Рассмотрим конкретный пример для Случая 1, когда f0>>1. принимаем исходные данные (см.рис.13.2): l1=50 м; l=2l1; h0=10 м; η*=0,5; æ*=5; K=10 мДа=1,02∙10-14м2; μ=1 мПа·с; ΔР=2∙106 Па; L=100 м; =0,5.
– определяем безразмерные параметры: l*=0,2; ρ=1,2.
– По соответствующим формулам находим: F(ρ, ,l*)=0,1033; Jтр=25,18; q0=0,44 м2/сут, следовательно дебит трещины составит Q=q0L=0,44·100=44 м3/сут; lnZ=1,0415; Jвк=38,543; q0г=0,174 м2/сут (в 2,5 раза меньше вертикальной трещины);
– Удельный дебит вертикальной трещины, рассчитанные по формуле (13.6.11) для полного вскрытия пласта, т.е. при =1, составляет q0=0,881 м2/сут или Q0=88,1 м3/сут, что превышает ровно в 2 раза дебит трещины при вскрытии =0,5.
Полученные результаты дают возможность построить графическую зависимость q0=q0( ), рис.13.12, позволяющий легко определять удельный расход трещины при любом вскрытии , не прибегая к утомительным расчетам.
Рис.13.12. графическая зависимость удельного расхода
Трещины q0 от степени вскрытия пласта
Заметим, в данном примере qmaх=0,88 м2/сут при =1. Однако относительно горизонтального ствола такое утверждение маловероятно и требует дополнительного исследования, и в особенности при наличии подошвенной воды с активным или пассивным напором.
Расчет, произведенный по формулам (13.6.31') и (13.6.34) при найденных коэффициентах добавочных фильтрационных сопротивлений:
и
определенных по таблицам (Прил.1) и формуле (9.5.6), дают Jвк=7,658 и q0г=0,337. Как видим, здесь превышение дебита трещины над дебитом горизонтального ствола составляет в 1,31 раза, что явно занижено.
5. Приток реального газа к вертикальной трещине грп и горизонтальному стволу по нелинейному закону фильтрации
Уравнение линии тока h(x) для I зоны (рис.13.13) будем аппроксимировать функцией, которая называется "Локоном Аньези" [38]:
,
(13.6.35)
где
;
Для нелинейного закона фильтрации известно уравнение:
,
(13.6.36)
где
ω – половина ширины трещины или радиус скважины rс;
м – коэффициент макрошероховатости;
r – плотность газа в пластовых условиях;
m=m(Р) – вязкость газа в пластовых условиях;
– скорость
фильтрации;
K – проницаемость пласта по горизонтали.
Выражая скорость фильтрации через расход, приведенный к поверхностным условиям, отнесенный к длине трещины 2L или ствола скважины
(13.6.37)
и учитывая уравнение состояния реального газа
,
(13.6.38)
из (13.6.36) получаем:
,
(13.6.39)
где
.
(13.6.40)
Рис. 13.13. Схема притока к несовершенной галерее (вертикальной трещине)
Интегрируя уравнение (13.6.39) по давлению от Ртр до Р0 и по х от w до х=h0 (см. рис.13.13, I зона) и замечая, что интеграл
(13.6.41)
представляет собой функцию Лейбензона, внося выражение (13.6.35) в уравнение (13.6.38), после ряда преобразований получаем
,
(13.6.42)
;
(13.6.43)
(13.6.44)
(13.6.45)
.
(13.6.46)
Используя
решение (13.6.20) и (13.6.27) и переходя от
давления к функции Лейбензона, для
притока реального газа в зоне II
при q(t)=q=const
и
получаем:
,
(13.6.47)
где
;
(13.6.48)
При
функция
т. е. процесс фильтрации становится
стационарным и будет описываться
формулой (13.6.47).
Решая совместно (13.6.42) и (13.6.47) и переходя к объемному расходу согласно соотношениям с двухсторонним симметричным притоком к вертикальной трещине (см. рис. 13.13)
,
(13.6.49)
получаем уравнение притока:
,
(13.6.50)
где
,
(13.6.51)
(13.6.52)
Для
определения дебита трещины ГРП надо
решить квадратное уравнение (13.6.50)
относительно Q
и
принять положительный его корень.
Коэффициенты гидравлического сопротивления
А
и
В можно
определить по индикаторным диаграммам,
построенным по результатам исследования
скважин на стационарных режимах, положив
в основу уравнение (13.6.50) при
=0.
Для приближенной интерпретации
результатов исследования в уравнении
(13.6.50) можно перейти к давлениям, усредняя
коэффициенты m(Р)
и Z(P).
Для интерпретации кривой восстановления давления (КВД) после остановки скважины уравнение (13.6.50) записывается в следующем виде:
,
(13.6.53)
где
– значение
функции Лейбензона на забое скважины
в момент ее остановки, Па/с;
a0 – отрезок,
отсекаемый экстраполяционной прямой
на оси ординат зависимости (13.6.53) в
координатах
;
b0 – угловой коэффициент прямой.
Если для конкретной скважины определены значения А и В при наличии вертикальной трещины ГРП по результатам исследования на установившихся режимах, то по отрезку a0 можно оценить длину трещины L.
Как видим, формула для a0 в выражениях (13.6.53) и (13.6.51) представляет собой квадратное уравнение относительно L-1, решение которого после ряда преобразований представляется в виде:
,
(13.6.54)
где
,
(13.6.55)
,
(13.6.56)
Определив L,
из формулы (13.6.52) можно определить
пьезопроводность пласта æ.
В приведенных формулах принимаются
размерности: [А]=[а]=
[В]=
[в]=
[Р
]=
[K]=м2;
[Q]=м3/с.
Если вертикальная трещина распространяется
по всей толщине пласта, что имеет место
практически при отсутствии непроницаемых
пропластков, тогда в приведенных формулах
следует принять С1=0
и С2=0.
Заметим, для
линейного закона фильтрации в формулах
(13.6.38), (13.6.42), (13.6.51) следует принять В=0
и
=0.
Выводы:
1. На основе теории потенциала получены строгие и приближенные аналитические зависимости, позволяющие решать основные промысловые задачи: определение дебитов или депрессий, распределение потенциала (давления) в однородно-анизотропном полсообразном пласте при нестационарном и установившемся притоках жидкости и газа, при линейном и нелинейном законах фильтрации к вертикальной трещине ГРП и горизонтальному стволу скважины при различных начальных и граничных условиях.
2. Если известны аналитические решения для трещины ГРП, то они непорсредственно могут быть преобразованы для горизонтального ствола введением дополнительного фильтрационного сопротивления, обусловленного конвергенцией линий тока.
3. Предлагаемые сложные формулы и уравнения затабулированы и достаточно доступны к использованию; в большинстве случаев не требуют сложной вычислительной техники и вполне пригодны для периодического контроля за работой скважин.