Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.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 рода; контакт между коронкой и массивом происходит

по контуру 123—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