Дифференциальное уравнение движения

Лекция 2

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

Для подземной гидромеханики характерно изотермическое изменение параметров вследствие значительных величин удельной поверхности коллекторов и их теплоёмкости. Т.о. для таких процессов можно не рассматривать уравнение энергии и ограничиваться уравнениями балланса массы (неразрывности) и движения.

Уравнение энергии необходимо рассматривать в локальных областях призабойной зоны из-за значительных перепадов давления, проявления дроссельного эффекта, а также при применении тепловых методов повышения нефте-газоотдачи.

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

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

Рассмотрим фильтрацию флюидов в пористых средах, принимая во внимание линейный закон Дарси.

Выделим два сечения – первое на расстоянии S от начала отсчета вдоль линии тока, второе – на расстоянии DS от первого (рис. 1).

Движение флюида происходи в направлении возрастания координаты S. В сечении с координатой S обозначим приведенное давление через p*(S, t), в сечении координат S + DS – через p*(S + DS,t), используя формулу ,

получаем

, (20) Рис. 1. Трубка тока

или перейдем к пределу при ,

, (21)

Знак (-) в правой части означает, что приведенное давление падает по движению жидкости, т.е. градиент приведенного давления отрицателен .

Формула (21) справедлива только для изотропной среды, для которой характерно постоянство проницаемости по всем направлениям в окрестности рассматриваемой точки. Однако с переходом от точки к точке пласта проницаемость может и изменяться, таким образом (модель изотропного неоднородного пласта).

Запишем уравнение (21) в проекциях на оси координат x, y, z. Если обозначить через ,, единичные векторы вдоль осей координат, вектор скорости фильтрации можно записать в виде

, (22)

, (23)

тогда

, (24)

или в проекциях на оси координат

, , , (25)

если ось z направлена вверх и дифференциальные уравнения движения примут вид

, , , (26)

в векторной форме . (27)

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

где S – координата, взятая вдоль линии тока по движению жидкости.

В векторной форме двучленный закон выведен из теории размерностей, в виде

(29)

В прекциях на оси координат имеем

, (30)

При фильтрации неньютоновских вязкопластичных жидкостей, а также при фильтрации с очень малыми скоростями имеет место закон фильтрации (5), который отличается от закона Дарси наличием предельного градиента , по достижении которого начинается движение. В векторной форме закон фильтрации с предельным градиентом выведен из теории размерностей и имеет вид . (31)

; (32)

в проекции на оси координат:

; (33)

;

Лекция 3.

Вывод дифференциального уравнения неустановившейся фильтрации однородного флюида по закону Дарси. Функция Л. С. Лейбензона.

Для вывода дифференциального уравнения неустановившейся фильтрации используем уравнение неразрывности

или

(34)

Сумма в скобках в левой части уравнения (34) представляет собой дивергенцию вектора скорости фильтрации и кратко записывается таким образом:

, (35)

поэтому уравнение (34) можно записать в виде:

. (36)

Уравнение (34) (или 36) справедливо только в том случае, если внутри объема нет источников или стоков, выделяющих или поглощающих флюид, не происходит химических реакций, фазовых превращений и т.д.

И уравнения движения

(37)

В уравнении (11) не будем учитывать силу тяжести.

Введем функцию (функцию Лейбензона), тогда дифференциал этой функции равен:

, (38)

тогда

, (39)

т. к. функция Лейбензона и давление зависит от координат x, y, z и времени t, то (38) можно записать в развернутом виде, используя понятие полного дифференциала функции от многих переменных:

Сравнивая коэффициенты при x, y, z получаем:

, , , (40)

Запишем выражение для составляющих массовой скорости фильтрации, умножив правую и левую части уравнения (37) на плотность и используя соотношения (40):

, (41)

Подставим выражение (41) в уравнение неразрывности (34), получим:

(42)

или

, (43)

где - оператор Лапласа от функции Лейбензона (39).

Уравнение (42) справедливо для неустановившегося движения однородного флюида в однородной пористой среде по закону Дарси.

При установившейся фильтрации и будет удовлетворяться уравнение Лапласа для функции Лейбензона:

(44)

При k = const, m = const, и , тогда можно ввести функцию Лейбензона в виде:

. (45)

Тогда дифференциальное уравнение неустановившейся фильтрации примет вид:

. (46)

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

Для несжимаемой жидкости r о = const, тогда

, (47)

т. е. функция Лейбензона пропорциональна давлению.

Для упругой жидкости:

, (48)

т. е. имеем тот же вид, что и для несжимаемой жидкости.

Для совершенного газа с уравнением состояния

, (49)

получаем

, (50)

т. е. функция Лейбензона пропорциональна квадрату давления.

Для реального газа с уравнением состояния

, (51)

тогда

, (52)

т. е. функция Лейбензона записывается в виде интеграла.

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

Лекция 4.

1. Установившиеся потоки флюида в пористой среде.

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

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

Примеры.

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

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

в) в лабораторных условиях при течении через цилиндрический керн или прямую трубу постоянного сечения, заполненную пористой средой.

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

Примеры.

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

б) Гидродинамически- несовершенная скважина - вблизи скважины линии тока искривляются и поток можно считать плоско-радиальным только при некотором удалении от скважины.

в) Круговая батарея эксплуатационных скважин - поток плоско-радиален на некотором удалении, т.к. жидкость движется как бы к укрупнённой скважине радиуса, равного радиусу окружности батареи.

3. Радиально-сферический поток. Траектории всех частиц жидкости - прямолинейные горизонтальные прямые, радиально сходящиеся к центру полусферического забоя; изотахи и эквипотенциальные поверхности перпендикулярны траекториям и образуют сферические поверхности. Скорость фильтрации в любой точке потока является функцией только расстояния этой точки от центра забоя. Следовательно, этот вид фильтрационного потока также является одномерным.

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

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

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

Задача исследования установившегося фильтрационного потока заключается в определении дебита (расхода), давления, градиента давления и скорости фильтрации в любой точке потока, а также в установлении закона движения частиц жидкости (или газа) вдоль их траекторий и в определении средневзвешенного по объёму порового пространства пластового давления.

При условии вытеснения флюида из пласта или его нагнетания в пласт через галерею или скважину условимся принимать за координату произвольной точки пласта расстояние r до этой точки от:

1) галереи (для прямолинейно- параллельного потока);

2) центра контура скважины в основной плоскости (плоскости подошвы пласта) фильтрации (для плоско-радиального потока);

3) центра полусферического забоя скважины (для сферически-радиального потока).

2.Характеристики одномерных фильтрационных потоков

жидкости и газов.

Для расчета перечисленных характеристик одномерных фильтрационных потоков жидкости и газа можно использовать два подхода. Первый из них – вывод дифференциальных уравнений и их решение отдельно для прямолинейно-параллельного, плоскорадиального и радиально-сферического потоков жидкости и газа. Второй – вывод обобщенного уравнения одномерного течения флюида в недеформируемой трубке тока переменного сечения с использованием функции Лейбензона и получение из него конкретных формул применительно к различным схемам фильтрационных потоков. Второй подход более эффективен, позволяет исходить из обобщенных характеристик течения, он был использован и при выводе следующих характеристик:

Прямолинейно – параллельный фильтрационный поток.

Площадь поперечного сечения ; на контуре питания x₁=0, P₁=P_K, на галерее x₂=L, P₂=P_Г;

Схема прямолинейно параллельного течения в пласте

(53, 54, 55)

Плоскорадиальный фильтрационный поток

От координаты S переходим к r, отсчитываемой от центра скважины. Для добывающей скважины , , площадь фильтрционной поверхности - боковая поверхность цилиндра; на контуре питания r 1 = R k,P₂=P_K на забое скважины r 2 = r c,P₂=P_C.

Схема плоско – радиального потока в круговом пласте.

(56)

, (57)

(58)

Радиально – сферический фильтрационный поток.

В этом случае для добывающей скважины с полусферическим забоем имеем: , , - площадь поверхности полусферы с радиусом r, r 1 = R k, P₁=P_K, r ₂= r _C, P₂=P_C.

(59)

; (60)

. (61)

3. Анализ одномерных потоков несжимаемой жидкости и газа.

Рассмотрим конкретные модели флюидов – несжимаемую жидкость и совершенный газ. Выпишем для них формулы для расчета основных характеристик одномерных фильтрационных потоков. Сопоставление этих формул позволит оценить эффект сжимаемости при прочих одинаковых условиях.

Прямолинейно – параллельный поток несжимаемой жидкости и совершенного газа

Подставим в основные расчетные формулы (53), (54), (55) выражение функции Лейбензона (для несжимаемой жидкости)и (для совершенного газа) , а также на контуре и на галерее (аналогично и для совершенного газа)

Расчетные формулы для прямолинейно-параллельного потока несжимаемой жидкости и совершенного газа

Характеристика	Несжимаемая жидкость	Совершенный газ
Функция Лейбензона
Распределение давления по пласту, 0 £ x £ L	(62)	(63)
Массовый расход Qm	(64)	(65)
Массовая скорость фильтрации	(66)	(67)
Объемный расход Q	(68)
Скорость Фильтрации (объемная)	(71)	(72)
Средневзвешенное давление	(73)	(74)
Время движения отмеченных частиц t	(75)	(76)
Время продви – жения до галереи Т	(77)	(78)

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

Плоскорадиальный фильтрационный поток

Модель флюида

Характеристика	Несжимаемая жидкость	Совершенный газ
Распределение давления по пласту	(79) (81)	(80) (82)
Массовый расход Q m	(83)	(84)
Массовая скорость фильтрации rW	(85)	(86)
Объемный расход Q	(формула Дюпюи) (87)	(88) (89)
Объемная скорость фильтрации	(90)	(91)
Средневзвешенное давление	(92)	(93)
Время движения отмеченных частиц	(94)	________________
Время движения частицы от контура до забоя Т	(95)	(96)

Рис.2 Кривые распределения давления в плоскорадиальном потоке:

1 – для жидкости, 2 – для газа.

Для несжимаемой жидкости давление меняется вдоль координаты r по логарифмическому закону (Рис. 2, кривая 1). Вращение кривой p (r) в пространстве вокруг оси скважины образует поверхность, называемую воронкой депрессии.

Зависимость дебита от перепада давления называется индикаторной линией. В потоке жидкости по закону Дарси индикаторная линя – прямая (Рис. 3).

Вид индикаторной линии не зависит от геометрии потока и определяется только законом фильтрации. Отношение массового дебита скважины Q m к перепаду давления D р называется коэффициентом продуктивности скважины k. Рис. 3.

Из (60) следует, (для жидкости):

, (97)

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

, (98)

, (99)

где - приведенный радиус скважины.

Приведенный радиус скважины – это радиус гидродинамически совершенной скважины, которая обеспечивает при равных прочих условиях такой же дебит, как гидродинамически несовершенная скважина,

Радиально – сферический фильтрационный поток несжимаемой жидкости и совершенного газа.

Модель флюида

Характеристика	Несжимаемая жидкость	Совершенный газ
Распределение давления р (r)	(100)	(101)
Массовый расход Q m	(102)	(103)
Массовая ско –рость фильтрации	(104)	(105)
Объемный расход	(106)	(107)
Объемная скорость фильтрации	(108)	(109)
Время движения частиц t	(110)	______________
Время движения от контура до забоя	(111)	(112) где

Лекция 5.

Фильтрация по степенному закону

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

При плоскорадиальном движении закон приобретает вид:

, , (113)

где с и n – константы, определяемые из опыта или по результатам исследования скважины.

Для вывода формул введем функцию давления для несжимаемой жидкости и совершенного газа соответственно получаем:

, (114)

. (115)

Расчетные формулы для плоскорадиального течения несжимаемой жидкости и газа по степенному закону

Модель флюида

Характеристика	Несжимаемая жидкость	Совершенный газ
Распределение давления р (r)	(116)	(117)
Массовый расход	(118)	(119)
Распределение давления р (r)	(120)	(121)
Массовый расход Q m	(122)	(123)

Массовый расход для жидкости пропорционален депрессии в степени 1/ n, поэтому индикаторная линия при 1< n < 2 будет иметь вид выпуклой к оси дебитов степенной кривой с дробным показателем меньшим 2. В случае фильтрации по закону Краснопольского, индикаторная линия является параболой второго порядка.

Рис. 4. Индикаторнаые линии,

соответствующие различным

законам фильтрации жидкости.

Рис. 4.

Фильтрация по двучленному закону.

Модель флюида

Характеристика	Несжимаемая жидкость	Совершенный газ
Функция Лейбензона
Распределение Давления	(124)	(125)
Уравнение притока к скважине	(126)	(127)

Из (126) и (127) видно, что индикаторная линия, построенная в координатах для жидкости идля газа, является параболой (Рис. 5, 6).

Рис. 5. Индикаторная линия Рис. 6. Индикаторная линия

при фильтрации жидкости при фильтрации газа по

по двучленному закону. двучленному закону.

Уравнение притока к скважине для несжимаемой жидкости имеет вид:

(128)

для газа

(129)

где

(130)

(131)

А, В, А1, В1, - коэффициенты фильтрационного сопротивления, являются постоянными для данной скважины.

Скважины исследуют на 5 – 6 режимах (однако ка показывают исследования и результаты обработки индикаторных линий этих замеров недостаточно, необходимо увеличить число замеров для более точного определения коэффициентов фильтрационного сопротивления. Кроме того можно упомянуть об аномальных видах индикаторной линии, о случаях кольматации и наоборот раскольматирования при высоких отборах).

Затем скважину закрывают и давление на забое остановленной скважины принимают за контурное давление рк.

Уравнения (128) и (129) можно представить соответственно к уравнению прямой:

(132)

(133)

Рис. 7. График зависимости от при фильтрации газа по двучленному закону

Коэффициент А – отрезок, отсекаемый на оси ординат, В – тангенс угла наклона прямой к оси абсцисс. (Рис. 7).

По значениям коэффициентов А и В определяют коллекторские свойства пласта: коэффициент проницаемости (эффективный), эффективную мощность пласта, коэффициент гидропроводности:

Для нефтяной скважины

(134)

Для газовой скважины

. (135)

Лекция 6.

Гидродинамические методы исследования скважин и пластов

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

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

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

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

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

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

Методики интерпретации КВД с учетом притока позволяют определять параметры призабойной и удаленной зон пласта. Решению данной задачи посвящены работы Щербакова, Чарного-Умрихина, Баренблатта и др, Борисова, Каменецкого, Чекалюка и др.

Для пласта неограниченных размеров при мгновенном прекращении притока жидкости к скважине после ее остановки (или после пуска с постоянным дебитом) повышение давления в скважине с высокой степень точности определяется по формуле Тсейса:

, (136)

Если мало и в связи с этим нельзя пренебречь суммой , практически , тогда

, (137)

где , - коэффициент упругоемкости пласта (приведенный коэффициент сжимаемости жидкости и породы) ; - коэффициент объемного сжатия жидкости , - коэффициент объемного сжатия пласта ; - изменение давления в скважине после ее остановки, или - изменение давления после пуска скважины.

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

, (138)

а по отрезку, отсекаемому этой прямой на оси ,

, (139)

определим пьезопроводности пласта

. (140)

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

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

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

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

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

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

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

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

Разработанные за последние время методы исследований скважин и пластов при установившихся режимах эксплуатации имеют ограниченность применения в части: определения фильтрационных характеристик пласта; оценки состояния ПЗП; оценки состояния фильтровой части скважины; выявления динамики фильтрационных характеристик пласта и продуктивности скважин; количественной оценки критических дебитов скважин в условиях выноса песка и пластовой воды, а также режимов энергосбережения.

Оценке состояния ПЗП посвящены работы С.М. Тверковкина, Г.А. Зотова, Гриценко А.И., В.К. Зинченко, О.М. Ермилова, С.Н. Бузинова, И.Д. Умрихина, А.И. Петрова, В.Н. Васильевского и др.

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

Анализ данных гидродинамических исследований в соответствие с предлагаемой инструкцией позволяет:

- определять работающие интервалы пласта;

- определять фильтрационные свойства пласта (коэффициент проницаемости);

- оценивать степень кольматации забойных фильтров и высоту песчаной пробки на забое скважины;

- оценивать изменение во времени фильтрационных характеристик призабойной зоны пласта;

- прогнозировать продуктивную характеристику эксплуатационных скважин;

- оценивать интенсивность выноса песка из пласта в скважину и предельно-допустимые дебиты скважин.

Методика обработки результатов исследований

Уравнение притока к забою скважины при нелинейном законе
фильтрации имеет вид:

для жидкости

; (141)

для газа

, (142)

где - пластовое давление; - давление на забое работающей скважины; Q - дебит жидкости (газа) при стандартных условиях; , -коэффициенты фильтрационного сопротивления.

Уравнения (141), (142) записывают также в виде:

для жидкости

; (143)

для газа

, (144)

где для жидкости и соответственно для газа

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

В промысловой практике форма индикаторной линии может существенно отличаться от стандартной в силу ряда причин геологического, технологического и технического характера.

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

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

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

Типовые индикаторные линии (для газа), выявленные по данным исследований, приведены на рисунке 8.

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

Методика проведения испытаний, замера дебита скважины и определения забойного давления должны соответствовать требованиям инструкции или руководству по исследованию скважины.

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

ТИПОВЫЕ ГИДРАВЛИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПЛАСТА

	Гидравлическая характеристика	Индикаторная линия	Возможные причины аномалии
а			Стандартная индикаторная линия 1. - 2. -
б			Занижено пластовое давление
в			Завышено пластовое давление
г			Кольматация каналов фильтрации, вынос песка
д			Включение новых каналов фильтрации, очистка ПЗП и скважины
е			Работает 2 пласта и более

Рис. 8.

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

1. Строится гидравлическая характеристика пласта в координатах

или и методом графической экстраполяции кривой до оси ординат определяется условное пластовое давление или . Данные исследования обрабатываются по уравнению:

для жидкости

; (145)

для газа

; (146)

2. Используя значение замеренного пластового давления, строят графическую зависимость для жидкости или для газа , и методом графической экстраполяции до оси ординат определяют отрезок . В дальнейшем при обработке данных исследований используют уравнение:

для жидкости

; (147)

для газа

, (148)

при этом условное пластовое давление равно

, (149)

или

. (150)

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

3. Условное пластовое давление рассчитывается с использованием метода наименьших квадратов.

Решается система линейных уравнений (на ЭВМ по стандартной программе):

- для жидкости , (151)

или

-для газа , (152)

где ; - число режимных точек, откуда определяются все искомые параметры , и .

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

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

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

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

Диагностика выполняется на основе анализа коэффициентов фильтрационного сопротивления.

Теоретические значения коэффициентов А и В для случая плоскорадиальной фильтрации нефти (газа) с учетом гидродинамического несовершенства скважины по степени и по характеру вскрытия пласта рассчитываются по формулам:

для жидкости

; (153)

; (154)

для газа

; (155)

, (156)

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