Диагностика газовой скважины по результатам гидродинамических исследований при установившейся фильтрации
.pdfСПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Уравнение притока (12) с экспериментально определен ными коэффициентами широко используется в расчетах при проектиро вании разработки месторождений. Кроме того, по значению А, найденному в результате исследования скважины, можно определить коллекторские свойства пласта, например коэффициент гидропроводности:
kh |
|
pатм |
ln RК |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(14) |
|
|
|
|
|
|
|
|
|
||||
|
|
A |
rс |
|
|
|
|
|
|
|
|
|
||||
Уравнение притока реального газа к скважине по двучленному закону |
||||||||||||||||
фильтрации имеет вид |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
RК |
|
|
|
|
|
|
|
|
|
|
p2 |
p2 |
|
zpатм |
ln |
Q |
атм |
|
атм zpатм |
Q2 |
|
||||||
|
kh |
r |
|
|
||||||||||||
к |
|
с |
|
|
|
2 |
2 |
r |
k |
атм |
|
|||||
|
|
|
|
|
|
|
C |
|
|
|
2 |
h |
|
(15) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
к ) / 2 |
|
|
|
||||
где |
z (zc zк ) / 2 |
; |
( c |
и являются константами. |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Отметим, что в реальных условиях нельзя считать, что во всем пласте - от стенки скважины до контура питаниясправедлив единый нелинейный закон фильтрации.
1.3 Приток газа к несовершенной скважине
Виды несовершенства скважин.
Скважина называется гидродинамически совершенной, если она вскрывает продуктивный пласт на всю толщину и забой скважины открытый, т. е. вся вскрытая поверхность забоя является фильтрующей.
Если скважина с открытым забоем вскрывает пласт не на всю толщину h, а только на некоторую глубину b, то ее называют гидродинамически
несовершенной по степени вскрытия пласта. При этом h b / h называется относительным вскрытием пласта.
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Если скважина вскрывает пласт до подошвы, но сообщение с пластом происходит только через специальные отверстия в обсадной колонне и цементном камне или через специальные фильтры, то такую скважину называют гидродинамически несовершенной по характеру вскрытия пласта.
Нередко встречаются скважины с двойным видом несовершенства-как по степени, так и по характеру вскрытия пласта.
Степень и характер вскрытия пласта имеют важное значение при разработке месторождений нефти и газа, так как они определяют фильтрационные сопротивления, возникающие в призабойной зоне, и, в конечном итоге, производительность скважин. Выбор степени и характера вскрытия осуществляется в зависимости от физических свойств пластов, их толщины, степени неоднородности, способа разра ботки и т. д. Несовершенство скважин по степени и характеру вскрытия приводит к таким деформациям линий тока, которые приводят к возникновению в призабойной зоне сложных неодномерных течений. В связи с этим рассмотрение особенностей притока к гидродинамически несовершенным скважинам имеет большое практическое значение.
Приток газа к несовершенным скважинам при двучленном законе фильтрации.
Несовершенство газовых скважин при выполнении закона Дарси
Q k (Pк Pс )
m r (16)
учитывается так же, как несовершенство нефтяных скважин, т. е. радиус скважины в формуле дебита заменяется приведенным радиусом:
rc rce (c1 c2 ) (17)
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Для расчета дебитов газовых скважин несовершенных по степени и по характеру вскрытия при нарушении закона Дарси может быть предложена следующая схема. Круговой пласт, в центре которого находится скважина, делится на три области (рис. 6).
Рис.6 - Схема притока газа к несовершенной по степени и характеру вскрытия скважине
Первая область имеет радиус R1=(2-3)rс, здесь из-за больших скоростей вблизи перфораци онных отверстий происходит нарушение закона Дарси, т. е. в основном проявляется несовершенство по характеру вскрытия. Линии тока пока заны на рис. 9.
Вторая область представляет собой кольцевое прост ранство R1< r < R2, R2≈h; здесь линии тока искривляются из-за несовершенства скважины по степени вскрытия, имеет место двучленный закон фильтрации.
В третьей области R1< r < Rк, действует закон Дарси, течение можно считать плоскорадиальным. Обозначив давления на границах областей через р1 и р2, запишем для третьей области в соответствии с формулами нахождения дебита скважины для плоскорадиальной фильтрации:
Qm 2 h (Pк Pс )
ln(Rк / rc ) (18)
Подставив (18) в (19), получим:
|
|
|
|
pк |
pdp |
|
|
|
2 hk атм |
|
|||
Qm |
|
|
|
|
|
|
p |
|
ln(R / r ) |
( p)z( p) |
(19) |
||
|
|
атм |
к c p |
|
||
|
|
|
|
c |
|
|
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Перейдём к дебиту, приведённому к атмосферному давлению:
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
|
|
2 hk |
|
pк |
pdp |
|
|
|
|
||||
Qатм |
|
|
|
|
|
|
|
|
|||||
p |
ln(R |
/ r ) |
( p)z( p) |
(20) |
|
|
|||||||
|
атм |
|
к |
c p |
|
|
|
|
|
|
|||
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
подстам в (20) |
|
получим: |
|
||||||||||
z и |
|
|
|
||||||||||
|
Qатм |
|
2 hk |
|
|
|
pк |
|
hk( p2 p2 ) |
|
|||
|
|
|
|
|
|
pdp |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
к с |
|
|
|
|
pатм z ln(Rк / rc ) p |
|
pатм z ln(Rк / rc ) |
(21) |
|||||||
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из уравнения (21) получим течение газа в третьей области
|
|
|
Rк |
|
|
( pк2 p22 ) |
Qатм pатм z |
ln |
|
||
|
R2 |
||||
|
hk |
|
(22) |
Во второй области примем, что толщина пласта переменна и изменяется по линейному закону от значения b при r = R1 до значения h при r = R2, т. е.
z(r) = α+βr (23)
где αиβопределяются из условий z=b при r=R1,z = h при r = R2. Чтобы получить закон движения в этой области,надо проинтегри ровать уравнение (2), предварительно подставив вместо постоянной толщины h переменную толщину по формуле (23).
|
|
|
|
|
R2 |
|
|
|
|
2 |
|
1 |
|
1 |
|
|
|
p2 |
p2 |
|
Qатм pатм z |
ln |
C |
|
атм pатм zQатм |
|
C' |
||||||||
|
|
|
|
|
|
|
|
|
|||||||||
2 |
1 |
|
hk |
R |
1 |
|
2 2h2 k |
|
R R |
1 |
|||||||
|
|
|
|
|
|
|
|
|
|
(24) |
|||||||
|
|
|
|
|
1 |
|
|
|
|
1 |
|
2 |
|
Здесь C1 и С1’ -коэффициенты, характеризующие несовершенство скважины по степени вскрытия.
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
h ln |
|
|
h |
|
|
|
|
|
|||||||||||||
C1 |
ln h |
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
(25) |
||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
h |
|
|
h |
|
|
|
|
R1 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
C1' |
|
|
1 |
1 |
|
|
1 |
|
|
|
|
|
b |
||||||||||
|
|
|
|
|
h |
||||||||||||||||||||
|
|
|
|
2 |
|
R1 |
, |
||||||||||||||||||
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
h |
|
|
|
|
|
|
|
|
h (26) |
Обе последние формулы-приближенные, они имеют место при значениях b » R1.
В первой области фильтрация происходит по двучленному закону, плоскорадиальное течение нарушается из-за перфорационных отверстий; несовершенство по характеру вскрытия учитывается коэффициентами С2 и C1:
|
|
|
|
|
R1 |
|
|
|
|
|
2 |
|
1 |
|
1 |
|
|
|
p2 |
p2 |
|
Qатм pатм z |
ln |
C |
|
|
атм pатм zQатм |
|
C' |
||||||||
|
|
|
2 |
|
|
|
|
|
|
|||||||||
1 |
2 |
|
hk |
r |
|
2 2h2 k |
|
r R |
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
(27) |
||||||||
|
|
|
|
|
c |
|
|
|
|
c |
1 |
|
Здесь С2 определяется по графикам В. И. Щурова, для С2’ предла гается приближенная формула
C2' h2
3N 2l'3 (28)
где N- суммарное число перфорационных отверстий; ℓ’-глубина проникновения перфорационной пули в пласт.
Складывая почленно уравнения (22), (24)и(27) и пренебрегая величиной 1/R2, получим уравнение притока газа к несовершенной скважине в виде
2 |
2 |
|
|
|
|
Rк |
|
|
|
|
|
2 |
|
' |
' |
|
|
|
Qатм pатм z |
|
|
атм pатм zQатм |
1 |
|
|||||||||||
pк |
pс |
|
|
|
ln |
|
|
C1 |
C2 |
|
|
|
rcС1 |
rcC2 |
|||
hk |
|
r |
2 2h2 kr |
|
|||||||||||||
|
|
|
|
|
|
|
c |
|
|
|
|
c |
|
|
|
|
(29) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Если записать уравнение (29) через коэффициенты фильтрационных сопротивлений А и B в виде (12), то для несовершенной скважины получим:
|
|
|
Rк |
|
|
|
|
|
|
' |
' |
|
|
|
pатм z |
|
|
атм pатм z |
1 |
|
|||||||
A |
|
ln |
|
|
C1 |
C2 |
; B |
|
rcС1 |
rcC2 |
|||
hk |
r |
2 2h2 kr |
|||||||||||
|
|
|
|
c |
|
|
|
|
c |
|
|
|
(30) |
|
|
|
|
|
|
|
|
|
|
|
где C1 и C1’ определяются по формулам (25) и (26), С'2-по формуле (28), а С2-по графикам В. И. Щурова(рис.10).
Рис. 10 - Графики В. И. Щурова для определения коэффициента С2 при ℓ= 0,5.
Номерам кривых соответствуют значения α: 1 -_0,02; 2 - 0,04; 3 - 0,06; 4 - 0,08; 5 - 0,1; 6 - 0.1; 7 - 0,14; 8 - 0,16; 9 - 0,18; 10 - 0,2
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
2.Расчётная часть
2.1 Определение коэффициента фильтрационного сопротивления по данным исследований.
В ходе проведения исследований были установлены следующие значения для забойного давления(pзаб) и дебита скважины(Qат),
Таблица 1
pзаб , МПа |
Qат, м³/сут |
|
|
1,5 |
124000 |
|
|
1,6 |
76000 |
|
|
1,6 |
36000 |
|
|
1,66 |
14000 |
|
|
Взяв за основу эти pзаб и переведя Q из куб. метров в сутки в куб.метры в секунду(таб.2), а также зная тот факт, что при Q=0 pзаб=pпл , то есть при дебите скважины равном 0 забойное давление равно пластовому, можем найти пластовое давление, построив график зависимости между забойным давлением и дебитом скважины(рис.11)
Таблица 2. Зависимость между пластовым давлением и дебитом скважины.
pзаб, МПа |
Qат, м³/с |
|
|
1,5 |
1,435185185 |
|
|
1,6 |
0,87962963 |
|
|
1,6 |
0,416666667 |
|
|
1,66 |
0,162037037 |
|
|
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Рис.11 - Зависимость между пластовым давлением и дебитом скважины
Из этого рисунка видно, что при исследованиях была допущена ошибка в измерении pзаб. и соответствующего ему дебита, а именно при pзаб =1,6 МПа дебит скважины, в данном случае, не равен Qат ≠0.416666667 м3/с.
Исключая это значение и продолжая график до пересечения с осью Y, когда Qат=0 построим новый график зависимости между забойным давлением и дебитом скважины (рис.12) и найдем из него pпл.
Рис.12 - Зависимость между квадратом пластового давления и дебитом скважины
Видим, что пластовое давление равно pпл 1,73 МПа
Теперь, зная, что пластовое давление =1,73 МПа и фильтрация
происходит по двучленному закону построим график зависимости ( pпл2 pзаб2 )/Qат от Qат для фильтрации газа(рис.13), взяв значения из Таблицы 3.
СПБГУАП группа 4736 https://new.guap.ru/i03/contacts
Таблица 3.
2 |
2 |
Qат , м3/с |
( pк |
pс )/Qат , МПа2*с/м3 |
|
0,515679 |
1,435185185 |
|
|
|
|
0,488636 |
0,87962963 |
|
|
|
|
Рис. 13 - График зависимости ( pпл2 pзаб2 )/Qат от Qат при фильтрации газа по двучленному закону
А и В – коэффициенты фильтрационного сопротивления. Коэффициент А находим, как расстояние между осью абсцисс и точкой
пересечения прямой с осью ординат, а коэффициент B, как тангенс угла наклона прямой к оси абсцисс, то есть B=tgβ.
Из теоремы о нахождении тангенса угла в прямоугольном треугольнике знаем, что он равен отношению противолежащего катета к
(0,515 |
0, 44) |
МПа2 с |
|
|
МПа2 |
с2 |
|||
|
м3 |
0,05 |
|||||||
прилежащему поэтому B tg |
|
|
|
м6 |
|
||||
|
1, 44м3 / с |
|
|
||||||
|
|
|
A=0,44 |
МПа2 с |
|
|
|
||
Коэффициент A в свою очередь равен: |
м3 . |
|
|
||||||
|
|
|
|