Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 100
Скачиваний: 0
В решениях (7.88), |
(7.89) через / и ср обозначено: |
|
|
/ = |
— |
[Я 0/ і ( а „ Д 0) - ■г0/ 0 (а„г0)], |
(7.91) |
|
ОСЛ |
|
|
ф = ч р |
[Я < Л (а«Я о) — i'fjjх (а „ /‘о)|. |
(7.92) |
|
|
и.п |
|
|
Выведенные выше |
аналитические зависимости (7.69), |
(7.70),. |
(7.88) позволяют получить распределение тепловых потоков между элементами пары трепня резец — массив в процессе бурения. Пе рейдем к рассмотрению температурного поля твердосплавной корон ки действительной конфигурации. Коронка представляет собой па раллелепипед, на верхнем основании которого действует постоянный плоский равномерно распределенный источник тепла. Имеющееся в литературе [16] решение задачи о температурном поле параллеле пипеда с плоским источником энергии, занимающим одну из граней, относится лишь к случаю, когда на боковых гранях поддерживается постоянная температура. В монографии [161 приведено также ре шение для стержня, в торец которого входит тепловой поток, пере дающийся кондуктивпо по стержню, рассеиваясь при этом с поверх
ности |
в окружающую среду. Температурное поле по |
сечению при |
|
этом, |
естественно, считается равномерным |
(градиент |
температуры |
в поперечном сечении равен нулю). |
|
|
|
Отметим, что приведенные в работе [16] |
аналитические зависимо |
сти для расчета температурного поля параллелепипеда включают в себя двойной бесконечный знакопеременный ряд. Расчеты по этим зависимостям, выполненные нами с помощью ЭВМ «Мир», дали ре зультаты, расходящиеся с экспериментальными данными на величину одного порядка, т. е- примерно в десять раз.
Наиболее приемлемым для нашего случая является решение, полученное в работе С. А. Попова и В. М. Давыдова применительно' к шлифованию методом врезания. В этой работе задача формулируется следующим образом. Торец полубесконечиого стержня без теплоизо ляции нагревается тепловым потоком заданной интенсивности д, действующим в течение времени т. Необходимо найти распределении температуры по длине стержня во времени при условии, что на торце- и боковой поверхности стержня осуществляются теплообмен с гра ничными условиями III рода со средой постоянной температуры Г 0. Учитывая малый поперечный размер твердосплавной коронки, при меняемой для оснащения бурильного инструмента, и равномерное распределение теплового источника на ее торце, принимаем эту за дачу одномерной, как это принято с теми же допущениями мно гими исследователями.
Теплоотдача с боковой поверхности учитывается в дифференци альном уравнении теплопроводности в качестве отрицательного'
источника тепла [24] |
|
U = a1[T(j/,t) — T0]-jr, |
(7.93) |
149
где а-у — коэффициент |
теплоотдачи |
боковой поверхности |
стержня: |
|
Т (у, т) — функция |
температурного |
состояния стержня; |
Р, F — |
|
соответственно периметр и площадь стержня. |
|
|||
Краевая задача для стержня записывается в виде |
|
|||
дТ{у, т) _ „ д*т (ю, т) |
ß [Т (уі, т) — Г01; |
(7.94) |
||
дх |
|
Ö1J- |
||
|
|
|
||
начальное условие |
|
Т(у,0) = |
Т0; |
(7.95) |
граничные условия |
|
|||
|
Т (оо, т) = Г0; |
(7.96) |
||
|
|
|||
А 9 |
T t j |
T) + g — |
(Q^ X) — Г0] = 0, |
(7.97) |
где А — теплопроводность; а — температуропроводность; а — коэф
фициент теплоотдачи с поверхности торца; ß = |
(с — теплоем |
кость; у — удельный вес). |
|
Решение задачп (7.94) — (7.97) операционным |
методом имеет |
вид |
|
У а |
|
|
|
|
|
|
|
|
|
ехр[ у т + |
( і ^ - р |
) т] |
егіс I |
2 |
!' |
+ |
У ах |
(7.98) |
|
|
а -а |
|
|
|
у ах |
А |
|
|
|
|
~Т~ |
|
|
|
|
|
|
|
|
■* ь |
|
приводится к виду |
~ |
Q |
|||||
При а ^ = 0 |
уравнение (7.98) |
||||||||
т {у> Т)- |
гр __ V |
erfc (— |
|
■exp |
/ |
а |
. а-а |
\ |
|
0“ 1Г |
\2 V a x |
|
( т у + — т) х |
|
|||||
|
X erfc |
|
■Уах |
|
|
|
(7.99) |
Температура на поверхности контакта определяется по формуле
(7.98) при у = 0
T (0 ,T ) - T 0 = J L [ l - e x p ( ^ J - ü ) er fc ( - J - y ^ F ) ] . (7.100)
По формулам (7.98)—(7.100) были вычислены значения темпера туры при трении твердосплавного образца (штабика) о горную по роду (песчаник). Для расчета принимались следующие исходные
150
дапі-ше: Т 0= 20° С; X = 58,2 Вт/м-°С; I = 0,005 м. Данные получены для одной скорости вращения, равной 2,0 м/с, когда расход мощно сти составлял 0,3 кВт.
Сопоставление результатов аналитического решения с данными электрического моделирования температурного поля и эксперимен тальными данными, полученными на физической модели, приведено' ниже.
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ В БУРИЛЬНОМ ИНСТРУМЕНТЕ
ИГОРНОМ М АССИВЕ
ВПРОЦЕССЕ БУРЕНИЯ
Изложенное выше аналитическое решение задачи о температурном поле бурильного инструмента и горного массива в процессе бурения получено, как указывалось, для стационарного случая при упрощен ной схеме взаимодействия термически сопряженных тел. Для не стационарного случая получить замкнутое аналитическое решение,, пригодное для выполнения расчетов, даже в такой постановке невоз можно. В связи с этим представляет
интерес |
математическое |
моделирова |
|
ние на і?-сетках интересующей нас |
|||
задачи в постановке, отвечающей |
|||
реальным условиям. |
|
|
|
Для |
электрического |
моделиро |
|
вания |
температурных |
полей |
пары |
трения |
резец — породный |
массив |
были избраны следующие |
варианты |
|
|
|||
постановки |
задачи: |
|
Рис. 34. |
Схема коммутации со |
||
1) боковые |
поверхности |
твердо |
||||
противлений при электромодели |
||||||
сплавной коронки теплоизолированы; |
||||||
ровании |
температурного поля |
|||||
контакт между коронкой и массивом |
па забое скважины по 3-му ва |
|||||
в процессе |
опыта происходит иа по |
|
рианту. |
|||
верхности z = |
0; |
|
|
|
2) на боковых поверхностях коронки осуществляются граничные условия теплообмена III рода; контакт — на поверхности z = 0; 3) на боковых поверхностях коронки осуществляются граничные условия III рода; контакт между коронкой и массивом происходит
по контуру 1—2—3—4 (рис. 34).
В первых двух вариантах модели не учитывается разрушение горного массива в процессе бурения (внедрение породоразрушающего инструмента в массив); в третьем варианте внедрение коронки в мас сив учитывается статически, т. е. для одного положения заглублен ной коронки (на глубину 5 мм).
При этом для всех вариантов постановки задачи были сделаны следующие допущения:
1)тепловой поток в горном массиве распространяется по нормали
кповерхности трения и параллельно этой поверхности;
151
2)поверхность массива z = 0 теплоизолирована, а на остальных границах модели массива температура конечна и равна начальной;
3)твердосплавная коронка представляет собой часть цилиндри
ческой стенки (сегмент), условно имеющую квадратное поперечное ■сечение со стороной I = R 0 — г0\
4)теплофизические свойства коронки и массива не изменяются во времени и пространстве;
5)температура среды, окружающей резец, в процессе опыта по
стоянна.
Исходя пз перечисленных допущений, были рассчитаны Д-сеткп
омических |
сопротивлений. При |
конструировании сеточной модели |
|||||||
(выбор |
ее |
размеров и |
шагов |
измерений) |
учитывались |
реальные |
|||
|
|
|
|
|
размеры коронки, различие в теп- |
||||
|
|
|
|
|
лофизнческих свойствах резца и мас |
||||
|
|
|
|
|
сива и влияние числа шагов на по |
||||
|
|
|
|
|
грешность опыта. Сеточная модель |
||||
|
|
|
h, |
|
№ 1 использовалась для моделиро |
||||
|
|
hrl |
hr, |
вания 1-го и 2-го, |
а модель № 2 — |
||||
|
|
|
для 3-го варианта задачи. Следует |
||||||
|
|
|
|
|
|||||
|
|
|
|
|
отметить, что практически модель №2 |
||||
|
|
|
|
|
была получена из модели № 1 соот |
||||
|
|
|
|
|
ветствующим |
изменением |
коммута |
||
|
|
|
|
|
ции в окрестности поверхности кон |
||||
Рпс. 35. Схема элементарного |
такта. |
|
|
|
|||||
объема тела при электромоделн- |
|
Для моделирования исследуемой |
|||||||
ровапші |
температурного поля на |
задачи, как указывалось выше, была |
|||||||
забое скважины. |
|
принята схема «узел внутри объема». |
|||||||
|
|
|
|
|
Схема линейных |
размеров |
и распо |
ложения элементарного объема модели по отношению к оси сим метрии показана на рис. 35.
Поскольку задача осесимметрична и моделируется в цилиндриче ских координатах, расчет сопротивлений і?-сеток выполнялся по
методике Л. А. Коздобы [26]. |
Расчетные зависимости для случая |
||||
«узел внутри объема» при hpj — hpо |
и |
|
|
2 имеют следующий |
|
вид: |
2hr |
|
|
|
|
R, |
|
|
(7.101) |
||
|
|
|
|
||
U z ( 1 + ^ r ) r |
|
|
|||
я.-, |
2hr |
|
|
(7.102) |
|
/ |
h~ \ |
|
’ |
||
|
|
||||
x,!z( 1- |
i r ) |
r |
|
|
|
D |
|
2hz |
_ |
|
(7.103) |
(1>a) |
khrr |
’ |
|
||
|
|
||||
я - |
860N |
|
|
(7.104) |
|
f |
; |
|
|
||
Rt |
cyhzhrr * |
|
|
(7.105) |
|
|
|
|
152