Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Нефтегазовая гидромеханика

..pdf
Скачиваний:
10
Добавлен:
12.11.2023
Размер:
4.17 Mб
Скачать
Рис. 2.62. КВД в координатах Р–t

Pc =

qμ

 

2,246χt

 

qμ

 

2,246χt

 

qμ

ln t (2.364)

 

ln

r2

 

=

 

ln

r2

+

 

4πkh

4πkh

4πkh

 

 

 

c

 

 

 

 

c

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pc = A + B ln t.

 

 

(2.365)

Формула (2.365) является в координатах ( P, ln t) уравне-

нием прямой линии с угловым коэффициентом В и отсекаемым на оси ординат отрезком А (рис. 2.63). На практике форму КВД искажают такие причины, как продолжающийся приток жидкости в скважину после ее ос-

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

 

Рис. 2.63. КВД в координатах

P, ln t

Обработка кривой восстановления давления производится сле-

дующим образом:

 

1.

Строится график КВД в координатах Р–t (см. рис. 2.62).

2.

Строится график КВД в координатах

P ln t (см. рис. 2.63).

 

 

171

3. На КВД в координатах ( P, ln t) выделяется прямолинейный

участок.

4. Определяется уклон выделенного прямолинейного участка (коэффициент В) по координатам точек, соответствующих началу и концу этого участка:

B = y2 y1 .

x2 x1

1.Определяется гидропроводность пласта:

ε= 4πqB .

2.Определяется проницаемость пласта:

k = ε hμ .

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

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

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

2.14.7. Неустановившееся движение газа в недеформируемой пористой среде

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

172

(ρν x )

+

(ρν y )

+

(ρν z )

= −m

∂ρ

.

x

y

z

t

 

 

 

 

Последовательно подставим в него уравнение состояния газа при изотермическом процессе

ρ= ρатм pатм pz ( p)

иуравнение линейного закона фильтрации

ρν = − μk ρ gradp.

После таких

подстановок

 

исходное

уравнение

приводится

в следующему виду:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

p

 

+

 

 

 

 

 

p

 

 

 

p

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

μ ( p) z ( p)

x

 

y

μ

( p) z ( p)

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

p

 

 

 

p

= −

m

 

p

 

 

 

 

(2.366)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

z

μ ( p) z ( p)

z

k t

z ( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.15.Фильтрация жидкости в трещинных

итрещинно-поровых коллекторах

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

Если матрица является плотной непроницаемой породой, то такая среда называется трещинной. Здесь трещины играют роль хорошо проводимых поровых каналов, а матрица – роль укрупненных зерен своеобразной пористой среды.

Если матрица представлена пористой породой, то среда называется трещинно-поровой.

173

Основные характеристики трещинной среды:

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

2.Густота трещин n* – количество трещин на единицу длины некоторого отрезка, перпендикулярного направлению распространения трещин:

n* = Nh ;

3. Раскрытость трещин δ – средний поперечный размер трещин (средняя ширина трещин).

Если трещинный пласт моделируется одной сеткой горизонтальных трещин, одинаково раскрытых и равно отстоящих друг от друга, то

m* = Vт

=

F N δ

= n* δ.

F h

V

 

 

п

 

 

 

Закон фильтрации – уравнение Буссинэ:

u = −

δ2

 

dP

,

 

 

(2.367)

12μ

dx

 

 

 

 

 

 

 

 

 

где u – скорость движения.

 

 

 

 

 

 

 

 

Переход к скорости фильтрации:

 

 

 

 

 

v = m* u = −m*

δ2

 

 

 

dP .

(2.368)

12μ

 

 

 

 

dx

 

Также по аналогии с пористой средой вводится понятие трещинной проницаемости:

v = − k*

 

dP

= −m*

δ2

 

dP .

dx

12μ

μ

 

 

 

dx

174

Тогда

k* = m* δ2 .

(2.369)

12

 

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

Если изменение раскрытости трещин при изменении пластового давления определяется упругими деформациями породы, то такое изменение описывается формулами:

δ= δ

0

δ= δ

0

1

β(Р

 

Р)

,

 

 

 

 

 

0

 

 

 

k

т

= k

1

β

(

Р Р

 

3

,

(2.370)

 

 

т0

 

 

0

)

 

 

где δ0 – раскрытость трещин при давлении Р0; β – комплексный параметр трещинной среды:

β= βт l ,

δ0

где

βт = 1E,

где βт – упругая константа; σ – коэффициент Пуассона; Е – модуль Юнга для породы; l – среднее расстояние между трещинами.

При малых изменениях давления зависимость коэффициента трещинной проницаемости от давления можно считать линейной:

k

т

= k

1 α

(

Р

Р

 

3

при α = 3β.

 

 

т0

0

 

)

 

 

175

Также существует экспоненциальная зависимость коэффициента проницаемости трещинного пласта от давления:

kт = kт0 eα(P0 P).

(2.371)

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

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

 

k (P) ρ(P)

 

 

P =

μ(P)

dP + C.

(2.372)

 

 

 

При неизменных свойствах жидкости функция Лейбензона примет вид

P = ρ k (P)dP + C.

(2.373)

μ

 

С учетом (2.373) функция Лейбензона имеет вид

 

P = μρ kт0 eα(P0 P)dP + C.

(2.374)

Значения функции Лейбензона на контуре питания и у стенки скважины:

P

=

ρkт0

 

еα(Р0 Рк )

,

к

 

μ

 

α

 

P

=

ρkт0

 

еα(Р0 Рс ) .

μ

с

 

 

α

 

176

Массовый дебит

 

 

 

kт0 hρ

 

 

 

 

α(Рк

Рс )

 

 

 

 

 

Q

=

 

1e

 

 

 

 

 

.

 

 

(2.375)

 

 

 

 

 

 

 

 

 

 

 

 

м

 

 

μα

 

 

 

 

lnRк

/ rc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Объемный дебит

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kт0

h

 

 

e

α(Рк

Рс )

 

 

 

 

 

 

Q =

 

1

 

 

 

 

 

.

 

 

 

(2.376)

 

 

μα

 

 

 

 

lnR

 

/ r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

c

 

 

 

 

 

 

При линейной зависимости трещинной проницаемости от дав-

ления

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Р

 

Р )

1α

(

Р

Р )

 

 

 

kт0 h

к

 

 

с

 

 

 

2

 

к

 

с

 

 

 

Q =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(2.377)

 

μ

 

 

 

 

 

 

 

lnR

 

/ r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

c

 

 

 

 

 

 

При кубической зависимости трещинной проницаемости от давления:

Q =

πkт0 h

1

1

β(Р

Р )

4

}

.

(2.378)

2μβlnR

/ r

 

 

{

 

к

с

 

 

 

 

к

c

 

 

 

 

 

 

 

 

Некоторые особенности фильтрации жидкости в трещинно-

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

177

в трещины постоянно поступает жидкость. Интенсивность перетока зависит от этого перепада давления ( p2 p1 ).

q = β** ρμ0 ( p2 p1 ),

где β** – параметр трещинной среды, зависящий от упругих свойств среды и от геометрии трещин.

β** = βср δl0 ,

где βср – коэффициент объемного сжатия матрицы трещинной сре-

ды; l – линейный размера матрицы (расстояние между двумя соседними трещинами); δ0 – раскрытость трещин при начальном пласто-

вом давлении.

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

щинно-поровом коллекторе. Уравнения неразрывности для трещинной и пористой сред с учетом перетока между ними записываются в следующем виде:

 

 

(ρν*x )

 

 

(ρν*y )

 

 

(ρν*z

)

 

 

(ρm* )

 

 

 

 

 

 

+

 

 

 

+

 

 

 

 

 

 

=

 

 

 

q,

 

x

 

 

 

 

y

 

z

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ρν

x

)

 

 

(ρν y )

 

(ρν

z

)

 

(ρm)

 

 

 

 

 

 

 

 

+

 

 

 

+

 

 

 

 

=

 

 

+ q.

 

 

x

 

 

 

 

y

 

z

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

k* ( p

)

ρ ( p )

p

 

 

k* ( p

)

ρ ( p )

p

 

 

 

 

1

 

1

 

+

 

 

1

 

1

 

+

 

 

 

 

 

 

 

x

μ

 

 

x

 

 

y

μ

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

178

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k*

( p )

ρ ( p )

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

μ

 

 

 

 

 

 

 

1

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

ρ

( p )m*

( p

)

− β** ρ0

f ( p f

( p

) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

μ

 

 

 

2

 

 

 

1

 

 

 

 

 

 

k* ( p

 

)

 

( p )

p

 

 

 

 

 

 

 

 

k* ( p

 

)

 

ρ ( p )

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

ρ

 

1

 

 

+

 

 

 

 

 

 

1

 

 

 

 

1

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

μ

 

 

 

 

 

 

1

 

x

 

 

 

 

 

 

 

 

μ

 

 

 

 

1

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k*

( p

)

ρ ( p )

p

 

 

 

ρ ( p )m* ( p )

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

μ

 

 

 

 

 

1

 

 

z

 

 

t

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− β** ρ0

f ( p f

 

( p ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ

 

 

 

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k* ( p

2

)

ρ ( p

)

 

p

 

 

 

 

 

 

 

k* ( p

 

)

ρ ( p )

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

+

 

 

 

 

 

 

 

2

 

 

 

2

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ

 

 

 

 

 

 

2

 

 

x

 

 

 

 

 

 

 

μ

 

 

 

 

1

y

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k* ( p

 

)

 

ρ ( p )

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

μ

 

 

 

 

 

 

 

1

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

ρ

( p )m* ( p

 

)

 

 

+ β** ρ0

 

f ( p

f

( p

)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

μ

 

 

 

 

 

 

2

 

 

1

 

 

 

Для жидкости

f

( p) = p,

для газа

f ( p)

=

p2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ат

 

 

 

 

 

 

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

p =

ρ k0*e−β** ( p0 p)dp + C =

ρ k0*

e

−β** ( p0

p)

 

 

+ C.

 

β**

 

 

μ

μ

 

 

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

 

 

ρ k*

e

−β** ( p0 p)

 

 

ρ k*

e

−β** ( p0

p)

p

=

 

+ C; p

=

 

 

+ C.

 

 

 

β**

 

k

 

μ 0 β**

c

 

μ 0

 

 

179

Дебит скважины в трещинном коллекторе

 

*

 

 

−β** ( p0

p)

 

Q =

2πk0 h

1

e

 

 

 

 

.

(2.379)

 

 

**

 

rк

 

 

 

 

 

 

 

 

 

 

μβ ln

r

 

 

 

 

 

 

 

 

 

c

 

 

 

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

 

 

pк pc

 

*

 

 

−β** ( p0

p)

 

Q =

2πkh

+

2πk0 h

1

e

 

 

 

 

.

(2.380)

 

 

 

 

 

rк

 

 

μ

 

rк

 

 

**

 

 

 

 

 

 

ln

 

 

 

 

μβ ln

 

 

 

 

 

 

 

r

 

r

 

 

 

 

 

 

c

 

 

 

 

c

 

 

 

180