Вытеснение нефти водой и газом из пласта

Рис. 3.1. Построение зависимости Q(Рc) по четырем фактическим точкам

1. Образованием вокруг  скважины области с двухфазной  фильтрацией при забойном давлении  ниже давления насыщения. Чем  больше эта разница, тем больше  радиус области двухфазной (нефть+газ)  фильтрации и, следовательно,  больше фильтрационное сопротивление. 

2. Изменением проницаемости  и раскрытости микротрещин в  породе при изменении внутрипластового  давления.

3. Превышением скоростей  движения жидкости в призабойной  зоне критических значений, при  которых линейный закон Дарси  нарушается.

       Искривление в сторону оси дебитов объясняется неодновременным вступлением в работу отдельных прослоев или пропластков и разными значениями в них пластовых давлений. Это можно пояснить на примере двухслойного пласта. Если пластовое давление в первом пропластке Рк1 а во втором Рк2 причем Рк1 > Рк2, то при всех забойных давлениях Рс, лежащих в пределах Рк1 > Рс > Рк2, приток будет только из первого пропластка. При снижении Рс до величины Рс < Рк2 < Рк1 будут работать оба пропластка, т. е. дебит непропорционально возрастет. Если бы в реальных скважинах процессы фильтрации были обратимы, т. е. расход при отборе был бы равен расходу при поглощении в условиях численного равенства депрессии и репрессии, то это не могло обусловить кривизну индикаторной линии, так как известно, что алгебраическая сумма линейных уравнений всегда дает результирующее линейное.

Однако по целому ряду причин процессы фильтрации необратимые, и  в частности потому, что на забое  всегда имеются взвесь, илистые и  глинистые осадки, которые при  возникновении репрессии закрывают  поры, т. е. работают как обратный клапан: выпускают жидкость из пласта в скважину и задерживают в обратном направлении. Наличие многих прослоев сглаживает ход результирующей индикаторной линии, плавно загибающейся в сторону оси дебитов.

Известны многочисленные факты, когда разница в пластовых  давлениях соседних пропластков, вскрываемых  общим фильтром, доходила до 3 МПа, а  также факты внутрипластовых  перетоков. Кроме того, нелинейный характер индикаторных линий находит объяснение в аномальных неньютоновских свойствах  некоторых пластовых жидкостей.

При любом виде искривления  индикаторной линии ее всегда можно  аппроксимировать уравнением

,                                                                                (3.11)

Это уравнение называют общим  уравнением притока.

При n = 1 уравнение (3.11) описывает прямолинейную индикаторную линию. При 1 > n > 1/2 - индикаторные линии с искривлением в сторону оси P, при n > 1/2 - индикаторные линии с искривлением в сторону оси Q.

 

Случай линейной фильтрации является аналогом ламинарному течению  жидкости в трубной гидравлике. Ламинарное течение с энергетической точки  зрения наиболее экономичное, поэтому  в общем уравнении притока n больше единицы быть не может.

При n = 1/2 приток имеет четко  выраженное турбулентное течение жидкости, когда коэффициент трения l не зависит  от числа Re. Таким образом, индикаторная линия с искривлением в сторону  оси дебитов (n > 1) физически невозможна. В противном случае это означало бы существование течения жидкости с энергетической точки зрения более  экономичного, чем в случае ламинарного  течения. Поэтому индикаторные линии  с искривлением в сторону оси  дебитов, когда n > 1, просто считали  дефектными.

С появлением скважинных дебитомеров  удалось установить истинные причины  искривления индикаторных линий. При n ¹ 1 коэффициент К в общем уравнении  притока теряет свой физический смысл  коэффициента продуктивности и превращается просто в коэффициент пропорциональности или в угловой коэффициент.

При прямой индикаторной линии  коэффициент продуктивности К может  быть найден по любым двум фактическим  точкам как 

,  (3.12)

Зная К, можно определить гидропроводность e = kh/m. Для этого  надо решить формулу (6.10) относительно e

,  (3.13)

Зная по геофизическим  данным или по результатам глубинной  дебитометрии h, а по лабораторным данным m, можно определить проницаемость k в районе данной скважины. Обычно вместо Rк берут половину среднего или  средневзвешенного по углу расстояния до соседних скважин. Для одиночно работающих скважин Rк принимают равным 250 - 400 м, исходя из физических представлений  о процессах фильтрации.

Если имеется ряд фактических  замеров дебитов Qi и соответствующих  этим дебитам замеров забойного  давления Pi, то по этим данным можно  определить все постоянные коэффициенты общего уравнения притока. Поскольку  их три (К, Pк, n), то нужно иметь по крайней мере замеры дебитов и  давлений при трех режимах эксплуатации. Полагая, что индикаторная линия  описывается уравнением вида (3.11), то для каждого режима будем иметь

(3.14)

 

Деля уравнение 1 на 2, сокращая на К и логарифмируя правую и левую  часть, найдем

     или                                 (3.15)

Решая 2 с 3 аналогично, найдем

                                                                               (3.16)

Деля (3.15) на (3.16) и сокращая на n, найдем:

                                                                 (3.17)

Левая часть, число N, известна, так как Q1, Q2, Q3 известны. Неизвестно Pк. Его аналитическое определение  представляет сложности, так как  уравнение (3.17) трансцендентное. Поэтому из (3.17) Pк может быть найдено графоаналитически. Задаваясь несколькими произвольными, но близкими к реальным значениям Pк, вычисляем отдельно М - правую часть равенства (3.17) и строим график зависимости М(Рк) (рис. 3.2). Находя на этом графике точку пересечения А прямой N = const с линией М(Pк), получаем абсциссу этой точки, т. е. такое значение Pк, при котором N = M. Дальнейшее определение коэффициентов К и n трудностей не представляет.

 

                        

               Рис. 3.2. Графоаналитический метод определения величины Рк.

Решая уравнение (3.15) или (3.16) относительно искомого n и подставляя в него уже найденное Pк, найдем

 

 

                                                                                (3.18)

И далее, из любого 1, 2, 3 уравнений, подставляя в них уже найденные Pк и n, определим К.

                                                                                 (3.19)

Поскольку точки Qi и Pi всегда дают разброс, то из уравнений системы (3.14) получится три разных значения К. Рекомендуется принять среднее арифметическое.

Аналогично получится  и с величиной n, из трех значений которых принимают средние арифметическое. Тогда уравнение притока в  данной скважине будет 

                                                                           (3.20)

Для вычисления коэффициентов  К, Pк и n можно применить другие методы, например метод наименьших квадратов.

В подземной гидравлике обосновывается описание процесса фильтрации так называемой двучленной формулой

                                                                        (3.21)

По существу это есть аппроксимация  индикаторной линии полиномом второй степени. При малых Q второе слагаемое  в (3.21) мало. При увеличении Q оно увеличивается пропорционально квадрату Q. Уравнение (3.21) преобразуется в прямую делением всех членов на Q:

                                                                            . (3.22)

                            

          Рис. 6.3. Индикаторная линия при двучленной формуле притока

Имея фактические точки, т. е. Qi, и Pci, и строя по ним график в координатах у = (Pк - Рc) / Q и x = Q, получаем прямую, не проходящую через начало координат и отсекающую на оси  у отрезок a, b - угловой коэффициент  этой линии (рис. 6.3):

                                                                                       (3.23)

При Q ® 0 у = а = (Рк - Рc) / Q, откуда

следовательно,                            (3.24)

Таким образом, в случае аппроксимации  криволинейной индикаторной линии  двучленной формулой необходимо по фактическим  точкам перестроить индикаторную линию  в координатах у = (Рк - Рc)/Q, х = Q, определить отрезок а и далее согласно формуле (3.24) найти искомые параметры: гидропроводность e или проницаемость k. В случае фильтрации в пласте газированной жидкости, т. е. при Рпл < Рнас, параметры пласта определяют по формуле установившегося радиального притока газированной жидкости, которая решается относительно искомой e.

Формула имеет вид 

                                                             (3.25)

где Нк и Нс - функции, учитывающие  изменение газонефтенасыщенности  и фазовой проницаемости для  жидкости при изменении давления. Остальные обозначения прежние.

Физически величина (Нк - Нс) представляет собой ту часть общей  депрессии, которая необходима для  проталкивания к забою скважины только жидкости.

Функция Н через известные  забойные давления вычисляется с  помощью специальных аппроксимирующих формул, различных для различных  интервалов давлений. Обычно методы определения  параметров пласта при установившихся отборах газированной жидкости описываются  в соответствующих инструкциях, как правило, не отличаются надежностью  и применяются редко.

В заключение необходимо заметить, что использование формул установившегося  радиального притока для определения  гидропроводности пласта e = kh/m дает значения этого параметра, характерные для  призабойной зоны пласта, так как  в этой зоне происходит наибольшее падение давления.

 

3. ИССЛЕДОВАНИЕ СКВАЖИН ПРИ НЕУСТАНОВИВШИХСЯ РЕЖИМАХ

Если давление на забое  Рс, а тем более пластовое Рк превышает давление насыщения Рнас, то предполагается, что перераспределение  давления в пласте после любых  возмущений происходит по законам упругого режима. В подземной гидродинамике  рассматривается задача притока  упругой жидкости к скважине в  бесконечном упругом пласте после  ее внезапного пуска или остановки. Решением этой задачи является формула 

                                                   (4.1)

Физическая интерпретация  этой формулы следующая: Dp{r,t) означает изменение давления в упругом  пласте в точке М, удаленной от точки возмущения - скважины на расстояние r через время t после начала возмущения. В данном случае под возмущением  понимается либо пуск скважины с дебитом Q, либо внезапная остановка скважины, работавшей перед этим длительное время, с дебитом Q (Q - расход при стандартных  условиях). При пуске скважины давление в точке М уменьшается на DP по сравнению с первоначальным, а  при внезапной остановке скважины, длительно работавшей с дебитом Q, DP - увеличение давления в точке  М по сравнению с первоначальным, Еi( - х), где x = r2/4c×t - специальная табулированная экспоненциальная функция, значения которой  можно найти в таблицах специальных  функций. Здесь ×c = k/mbx - пьезопроводность, причем bx - приведенный объемный коэффициент  упругости среды (вода, нефть, порода), t - время с момента пуска или  остановки скважины.

Решение (4.1) является строго аналитическим, поэтому оно справедливо для любых радиусов и в частности для радиуса r, равного радиусу скважины rс. В этом случае формула (4.1) будет описывать закон изменения давления на стенки самой скважины и является характеристикой процесса «самопрослушивания» скважины. Таким образом, если остановить скважину и зарегистрировать изменение во времени давления на забое скважины, можно будет найти те параметры пласта, при которых закон изменения DP(t) совпадет с фактически зарегистрированным. Для практического использования формулу (4.1) несколько упрощают. Дело в том, что при исследовании скважин на неустановившихся режимах, т. е. при самопрослушивании, приходится иметь дело с малыми значениями аргумента x = rc2/4ct, так как rс - радиус скважины мал, a t составляет сотни и тысячи секунд.

При малых х экспоненциальная функция 

хорошо аппроксимируется логарифмической функцией Ei( - х) = Ln (х) +0,5772, где 0,5772 - постоянная Эйлера. Поэтому  формулу (6.26) можно переписать следующим  образом:

                                                    (4.2)

Вводя знак минус в скобки и учитывая, что Ln (e) = 1, можем записать:

. Но e0,5772 = 1,781.

Следовательно,

    или      (4.3)

 

Обычно числовой коэффициент  под логарифмом округляют, так что 2,24587 = 2,25. Итак, если остановить скважину, работавшую с дебитом Q, то на ее забое  давление начнет повышаться в зависимости  от времени t согласно формуле (4.3). При этом предполагают, что режим упругий и давление на забое больше давления насыщения.

На формуле (4.3) основана методика исследования скважины при неустаповившихся режимах. Следует отметить, что формула (4.3) предполагает мгновенную остановку скважины (при t = 0, Q = 0). Это равносильно срабатыванию крана или клапана непосредственно на забое скважины. В действительности остановка, например, фонтанной скважины производится на устье путем закрытия задвижки. В НКТ находится газожидкостная смесь, которая после остановки начнет сжиматься под действием возрастающего забойного давления. В затрубном пространстве также произойдет рост давления и сжатие газовой шапки

Мгновенной остановки  скважины не произойдет, а будет  продолжающийся последующий затухающий приток жидкости из пласта в скважину, чего формула (4.3) не предусматривает. Поэтому последующий приток является источником некоторых погрешностей, которые возможно исключить путем специальной обработки фактических данных.

Возвращаясь к формуле (4.3), перепишем ее так, чтобы время t было выделено, а именно

                                . (4.4)

Обозначим:

                                                                         , (4.5)