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

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

..pdf
Скачиваний:
10
Добавлен:
12.11.2023
Размер:
4.17 Mб
Скачать

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

Рис. 2.19. Индикаторные диаграммы при нелинейных законах фильтрации

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

P = A Q + B Q2 .

(2.59)

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

Коэффициент В зависит от конструкции забоя скважины. Для его определения И.А. Чарным приводятся следующие формулы:

– для совершенной скважины:

B =

120

dэф 2

r

 

ρ

 

(2.60)

 

 

 

 

 

c

 

 

,

m

 

k

f

2

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

61

где dэф – эффективный диаметр зерен горной породы; k – коэффициент проницаемости, мкм2; ρ – плотность жидкости, г/см3; f – площадь вскрытия забоя;

– для скважины, несовершенной по характеру вскрытия:

B =

40ε dэф 2

D

 

ρ

 

(2.61)

 

 

 

 

 

 

 

 

,

m

 

k

f

2

 

 

k

 

 

 

 

 

где f – суммарная площадь перфорационных отверстий; D – диаметр перфорационных отверстий; ε – коэффициент, зависящий от проникновения пуль в породу, изменяется в пределах 0,15–0,4. Значение

ε= 4 соответствует нулевому внедрению пуль в породу;

для скважин с щелевым фильтром:

B =

120ε dэф 2

t

 

ρ

 

(2.62)

 

 

 

 

 

 

 

 

,

m

 

k

f

2

 

 

k

 

 

 

 

 

где f – суммарная площадь щелей; t – ширина щелей; 0,25 < ε <0,5. Двухчленная формула притока является физически наиболее

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

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

62

2.5.5. Исследование скважин методом установившихся отборов (при установившихся режимах)

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

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

Обработка прямолинейной индикаторной диаграммы основана на использовании формулы Дюпюи. На прямой линии выбирается

произвольная точка с координатами Qи P,

по которым опреде-

ляется коэффициент продуктивности Kпрод = Q/

P.

По коэффициенту продуктивности в соответствии с формулой (2.43) можно вычислить значение коэффициента проницаемости.

 

Kпрод μ ln

rк

 

 

 

k =

rc

.

(2.63)

 

 

 

 

2πh

 

 

 

 

 

 

63

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

 

 

 

P

= A + B Q.

(2.64)

 

 

 

Q

 

 

 

 

 

Индикаторная диаграмма (рис. 2.20) перестраивается в коорди-

натах

P

,Q

. Прямая линия отсекает на оси ординат отрезок А:

Q

 

 

 

 

 

 

 

 

 

 

 

A =

μ

ln

rк

.

(2.65)

 

 

 

2πkh

 

 

 

 

 

 

r

 

 

 

 

 

 

 

c

 

По углу наклона прямой определяют коэффициент B = tga.

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

 

сти Q Р это условие не вы-

 

полняется.

Для

обоснованного

 

определения

коэффициента

 

продуктивности

необходимо

 

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

 

фильтрации жидкости, для чего

 

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

 

мероприятие, как исследование

Рис. 2.20. Обработка индикаторной

скважины.

Разделив фактиче-

диаграммы по двухчленной формуле

ский дебит

скважины на де-

 

прессию, можно получить вели-

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

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

64

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

Также существуют некоторые разногласия при оценке части пласта, характеризуемой определенным при обработке индикаторной диаграммы коэффициентом продуктивности. Некоторые специалисты считают, что этот коэффициент характеризует так называемую околоскважинную зону пласта (ОЗП), основываясь на том, что основной перепад давления приходится как раз на эту зону. Данное утверждение теоретически не обоснованно, так как в формулах, по которым и определяется коэффициент проницаемости (2.63), четко указаны границы, для которых он определяется: rc rк. Таким образом, теоретически более обоснованно считать, что коэффициент проницаемости, полученный при обработке индикаторных диаграмм, характеризует в среднем всю зону дренирования пласта скважиной.

2.6. Установившееся безнапорное движение жидкости

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

Рис. 2.21. Безнапорное движение жидкости

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

65

ные потоки имеют переменную высоту живого сечения z; высота живого сечения максимальна на контуре питания и соответствует напору на контуре питания hк; минимальный напор соответствует

живому сечению в области стока hг ,hс. Положение уровня свобод-

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

разница между

этими уровнями называется понижением уровня

S = S0 Sд. На

практике безнапорные потоки встречаются при

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

Дренажный грунтовый поток. Движение жидкости осуществляется под действием перепада напоров между контуром питания hк

и прямолинейной галереей hг , протяженность потока есть расстояние между контуром питания и галереей Lк, ширина потока a. Площадь живого сечения потока является величиной переменной в свя-

зи с переменной отметкой свободной поверхности z: При линейном законе фильтрации

ν =

Q

= k γ

 

dz

 

,

F (x)

 

 

μ

 

dx

 

 

h

 

 

Q μ

0 dx.

 

 

к zdz = −

 

 

a k γ

 

 

h

 

 

L

 

 

г

 

 

 

 

к

 

 

Расход жидкости в потоке

Q = a k γ hк2 hг2 . 2μ Lк

F ( x) = a x.

(2.66)

(2.67)

(2.68)

66

Распределение напора и форма свободной поверхности жидкости определяется из следующего уравнения:

z2 = hк2 hк2 hг2 x. Lк

При нелинейном законе фильтрации:

 

 

 

Q

 

 

 

 

 

dz

 

1

 

 

 

 

 

 

 

 

 

n

 

ν =

 

 

 

 

= a z

c

 

 

x

,

 

F

(x)

d

 

 

 

 

 

 

 

 

 

 

hк

 

n

 

 

 

Q

n

 

0

 

 

 

z

 

dz = −

 

 

 

 

 

dx,

 

 

 

 

 

 

 

h

 

 

 

 

a c

L

 

 

 

г

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

hn+1

hn+1

n

,

 

Q = a c

к

 

 

г

 

 

 

 

 

 

 

 

 

 

(n + 1) L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zn+1

= hn+1

hn+1 hn+1

x.

 

к

 

г

 

 

 

 

 

 

 

к

 

 

Lк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.69)

(2.70)

(2.71)

(2.72)

(2.73)

Приток жидкости к грунтовому колодцу. Движение жидко-

сти к грунтовому колодцу происходит под действием перепада напоров на контуре питания hк и на забое колодца hс. Напор на забое

скважины соответствует высоте столба жидкости в колодце. Площадь живого сечения: F (r ) = 2π r z.

При линейном законе фильтрации:

ν =

Q

 

= k γ

dz

 

,

F (r )

 

μ

 

dr

 

 

h

 

 

Q μ

 

r

dr .

 

кzdz =

 

 

к

 

 

 

 

 

h

 

2π k γ r

r

 

 

c

 

 

 

 

c

 

 

 

(2.74)

(2.75)

67

ln rк rc

Дебит грунтового колодца

Q = π k h hк2 hc2 .

μ

Уравнение свободной поверхности

 

 

z2 = h2

 

h2

h2

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

к

 

 

rк

c ln

 

 

 

к

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

ln

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При нелинейном законе фильтрации

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

dz

 

1

 

 

 

 

ν =

 

 

 

=

 

 

 

 

 

 

 

 

= c

n

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(r)

 

2π r z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hкn+1

 

hcn+1

 

 

 

 

 

 

Q = 2π c

n

1

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 1

 

1

 

n1

 

 

1

 

 

n1

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rc

 

 

 

 

 

rк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zn+1 = hn+1

 

 

 

 

 

hn+1 hn+1

 

 

 

1

 

1

 

.

 

 

 

 

 

к

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

 

 

 

n1

n1

к

 

 

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

rк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.76)

(2.77)

(2.78)

(2.79)

(2.80)

Безнапорное движение жидкости к прямолинейной галерее.

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

ν = kγ

 

dz

,

(2.81)

μ

 

dx

 

 

68

Q = ν F = a z ν = a z kγ

dz

,

(2.82)

 

 

μ

dx

 

 

zdz =

Q μ

dx.

 

 

(2.83)

a k γ

 

 

 

 

 

 

 

Рис. 2.22. Безнапорное движение жидкости к прямолинейной галерее

Для определения дебита проинтегрируем это уравнение по z в пределах от hк до hc, и по x в пределах от 0 до L.

Q =

a k h

(hк2

hг2 )

.

(2.84)

2μL

 

 

 

 

 

Безнапорное движение жидкости к скважине. Допустим,

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

(рис. 2.23).

69

 

 

ст

 

 

H

d

H

 

 

 

 

 

 

 

Рис. 2.23. Безнапорное движение жидкости к скважине

 

ν = kγ

dz

,

 

 

 

 

 

 

 

 

(2.85)

 

 

 

 

 

μ

dr

 

 

 

 

 

 

 

 

 

 

 

Q = 2π r z ν = 2π r z

kγ

 

dz

,

(2.86)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ

 

dr

 

 

zdz =

 

Q μ

 

dr

,

 

 

 

 

 

(2.87)

2π k γ

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

h

 

 

 

 

Q μ

 

 

R

dr .

 

 

 

 

к zdz =

 

 

 

к

 

 

 

(2.88)

 

 

2π k γ

 

 

 

z

 

 

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2 = h2

=

Q μ

ln

Rк

.

 

 

 

(2.89)

 

 

 

 

 

 

 

 

 

к

 

 

 

π k γ

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в эту формулу

z = hc

 

при r = rc , получим следую-

щую формулу дебита скважины:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q =

π k γ (hк2 hc2 )

.

 

 

 

 

(2.90)

 

 

 

 

μ ln

 

Rк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

70