08-RUS Permeability PDF
08-RUS Permeability PDF
добычи
NExT Training
E. Standen
Оценка добычи
• Дебиты рассчитываются по закону Дарси
– Чтобы упросить ввод исходных данных, для
большинства уравнений принимаются, допущения
• По каждому слою, который мы предполагаем
однородным на некотором удалении от ствола
скважины, нам нужен замер проницаемости
– Прямой замер по керну
– Оценка по данным ГИС
– Замер подвижности флюидов по данным ГДИ
• Пористость – это емкость пласта, которая
может иметь или не иметь отношения к
проводимости
2
Сопоставление керна
и каротажной диаграммы
3
Замеры пористости песчаника
различными методами
Общая пористость – нейтронный каротаж
Скелет
Пористость по керну, высушенному в печи
породы
d
n
a
S P
e
n
to
s o
yrsit s
a
e
M d
re
u
b rio
a
Vy c
e
sT
u iq
sn
e
u
h
T
T o
rsiotyalPrsi
lP
yta tr-tyLogDensiNur-
u
siN
n
e
D
g
o
yL
xxR
triRockMatri
a
ckM
o riOven-Dri
-D
n
ve
O a
n
A
re
o
C
d
e
a
n
rA
o
C
d lu
b
A
rsite
o
yP
u
ts
l
b
A
rie
o
P
y rsrTotalPs
lP
ta
o
T tiyity
rtzQuartz
a
u
Q la
C
ly
V
aV
C
y seShalaySurfces
rfce
S lyu
a
h lL
a
m
S
C
H
re
o
H
re
L
lS
a
mo
C itP
m
y-D
ld
a ,Icrge,IntcasydiPol-DmuA
rsn
o
e
A
u g tyisn
cd
e
o
osy
td
c
e
ni d
te
la d
te
Isola
Iso
(
( k)F
Fo
w
e
m
ra
k)
o
w
e
m
ra rs&
ye
a
L
&
rs
e
y
a
L tla
rsIn
ye
a
tl
In
rs
e
y syC
sae
o
rP
il
p
a
ly
re
o
Pi
p
C re
o
P
re
o
P ss sPores
re
o
P
t(O
a
)W
a
c c re
ru
tru
S
HS
H
(O
W
-)llB
H
H
B
e
r iyyoundW
ro
te
a
W
d
n
u ri
te
a rW
te
a
W
r
te
a H
ca
yd
V
re
o
P
H
a
c
d
y
V
re
o
P n
o
b
e
m
u
lrumebonlr
Пористость по керну,
Vглин Высушенному с охранением связанной воды
Структурная вода
вода- (OH Гидратная или Эффективное
связанная вода поровое пространство
Irreducible or
Immobile Water
(modified from Eslinger and Pevear, 1988)
4
Насыщенность пор горных пород
Pc
5
Sw
Зависимость проницаемости и
пористости
permeability vs porosity
10000
1000
Permiability
100
10
0.1
0 0.1 0.2 0.3 0.4
porosity
Повышение
зернистости
и
ст
но
тн ие
ен
ло н
у п че
и/ и
ци ел
Увеличение
та Ув
видимой
ен
пористости
Проницаемость
ем
Ц
ы а
овк
ли н
тир
сор
ая
ой г
нн
зер sing
е
учш
Ул
нов
Меж Increa
Увеличение микропористости
Пористость
(Modified from Evans and Others, 1997) 7
• Напредыдущем графике показано
влияние различных факторов на
поровые системы обломочных
осадочных пород.
• Стрелками обозначены
направления лишь качественных
изменений, но не количественных
показателей. Значение имеет не
место расположения стрелок, а
лишь их направление.
8
Проницаемость по данным ГИС
Первая попытка определить проницаемость по каротажным
данным была предпринята Вилли и Роузом в 1950 году,
когда они в своей работе привели следующее уравнение :
где:
C и C’ – константы,
a ≈ 2.
(C действительно связана с пористостью)
10
Дополнительные попытки определения проницаемости по данным ГИС
• Tixier;
• K^1/2 = 250 Φ^3 / Swi = 1,139 md
• Timur;
• K^1/2 = 100 Φ^2.25 / Swi = 1,109 md
• Коутс – Дюмануар;
• K^1/2 = 300 Φ^w / w^4 Swi^w = 1,779 md
• Коутс (с изменениями)
• K^1/2 = 70 Φ^2 (1 – Swi) / Swi = 635 md
13
Проницаемость по Сарабанду
15
Everyone assumes this is the Wyllie and Rose
permeability chart, however it is the Timure
permeability chart.
The chart book says the formula presented is the
empirical relation for one study the chart actually
correlates to a relation of K = (100 Φ^2.25 / Swi)^2.
In the 1969 chart book the formula for an approach
to Wyllie & Rose permeability was K^1/2 = 250
Φ^3/Swirr for light oil and K^1/2 = 79 Φ^3/Swirr
for dry gas which is the Tixier equation.
16
Палетка проницаемости K4
– 2000 Schlumberger Chart book (Коутс с изменениями)
K= {70* Φe**2[(1-Swi)/Swi]}**2
17
Оценка дебита
• Зависит от ухудшения коллекторских свойств
в процессе бурения (положительный скин).
• Допущения к уравнению радиального потока
Дарси.
• В большинстве случаев предполагается
абсолютно свободное фонтанирование.
18
Хорошая
газовая
скважина?
В Западной Канаде
неглубокие газовые
скважины добывают из
пласта Милки Ривер по 100
тысяч кубофутов в сутки
при перепаде давления в 300
psi.
20
Упрощенная форма для
негидродинамиков
Qoil = k h P / 1000 µ
Примечание: это дебиты при абсолютно
свободном фонтанировании без учета скин-
эффекта и многофазного потока. Поэтому эти
оценки будут выше, чем результаты
опробований скважин.
21
Φ = 12%
Φ= 6 to 15%
Φ = 19%
Φ = 19%
Sw=100% Φ = 18%
Sw=100% Φ = 19%
22
Зависимость общего объема воды
BVW от капиллярного давления
Swirr x Porosity = BVW
(Остаточная водонасыщенность
умноженная на пористость дает
общий объем воды) в условиях
остаточной насыщенности. Это
значит, что когда BVW
приближается к низкому
Давление постоянному значению для пласта,
то при значениях выше этой точки
или пласт будет давать безводную
глубина продукцию.
Выше точки Swirr изменения в
будут отражать изменения в
размерах пор (зернистости) или
же изменения в содержании
углеводородных флюидов.
Помните, что у каждого
литологического комплекса своя
0 Swirr Swirr 100 остаточная водонасыщенность
Sw Swirr.
23
Низкий BVW Высокий BVW
Зависимость общего объема воды от проницаемости – палетка K4
“График кривых”
K= {70* Φe**2[(1-Swi)/Swi]}**2
Литотип 1
Литотип 2
Водонас.
зона
Переходная зона
24
Пример по м. Базальный Кварц
• Нефть
• Q = .00708 k h (Pe – Pw)/ µ ln(Re/Rw)
• Простой q = k h P / 1000 µ
• Q = 1,200 барр/сут (высокий при наличии газа)
25
Пористость и проницаемость по данным ЯМК
Многоэкспонентное затухание
1.0
Амплитуда эха
0.6
Поры трех размеров
3 sizes of pores 0.4
80
0.2 20
0.0
T = 20 ms T = 5 ms
2 2
-0.2
30% порового объема
0 20 40 60 80 100
t (ms) 950505-01c
T = 5 ms
2 Распределение T2
30% порового объема 0.5
0.4
T = 80 ms 960909-01
2
0.3
40% порового объема
0.2
0.1
0.0
1 10 100 1000
T (ms) 950505-02a
2
26
It must be understood that while effective rate of decay of any given fluid
component in any given pore is a sum of the rates of decay due to surface, bulk
and diffusion effects, the signal we measure is a further sum of all the components
within the pore space. Hence, the decay we measure is no longer a simple
exponential decay. It becomes a multi-exponential decay. To explain this, let us
take a simple case of a triple porosity system in a rock.
Let’s say that 30% of the porosity exists in very fine pores having an effective T2
of 5ms; a further 30% of the porosity is in larger pores with a T2 of 20ms; and the
remaining 40% of the pores are very large with a T2 of 80ms.
The contribution to the echo train from each component is shown here on the top
right. The small pores start at 0.3 representing 30% of the porosity and decay
rapidly. The medium pores also start at 0.3 and decay slightly slower, while the
large pores start at 0.4 (representing 40%) and decay very slowly. The net rate of
decay that we measure is shown in black dots at the top and we can see that we
start at 1.0 (representing all of the porosity) and decay in a multi-exponential
manner.
27
The NMR inversion software attempts to extract the 3
individual rates of decay from the measured echo train by
solving for a given number of decay components distributed
logarithmically from T2 Minimum to T2 Maximum. In this
example, we have used 50 components from 1ms to 1000ms.
On the downhole log we typically solve for 30 components
from 0.3ms to 3000ms.
The inversion appears to have correctly extracted the 3
components of our example with the peaks located at the
appropriate T2 values and having the correct amplitudes. We
see now how the measurement can give us the porosity
distribution in terms of T2. If surface effects predominate,
then the T2 distribution can be interpreted as the pore size
distribution.
28
Зависимость между T2 и размером пор
(V/S) pore (µm)
-3 -2 -1 0 1
Форма пор V/S
10 10 10 10 10
Диаметр сферы = d d/6
Диаметр трубки = d d/4
Berea 100 Sandstone
ρ = 5µm/s Ширина покрова = d d/2
Population
29
• Provided surface effects predominate, we can transform the
T2 distribution into pore size distribution through a
knowledge of the surface relaxivity of the rock.
• Regular geometrical shapes have a fixed relation between
the volume to surface ratio and their size. Can you derive
the relationship for a sphere?
30
Зависимость между свободной жидкостью по
данным ЯМК и пористостью по методу
центрифугирования
20
Well A
Well B
15
φff (3 3) 10
0
0 5 10 15 20
Φ c e n trifu g e
31
NMR Free Fluid Porosity vs Centrifuge Porosity
Free fluid porosities calculated from NMR laboratory spectrometer measurements
are compared here with centrifuge porosity measurements on sandstone core
samples from two wells.
Free fluid porosity is calculated using a 33-msec cutoff on T2 distribution
(above 33 msec is considered to be free fluid porosity, and below 33 msec is
bound fluid porosity). The centrifuge porosity is determined by spinning the
cores at 100 psi and measuring the amount of fluid extracted.
The excellent correlation between the calculated free fluid porosity and the
centrifuge porosity demonstrates that a good free fluid porosity can be
determined by analyzing the T2 distribution or the pore size distribution of the
rock.
32
Увязка ЯМК с проницаемостью
• Метод 1: Формула Козени-Кеньона
(проницаемость SDR) a = prefactor, ~ 4 для песчаников
K SDR = a ⋅ φ ⋅ (T2LM )
b c b = пористость как экспонента, ~ 4
c = экспонента T2 , ~ 2
T2LM = логарифмическое среднее T2
35
Граничные значения для уравнений проницаемости
K SDR = a ⋅ φ ⋅ (T2LM )
b c
Pore diameter (microns)
0.01 0.1 1 10
T2 original
T2 spun sample
1 10 100 1000
Time (T2 msec)
c′
b′ FFV
Clay- Capillary-
′ 4
K Timur = a ⋅ 10 ⋅ φ ⋅
bound bo un d
Producible
fluids
BFV
water water
φ = FFV + BFV
Total CMR porosity
Эмпирически граничное
Free fluid porosity
значение (φff) для песчаников
CMR 3ms porosity
было определено равным 33
мсек. Если его изменить, то
изменится соотношение36
FFV/BFV для формулы Тимура.
Зависимость между проницаемостью по данным ЯМК
и проницаемостью керна для модели пластовой воды
10 0
W e ll A , C = 2 .8
W e ll B , C = 3 .4
k b rin e
k N MR = C ( φ N MR ) 4 (T2 ,lo g ) 2
0 .01
0 .01 10 0
k N MR
37
• NMR Permeability vs Core Brine Permeability
• Permeabilities calculated from NMR laboratory
spectrometer measurements are compared with
conventional core brine permeabilities of sandstone core
samples from two wells.
• NMR permeabilities are calculated from NMR
porosity and the logarithmic mean of T2 relaxation times.
• The excellent correlation between the conventional
brine permeability measurements and the NMR-derived
permeabilities shows that good permeability values can be
computed from NMR measurements.
38
10 3
10 2 ++
++++
++++ ++
+
Проницаемость
++++
++ +
++ ++
+ + + ++ +++ ++
+ ++
по данным ЯМК
1 ++++++++
k (m d ) 10 +++++ ++ +
+ + + ++ +
+++ ++ При наличии хорошей
+ + зависимости между размерами
0
+ + ++
10
пор и размерами поровых
+
+ связок мы получим хорошее
10 –1 + соотношение между
10 –1
10 0
101 102 103 проницаемостью по данным
ke s tim a te (m d ) ЯМК и проницаемостью по
3
керну. Обычно такое случается
10
в хорошо отсортированных
+ обломочных породах.
2
+ +++++
10 ++
+ + +++ +++
+
++
+ ++ + ++++ + Если же зависимости между
+ ++++ + +
+ +++ +++ + +++ + размерами пор и размерами
+ + +
k (m d ) 10 1
+
+
+
+++ ++++ ++ + поровых связок не наблюдается
+ + + + + +++
+ (например, в карбонатах), тогда
+ ++ + проницаемость, выведенная по
0
+
10
данным ГИС, будет неверной,
+ как бы мы не подгоняли
+
10 –1
+ пористость, потому что
5 10 15 20 25 пористость в уравнении
φ (p . u . ) возводится в четвертую степень.
39
• NMR Permeability and Crossplot Porosity vs Conventional Core Brine
Permeability
• Here, we compare two methods of obtaining continuous permeability
information in a well.
• The traditional method derives a transform between porosity and
permeability from core brine measurements. The lower plot shows
how scattered a simple permeability/porosity plot can be, and it is clear
that a simple transform cannot lead to quantitatively reliable
permeabilities.
• Permeabilities derived from NMR measurements incorporating
both porosity and pore size information (upper plot) are much more
coherent with the core brine measurements.
• The CMR tool can give a continuous NMR measurement of the
formation, allowing continuous permeability calculations. This has two
clear advantages over the traditional porosity/permeability plot
approach:
• •
40
• The permeability values derived are much more
precise because pore size information is included.
• • Permeability measurements are not dependent
on extensive core sampling.
• This slide also gives us a feel for the real
improvement in permeability estimation delivered
by the CMR tool. The poro-perm method typically
has an uncertainty of nearly 2 orders of
magnitude. The NMR estimate, on the other hand,
has an uncertainty of just about 1/2 an order of
magnitude.
41
Пористость и проницаемость на скважине
Примечание: в данном примере пористость по ЯМК ниже, чем по керну, из-за недостаточного времени
ожидания углеводородов, а проницаемость по Тимуру ближе к проницаемости по керну, потому что по
42
результатам анализа керновых данных граничное значение Т2 в данных ЯМК были смещены.
В примере от Австралазии мы строим график зависимости
пористости и проницаемости, определенных на скважине, от
пористости и проницаемости по керновым данным.
43
This zone is known to contain very light oil of downhole viscosity of 0.1Cp. The
log was acquired with a WT time of 6s, while the oil has a T1 of 13.3s. Hence, we
are underpolarized for the oil signal. This explains the lower TCMR porosity
compared to the core porosity. In the two shale44 streaks at x050’ and x140’, the
TCMR matches perfectly with core porosity. Inspite of the mismatch in the
absolute porosity in the sands, the porosity profile is faithfully reproduced.
We can also see that the estimated permeabilities underestimate the core
permeability. However, here again the permeability profile is very well matched.
A polarization corrected porosity and calibrated permeability coefficients will
permit us to match core permeability in absolute value also. Even the uncalibrated
permeability is useful in that we can still compare sand quality using the wellsite
permeability log and locate the sweetest spots for planning formation tester
surveys.
44
Классификация карбонатных пород по Данему
Осадочная структура различима Осадочная
Первичные компоненты не были связаны при отложении
Первичные компоненты структура
Содержит частицы Глинистые,
во время отложения неразличима
Глинистой, илистой и алевритовой размерности были скреплены
Частицы
отсутствуют
Скрепл. Илист. частицами Порода
скреплена Порода зернистая
Зерен >10 %
зернами
u h
D m
a
n C
n te
o
rb
a o
R
a cs
la cfia
C
k ts n
o
o
p
e
D so
tin
x
e
lT
a zre
iu
n
g
o
c
R a
e
l
b
N
ts
e
n
p
m
o
C n
e
g
T
o
iB
ta
d
u
M
s n
it
ru
h
e
D
g o
tin
sp n
g
iB
r
O
m
lC
a
T
d
u
o D
tpherns
o
N o
tin
T se
lp
a
zb
r
u
x
g
c
R
a
p
S
d
u
lM
(c
y tz
le
id
rs
a
p
o
n e
n
c
ls
ra
iG
t n
i
,G
M
s
k
c
a
L
d
rte
o
p
u
-S n
ri
u
o
p
e
D
g o
tin
s
n
i
ra
G
%
0
1
<
s n
i>
ra
G
%
0
1
s d
rte
o
p
u
S C ry
e
lin
ta
s
M u
a
W
e
n
to
s
d c
P
n
to
s
e
k G
e
n
to
s
k
c
a e
to
s
in
a
r B
s
d
n
u
o tC
e
n
o a
te
n
o
rb
<10 % Зернистый
Grains Кристаллический
известняк
Мадстоун Вакстоун Пакстоун Грейнстоун Баундстоун Карбонат
По-видимому, диапазон
зернистости от 100 до 20
Rock type 1
100µm микрон разделяет карбонаты
2 категории пород – с
Rock type 2 высокой и низкой
20µm
проницаемостью. Как ни
Rock type 3 странно, зависимость между
проницаемостью и
пористостью легче
Rock type 4 спрогнозировать для пород с
меньшей зернистостью.