Другие аналитические зависимости для определения температуры бурового раствора
В дальнейшем появились иные аналитические зависимости, также построенные по принципу теплового баланса. Но при этом использовались другие составляющие тепла, иначе учитывалась зависимость пластовой температуры от глубины и т. п. Уравнения теплового баланса, получившие наибольшую известность, записаны в обозначениях авторов в табл. 22. Там же для ясности даны — стрелки, которые показывают, что речь идет о температуре восходящего ( f ) или нисходящего (| ) потока. Буквы «н», «к», «п» означают «начало», «конец», «пласт». Из рассмотрения табл* 22 вытекает, что в принципе все приведенные в ней выражения в конечном счете являются как бы усовершенствованными зависимостями (VIII.7) и (VIII.8).
Различные авторы неодинаково решают и сами уравнения теплового баланса относительно температуры восходящего или нисходящего потока. Так, одни авторы от уравнений теплового баланса переходят к дифференциальному уравнению второго порядка согласно приему, аналогичному тому, который использован при составлении дифференциального уравнения (VIII.14). Такие уравнения получены А. А. Афанасьевым, И. А. Чарным, Б. Б. Кудряшовым, А. Н. Щербанем и В. П. Черняком.
• Другая группа авторов — А. М. Погорельский и И. А. Кулиев, Н. Р. Акопян, Г. А. Обабко, Ю. М. Проселков, А. Г. Потапов, А. Н. Щербань и В. П. Черняк (упрощенное решение) получают зависимости для определения температуры, используя иные пути решения. Так, например, Н. Р. Акопян, Г. А. Обабко, Ю. М. Проселков вместо дифференциального уравнения второго порядка получают обыкновенное дифференциальное уравнение первого порядка, принимая условие
. — dti — (It, .
где dx— приращение температурного напора со стороны жидкости в кольцевом пространстве к нисходящему потоку в трубах, причем
12у — z)
Я |
г =
Авторы другой работы [69] предложили упрощенный вариант, при котором уравнения теплового баланса записываются в конечно-разностном виде, соответственно для бурильной колонны, при-
Автор |
Нисходящий поток |
Восходящий поток |
Б. И. Есьман |
t 1 fcjjidj (t2 — t{) dz = QiCiYldtl — f , QiVi’i. + c dz E |
k2nd2 (t„ — t2) dz =—Q2c2y2dt2A- dz + QiCiYidti +—(Q2y2i2+QtYiii> |
А. М. Погорельский И. А. Кулиев |
Q = k^HAE, Q + ZAQi = Gc x | к 4 H X ( ^2 ^2) lAQi — тепло внутренних источников |
|
А. А. Афанасьев |
‘ -1 4 6Q* = kndjh (т — t) dz |
* — t П 6Q2 = а2л^2 (л — 0) dz + tf kiJidjn (t — t) dz |
И. А. Чарный |
4 / dT, cGT! — cG Ti— — dx j = Л ) 4 t = ki2ndl{Tl — T2) dx |
f t dT2 cG ( T2 + ■ ^ dx) cGT2 — f — 4 + + ^12^1 (T i — T2jdx = t П = k2giid2 [T2 T3 (x)J dx |
Н. Р. Акопян, Г. А. Обабко, Ю. М. Проселков |
— |
annD (tn —12) dz = Gcdt2 — . G (4 + /2) . — Gcdt^ + J dz |
Б. Б. Кудряшов |
t — . 4 ‘ dQi = kn (t2 — <i) dh — j — AGiidh; dQ = Gcdti |
tn t dQ2 — kxnD (/„ + an —■ t2) dh + -)- AGi2dh — kn (t2 4) dh, dQ2 = — Gcdt2 |
А. Г. Потапов |
t 4 AQTn = k^di (tKn — tTn) dz |
AQ = dqH — dQ3 hndj X t 4 X (tкп — tTn) dz qn — тепло, полученное в нижней части скважины qB — тепло, отданное в верхней части скважины |
Восходящий поток |
Автор |
Нисходящий поток
A. Н. Щербань B. П. Черняк |
Gcp (^2 — ^i)- —^ ^тр^тр Ei) |
)к I н |
f U + h X (1 — 62) ^пср „ ^ |
-f — kTpFтр (1 ej) X |
f k^iFjp |
-2~r-i~) |
GcP (*3 К — *4 Н) = kxFci Х |
забойной зоны и кольцевого (затрубного) пространства. Своеобразное решение предложено А. Г. Потаповым [50]; ряд специфических допущений делают А. М. Погорельский, И. А. Кулиев и т. д. . ‘
В литературе можно встретить ссылки на работы Егера, Реми, Реймонда, в которых также рассматриваются аналитические методы определения температуры в бурящихся скважинах.
Каждая. из указанных зависимостей может быть использована для определения температуры потока бурового раствора, но для того, чтобы выбрать одну из них в качестве расчетной, необходимо учесть следующие обстоятельства.
Первые аналитические формулы (А. X. Мирзаджанзаде,
А. А. Мирзоян, Б. И. Есьман, М. А. Абдинов оказались ошибочными, так как, копируя задачу И. А. Парного о подогреве трубопровода, они не учитывали специфику бурящихся скважин, т. е. условие равенства температры нисходящего и восходящего потоков на забое. В дальнейшем, отмеченный недостаток был устранен (Б. И. Есьман, Г. Я. Дедусенко, Е. А. Яишниковой), но оказалось, что необходимо еще учитывать подъем температуры на забое за счет работы долота.
Указанный скачок температуры и, кроме того, приток тепла за счет работы забойного двигателя, а также тепловыделения в то — коподводе (при бурении электробурами) впервые были учтены в зависимости, предложенной А. М. Погорельским и И. А. Кулиевым. Однако поскольку они прибегли к ряду допущений и отошли от строгого решения, то предложенная ими зависимость перестала быть чисто аналитической и скорее должна рассматриваться как полуэмпирическая. Сразу же отметим, что аналогичные замечания относятся и к формулам, предложенным Н. Р. Акопяном, Г. А. Обабко, Ю. М. Проселковым.
Более строгое решение этой же задачи принадлежит А. А. Афанасьеву. Однако ни в зависимости А. А. Афанасьева, ни в зависимостях многих других авторов не учитывался фактор времени,
т. е. влияние на теплообмен продолжительности циркуляции бурового раствора.
Попытка учесть это обстоятельство была предпринята И. А. Чар — ным, который использовал приближенный метод, основанный на схеме последовательной смены квазистационарных состояний, и ввел в рассмотрение радиус теплового влияния, величина которого является функцией времени. Вместе с тем он ввел условие, согласно которому время промывки скважины определяется через коммерческую скорость бурения. Но, как справедливо заметил Б. Б. Кудряшов, учитывать время промывки на основе коммерческой скорости бурения неверно, так как процессы теплообмена при движении бурового раствора в период бурения и в период простоя совершенно различны. Он предложил свою зависимость, в которой учет изменения интенсивности теплообмена между циркулирующим буровым раствором и стенками скважины выполнен ’ за счет использования коэффициента (нестационарного теплообмена kx:
tt = m#r’h + щег*н + b + t0+Th, (VIII.24)
где
TOC o "1-5" h z • С1Г1еГ>н + С2 — Сгг2 eriH + С2 —
— г2 т 1
«3 о ———————————————————— !——————————————————————————— /7о — . ;■ •
2 В1 В1
Причем ■ ‘
ktt
Ci = tlH — t0 + a — b-, С2 = (а — А^з); В1 = г#?1*1 — г2ег*я;
(JC
Gc Г А л G
а==~ы { ~01 ь = А~к^о^ + ^’ G = Qy —
Здесь ts — температура восходящего потока бурового раствора в затрубном пространстве; h — глубина (текущая координата); t0 температура пород на поверхности; Г — геотермический градиент; tн — начальная (устьевая) температура нисходящего потока бурового раствора в бурильных трубах; k — коэффициент теплопередачи через стенку бурильной трубы, отнесенный к 1 м длины трубы; G — массовый расход бурового раствора (кг/ч); с—: удельная весовая теплоемкость бурового раствора; At3 — прирост температуры бурового раствора за счет генерируемого на забое тепла
— N
А^3 ЕЕ 860 —— , ис
где N — мощность, затрачиваемая на забое; г и г2 — корни харак-
л (kxD теристического уравнения rt = — g^- ( —+ V J ;
Л Т л / ‘
гг=~о~ ~ V = y ■-r — + krkD. (VIII.25)
kx—коэффициент нестационарного теплообмена; D — диаметр скважины; Л =1/427; Я — глубина забоя.
Аналогично для температуры жидкости, движущейся в бурильных трубах,
t1 = mi&rih — f — n1eraft—a — f — b + t0 + Г/г, (VIII.26> »
где
—C2 C1rleriH + C2
Щ = — gi ; rif = gi •
Приблизительно одновременно с Б. Б. Кудряшовым, но независимо от него аналогичную зависимость предложили А. Н. Щер — бань и В. П. Черняк [69]. В отличие от формулы Б. Б. Кудряшова в этой зависимости помимо коэффициента нестационарного теплообмена учитывается еще эксцентричное расположение колонны бурильных труб в скважине, изменение диаметра колонн с глубиной, схема циркуляции бурового раствора (прямая или обратная),, а также нелинейность изменения температуры горных пород.
Таким образом, формула А. Н. Щербаня и В. П. Черняка наиболее универсальна из всех известных формул и в принципе может учитывать почти все нюансы теплообмена, происходящего в бурящейся скважине. Однако для практических расчетов учитывать все эти тонкости обычно нет никакой необходимости, а в таком случае весьма громоздкая формула А. Н. Щербаня и
В. П. Черняка автоматически переходит в формулу Б. Б. Кудряшова.
И. М. Астрахан и В. И. Марон составили дифференциальные уравнения, характеризующие нестационарные условия теплообмена в бурящихся скважинах, и, используя преобразование Лапласа— Карсона, в результате достаточно сложных приемов нашли способ для обращения изображений функций A, t2 и Tn(z). Однако определение температуры бурового раствора по методике возможно только при наличии специальной программы на ЭЦВМ, что значительно снижает вероятность использования этой методики для практических расчетов.
На основании сказанного представляется целесообразным в качестве основной расчетной формулы принять формулу Б. Б. Кудряшова, так как она является наиболее простой и в то же время правильно отображает физическую сущность рассматриваемого’ явления. При более точных расчетах следует пользоваться формулами А. Н. Щербаня — В. П. Черняка.
Для выяснения распределения температуры по стволу бурящейся скважины может быть использован и экспресс-метод, предложенный Ю. М. Проселковым [53]. В основе метода лежит утверждение, что наибольшая темпер. атура в кольцевом пространстве отличается на расстоянии около 2/3Я от забоя (Н — глубина скважины). Исходя из такой закономерности легко графически (рис. 45) найти эту точку (bi) на геотерме tobbtm а затем используя зависимости.
hy = А + А — Ai; tig — t2y — ДА; A == tn + A — t2y
{где tBi, t2y, to, in, ty, t3—температуры соответственно в точке bx выходящего раствора; нейтрального слоя; пластовая на забое; входящего раствора; на забое в процессе циркуляции; At-—температура охлаждения раствора в желобах; зимой At = 8—10 °С; летом— 4—5°С), найти точки-а, аь Ь2, соединяя которые между собой, можно получить наглядное представление о распределении температуры в стволе бурящейся скважины. При этом отрезок ab2 характеризует распределение температуры в бурильных трубах, а ломаная Ь2Ьау— в кольцевом пространстве. По утверждению автора работы [53], экспресс — метод дает погрешность не выше 10 %. .
Рис. 45. Приближенное графическое представление распределения температуры в бурящейся скважине |
Ниже приводится формула, позволяющая судить о характере распределения температуры в восходящем потоке бурового раствора. Для получения этой формулы достаточно рассмотреть рис. 20 и убедиться, что температура в точке М, расположенной на произвольной глубине h в скважине с забоем Н, может быть найдена из выражения
*ь. н = — + {у н, ■ (V.30a)
которое аналогично выражению (V.30). 4
С учетом же зависимости (V.31) и того, что 1 [с = а, выражению (V.30a) можно придать более удобный вид
. — thH = ah + bH + t0. (VIII.27)
. Коэффициенты а и b находятся по наблюдениям в ряде скважин (см. табл. 3).
В случае, если Ъ = Н, то зависимость (У. ЗОа) автоматически
превращается в зависимость (V.30), а выражение (VIII.27) переходит в формулу (V.32), причем становится справедливым условие, что а + Ь = ср.