книги из ГПНТБ / Ониани, Ш. И. Тепловой режим глубоких шахт при гидравлической закладке выработанного пространства и сложном рельефе поверхности
.pdfГЛАВА 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 |