Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Ониани, Ш. И. Тепловой режим глубоких шахт при гидравлической закладке выработанного пространства и сложном рельефе поверхности

.pdf
Скачиваний:
2
Добавлен:
22.10.2023
Размер:
12.01 Mб
Скачать

ГЛАВА V

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ЗАДАЧИ ПОСТРОЕНИЯ ТЕМПЕРАТУРНОГО ПОЛЯ СИСТЕМЫ УГОЛЬ—ЗА­ КЛАДКА—ПОРОДА ПРИ ЭКСПОНЕНЦИАЛЬНОМ НАЧАЛЬНОМ Р А С П Р Е Д Е Л Е Н И И ТЕМПЕРАТУРЫ В [БОКОВЫХ ТЕЛАХ

;§'1. Постановка задачи

Как было отмечено выше, начальная температура заклад­ ки не превышает зимой (январь, февраль) 6, а летом (июль, август) 20°С. Естественная температура горных пород на пло­ скости -контакта угольной толщи с нижними песчаниками для шахт «Западная-2» и «Комсомольская» достигает 37—38°С, а для запроектированной ша.хты «Шаори № 1» она доходит до 48—49°С.

Таким образом, температурный напор от неохлажденногоугольного и народного массивов к закладке на шахтах «Ком­ сомольская» и «Западная-2» может достичь в зимний период 32, а в летний — 18 град. На нижних горизонтах шахты «Ша­ ори № 1» температурный перепад между закладкой и окру­ жающей чредой существенно увеличится. Это обстоятель­ ство при управлении кровлей полной гидравлической заклад ­ кой выработанного пространства является очень важным для рассматривааМ'Ого месторождения.

В результате столь значительной начальной разности температур между закладкой и прилегающими к ней масси­ вами заложенный материал окажет существенное влияние на распределение температуры вокруг очистной выработки со­ седнего с закладкой угольного слоя.

Определенное охлаждающее влияние на тепловые атмос­ ферные условия очистной выработки гидравлическая заклад­ ка М'Ожет оказать и при выемке первого слоя. Это прежде все-

120

mo .-происходит потому, что боковые стенки нейтрального брем­ сберга представлены заложенным материалом, а в очистной выработке (при выполнении первого цикла отбойки угля) за­ кладка занимает примерно 25% теплоотдающей поверхности лавы (рис. 10). Кроме того, .некоторое количество тепла от вентиляционной струи передается отработанной воде заклад­ ки через стенки трубопровода, проложенного по откаточным выработкам. Все эти факторы, несомненно влияющие на теп­ ловой режим лавы, при выемке первой горизонтальной поло­ сы первого на;кло.нно1по слоя отсутствуют. Затем они вступа­ ют в силу и с увеличением .количества отработанных полос их роль в формировании микроклимата очистной выработки по­ степенно возр астаег.

Выше было сказано, что ткибулыские угли характеризу­ ются высокой газоносностью и газопроницаемостью. В ре­ зультате этого при выемке первого слоя и 'особенно первой горизонтальной полосы угольной толщи, интенсивная дега­ зация угольного массива и частичная — вмещающих пород, предопределяет характер распределения температуры вокруг очистной выработки. С увеличением глубины разработки изза повышения газоносности пластов будет иметь место рост степени охлаждения и толщины охлажденной зоны угля и по­ роды около поверхности обнажения, вызванный охлаждаю­ щим эффектом десорбции и расширения газа. Поэтому вли­ яние гидравлической закладки на формирование теплового режима лавы в данном случае имеет подчиненное значение.

Таким образом, при выемке первого слоя гидравлическая закладка воздействует на температуру вентиляционной струидо очистного забоя, но в лаве ее влияние практически отсут­ ствует. Температурное поле вокруг очистной выработки и тем­ пература обнаженной поверхности полностью определяются'

естественным тепловым

состоянием

массива и условиями и:

и нтеноивностьто дегаз ации.

 

Совершенно иное положение следует ожидать три выемке-

второго, прилегающего

к закладке,

слоя. Второй слой угля

в значительной .степени дегазируется при выемке первого. По­ этому при его отработке дегазация не может играть сущест­ венной роли в перераспределении температуры в, массиве

121

•-около плоокости обнажения в очистной .выработке. Однако, закладка .в .выработанном .пространстве первого слоя со зна­ чительно низкой начальной температурой передает большую часть аккумулированного холода прилегающей среде и к моменту отработки второго слоя (т. е. спустя 4—6 месяцев) создает совершенно новое распределение температуры вокруг очистной выработки, существенно снижая температуру по­ верхности обнажения в лаве.

•Охлажденная зона, обусловленная теплообменом с руд­ ничным воздухом, в очистных выработках настолько мала, •что практически не оказывает влияния на интенсивность пе­ реноса тепла от гарного массива к воздуху. Поэтому в ла­ вах глубоких шахт повсеместно наблюдается максимальное приращение температуры рудничного воздуха, достигающее нескольких (5—8) градусов. Без знания распределения тем­ пературы вокруг очистной выработки немыслим тепловой ра­

счет и,

следовательно, решение задачи

прогноза и регул иро-

• зания

теплового режима глубоких шахт с необходимой точ­

ностью.

Таким образом, при слоевой

разработке мощных

угольных пластов глубоких шахт и полной гидравлической закладке выработанного пространства решение вопроса по­ строения температурного ноля системы уголь—закладка—по­ рода является жизненно важной и актуальной задачей.

•Строго говоря, поставленная задача трехмерная, так как аккумулированный в закладочном материале холод переда­

ется

не

только

углю и породе,

расположенным

соответст­

венно

в

кровле

и почве закладки,

но

и в какой-то степени

угольному массиву через торцевые поверхности

раздела

(по

нижней и одной из

баковым границам выемочного блока,

рис.

10). Неамотря

на

это, без снижения

точности

инженерных

расчетов теплообмен между закладкой и прилегающими мас­ сивами можно расематрмюагсь как одномерную задачу по .сле­ дующим соображениям. Во-первых, теплоот.да.ющ.ая поверх­ ность по границам выемочного блока в (несколько десятков раз меньше по сравнению с твплоотдающими поверхностями в кровле и почве, н, во-вторых, при закладке второй и .каж­

дой

последующей горизонтальной полосы

отработанная во­

да

закладки попадает в перфорированный

трубопровод после

122

инфильтрации ранее заложенного материала и отвода зна­ чительного .количества тепла от последнего. Очевидно, унос тепла этой водой с избытком будет .комиенсировать приток

тепла через боковые границы

выемочного блока.

Фильтра­

ции

воды низкой температуры

через

старую закладку и

вдоль

торцевых поверхностей (раздела

практически

исключа­

ет возможность более быстрого повышения температуры на

границе

блока по

сравнению

с его центральной

частью.

Шаг по

времени между

закладкой соседних полос

настоль­

ко незначительный

(3—4

сутки), что температура

той

или

виной полосы практически

не отличается от температуры

со­

седних.

 

 

 

 

 

 

 

Таким образом, размеры горизонтальных полос .и выемоч­

ного блока

и технология

закладки дают право считать, что

на границе

блока

потери

холода

отсутствуют, т. е. закладку

можно рассматривать .как неограниченную пластину, поме­ щенную между угольным и породным массивами. Более того, инфильтрация отработанной воды должна сдерживать и на­

зревание закладки за счет тепла

угольного массива

в кровле

и породного

в почве.

 

.

В главе

IV

было показано,

что распределение

темпера­

туры в угольном и породном массивах перед закладкой вы­ работанного пространства первого слоя принимает экспонен­

циальный

характер, начиная от

поверхности обнажения

до

некоторой

глубины

х'у_п.

Затем

(для больших расстояний

от поверхности обнажения)

распределение

температуры

в

массиве сохраняет

естественный

(линейный)

характер.

 

Температура заложенного материала сразу после окон­ чания закладочных работ по всей толщине одинакова.

Рассмотрим задачу аналитического построения темпера­ турного поля вынимаемого второго слоя при заложенном вы­ работанном пространстве первого и восходящей последова­ тельности отработки слоев без деления угольной толщи на отдельные пачки.

Принимая во внимание вышеизложенное, постановку за­ дачи М'ожно сформулировать следующим образом. Дана не­ ограниченная пластина с равномерной начальной температу-

123

I

рой

Ц, состоящая из

заложенного

материала с

теплофизи-

ческими свойствами а3,

Х3 , с 3 . В начальный момент времени

она

помещается между

массивами

угля и породы

(лолуопра-

ниченные тела), температура которых, начиная от поварьхности раздела, с увеличением координаты повышается экспонен­ циально и по величине на 16 — 30 град выше /„. Теплофизичешие свойства угля, закладки и породы существенно .различ­

ны. Необходимо найти распределение температуры

в системе

в любой заранее заданный момент времени т.

 

•Подобные задачи в теории теплопроводности получили

название краевых задач с граничными условиями

четвертого

рода или задач с разрывными коэффициентами, ибо теплофизичеокие свойства составляющих тел являются кусочнопостоянными функциями координаты х.

Аналитическому построению температурного поля сос­ тавным неоднородных тел теория теплопроводности придает большое значение. Известно несколько аналитических реше­ ний. Но в одних задача симметрична с равномерным началь­ ным распределением температуры во всех телах [91], другие рассматривают симметричную систему трех неограниченных

пластин с одинаковой

и равномерной начальной температу­

рой при теплообмене

между поверхностями боковых пластин

и окружающей

средой по закону свободной конвекции [127].

В специальных

монографиях приводятся общие решения для

многослойной системы при теплообмене на свободных поверх­ ностях по закону конвекции и граничных условиях первого,

второго и третьего рода

[92]. В работе [64] рассматривается

нестационарная

теплопроводность неоднородных слоистых

тел (пластина,

цилиндр,

шар), вызванная действием (распре­

деленных источников внутри и на границе тел, опять-таки при теплообмене на .свободных поверхностях системы по закону конвекции. В монографии [61] даются решения для полуограниченных и конечных составных тел при равномерном на­ чальном распределении температуры и разных граничных ус­ ловиях. Таким образом, известные аналитические решения для рассматриваемой задачи неприменимы.

Уголь и порода (нижние песчаники) вследствие разных теплофизических свойств будут оказывать существенно раз-

124

личное влияние на распределение температуры в пластине. Но внутри пластины всегда будет существовать некоторая параллельная ограничивающим поверхностям плоскость, на /которую тепловое воздействие с обеих сторон будет совер-

 

 

*

*

 

 

 

 

УЩ

 

 

 

•0.'

•У,:::.'.

 

 

 

 

 

 

У : . : :

 

- \ у

;

 

 

Порода

 

 

'. Л„,а.пп

•. .

 

 

о

 

 

 

 

 

 

 

 

 

 

о"

 

 

 

 

 

о

 

 

 

 

 

 

«

 

 

« ,

. *

'

 

 

о

• О

 

 

 

 

 

 

 

' о

 

 

 

 

 

 

• о

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 4 5 .

Система

трех тел

при граничных условиях четвертого рода

шанню одинаковым. Расстояния от этой плоскости

(назовем

ее нейтральной)

до поверхностей

раздела

различны

(несим­

метричная

задача).

 

 

 

 

 

 

 

 

Расположим

начало

координат

на

нейтральной

плоскости.

Расстояние

от нее

до угольного массива /х ,

до песчаников /2 ,

причем / =

tx-\-/2—толщина

пластины (рис.

45). Тогда

матема­

тическая

запись

поставленной

задачи

принимает

следующий

вид [ Ш ] :

 

 

 

 

 

 

 

 

 

 

125

для системы'тел уголь—верхний слой закладки толщиной' lv

д t3 (х,

х)

д%

(х, х)

 

 

 

 

 

 

 

(5.1).

 

дх

=

а,

 

 

(т >

0,

0 <

х <

 

 

 

 

 

д х2

 

 

 

 

 

 

 

 

 

д ty(x,

х)

дУу(х,

х)

( т > 0 , / 1 < л < о о ) ;

 

(5.2)-

 

дх

 

 

дх*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ир а евые

условия

 

 

 

 

 

 

 

 

 

 

 

/э(х,

0) = If, ty(x, 0) =

t l -

 

Ъх

 

 

,r°il*-ix)-?y(x-lj,

 

(5.3) '

 

 

 

1

dty(lltx)

 

 

 

 

 

 

(5.4) »

 

 

 

1 L

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.5)

(3(lv

x) = ty(lv

x),

 

 

 

 

 

 

 

 

 

dt3(0,x)

 

dty(<x>,T)

_

 

 

 

 

 

(5.6);

 

 

 

Ol

J ..

-

1

*• yi

 

 

 

для системы тел порода — нижний слой

закладки

толщиной /,

д ts (х,

х)

д2t3

(х,

х)

( т > 0 , 0 < | х | < / 8 ) ,

 

(5-7)

 

дх

= а,

дх2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д t„ (х, х)

 

 

 

 

( т > 0 ,

/ 2

< Ы < оо);

(5.8)

 

 

 

 

дх2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

краевые

условия

 

 

 

 

 

 

 

 

 

 

 

t n ( - x ,

0) = *s _ вя в-в »(и-'Я )

+ г п д о

-

д ,

 

 

(5.9)

^ з ( ^ , т )

1

 

 

 

 

 

 

 

 

 

(5.10)'

 

дх

 

к .

 

 

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t) = /„(/„, т),

 

 

 

 

 

 

 

 

 

(5.11)

dt3(0,x)

 

df„(-oo,T)

 

 

 

 

 

 

где ^ и ^—естественная

невозмущенная начальная

температу­

 

 

ра на поверхности

раздела (х = 1г и —л; = /2 ) уголь­

 

 

ного и породного

массивов

соответственно;

Гу , Г„ —естественный

(геотермический)

градиент

изменения

 

 

температуры

с глубиной в массиве угля

 

и породы

 

 

соответственно (в данном

случае для песчаников Г

• 126

величина положительная, для углей—отрицатель­ ная);

bv 52, ох , а2 —некоторые постоянные эмпирические величины, за­ висящие от начального распределения температуры в охлажденной зоне угольного и породного л-асси-. вов соответственно;

 

 

Лу

 

Л п

 

 

 

 

Строго говоря, при х-+со

геотермический

градиент не может-

быть равным

Г у

или

Г п ,

ибо толщина

песчаников

(нес­

колько сот метров) и тем более мощность

 

угольной

толщи

(околю 30 :м)

величины

конечные. Но в

данном случае

это

не .имеет значения,

так

как .возмущающее

 

влияние закладки

на температурное поле недр не может распространиться

на

(бесконечное расстояние

и при конечных

значениях

времени

(до 5 лет)

оно не будет проникать дальше 20—30 м. Поэтому

вполне допустимо

считать,

что к неограниченной

пластине-

закладки

с одной

стороны

прилегает полуограниченное тело

из угля, а

с другой — такое же тело, состоящее из

песчани­

ков. Такому физическому представлению задачи и соответст­

вуют граничные условия (5.6)

и (5.12) при х - э - о э .

§ 2. решение

задачи

Существуют три основных метода решения дифференци­ альных уравнений нестационарной теплопроводности при из­ вестных краевых условиях. Первым из них следует считать наиболее распространенный классический метод, разработан­ ный основателем аналитической теории теплопроводности Фурье. Главной особенностью этого метода является пред­ положение о том, что решение данного уравнения может быть представлено в виде произведения двух функций, одна из ко­ торых зависит только от «ременной независимой переменной, а другая — только от пространственной независимой пере­ меннойРешение получается в виде бесконечного ряда. При •точных 'расчетах :беретоя бесконечное число членов ряда, а при приближенных — сумма конечного числа членов. Этим методом получены удовлетворительные решения многих тех-

127

. шичеших

задач,

касающихся

в основном

тел простейших

'форм. До

40-х годов при решении

задач

нестационарной

теплопроводности

преимущественно

пользовались

методом

Фурье. В

горной

теплофизике

классический

метод

успешно

• применил А.

Ф. Воропаев для

описания

процесса

тепло- и

маооообмена

между

горным

массивом

и вентиляционной

струей, протекающей

по выработкам вентиляционной сети

[15].

 

 

 

 

 

Метод операционного .исчисления дли решения дифферен­

циальных уравнений в частных производных второго

порядка

в математической физике широкое распространение

получил

в последние годы- А. В. Лыков в своих работах изложил ос­ новные принципы практического приложения метода преоб­

разования

Лапласа" и успешно применил

его для

решения из­

вестных

и

новых

задач

 

нестационарной

тепло- и масеопро-

водности

[91, 92,

93]. В

теории теплопроводности

этим мето­

дом

широко пользуются

Г. Карелоу и Д. Егер [61]. В обла­

сти

гарной

теплофизики

примером

успешного

применения

метода операционного

исчисления

служат известные работы

А. Н. Щербаня и

О. А.

Крайнева [147,

148, 152].

 

Несмотря на свою распространенность, два точных ана­ литических метода охватывают далеко не все технически интересные задачи нестационарной теплопроводности. Кроме того, точные аналитические решения иногда настолько .слож­ ны и громоздки, что теряют практическое значение. Вследст­ вие этого в последнее время наблюдается стремление к раз­ работке и применению приближенных методов решения крае­ вых задач теплопроводности. В большинстве случаев такие решения при достаточной точности приближения просты и удобны для ведения инженерных расчетов. Среди приближен­ ных методов наибольшее (распространение получил метод ко­ нечных разностей, или так называемый метод сеток, сущность которого заключается в разбивке пространственной перемен­ ной на отдельные отрезки, а временной переменной на от­ дельные периоды. Искомая температура определяется для каждого отрезка и периода времени последовательно. Даль­ нейшим развитием метода конечных разностей является метод элементарных балансов, примененный для решения прост-.

128

ранственных нестационарных -задач теплопроводности яри гра­ ничных условиях первого, второго и третьего рода [9]. Для решения многих технических задач нестационарной теплопро­ водности с граничными условиями всех родов используется и •метод исключения дифференциального уравнения теплового баланса одной из нескольких независимых пространственных переменных. При этом действительная температурная кривая в теле заменяется параболой n-vo порядна [10]. В области гор­ ной теплофизики приближенный метод (метод конечных раз­ ностей) для решения нестационарных задач то прогнозу и регулированию теплового режима глубоких шахт был при­ менен А. Ф. Воропаевым [14].

Классический .метод точного аналитического решения диф­ ференциального уравнения теплопроводности и метод конеч­ ных разностей не дают возможности получить удовлетвори­ тельное решение (поставленной перед нами задачи, а прибли­ женный метод А. И. Вейника [10] приводит к более сложным зависимостям, нежели метод точного аналитического реше­ ния с применением интегрального преобразования Лапласа. Поэтому для точного аналитического решения задачи неста­

ционарного теплообмена между

закладкой и прилегающими

к ней массивами угля и породы

мы пользуемся методам 'Опе­

рационного исчисления, который согласно современной тео­ рии теплопроводности [61, 91], является наиболее удобным методом для решения задач нестационарной теплопроводно­

сти

в

составных

телах-

 

 

 

 

Применим преобразование

Лапласа

к уравнению

(5.1):

 

 

L

 

 

 

L

д \ (х,

х)

(5.13)

 

 

 

дх

 

дх2

 

 

 

 

 

 

 

 

 

 

 

L [t (х,

т)] =

j

t(x, x)e~Pxdx

= t(x,

p),

(5.14)

 

 

 

 

 

0

 

 

 

 

 

 

L

Ma (x,

x)

=

pt3(x,p)-f{0),

(5.15)

 

 

 

дх

 

 

 

 

 

 

 

 

 

 

где

/(0) = fu

 

согласно

уравнению (5.3).

 

 

 

•С учетом (5.15) выражение (5.13) перепишется следую­

щим

образам:

 

 

 

 

 

9. Ш. Ониани

129

Соседние файлы в папке книги из ГПНТБ