Задача о притоке реального газа к несовершенной скважине при нелинейном законе фильтрации является весьма сложной и до сих пор не получила точного аналитического решения. В работах Е.М. Минского, А.Е. Хейна, Г.А. Зотова, С.М. Тверковкина и др. рассматривалась данная задача в приближенной постановке. Здесь рассматривается задача о притоке реального газа к несовершенной скважине в однородно-анизотропном пласте, т. е. с учетом анизотропии, а также предлагается несколько иной подход к расчету фильтрационных сопротивлений, обусловленных несовершенством скважины по степени вскрытия.
В работе Е.М. Минского [1] показано, что коэффициент фильтрационного сопротивления как при линейном, так и при квадратичном законе фильтрации зависит только от геометрии потока. В связи с этим к выводу уравнения притока газа можно подойти следующим образом. Для нелинейного закона фильтрации имеем уравнение
(1)
где
– плотность газа.
Геометрия потока, очевидно, будет определяться функцией в области пространственного потока
(рисунок 1). Вся трудность решения состоит в нахождении уравнения кривой
, ограничивающей область потока, или, другими словами, уравнения линии тока. Размер зоны пространственного движения будет зависеть от многих факторов, например, не только от геометрии пласта (
,
,
), но и от анизотропии пласта æ*, дебита
, градиента давления (
) и т. д. Следуя И.А. Чарному [5], примем радиус зоны пространственного притока
. Будем аппроксимировать упомянутую линию тока уравнением вида [4]
(2)
где
– некоторая функция, зависящая от несовершенства скважины по степени вскрытия, геометрии пласта и скважины, анизотропии пласта.
Рис.1. Двухзонная схема притока газа к несовершенной скважине при нелинейном законе фильтрации
Умножая левую и правую части уравнения (1) на , применяя двухзонную схему притока (см. рисунок 1), учитывая уравнение состояния реального газа, уравнение (2) и интегрируя в соответствующих пределах по давлению и радиусу, после некоторых преобразований получаем известную двучленную формулу притока:
(3)
где
(4)
(5)
;
(6)
;
. (7)
Выражение для представляется сложной функцией, выраженной суммой рядов от 1 до
и зависящей от параметров
и
[4], где
;
;
. (8)
Коэффициент фильтрационного сопротивления , обусловленный относительным вскрытием пласта
, определяется формулой [2]
(9)
где
– некоторая функция, связанная с распределением потенциала скорости фильтрации, вызванного работой несовершенной скважины, рассчитана на ЭВМ, затабулирована и представлена графиками (рисунок 2).
Из совместного решения (6) и (9) определена функция , рассчитана на ЭВМ, затабулирована в широком диапазоне параметров и представлена графическими зависимостями (рисунок 2). При найденных значениях
функция
, также была рассчитана на ЭВМ в широком диапазоне параметров и затабулирована.
Рис.2. Зависимость функции , связанной со средним значением потенциала скважины, от относительного вскрытия пласта
Методика расчета безводных дебитов в случае притока реального газа к несовершенной скважине по нелинейному закону фильтрации при наличии подошвенной воды рассматривалась в работе [6]. Однако сама задача о предельных безводных дебитах в точной постановке не решена из-за того, что не известно уравнение границы раздела двух жидкостей при наличии конуса подошвенной воды.
Приведем приближенное решение этой задачи, используя приближенное уравнение границы раздела, когда конус воды находится в предельно-устойчивом положении [2]:
;
;
(10)
где
– безразмерный предельный безводный дебит, определяемый по известным формулам, графикам или таблицам [2].
Переменная толщина пласта, ограничивающая область пространственного притока, как это следует из (10), выражается формулой вида
. (11)
Умножая левую и правую части уравнения (1) на , учитывая, что объемный предельный дебит Q=Q0q(ρ0,
) и (11), получаем формулу притока, выраженную через предельный безразмерный дебит, из которой легко определить предельную депрессию
(12)
После интегрирования и некоторых преобразований уравнение притока примет вид (3), где A и В определяются по формулам (4) и (5), в которых коэффициенты представляются выражениями:
;
; (13)
; (14)
где
– коэффициент за счет перфорации,
– интегральный логарифм, который связан с интегральной показательной функцией
зависимостью
. (15)
При >1 интеграл (15) расходится в точке
=1. В этом случае
представляет значение несобственного интеграла.
Заметим, что из полученных формул притока для несовершенной скважины как частный случай вытекают формулы для линейного закона фильтрации.
Пример: Скважина работает при наличии устойчивого положения границ раздела. Исходные данные R0=1000 м; æ*=1; rс=0,1 м; h0=10 м; b=4 м. Тогда имеем следующие безразмерные параметры: ρ0=100; =0,4;
=104. Требуется определить фильтрационные сопротивления, обусловленные наличием конуса воды.
По графикам [3] определяем безразмерный предельный безводный дебит: q(ρ0,)=q(100;0,4)≈0,4. По формуле (13) находим
≈2,3 и по формуле (15) подсчитываем С2≈2,9. Как видим, значения
и С2 оказались сравнительно небольшими, это объясняется тем, что предельный безводный дебит очень мал, а он и определяет геометрию потока, т. е. форму конуса подошвенной воды. После чего предельная депрессия
легко подсчитывается по формуле (12).
Рецензенты:
Грачев С.И., д.т.н., профессор, заведующий кафедрой «Разработка и эксплуатация нефтяных и газовых месторождений», Институт геологии и нефтегазодобычи, ФГБОУ ВПО ТюмГНГУ, г. Тюмень;
Леонтьев С.А., д.т.н., профессор, профессор кафедры «Разработка и эксплуатация нефтяных и газовых месторождений», Институт геологии и нефтегазодобычи, ФГБОУ ВПО ТюмГНГУ, г. Тюмень.