Осипов П.Ф. Гидроаэромеханика бурения и крепления скважин
Подождите немного. Документ загружается.
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
111
2
R
r
d
ur
τ
τ
ττ
=−τ
µ
τ
∫
. (14.12)
Интегрирование дает результат:
ln
R
R
r
u
τ
τ
=−
µ
τ
. (14.13)
Из (14.8) следует, что
τ
R
/
τ
= y/R и
τ
R
/
τ
r
= r/R . Учитывая, что на стенке тру-
бы
u = u
т
, а на стенке скважины u = 0,
получаем зависимость
u(y):
ln ln
ln ln
т
y
R
uu
rR
−
=
−
. (14.14)
Чтобы найти
Q
ск
, нужно просумми-
ровать расходы
dQ
ск
элементарных стру-
ек толщиной
dy на радиусе y при
скорости струй (рис. 14.6)
u:
dQ
ск
= 2
π
y
.
dy
.
u.
После интегрирования этого выра-
жения в пределах от
r до R и подста-
новки результата в (14.6) получаем:
2
2
11
2ln
1
ск
K
R
R
r
r
=−
−
.
Аппроксимация этой зависимости степенной функцией дает результат:
K
ск
=0,5α
0,38
, (14.15)
где
α
= r/R=d
н
/D.
Формула (14.15) получена в предположении, что течение вязкой жидкости ла-
минарное. Однако она вполне применима и при турбулентном режиме течения.
14.2.3. Расчет коэффициента К
ск
для случая, когда в скважине
вязкопластичная (бингамовская) жидкость
Картина движения вязкопластичной жидкости будет существенно отли-
чаться от описанной для вязкой жидкости.
Предположим, что при скорости
u
т1
(рис. 14.7) в кольцевом пространстве, за-
полненном ВПЖ, сформировалась эпюра
τ
R1
-
τ
r1
. И если при этом динамическое
напряжение жидкости равно
τ
o
, то становится очевидным, что при y > r
o
жид-
кость останется без движения (поскольку
τ
будет меньше
τ
o
), образовав “струк-
турную” оболочку из неподвижного бурового раствора на стенке скважины.
Рис. 14.6. Эпюра скоростей
спутного потока
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
112
На рис. 14.7 показаны три эпюры напряжений, соответствующие разным
скоростям движения колонны (
u
т1
<u
т2
< u
т3
). При скорости u
т2
напряжения
на стенке скважины равны
τ
o
, что соответ-
ствует случаю, когда
r = r
o
.
Воспользуемся уравнением (14.10), под-
ставив в него реологическую функцию
ϕ
(
τ
)=(
τ
–
τ
o
)/
η
. (14.16)
Интегрирование полученного выраже-
ния в пределах от
τ
до
τ
o
(до поверхности
“структурной” стенки, где скорость равна
нулю), дает результат:
ln 1
oo
r
r
u
ττ
τ
⎛⎞
=
−−
⎜⎟
ητ τ
⎝⎠
. (14.17)
На поверхности трубы u=u
т
, а
τ
=
τ
r
.
Следовательно,
ln 1
oo
r
т
rr
r
u
⎛⎞
ττ
τ
=
−−
⎜⎟
ητ τ
⎝⎠
. (14.18)
Нетрудно показать, воспользовавшись
формулой (14.7), что
τ
o
/
τ
r
=r/r
o
, а
τ
o
/
τ
=y/r
o
. Тогда уравнения (14.17) и
(14.18) примут соответственно вид:
ln 1
r
oo
r
yy
u
rr
⎛
⎞
τ
=
−−
⎜
⎟
η
⎝
⎠
, (14.19)
ln 1
r
т
oo
r
rr
u
rr
⎛
⎞
τ
=
−−
⎜
⎟
η
⎝
⎠
. (14.20)
Разделив (14.19) на (14.20), получаем:
1ln
1ln
oo
m
oo
yy
rr
uu
rr
rr
+
−
=
+
−
. (14.21)
Это уравнение справедливо для интервала изменения
y:
r
o
≥
y
≥
r.
В “структурной” оболочке (y
≥
r
о
) скорость u=0.
Расход спутного потока найдем тем же методом суммирования расходов
элементарных кольцевых потоков
dQ
ск
, как и для вязких жидкостей:
Рис. 14.7. Образование “
структурной” оболочки
на стенке скважины
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
113
2
oo
rr
ск ск
rr
QdQ ydyu==π
∫
∫
.
Интегрирование после подстановки
u по (14.21) дает результат:
Q
ск
=
ψ
[(r
o
3
- r
3
)/3 - r
2
ln(r
o
/r)/ 2 - (r
o
2
- r
2
)/4 ] , (14.22)
где
2
ln 1
т
оо
u
rr
rr
π
ψ=
−
−
.
Чтобы найти искомую величину K
ск
, обратимся к формуле (14.6). В резуль-
тате получим:
(
)
(
)
()
33 22
2
0
22
2ln
32 4
ln 1
o о
о
ск
оо
rr rr
r
r
rr
K
rr
Rr
rr
⎡⎤
−−
⎢−−⎥
⎢⎥
⎣⎦
=
⎛⎞
−− −
⎜⎟
⎝⎠
. (14.23)
Для того чтобы воспользоваться этой формулой, нужно прежде определить
величину
r
o
. Такую возможность предоставляет уравнение (14.18). Если задать-
ся
u
т
, то величину
τ
r
можно найти численным методом, а затем вычислить r
o
по формуле:
r
o
=
τ
r
r/
τ
o
. (14.24)
В общем случае расчетный внутренний радиус “структурной” оболочки
r
o
может оказаться как меньше, так и больше
R. Описанная выше методика,
формулы (14.18), (14.21), (14.22) и (14.23), справедлива только для этого
случая, когда
r
o
<R.
На рис. 14.8 приведен сравнительный пример распределения скоростей ув-
лекаемого трубами потока в заколонном пространстве при использовании раз-
личных жидкостей при одной и той же скорости движения колонны труб. Вид-
но, что количество жидкости, вовлекаемой трубами в движение существенно
меньше, когда в скважине вязкопластичная жидкость.
С увеличением скорости движения колонны или
с уменьшением
τ
o
толщина
неподвижной жидкостной оболочки уменьшается. Оболочка исчезает при
r
o
>R.
Этот случай требует специального рассмотрения.
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
114
0,0
0,5
1,0
1,5
2,0
00,020,040,060,080,10,120,140,16
расстояние от центра скважины, м
скорость жидкости, м/с
ВПЖ
ДНС=4 Па
ВЖ
r
R
y
r
o
Рис. 14.8. Эпюры скоростей увлекаемого трубой потока вязкопластичной и
вязкой жидкостей в заколонном пространстве при спуске бурильной
колонны в скважину диаметром 300 мм со скоростью 2 м/с
Если уравнение (14.10) с учетом реологической функции (14.16) проинтег-
рировать в пределах от
τ
r
до
τ
R
, то получим:
- уравнение для определения скорости спутного потока:
ln
ooo
r
RR
r
u
⎛⎞
τ
ττ
τ
=−−
⎜⎟
η
τττ
⎝⎠
; (14.25)
- уравнение для определения напряжения на стенке трубы при заданной
скорости движения (после замены отношений напряжений на отношение ра-
диусов):
ln
oo
r
т
Rr
R
r
r
u
Rr
⎛⎞
ττ
τ
=−−
⎜⎟
ητ τ
⎝⎠
; (14.26)
- уравнение для определения скорости потока на любом расстоянии
y от оси
скважины при заданной скорости движения колонны
u
т
:
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
115
ln
ln
оо
rr
т
o о
rr
R
y
y
rRr
uu
R
r
rR
τ
τ
+−
τ
τ
=
τ
τ
+−
τ
τ
; (14.27)
- искомую формулу для определения
К
ск
:
()
()
22 33
112
22
3
21 11
ln
42 42 2 3
ск
AAA
R
KRr rR
r
AR r
⎡⎤
⎛⎞⎛ ⎞
=−−+−−−
⎜⎟⎜ ⎟
⎢⎥
−
⎝⎠⎝ ⎠
⎣⎦
, (14.28)
где
A
1
=
τ
o
R/(
τ
r
r); A
2
=
τ
o
/(
τ
r
r);
3
ln
oo
rr
R
r
A
Rr
τ
τ
=− −
τ
τ
.
0
0,05
0,1
0,15
0,2
0,25
0,3
0,35
0,4
0,45
0 2 4 6 8 101214161820
Динамическое напряжение сдвига, Па
Относительный расход спутного
потока
r
o
< R
r
o
> R
Рис. 14.9. Изменение коэффициента К
ск
в зависимости от динамического
напряжения сдвига бурового раствора: диаметр долота – 216 мм, диаметр
труб – 127 мм, структурная вязкость буровых растворов − 0,02 Па
*
с
Величина
τ
r
, как ранее отмечалось, определяется численным методом из
уравнения (14.26).
На рис. 14.9 приведен пример зависимости
К
ск
от
τ
o
для некоторого
конкретного соотношения диаметров скважины и трубы.
Расчеты показали, что:
- при
τ
o
→ 0 К
ск
приближается к значению, равному К
ск
для вязких жидко-
стей, рассчитанному по формуле (14.25);
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
116
- при малых
τ
o
(ориентировочно до 3 Па), когда на стенке скважины не
формируется “структурная” оболочка, расчет ведется для случая, когда
R<r
o
;
- при
τ
o
>3 Па имеет место существенное уменьшения спутного
(увлекаемого трубами) потока и величины
К
ск
;
- применение буровых растворов, обладающих структурной прочностью,
сопровождается значительным, подчас кратным, уменьшением
Q
ск
, что обес-
печивает уменьшение величины
Q
экв
.
14.2.4. Расчет коэффициента К
ск
для случая, когда в скважине
псевдопластичная (степенная) жидкость
В последние годы, в связи с широким применением полимер-глинистых
растворов с малым содержанием твердой фазы, все большее применение нахо-
дят псевдопластичные жидкости, тяготеющие к степенным реологическим мо-
делям с реологической функцией
ϕ
(
τ
)=(
τ
/K)
1/n
. (14.29)
Если подставить эту функцию, как в предыдущих случаях с вязкопластич-
ными жидкостями, в уравнение (14.10) и затем проинтегрировать последнее в
пределах от
τ
до
τ
R
, то с учетом соотношения
τ
R
=
τ
r
r/R получим:
1
1
1/
(1)
rr
n
rr
u
KR
α+
α+
⎡
⎤
ττ
⎛⎞
=τ−
⎢
⎥
⎜⎟
α+
⎝⎠
⎢
⎥
⎣
⎦
, (14.30)
где
α
= 1/n - 2.
На стенке трубы u = u
т
, а
τ
=
τ
r
. После несложных преобразований имеем:
1
2
1/
1
(1)
r
т
n
r
r
u
KR
α+
α+
⎡
⎤
τ
⎛⎞
=−
⎢
⎥
⎜⎟
α+
⎝⎠
⎢
⎥
⎣
⎦
. (14.31)
В отличие от предыдущих случаев в данном уравнении
τ
r
можно выразить в
явном виде:
1
2
1/
1
(1)
1
n
т
r
uK
r
r
R
α
+
α+
⎧⎫
⎪⎪
α+
⎪⎪
τ=
⎨⎬
⎡⎤
⎛⎞
⎪⎪
−
⎢⎥
⎜⎟
⎪⎪
⎝⎠
⎢⎥
⎣⎦
⎩⎭
. (14.32)
Заменяя в (14.30) отношения напряжений на отношения радиусов, приходим
к результату:
1
1
2
1/
(1)
r
n
r
rr
u
KyR
α+
α+
α+
⎡
⎤
⎛⎞
τ
⎛⎞
=−
⎢
⎥
⎜⎟
⎜⎟
α+
⎝⎠
⎝⎠
⎢
⎥
⎣
⎦
. (14.33)
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
117
Чтобы найти
Q
ск
, нужно просуммировать расходы dQ
ск
элементарных стру-
ек толщиной
dy на радиусе y при скорости струй u:
dQ
ск
= 2
π
y
.
dy
.
u. (14.34)
Подставив (14.33) в это уравнение и интегрируя в пределах от
r до R, по-
сле некоторых преобразований получаем:
()
22
13
12
22
11
22
ск
BB
BB
QrR
rR
β− β−
⎛⎞
=−+ −
⎜⎟
β−
⎝⎠
, (14.35)
где
()
2
1
1/
2
1
r
n
r
B
K
α+
πτ
=
α+
,
1
2
r
B
R
α
+
⎛⎞
=
⎜⎟
⎝⎠
, B
3
=r
α
+1
,
β
=
α
+1.
Искомую величину коэффициента K
ск
вычислим по формуле:
22
()
ск
ск
т
Q
K
Rru
=
π−
. (14.36)
Анализ показал, что эпюра скоростей, а следовательно, и величина
К
ск
не за-
висят от индекса консистентности
К, а зависят только от показателя реологиче-
ского поведения жидкости
n. Установлено также, что K
ск
не зависит от скоро-
сти движения колонны.
На рис. 14.10 показаны зависимости
u(y) при различных значениях n. Вид-
но, что эпюра скоростей видоизменяется не столь сильно, как при вязкопла-
стичных жидкостях с явно выраженной прочностью структуры.
Рис. 14.11 иллюстрирует зависимость
К
ск
(n) для тех же расчетных условий,
причем
n=0,99 имитирует вязкую жидкость. Из анализа приведенной зависи-
мости можно заключить, что влияние
n на К
ск
весьма существенно, и им пре-
небрегать недопустимо.
Итак, мы рассмотрели три случая, отличающиеся реологическими моделями
жидкостей, которыми заполнена скважина. Установили, что при вязкой жидко-
сти эпюра скоростей в заколонном пространстве,
а также величина К
ск
зависят
только от соотношения диаметров скважины и труб.
Наиболее сложным является вариант применения вязкопластичной (бинга-
мовской) жидкости, состоящий из двух подвариантов в зависимости от соот-
ношения
r
o
и R, когда вначале нужно определить r
o
, чтобы идентифицировать
конкретный подвариант решения, затем вычислить
τ
r
на стенке трубы и, нако-
нец, определить искомую величину
К
ск
. Установили также, что при использо-
вании степенной жидкости
К
ск
зависит от реологических параметров, но только
от показателя реологического поведения
n.
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
118
0
0,2
0,4
0,6
0,8
1
1,2
1,4
1,6
1,8
2
020406080100
Расстояние от оси скважины, мм
Скорость жидкости, м/с
n=0,99
n=0,6
n=0,5
n=0,4
n=0,3
R
r
u
Рис. 14.10. Распределение скоростей жидкости в заколонном пространстве
при степенной жидкости: диаметр скважины 215,9 мм, диаметр труб – 127 мм
0,3
0,32
0,34
0,36
0,38
0,4
0,42
0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
Показатель реологического поведения
n
Относительная величина спутного
потока (коэффициент
К
ск
)
Рис. 14.11. Зависимость коэффициента К
ск
от показателя реологического
поведения n: условия расчета те же, что на рис. 14.10
Часть II. Раздел 14. Гидродинамические давления, возникающие при движении
колонны труб в скважине
119
14.3. Методика расчета допустимой скорости спуска (подъема)
"закрытой" колонны в скважине
Данная задача является обратной по отношению к описанным выше и пред-
полагает введение технологических ограничений на величину либо гидродина-
мического давления, либо общего давления на “слабый” пласт
р
сл
(с добавлени-
ем гидростатического).
Решение этой задачи сводится к определению такой скорости спуска
u
т.доп
колонны, при которой выполняется равенство
р
сл
=ρgL
сл
+р
гд
.
Задача является типично компьютерной, поскольку требуется организация
итерационного цикла по определению искомого
р
гд
путем изменения u
т
.
Задача определения
р
гд
, р
гд.доп
, u
т
и u
т.доп
существенно осложняется тем,
что не известен расход жидкости, устремляющейся, например, при спуске ко-
лонны в трубы. Решение задачи приведено в разделе 15.
Часть III. Раздел 15. Расчет гидродинамических давлений при равномерном
движении “открытых” трубных колонн в скважине
120
Часть III. Проектирование и оптимизация гидравлических
программ буровых процессов
15. Расчет гидродинамических давлений при равномерном
движении “открытых” трубных колонн в скважине
В разделе 14 изложена (с обоснованием) методика расчета гидродинамиче-
ского давления, возникающего при движении “закрытой” колонны, когда жид-
кость не может поступать внутрь труб при их спуске и вытекать из труб
при
подъеме. Нет нужды доказывать, что “открытые” колонны труб (если долотные
отверстия не забиты или в колонне нет обратного клапана, а также, например,
нет пластоиспытателя) создают меньшие по величине гидродинамические дав-
ления. Трубное пространство, куда при спуске устремляется (под действием
возникающего давления) часть вытесняемой жидкости, является как бы кана-
лом “утечки
” давления.
15.1. Методика расчета гидродинамического давления
при равномерном движении “открытой” колонны труб
Представим себе двухразмерную бурильную колонну (рис. 15.1), составлен-
ную из УБТ диаметром
d
1
и бурильных труб диаметром d
2
в скважине диамет-
ром
D. Пусть колонна заканчивается долотом с заданной площадью сечения
промывочных каналов (на рис. 15.1 промывочный узел представлен условно
одним отверстием).
При равномерном движении вниз в колонну устремляется жидкость с неко-
торым установившимся расходом
Q
в
, являющимся частью расхода вытеснения
Q
вт
:
2
1
1
4
т
вт
du
Q
π
=
или
2
2
2
4
т
вт
du
Q
π
=
.
Иначе говоря, если часть потока “уходит” с расходом
Q
в
в трубы, то ровно
на столько убывает расход вытеснения как за УБТ, так и за бурильными труба-
ми. А поскольку эквивалентный расход
Q
экв
состоит из Q
вт
и Q
ск
(см. часть II,
раздел 14), то и
Q
экв
уменьшается на каждом участке на Q
в
. Итак, для “откры-
той” колонны эквивалентный расход равен
Q
экв
, за трубой диаметром d
i
равен
Q
экв
для “закрытой” колонны за вычетом Q
в
:
Q
экв.i
= Q
cк.i
+ Q
вт.i
- Q
в
или более конкретно (с заменой d
i
на d
нi
):