Мангазеев П.В., Панков М.В., Кулагина Т.Е., Камартдинов М.Р. Гидродинамические исследования эксплуатационных и нагнетательных скважин. 2003 год
Подождите немного. Документ загружается.
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
375
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
376
Глава 6
Упражнение 1
0.1
1
10
100
0.001 0.01 0.1 1 10 100
Эффек т
ВСС
Ради альн ый
приток в
бесконечном
пласте
Переходная
зона
A
B
C
D
Диагностический график
t
p и p’
◊Давление
□ Производная
давления
Глава 6
Упражнение 1
• По фактическим данным был построен билогарифмический
диагностический график, который показал типовую кривую производной
давления для скважины в бесконечном пласте. По характеристическим
признакам кривой производной давления можно идентифицировать
несколько периодов в работе скважины:
– период влияния ствола скважины (в интервале между точками А и В),
который заканчивается в момент времени t
B
= 0.03 ч
– переходный период (в интервале между точками В и С), когда влияние
околоскважиной зоны постепенно уменьшается, и перепад давлений
распространяется дальше в пласт;
– период радиального притока к скважине (в интервале между точками С
и В), характеризующий свойства пласта. Примем время начала
радиального притока t
C
= 3.5 часа. Тогда начало прямолинейного
участка КПД в полулогарифмических координатах соответствует
log t
C
= 0.54.
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
377
Глава 6
Упражнение 1
245
247
249
251
253
255
257
259
261
263
265
0.001 0.01 0.1 1 10 100
A
B
C
D
p
wf
(t=1)
График в полулогарифмических координатах
log t
p
wf
t=0
m
log
= 2.0623
атм/лог. цикл
Глава 6
Упражнение 1
• На графике в полулогарифмических координатах, начиная с момента
времени t
C
= 3.5 часа, достаточно четко выделяется прямолинейный
участок, который можно определить по двум параметрам:
– наклон m
log
= 2.0623 атм/лог. цикл
– отрезок, отсекаемый прямолинейным (экстраполированным) участком
на оси ординат при t = 1, т.е. log t = 0, p
wf
(t=1) = 248.84 атм.
• Таким образом мы можем определить параметры пласта, используя
исходные данные по скважине и параметры m
log
и p
wf
(t=1).
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
378
Глава 6
Упражнение 1
мД
hm
Bq
k
s
1.14195.21
log
==
µ
7.309232.3log1513.1
2
log
1
=
+−
−
=
=
wt
ti
rc
k
m
pp
S
φµ
245
247
249
251
253
255
257
259
261
263
265
-3 -2 -1 0 1 2
Определение параметров пласта
Глава 6
Упражнение 1
• Среднее значение проницаемости оценивается с помощью уравнения:
• Величина скин-фактора определяется по выражению:
[
]
[
]
[][]
мД
мцикллогатм
спзсутм
hm
Bq
k
s
1.14
80.0623.2
11110
195.21195.21
3
log
=
⋅
⋅⋅
==
µ
[] []
[]
[]
[]
[]
[]
7.309232.3
08.01102.212.0
1.14
log
.0623.2
84.248265
1513.1
09232.3log1513.1
224
2
log
1
=
+
⋅×⋅⋅
−
−
−
=
=
+−
−
=
−
=
матмспз
мД
цикллогатм
атматм
rc
k
m
pp
S
wt
ti
φµ
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
379
Глава 6
Упражнение 2
0
1
10
100
0.001 0.01 0.1 1 10 100
Диагностический график
∆t
p и p’
◊Давление
□ Производная
давления
Глава 6
Упражнение 2
• Производная по данным восстановления давления имеет тенденцию к
снижению на поздних периодах исследования, даже если при этом приток
является радиальным в бесконечном пласте. Это объясняется тем, что
давление в закрытой скважине в конечном счете стабилизируется на
определенном значении, т.е. производная давления стремиться к нулю.
• Очень важно не спутать этот эффект с проявлением каких-либо других
свойств пласта, например, граница постоянного давления.
• С другой стороны, КВД можно интерпретировать таким же образом, что и
КПД. В случае, если пласт можно считать бесконечно-действующим, и,
если использовать соответствующую временную переменную (время
Хорнера, эквивалентное время Агарвала или временную функцию
суперпозиции). В этом случае производная принимает традиционный вид с
горизонтальным участком производной, характеризующим радиальное
течение в бесконечном пласте.
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
380
Глава 6
Упражнение 2
0
1
10
100
0.001 0.01 0.1 1 10 100
∆t
p и p’
Диагностический график
Эффек т
ВСС
Ради альн ый
приток в
бесконечном
пласте
Переходная
зона
A
B
C
D
◊Давление
□ Производная
давления
Глава 6
Упражнение 2
• При анализе билогарифмического диагностического графика (производная
вычислена по времени Хорнера), выделили следующие периоды работы
скважины:
– период влияния ствола скважины (в интервале между точками А и В),
который заканчивается в момент времени t
B
= 0.032 ч
– переходный период (в интервале между точками В и С), когда влияние
околоскважиной зоны постепенно уменьшается, и перепад давлений
распространяется дальше в пласт;
– период радиального притока к скважине (в интервале между точками С
и В), характеризующий свойства пласта. Примем время начала
радиального притока t
C
= 6 часов. Тогда начало прямолинейного
участка КВД на графике Хорнера будет составлять
log [(t
p
+ ∆t)/∆t] = 0.95.
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
381
Глава 6
Упражнение 2
245
247
249
251
253
255
257
259
261
263
265
110100100010000100000
График Хорнера
log [(t
p
+ ∆t)/∆t]
p
ws
A
B
C
D
p
ws
(∆t=1)
m
log
= 2.0442
атм/лог. цикл
t
p
+ 1
Глава 6
Упражнение 2
• На графике Хорнера, между точками С и D(определены по
диагностическому графику) можно выделить прямолинейный участок,
характеризующийся двумя параметрами:
– наклон m
log
= 2.0442 атм/лог. цикл
– отрезок, отсекаемый прямолинейным (экстраполированным) участком
на оси ординат при ∆t = 1, т.е.log (t
p
+ 1) = 1.7, p
ws
(∆t=1) = 261.6 атм.
• Таким образом мы можем определить параметры пласта, используя
исходные данные по скважине и параметры m
log
и p
ws
(∆t=1).
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
382
Глава 6
Упражнение 2
мД
hm
Bq
k
s
3.14195.21
log
==
µ
()
8.309232.3log
1
log
1
1513.1
2
log
=
+−
+
+
−=∆
=
wtp
ppwfws
rc
k
t
t
m
tptp
S
φµ
245
247
249
251
253
255
257
259
261
263
265
012345
Определение параметров пласта
Глава 6
Упражнение 2
• По наклону m прямолинейного участка на графике Хорнера определяем
проницаемость по формуле:
• По p
ws
(∆t=1) определяем величину скин-фактора:
• Так как в данном примере предполагается, что скважина находится в
бесконечном пласте, следовательно, экстраполированное давление (при
∆t→∞) соответствует начальному пластовому давлению p
i
= 265 атм.
[
]
[][]
мД
мцикллогатм
спзсутм
hm
Bq
k
s
3.14
80.0442.2
11110
195.21195.21
3
log
=
⋅
⋅⋅
==
µ
()
()
[] []
[]
[]
[]
[]
[]
[]
[]
8.309232.3
08.01102.212.0
3.14
log
48
148
log
./04423.2
4.2456.261
1513.1
09232.3log
1
log
1
1513.1
224
2
log
=
+
⋅×⋅⋅
−
−
+
+
−
=
=
+−
+
+
−=∆
=
−
матмспз
мД
час
час
цикллогатм
атматм
rc
k
t
t
m
tptp
S
wtp
ppwfws
φµ
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
383
Глава 6
Упражнение 3
0
2
4
6
8
10
12
14
16
18
20
0.001 0.01 0.1 1 10 100
MDH график
log ∆t
∆p
MDH
∆p
MDH
(∆t=1)
m
log
= 1.7176
атм/лог.цикл
Глава 6
Упражнение 3
• Для интерпретации данных КВД с помощью MDH метода строится график
зависимости ∆p
MDH
от log ∆t. На графике, когда влияние ствола скважины
станет незначительным, можно выделить прямолинейный участок и по
параметрам:
– наклон прямолинейного участка кривой m
log
= 1.7176 атм/лог.цикл
–p
ws
(∆t=1) – p
wf
(t
p
) = ∆p
MDH
(∆t=1) = 16.278 атм
определиться параметры пласта.
Центр Профессиональной Переподготовки Специалистов Нефтегазового Дела
384
Глава 6
Упражнение 3
0
2
4
6
8
10
12
14
16
18
20
-3 -2 -1 0 1 2
мД
hm
Bq
k
s
17195.21
log
==
µ
()
5.509232.3log
1
1513.1
2
log
=
+−
=∆∆
=
wt
MDH
rc
k
m
tp
S
φµ
Определение параметров пласта
Глава 6
Упражнение 3
• По наклону m прямолинейного участка на MDH графике определяем
проницаемость по формуле:
• По ∆p
MDH
(∆t=1) определяем величину скин-фактора:
• Скин-фактор и проницаемость, определенные по MDH методу получились
выше, чем соответствующие параметры, оцененные с помощью других
методов. Причиной этого является, то что не удовлетворяется условие
применимости MDH метода t
p
>> ∆t, так как t
p
= 48 часов, а∆t= 72 часа.
[
]
[][]
мД
мцикллогатм
спзсутм
hm
Bq
k
s
17
80.7176.1
11110
195.21195.21
3
log
=
⋅
⋅⋅
==
µ
()
[]
[]
[]
[]
[]
[]
5.509232.3
08.01102.212.0
17
log
./7176.1
2782.16
1513.1
09232.3log
1
1513.1
224
2
log
=
+
⋅×⋅⋅
−
−
=
=
+−
=∆∆
=
−
матмспз
мД
цикллогатм
атм
rc
k
m
tp
S
wt
MDH
φµ