X. движение границы раздела двух жидкостей в пористой среде
Вытеснение нефти водой
При проектировании разработки нефтяных месторождений в условиях водонапорного режима, когда нефть вытесняется в скважины напором краевых вод, необходимо учесть стягивание контура нефтеносности.
С вытеснением нефти водой приходится встречаться и при расчетах деформации водонефтяного контакта. Аналогичные задачи возникают и при эксплуатации газовых месторождении с краевой или подошвенной водой.
Предполагается, что вытеснение «поршневое» и граница раздела двух жидкостей является некоторой поверхностью. При решении задач о вытеснении учитывается различие в вязкостях нефти и воды. Плотности нефти и воды считаются одинаковыми. Это дает возможность рассматривать границу раздела двух жидкостей вертикальной. В общем случае на границе раздела двух жидкостей с различными физическими свойствами
происходит преломление линий тока. Учет этого преломления и составляет главную трудность в точном решении задачи о вытеснении нефти водой (или газа водой). Линии тока не преломляются при прямолинейно-поступательном и радиальном движениях, когда в начальный момент времени они перпендикулярны границе раздела. В этих случаях получены точные решения, в которых жидкости (нефть, вода) принимаются несжимаемыми, пласт — горизонтальным, режим пласта — водонапорным, фильтрация — происходящей по линейному закону.
При прямолинейном движении границы раздела (рис. 66), когда в начальном положении она параллельна галерее, в пласте с постоянными мощностью, пористостью и проницаемостью формула для дебита галереи имеет вид
(X.1)
где l - длина пласта; s — расстояние от контура питания до водонефтяного контакта.
Из приведенной формулы видно, что дебит нефти при заданных постоянных значениях рн и ргвозрастает при продвижении границы раздела, если mн>mв.
Время вытеснения нефти водой в случае прямолинейно-поступательного движения границы раздела подсчитывается по формуле
(X.2)
где s0 — координата, определяющая положение границы раздела в начальный момент времени.
Чтобы найти время полного вытеснения нефти, нужно в формуле (Х.2) положить s = l.
Аналогичная картина наблюдается и в условиях плоскорадиальной фильтрации (рис. 67). В этом случае дебит определяется по формуле
(X.3)
где r — координата, определяющая положение границы раздела нефть — вода в момент t.
Время радиального перемещения границы от начального положения r = ro (при t = 0) до r находится по формуле
(X.4)
Различие вязкости нефти и воды существенно влияет как на время извлечения нефти (газа) из пласта, так и на характер продвижения контура водоносности.
Допустим, что первоначальное положение водонефтяного контакта в пласте А В не параллельно галерее (рис. 68). Для решения задачи о продвижении водонефтяного контакта в указанных условиях используют приближенный метод «полосок», предложенный В. Н. Щелкачeвым. Рассматривается послойное движение частиц. Выделяют несколько узких полосок, и в пределах каждой полоски рассматривают вытеснение как поршневое с контуром водоносности, параллельным галерее. При условии mн > mв скорость точки В больше, чем скорость точки А, отсюда можно сделать вывод, что скорость движения «водяного языка» в наиболее вытянутой точке по мере его движения к галерее (или прямолинейной цепочке скважин) растет быстрее, чем скорость его основания и остальной части контура водоносности.
§ 2. Конус подошвенной воды. Определение предельного безводного дебита скважины
При отборе нефти (газа) из гидродинамически несовершенной по степени вскрытия скважины в пласте с подошвенной водой происходит деформация границы водонефтяного контакта. Образующееся повышение уровня воды называется конусом подошвенной воды (рис. 69). При увеличении дебита конус поднимается, и при некотором предельном значении Q = Qпред происходит прорыв подошвенной воды в скважину. Условием стабильности конуса является равенство градиента давления на вершине конуса удельному весу воды:
(X.5)
Методы расчета предельных безводных дебитов были предложены И. А. Чарным, II. Ф. Ивановым, Н. С. Пискуновым, Д. А. Эфросом, Г. Дж. Мейером, О. А. Гайдаром и др.
Н. А. Чариый, сопоставляя движение нефти при наличии
конуса подошвенной поды с напорным равнодебитным движением нефти в пласте постоянной мощности h(Rо)=hо и используя условие стабильности конуса (Х.5), получил формулу для верхнего значения предельного безводного дебита в однородно-анизотропном пласте, в каждой точке которого значение коэффициента проницаемости в горизонтальном направлении kгор резко отличается от значения коэффициента проницаемости в вертикальном направлении kнерт, в виде:
(X.6)
где = b/hо; q( ) —безразмерный дебит.
Кривые q( ) для различных значений r = R0/xh0 показаны на рис. 70. Здесь х = - коэффициент, учитывающий анизотропию пласта.
На рис. 70 приведены также графики для расчета высоты подъема конуса умах, соответствующей Q1.
Рассматривая предельный случай, в котором вершина водяного конуса находится у забоя скважины, П. Ф. Иванов вывел приближенную формулу для предельного безводного дебита скважины, аналогичную формуле (VI.5) дебита скважины при безнапорном движении
(X.7)
Задача 95
В полосообразном пласте имеет место поршневое вытеснение нефти водой. Первоначальная граница раздела вертикальна и параллельна галерее. Длина пласта Lк = 5 км, длина зоны, занятой нефтью в начальный момент, — 1 км. Динамические коэффициенты вязкости нефти mн = 4 сП, воды mв = 1 сП. Найти отношение дебита галереи в начальный момент эксплуатации и дебита той же галереи, когда весь пласт заполнен нефтью. Определить отношение времени вытеснения нефти водой и нефти нефтью.
Ответ: Qн-в/Qн-н = 2,5; Тн-в/Тн-н =0,325
Задача 96
Определить время продвижения нефти от контура водоносности до скважины в случае плоскорадиального движения по закону Дарси и сопоставить его со временем прохождения того же пути водой. Определить дебит скважины в начальный момент времени и в момент обводнения. Расстояние до контура питания Rк = 10 км, первоначальный радиус водонефтяного контакта rо = 450 м, мощность пласта h = 10 м, пористость пласта m = 20%, коэффициент проницаемости пласта k = 0,2 Д, коэффициенты вязкости нефти mн = 5 мПа×с, воды mв = l мПа×с, давление на контуре питания рк=9,8 МПа (100 кгс/см2), давление на забое скважины рс = б,86 МПа (70 кгс/см2), радиус скважины rс = 0,1 м.
Ответ:Т = 46,2 лет; Тв = 12,5 лет; Qнач = 72,2 м3/сут; Qкон = 283 м3/сут.
Задача 97
Положение водонефтяного контакта в пористом пласте, изображенном в плане на рис. 71, в начальный момент времени показано линией ab, не параллельной галерее. Найти скорость фильтрации в точках а и b.
Определить положение точки а, когда точка b достигнет галереи. Расстояние от галереи до контура питания Lк =10 км, расстояние от контура питания до точки а равно ха = 9200 м, расстояние до точки b xb = 9500 м, коэффициенты вязкости нефти mн = 6 сП, воды mв = 1 сП, коэффициент проницаемости пласта k = 1Д, коэффициент пористости пласта т = 20%, давление на контуре питания рк = 9,8 МПа (100 кгс/см2), давление на забое галереи рг= 6,86 МПа (70 кгс/см2).
Решение.Задачу будем решать приближенным методом полосок, предложенным В. Н. Щелкачевым. Выделим впласте две узкие полоски в окрестностях точек а и b и будем считать, что в каждой из них граница раздела нефть — вода вертикальна и параллельна галерее. В каждой полоске перемещение границы раздела будем рассчитывать по формулам для поршневого прямолинейно-параллельного вытеснения.
Найдем скорости фильтрации в точках а и b.
м/с,
м/с.
Определим время, за которое точка b достигнет галереи:
Найдем положение точки а, когда точка b достигнет галереи:
т. е. точка а будет отстоять от галереи на 360 м и граница раздела нефть—вода примет положение a'b'.
Задача 98
Определить предельный безводный дебит скважины, вскрывшей нефтяной пласт с подошвенной водой, если Rк = 200 м, радиус скважины rс= 10 см, нефтенасыщенная мощность пласта ho =12 м, разность плотностей воды и нефти rв — rн = 0,398 г/см3, динамический коэффициент вязкости нефти mн = 2,54 сП. Пласт считать однородным по проницаемости (х = 1), k = 1 Д.
Задачу решить по формуле Н. Ф. Иванова и по методу, предложенному И. А. Чарным при мощности вскрытой части пласта b, равной 6 м и 2 м.
Решение.Определим предельный безводный дебит по приближенной формуле Н. Ф.Иванова
1)
По графикам И. А. Чарного (см. рис. 70) найдем (r, h) = Qпр/Q0, где
q(16.6; 05) = 0.097, откуда Qпр = 0,097×123 = 11,95 м3/сут
2)
q(16.6; 0,166) – 0,14, Qпр=0,14×123 = 17,2 м3/сут
Как видно из расчетов, формула II. Ф. Иванова дает резко заниженный предельный безводный дебит по сравнению с предельным безводным дебитом по методу И. А. Чарного.
Задача 99
По данным предыдущей задачи определить высоту подъема конуса подошвенной воды по методу И. А. Чарного.
Решение.
1.Определим по графикам И. А. Чарного hmax = ymax/(ho—b} в зависимости от р = R0/xh0 = 16,6 и
;
hmax = 0,81, откуда высота подъема вершины конуса
Уmax = 0,81(12-6)=4,86 м.
2)hmax(16,6; 0,167) = 0,7,
Уmax = 0,7(12-2) = 7 м.
Задача 100
Определить предельно допустимую депрессию при отборе нефти из скважины, вскрывающей пласт с подошвенной водой на глубину b = 12,5 м. Мощность нефтеносной части пласта в отдалении от скважины hо = 50 м, проницаемость пласта k = 0,5Д, плотность воды rв = 1 г/см3, плотность нефти rн = 0,7 г/см3, динамический коэффициент вязкости нефти mн = 2 сП, расстояние до контура питания Rк = 200м, диаметр скважины dc = 21,9 см, пласт считать изотропным (х = kгор/kверт = 1).
Решение.По методу И. А. Чарного определим приближенное значение предельного безводного дебита нефти
По графику зависимости q от р и (см. рис. 70) при значении р = 4 и = 0,25 получаем
q(0.25; 4) = 0,173
Q1 = 1.175×10-2×0.173 = 2.04×10-3 м3/с.
Предельно допустимую депрессию найдем из решения Маскета о притоке к скважине гидродинамически несовершенной по степени вскрытия
здесь значение функции (0,25) = 4,6 (см. рис. 34).