
- •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)
Воды в условиях напорного притока к несовершенной скважине
Таблица 11.3
Значения ординаты x*=f(a, b) нейтральной линии тока
|
a |
|||||||
b |
0.1 |
0,2 |
0,3 |
0,4 |
0,5 |
0,6 |
0.7 |
0,8 |
0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 |
0,15 0,18 0,23 0,27 0,32 0,37 0,43 0,50 |
— 0,25 0,29 0,34 0,38 0,44 0,50 0,57 |
— — 0,35 0,40 0,45 0,50 0,56 0,63 |
— — — 0,45 0,50 0,55 0,62 0,68 |
— — — — 0,55 0,62 0,68 0,74 |
— — — — — 0,66 0,73 0,80 |
— — — — — 0,76 0,84 |
— — — — — — 0,87 |
В соответствии с
формулой (11.1) для удельного расхода
q0=Q/h
по верхней и нижней частям пласта (см.
рис. 11.7) можно записать следующие
соотношения [9, 16, 26, 27] при
:
;
(11.21)
где
. (11.22)
.
С учетом (11.22) формулы (11.21) принимают следующий вид
(11.23)
Рис. 11.8. К определению ординаты нейтральной линии тока x* в
Зависимости от расположения интервала вскрытия пласта
Безразмерные
предельные дебиты
определяются по табл. 11.1. Чтобы дебит
был одновременно безводным и безгазовым,
необходимо выбрать наименьший расход,
т. е. принять q0=min{q01,
q02}.
Тогда предельный расход нефти через
скважину будет выражаться как
.
(11.24)
Очевидно, этот дебит в общем случае является предельным либо для конуса воды (и меньше предельного для конуса газа), либо для конуса газа (и меньше предельного для конуса воды).
Выражения в правых частях формул (11.23)
,
(11.25)
,
(11.26)
представляют собой соответственно безразмерные предельные безгазовые и безводные плотности расходов. С учетом (11.25) и (11.26) формулы (11.23) принимают вид:
.
(11.27)
Для каждой пары
значений a
и b
и соответствующих им значений ординат
нейтральной линии тока (см. табл. 11.3) по
формулам (11.22) подсчитаны величины
относительных вскрытий
в зависимости от параметров a
и b
и значения параметров r01
и r02
.Затем, с помощью таблицы (см. табл. 11.1)
для предельных дебитов определялись
q1(a,b,r0)
и q2(a,b,r0),
а затем по формулам (11.25), (11.26) рассчитывались
плотности расходов q1
и q2.
Результаты расчетов сведены в табл.
(Прил. 3), которая охватывает все практически
интересные значения параметров a,
b
и r0
[16]. В силу симметрии каждая строка
таблицы дает одновременно значения
безразмерных предельных плотностей
расходов q1
и q2
для соответствующих значений a
и b,
т. е. q1,2(a,b)=q1,2(1–a,
1–b).
По данным таблицы нетрудно построить
сетку кривых зависимостей q1,2=q1,2(r0)
для фиксированных значений пары
параметров a
и b,
т. е. для заданного интервала вскрытия
(b–a).
При конкретных расчетах предельных безводных и безгазовых дебитов поступают следующим образом. По известным параметрам a, b и r0 из таблицы или графиков находят плотности расходов q1 и q2, затем по формулам (11.27) подсчитывают удельные расходы q01 и q02, из которых выбирают наименьшее значение q0=min{q01;q02}, и по формуле (11.24) подсчитывают искомый предельный дебит. Покажем применение метода на конкретных примерах.
Пример 2. Имеется подгазовая нефтяная залежь, подстилающаяся подошвенной водой. Исходные параметры: R0=200 м; h=25 м; Dr1=870 кг/м3; Dr2=200 кг/м3 (в пластовых условиях); mн=2,5 мПа×с; Кr=0,5×1,02×10—12 м2; Dr2=200 кг/м3 (в пластовых условиях); mн=2,5 мПа×с; Kr=0,5×1,02×10—12 м2; æ*=12. Требуется определить безводный и безгазовый дебит при безразмерных параметрах вскрытия: a=0,2; b=0,7 и a=0,2; b=0,5.
1. Определяем значение r0=R0/æ*h=0,66.
2. Из табл. (Прил. 3) находим плотности q1=0,145 и q2=0,290 при a=0,2 и b=0,7.
3. По формулам (11.27) находим удельные расходы:
q01=0,145×870eh=126,15eh;
q02=0,290×200eh=58eh.
4. Так как q02<q01, то выбираем q02 . По формуле (11.24) определяем Q=19,4 м3/сут.
5. Из табл. (Прил. 3) при a=0,2 и b=0,5 находим плотности q1=0,165 и q2=1,0.
6. Удельные расходы составят соответственно:
q01=0,165×870eh=143,55eh;
q02=1.0×200eh=200eh.
7. В этом случае принимаем q01=143,55eh. Тогда расход по формуле (11.24) составит Q»29,2 м3/сут.
Как видим, в этом случае предельный дебит оказался в 1,5 раза больше предыдущего. Таким образом, наибольший дебит зависит от положения интервала вскрытия.
Пример 3. Исходные параметры принимаются для Примера 1, интервал вскрытия, в котором определяемый ординатами b=14,84 м и а=2,34 м, соответствуют безразмерным ординатам: b=b/h=14,84/25»0,60 и a=а/h=2,34/25»0,1.
1. По табл. (Прил. 3) для параметров a»0,1, b»0,60 и r0=200/25=8 при æ*=1 определяем плотности q1»0,02, и q2»0,19.
2. По формулам (11.27) находим удельные расходы:
q01=0,02×870eh=17,4eh;
q02=0,19×200eh=38eh.
3. Выбираем наименьшую плотность q01. По формуле (11.23) находим предельный дебит Q=5,9 м3/сут. Сравнивая его значение с дебитом Q=9,87 м3 / сут., рассчитанное по приближенной методике (см. Пример 1), видим, что последний завышает в данном конкретном примере предельный дебит в 1,66 раза.
4. Для сравнения произведем расчет предельного дебита при тех же исходных данных по методике Курбанова-Садчикова, для чего пересчитаем параметры в обозначениях авторов [20]. Получаем:
Dr1/Dr2=870/200=4.35;
;
.
По графикам [20]
находим q»0,47
и
или h»0,095×25»2,38
м. Предельный дебит по формуле [20]
составляет
.
Завышение предельного дебита по сравнению с расчетным, учитывающим нейтральную линию тока, в данном случае составляет в 1,72 раза.
Пример 4. Принимаются исходные данные, для которых построены графические зависимости размерного предельного безводного и безгазового дебита, рассчитанные потенциометрическим методом [19] и приведенные на рис.14.6д: R0=1000 футов»305 м; h=100 футов»30.5 м; Dr1=500 кг/м3; Dr2=300 кг/м3; Kr=1 д=1 мкм2; mн=1 мПа×с и æ*=1.
Если принять интервал вскрытия l=20 футов»6,1 м, то по графику рис. 14.6д точка пересечения кривых В и b дает Qпр=750 барелей/сут»119 м3/сут и местоположение интервала перфорации а»30 футов»9,15 м (см. рис. 11.7). Следовательно, b=l+а=15,25 м или в безразмерном виде a=0,3 и b=0,5. Параметр r0=R0/æ*h=10. Определим Qпр по уточненному методу. По табл. (Прил. 3) находим плотности расходов q1(a,b,r0)=q1(0,3;0,5;10)»0,18 и q2(a,b,r0)=q2(0,3;0,5;10)»0,45. Затем, по формулам (11.27) определяем удельные расходы: q01=0,18×600eh=108eh и q02=0,45×300eh=135eh. Для наименьшего удельного расхода q02 по формуле (11.24) находим Qпр»109 м3/сут. В данном случае расхождение между двумя методами несущественное и составляет 8,4%.
Пример 5. За исходные примем данные в примере Курбанова-Садчикова [20]: R0=200 м; h=10 м; Dr1=700 кг/м3; Dr2=300 кг/м3; mн=2 мПа×с; Kr=0,5×1,02×10-12 м2; æ*=5; b–a=2 м; d=3,9 м (см. рис. 11.7). Из условия задачи имеем численные значения параметров a»0,3 и b»0,5 и r0=4. По табл. (Прил. 3) определяем безразмерные плотности расходов: q1»0,213 и q2»0,557. Удельные расходы составляют: q01»0,149eh и q02»0,167eh. Подсчитывая предельный дебит по формуле (11.24) по наименьшему удельному расходу q01, получаем Q»6,1 м3/сут.
По расчетам авторов [20, 21] этот дебит равен Q»4,33 м3/сут., т.е. отклонение составляет порядка 40%. такое расхождение, очевидно, объясняется тем, что авторы при решении задачи делают допущение, что нейтральная линия тока проходит через середину интервала вскрытия (см. рис. 11.7 и 11.9) при любом его положении, тогда как уточненная методика определяет положение нейтральной линии тока x* в зависимости от положения интервала вскрытия a и b. Заметим, что в своей предпосылке при решении задачи несовершенная скважина считалась линией стоков с постоянным удельным расходом. В действительности на скважине должен быть постоянным потенциал. Физически ясно, что картины линий тока будут отличаться несущественно, а, следовательно, положения горизонтальных линий тока будут близки друг к другу [9].
Метод Курбанова-Садчикова и предлагаемый уточненный метод решения задачи конусообразования имеют следующие преимущества перед потенциометрическим и другими существующими методами: они универсальные, т. е. расчетные зависимости представлены в безразмерном виде и применимы как для однородных, так и для однородно-анизотропных пластов; графические решения даны в широком диапазоне безразмерных параметров вскрытия (a, b) и радиуса контура питания (R0) и охватывают все практически интересные случаи; технически удобны и просты, не требуют сложной вычислительной техники.
11.3.5. Методика расчета предельной депрессии, обеспечивающей предельный безводный и безгазовый дебит. При вскрытии нефтяной оторочки в интервале (b–a) средние значения потенциала (давления) для каждой части пласта (см. рис. 11.7) будут равными. Вдоль всей вскрытой продуктивной толщи среднее значение потенциала можно определить по формуле
,
(11.28)
где
,
– средние
значения потенциалов вдоль скважины
по вскрытым толщинам (d–a)
и (b–d)
соответственно (см. рис. 11.7). В безразмерных
параметрах формула (11.28) представляется
в виде
.
(11.29)
Выразим разность потенциалов для каждой из частей пласта, принимая в качестве расхода предельные дебиты:
,
(11.30)
где в соответствии с формулами (11.27) имеем
.
(11.31)
Здесь q1 и q2 – безразмерные плотности расхода, рассчитанные по формулам (11.25) и (11.26) и затабулированные (см. табл. Прил. 3).
С учетом (11.30) и (11.2) имеем
;
(11.32)
.
(11.33)
Внося (11.32) и (11.33) в формулу (11.29), переходя к давлению и учитывая (11.31), получаем формулу для предельной депрессии DРпр=Р0–Рс в следующем виде:
´
(11.34)
где
,
,
r1,
r2
подсчитываются по формулам (11.22). Функция
затабулирована (см. табл. 11.2).
Пример 6. Рассчитать предельную депрессию для исходных данных Примера 2 при параметрах вскрытия a=0,2 и b=0,5.
1. По табл. 11.3
определяем ординату нейтральной линии
тока:
.
2. По формулам
(11.22) находим:
.
3. По таблице
(11.2) определяем значение функции:
=
=
.
4. Так как
(см. Пример 2), следовательно,
.
5. Подсчитываем DРпр по формуле (11.34): DРпр»0,243 МПа
Предельная депрессия может быть подсчитана более точно, если учесть добавочные фильтрационные сопротивления, обусловленные не только частичным вскрытием С1, но и нарушением линейного закона С2 и перфорацией колонны С0. Тогда определяя разность потенциалов по формулам (11.32) и (11.33) с учетом (11.4) и (11.22) для каждой части пласта и внося их выражения в формулу (11.29), после ряда преобразований получаем:
´
,
(11.35)
где
(11.36)
.
(11.37)
С1, С2, С0 — добавочные фильтрационные сопротивления в формулах (11.36) определяются по соответствующим формулам или таблицам [10, 12, 18, 28], приведенным в Гл.9 и Прил.1.
Пример 7. Применяются исходные данные Примера 6. Добавочные данные: фильтрация происходит по закону Дарси (С2=0); скважина перфорирована с плотностью m=3 отв/пог. м, глубина пулевого канала l0=0,3456 м, диаметр пулевого канала 2r=0,0127 м. Требуется определить предельную депрессию.
1. По формулам
(11.37) находим:
»7,84;
»18,13.
2. По формуле (9.5.6) находим: С01=60,07; С02=6,00.
3. По формуле (11.36) находим: S1=67,92; S2=24.13.
4. По формуле (11.35) имеем DРпр=0,653 МПа.
Как видим, предельная депрессия с учетом добавочных фильтрационных сопротивлений, обусловленных перфорацией, оказалась больше в 0,653/0,243»2,7 раза. Следует заметить, что при известных условиях (большой диаметр пулевых каналов, большая глубина проникновения и оптимальная плотность перфорации) значения С01 и С02 могут оказаться отрицательными, что влечет к уменьшению расчетных значений DРпр.