Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 103
Скачиваний: 0
производить набор сопротивлений с четырьмя знаками. Конструк ция его потенциометров обеспечивает плавное прохождение по всему диапазону от 0 до 100% без провалов напряжения и зада ние граничных условий с точностью до четвертого знака. Заданные напряжения контролируются по измерительному устройству.
При электромоделированцц расчетная часть бурильной колонны разбивалась семейством лпннй, параллельных координатным осям г
п |
ер с равномерным шагом разбивки вдоль радиуса |
h — 2,5 мм, |
|
а |
вдоль окружности — с неравномерным шагом Дф; |
в том числе: |
|
трп шага по 10°, |
три шага по 20° и три шага по 30°. Наружный диа |
||
метр бурильной |
колонны был принят равным 146 мм. При этом |
Рис. 28. Схема элементарных объемов |
при |
||
электромоделпровашш |
температурных |
полей |
|
прп роторном бурении. |
|
|
|
а — схема коммутации сопротивлений; б — схема |
|||
лпнейпых |
размеров. |
|
|
координаты г и ф изменяются в пределах 63 |
г ^ |
73 мм (толщина |
|
стенки трубы 10 мм), 0 ^ ф ^ |
180°. |
|
|
Схема элементарного объема /2-сетки для моделируемых тел показана на рис. 28,а.
Зависимости для расчета сопротивлений, составляющих сетку,
можно получить из |
соответствующих |
уравнений теплопроводности |
|||
с учетом граничных условий. |
|
|
|
||
Уравнение теплопроводности |
(6.3) |
переписывалось для |
этого |
||
в виде |
ди , |
д X |
ди \ |
ди |
|
д_ |
|
||||
~дг |
Хг ~г |
Т ’Щ ) ~ ~ сосуг |
(6.43) |
После этого конечно-разностное уравиеиие для любой внутрен ней точки данной расчетной области можно записать следующим образом:
Mr, |
ho |
Афі + Афз |
|
~2 |
2/22 |
||
|
|||
|
|
h o “Ъ h i ] |
|
|
|
r0 2 Дфі |
|
|
|
( h o |
|
|
|
|
|
Афі і-Дсрз ( и , - и 0) |
|
|
|
|
A |
h‘2+ hj\ |
2/14 |
|
{ U x - |
u |
о) |
(U3 — U0) = |
|
||
1 |
” |
и' ' |
Г 0 |
2 Дфз |
|
|
+ І 14) (Дфі + Д ф з ) |
_____ |
( 9и \ |
(6.44) |
|||
|
4 |
|
|
mcT,4 W V |
||
|
|
|
|
132
где Uх, Uо, |
U3, Uit U0 — значения температуры в рассматриваемом |
и четырех |
соседних узлах сетки. |
Из уравнения (6.44) получаем формулы для расчета сопротивле ний сетки:
2h,. Кг
А, ^го д- "7Г") (Афі + Афз)
огАфіг-А'! .
1 ^ (/г 3 + / г 4) »
о |
2ДфзrKj , |
3 |
А(Л2 + /м) ’ |
________2hjKx_______
В і
^ ( Г°— jp)
(6.45)
(6.46)
(6.47)
(6.48)
Тепловой поток, обусловленный трением бурильной колонны об обсадную, рассматривался как источник, действующий в зоне
контакта. |
|
|
|
Сопротивление |
источника: |
|
|
|
Д„ = (Л2+ |
и - Кі - Ко |
(6.49) |
|
Л4) (Дфх+Дфз) шсуго |
||
Напряжения, |
задаваемые |
на свободные концы Д н, |
|
г - - (£).*•• (6.50)
где коэффициент К х выбран из условия наилучшего использования магазина сопротивлений (Кх = 50-103), а коэффициент К 2 выбран из условия максимального увеличения соотношения между сопро тивлениями сетки и источников.
Для реализации граничного условия (6.11) запишем конечно
разностное уравнение для |
граничных точек |
области: |
|||||||
|
|
|
|
|
|
(Афі -j- Дфз) а 2 |
[f2 — С/„] + |
||
( г0 + X ) { - У \Х АЪ ) (С/з - ^ + го |
|
|
|
||||||
X h2 J-h± |
^ - t/o ) + ^ ^ ( t / 3 - C / o ) |
= |
|||||||
+ - 7 |
T |
^ |
|||||||
= |
_ |
(fea + |
ftj) |
(Дфі + |
Дфз) |
V |
( |
dU \ |
(6.51) |
|
|
|
|
4 |
|
\ |
дф / о ’ |
|
|
откуда определим |
величину |
і?вн |
внешнего |
сопротивления |
|||||
|
|
Р |
_ |
Жг |
1 1 |
|
|
(6.52) |
|
|
|
|
|
Дфі + |
Дфз |
а 2 г0 |
|
|
|
|
|
|
|
|
|
|
|||
Величина внешнего напряжения определяется равенством |
|||||||||
|
|
|
|
Vm = t2. |
|
|
|
(6.53) |
133
Уравнение теплопроводности (6.6) для обсадной колонны и мас сива перепишем в виде
|
д_ |
|
|
|
|
|
(6.54) |
|
дг |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Уравнение |
(6.54) в конечных |
разностях |
для |
любой внутренней |
|||
точки области имеет вид |
|
|
|
|
|
|
|
|
Лф12| 2Аф3- V* - Uo) + |
|
|
Дфі 4- Лфз |
( и * - и 6)- |
||
( Ч + х ) |
Ч |
“ Ч ) |
2/01 |
||||
|
/іо-)-/ij |
|
/?2 "Ь/і-І |
(U3- U 0) = |
о, |
(6.55) |
|
|
2Дфі " |
|
2 Дфз |
откуда формулы для расчета сопротивлений сетки будут иметь вид:
р |
2 Д ф і/?оЛ г . |
|
(6-56) |
1 |
(Й2 + Л4)А ’ |
|
|
|
|
||
1? ,= |
2/(2Л'і |
|
(6.57) |
|
|
||
4 ° |
А |
|
|
п |
2Лфз/?оА 'і _ |
|
(6.58) |
|
h2 + hi |
|
|
/?4= |
2/1,/щ |
|
(6.59) |
|
|
||
( /7 ° — I 1 ) (Д ф іТ Л ф з ) А |
|
|
|
где Л’, = 50-103, коэффициенты теплопроводности |
принимались: |
||
для труб А = 58,2 Вт/м-°С, |
для массива А -- 1,75 |
Вт/м-°С. |
|
Коэффициент теплоотдачи на стенках труб принимался равным |
|||
1160 Вт/м2-СС, температура промывочной жидкости |
80° С в затруб- |
||
иом пространстве и 70° С внутри бурильной колонны. |
Температура |
горных пород на расстоянии 115 мм от стенки скважины принима лась равной 110 °С.
Горизонтальные сопротивления сетки на границе раздела обсад
ной трубы и массива рассчитывались по формуле |
|
||
Ді = |
_______ 2/(]/рА х_______ |
(6.60) |
|
(Ло+/(.,) ( ^ Р + АГрУ„Та ^ |
|||
|
Процесс моделирования заключался в следующем. На сеточном поле интегратора замерялось распределение потенциалов исследуе мой области при заданных граничных условиях и со — 0. По рас пределению потенциалов на сетке, моделирующей бурильную ко лонну, определялась производная du1/ckp1 для равномерного шага разбивки сетки Асрх = Дср3 = Дер по формуле (рис. 28, б)
( ди\ \ _ |
Фі — Фз |
(6.61) |
I öcp Jo |
2 Дер |
|
134
и для неравномерного шага
( |
диі \ |
_ Аф§ (Фі— Фо) - Афі(Фз— Фо) |
|
г9) |
\ |
д(р J o |
Дфі Дфз (Дфх + Дфз) |
' |
I • ^ |
По значениям d u jd y определялись напряжения истоков согласно (6.49), которые задавались на свободные концы RUCT. В результате
Рис. 29. Изменение температуры в по перечном сечении бурящейся скважины с увеличением расстояния I от лішші кон такта бурильной и обсадной труб по дан ным электромоделпрования при роторном
бурении:
1 — на наружной поверхности бурильной трубы;
2 — на поверхности обсадной трубы; з — на гра нице горного массива н обсадной трубы; 4 —
температура промывочной жидкости.
моделирования было найдено распределение потенциалов на сетках, или, что то же, изменение температуры в окружном и радиальном направлениях в бурильной и обсадной колоннах и массиве при угло вых скоростях вращения бурильной колонны, равных соответ ственно 10, 100, 1000 рад/ч, исходя из которых построены темпера турные графики (рис. 29, а—в). На тех же рисунках нанесена рас четная кривая, (штриховая), полученная по формулам (6.21), (6.22).
Полученные результаты позволяют сделать следующие выводы.
1. При роторном бурении на поверхности контакта бурильной колонны и стенки скважины быстро наступает квазистационариый тепловой режим, характеристики которого зависят в основном от силы прижатия бурильной колонны к стенке скважины.
135
2. Повышение угловой скорости при квазистационарном режиме практически не влияет на вид температурного поля, что свидетель ствует об интенсивном отводе тепла промывочной жидкостью.
3. В теле бурильной колонны имеет место высокий температур ный перепад вдоль окружности между зоной контакта бурильной и обсадной колонн и непосредственно примыкающими к этой зоне участками колонн, подвергаемыми интенсивному охлаждению. На личие такого перепада, который может вызывать усталостные напря жения в теле обсадных труб и служить причиной их аварийного износа, необходимо учитывать при выборе конструкции, способа бурения и центровании бурильной колонны в стволе скважины.
Г Л А В А 7
ТЕМПЕРАТУРНОЕ ПОЛЕ РАБОЧЕЙ ЧАСТИ БУРИЛЬНОГО
ИНСТРУМЕНТА И ГОРНОГО МАССИВА В ПРОЦЕССЕ БУРЕНИЯ
СОСТОЯНИЕ ИССЛЕДОВАНИЙ ТЕМПЕРАТУРНЫХ ПОЛЕЙ БУРИЛЬНОГО ИНСТРУМЕНТА
ИМ АССИВА В ПРОЦЕССЕ БУРЕНИЯ
Вгл. 1 отмечалось существенное влияние температурного фак тора на эффективность работы бурильного инструмента.
Механизм разрушения забоя скважины рабочими органами бу ровых долот заключается в истирании и частичном резании горных пород. При этом согласно первому закону термодинамики количе ство выделившегося тепла эквивалентно затраченной работе, т. е. повышение температуры на поверхности контакта инструмента с за боем обусловливается превращением затрачиваемой на процесс раз рушения механической энергии в тепловую. Полная работа, реализуе мая долотом на разрушение забоя скважины, состоит из работы, за трачиваемой на пластическую и упругую деформацию породы, на трение рабочей поверхности долота о породу забоя, на резание — истирание и измельчение продуктов разрушения породы, на преодо ление сопротивления окружающей среды и др., причем главным источником тепла является теплота трения — скольжения рабочей поверхности долота о породу.
А. М. Абдулзаде рассмотрены тепловые явления при бурении алмазными долотами, имеющими произвольное геометрическое по
строение секторов, в предположении, что все тепло, поступающее в долото, передается промывочной жидкости. Полученное в этой работе выражение для определения общего количества тепла, вы деляемого при бурении, неудобно для практических расчетов. В ра боте С. М. Кулиева и сотрудников для определения температуры
137
нагрева элементарного объема лопасти долота турбобура в процессе бурения предлагается пользоваться уравнением
(7.1)
где /31і — температура окружающей среды, равная температуре при установившемся тепловом режиме промывки; qv — количество те пла, выделяемое единицей объема долота в единицу времени; к — коэффициент, учитывающий неравномерность выделения тепла и ог раниченные размеры лопасти; б — толщина лопасти; А. — теплопро водность материала долота; а — коэффициент теплоотдачи от лопа сти к промывочной жидкости; х — текущая координата. При х = = 0 выражение (7.1) определяет температуру в центре объема ло пасти, а при X — 8 — на ее поверхности.
Недостатком этой методики является трудность определения ве личин qv , к, которые в каждом отдельном случае необходимо за имствовать из опыта.
Ряд опубликованных за последние годы работ по изучению тем пературного поля при шлифовании и сверлении осиоваи на приме нении теории мгновенных тепловых источников, Так, например, в работе В. Ф. Совкниа и М. П. Шатунова температура рассчиты вается непосредственно по готовой функции влияния плоского ис точника, причем входящее в эту функцию количество тепла, посту пающее в инструмент, предполагается заранее известным. Как ука зывают те же авторы, поскольку функция мгновенного источника обращается в бесконечность па поверхности инструмента, то рассчи танная по ней температура может быть сколь угодно большой, что не соответствует действительности.
Кроме того, в применении к шлифованию метод источников обла дает следующим общим недостатком: функция источника сама по себе не отражает механизма распределения тепла между всеми взаимодей ствующими телами (коронка, массив). Чтобы учесть это взаимодей ствие, приходится делать дополнительные предположения.
Ученые Ленинградского горного |
института |
(Б. Б. Кудряшов, |
Н. Д. Михайлова, А. В. Касаточкин) |
выполнили |
комплекс работ |
по исследованию температурного режима алмазных буровых коронок и шарошечных, долот при бурении с продувкой. В результате этих работ получены ценные данные об оптимальных температурных ре жимах продувки в сочетании с рациональной конструкцией долот и подбором крупности алмазных зерен, позволяющих бурить разве дочные скважины алмазными коронками. Некоторым общим недостат ком этих работ является отсутствие в решаемых системах уравнений теплового баланса уравнения теплопроводности разбуриваемого массива, распределение температур в котором будет оказывать опре
деленное влияние |
и на температурное |
поле рабочей |
части долот. |
Г. А. Тирский |
рассматривал задачу |
о нагревании |
при помощи |
трения двух теплопроводящих полупространств, из которых одно покоится, а другое движется с постоянной скоростью вдоль линии
138