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

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

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

так как

'fQ "— — —

— .

==

1

'

2\f

ClyX

Т-э-оо

'

Таким образом, когда х ->• со температура любой точки под­ системы уголь—-закладка толщиной 1г становится равной первона­ чальной температуре tl Ту(х— /х ), а температура любой точки

подсистемы порода—закладка—температуре tn0

- f Гп (]х| — /.,). Ес­

тественное температурное поле любого месторождения

таково, что

всегда можно принять tl«

Следовательно,

мы

вправе

сде­

лать следующее заключение: температурное поле системы,

воз­

мущенное энергией холода

заложенного материала

с равномер­

ной начальной температурой tl, после истечения бесконечно дол­ гого времени t - v c o полностью восстанавливает свое первона­ чальное состояние полуограниченных тел и распределение темпе­ ратуры в них принимает естественный линейный характер.

Г Л А В А V I

АНАЛИТИЧЕСКОЕ ИССЛЕДОВАНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ СИСТЕМЫ

УГОЛЬ—ЗАКЛАДКА—ПОРОДА ПРИ ЛИНЕЙНО М И ПАРАБОЛИЧЕСКОМ НАЧАЛЬНОМ РАСПРЕДЕЛЕНИИ ТЕМПЕРАТУРЫ

,В БОКОВЫХ ТЕЛАХ

§ 1. Начальное распределение температуры — линейная функция координаты

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

случае

принимают

вид:

 

для

угол иного

массива

 

 

 

ф,0)=Ц-Ту(х-11);

(6.1)

для

породы

 

 

 

 

М - * . °) = + Г В ( М -1г),

(6.2)

т. е. о1 1 2 = 0 и сг1 ! 2 = 0.

Таким образом, если в решения (5.60) •— (5.63) подста­

вим нулевые значения постоянных

б и а, получим

новые ре­

шения,

описывающие температурное поле системы

при ли­

нейном

начальном распределении

температуры

в

угольном

и породном массивах [112].

 

со

 

Цх, т ) - П + j t l - t l )

^ ( - / ч ) " 1 - 1

X

1

m=l

151

X

 

erfc (2m-\)l1-x

 

 

erfc

( 2 m - 1

 

 

 

 

(6.3)

 

 

 

2]/a~x

 

 

 

 

 

 

2Va3x

 

 

 

 

 

 

ty(x,

t) =

 

tl-Tyix-lJ-W-Q

 

 

 

erfc

x—lx

 

 

 

 

 

 

2/a y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 +

^

У ^ - '

Ч

Г

x—l1+2??Hlx1

 

 

 

 

 

 

-

1 er/c

 

 

 

(6.4>

 

 

 

 

 

 

 

 

 

 

 

 

2Уау с

 

 

 

 

 

 

Ц-х,

x)=

« + ( « - « ) •

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

erfc

(2m—1)/3 -|JC|

- j

-

er/c ( 2 m - l t t . + | * |

 

(6.5)

*„(-*,

т) =

П - f Г п (|л:| - / . 2 ) - (^ - /2)

y q : tf8

 

 

2"|/ftT

 

 

 

 

 

 

GO

 

 

 

 

 

 

/а

+

2 т / а х а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.6)

 

 

 

 

 

 

 

 

 

 

 

 

2

yagt

 

 

 

 

 

 

 

m=l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

теперь допустим, что боковые тела

рассматривае­

мой

системы

имеют

одинаковые

 

теплофизические

свойства

(т. е. если они состоят из одного

и того же вещества) и

рав­

номерное

начальное

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

температуры (т. е. Г у

=

Гп

= 0),

 

то

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

ре­

шение для системы трех тел три

граничных

условиях

чет­

вертого

рода

 

(симметричная

 

задача) [91]. Это

свидетель­

ствует

о

том, что

известное

 

из

теории

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

решение

и решения

(6.3)

 

(6.6)

являются

частными

слу­

чаями

более

общей

задачи

рассмотренной

нами в

предыду­

щей главе.

Очевидно, решения (6.3) — (6.6) аналогично решениям; (5.60) — (5.63) полностью удовлетворяют исходным диффе­ ренциальным уравнениям и краевым условиям. Влияние на­ чального распределения температуры во всей системе (речь вдет о начальной разности температур неограниченной пла-

152

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

Решения (6.3) — (6.6) при известных значениях 1\ и U. позволяют произвести построение температурного поля всей

рассматриваемой системы в

любой

момент

времени.

Темпе­

ратура заложенного материала

на

границе

раздела

зон

1\

и 1% т. е. на нейтральной

плоскости, с обеих сторон оди­

накова в любой .момент времени

 

одно

и тоже

время в

любой точке

системы

двух

разных

 

значений

температуры-

быть не может). Поэтому, воспользовавшись решениями

(6.3)

и

(6.6), в

силу равенства (5.64)

получаем уравнение

для

оп­

ределения

зон /] и h теплового влияния угольного и породного-

массивов

(то же выражение можно

получить

из уравнения-

(5.65) путем подстановки б]=0 и сч=0):

 

 

 

 

х-

VI

vfc

(2m-1)/,

 

 

 

 

 

-о,

/

j ( - А х Г - 1

2

V

^

X

1

+KX

 

 

 

m=l

 

 

 

 

 

 

 

(6.7)

 

V

I

 

(2777- 1)/,

~

1 + ^ 2

 

 

 

 

 

 

( - АГ - 1

 

 

 

 

 

(tl-Qjj

 

erfc

2

y

~

 

 

 

 

 

В данном случае вторым уравнением для определения /] и /2 - служит также равенство (5-66).

При естественном распределении температуры в среде

можно допустить, что t\m t"0. Тогда зоны теплового влияния примыкающих к закладке массивов будут зависеть только от их тепловых активностей и времени т. На рис. 46 дается гра­ фик, построенный по результатам решений систем уравнений. (6.7) и (5.66) для разных значений т.

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

153.

м

И».

Рис. 46. График зависимости l l t 2 = fb)

середине

его

толщины (/i = /

2

= 0,5/ — симметричная

зада­

ча). В

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

активность

угля

значительно

меньше тепловой

 

активности

песчаников и

поэ­

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

смещена

в

сторону

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

{l\<h,

рис. 46).

 

 

 

 

Если

допустить,

что тепловая

активность

угля бесконеч­

но мала по сравнению с тепловой

активностью

песчаников,

то в этом случае нейтральная плоскость

совмещается

с по­

верхностью равдела

между

углем

и закладкой

(/i=0,

/2 =

/ ) ,

т. е. неограниченная

пластина будет нагреваться

только

за

счет притока тепла

из породного

массива,

и гидравлическая

закладка не окажет влияния на тепловое состояние угольной толщи.

При z

= 0 функция

ljli

=

f{t)

(выражение (6.7)) не

опре­

делена,

так

как,

согласно

условиям

задачи, она

имеет

смысл

только

при

т > 0 .

При

т - > о э

и t\

=

/J выражение

(6.7)

прини­

мает вид

2

V I

,

(2m-1)/,

л т к х

2 а М

ч )

e r f c ~ W ^ r

 

 

 

(2m-1U,

e r f c

erfc-i2m-l)k

m=l

(6.8)

так как сумма ряда

( - Л Г - 1 » 2

2

m=l

откуда следует, что при очень большой продолжительности

времени и одинаковых начальных температурах

угля и

по­

роды зоны

теплового влияния

одинаковы

и не

завиоят

от

тепловых свойств

.сред, которые

воздействуют на пластину

(/i = / 2 ) . Очевидно,

что причиной

этого является

замедление

процесса теплообмена между закладкой и

прилегающими

массивами

с увеличением времени. При больших

значениях

155

времени (т>1 100000 час) процесс становится «вазистационарным, а влияние тепловых активностей тел — незначи­ тельным. Роль тепловой активности в распределении темпе­ ратуры снижается с уменьшением скорости течения процес­ са.

Из

уравнения

(6.7) и графика зависимости

/ 1 5 3

= Дт)

следует,

что при

одинаковых

невозмущенных температурах

угля и породы функция IJl2

= /(т) в начале

процесса

умень­

шается,

достигает

некоторого минимума, а

затем

'начинает

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

.малой, что инерционность системы теряет влияние на про­ цесс, и /[ 12 оказываются равными величинами. То же са­

мое

следует и из выражения (6.8).

Так, например, при % —

= 100000 час /1 = 1,24 м и /2 =1,26 м

(рис. 46).

 

 

Довольно быстрое перемещение нейтральной плоскости

при

малых

т становится

причиной

некоторой

погрешности

при

определении искомой

температуры, тане

.как приведен­

ные

выше

решения не учитывают

влияние перемещаемости

плоскости отсчета. Поэтому при строгой постановке вопроса для каждого момента времени в решении подставляются средневзвешенные значения Ц и /2 , учитывающие перемещае­ мость плоскости отсчета, а не результаты непосредственного решения уравнений (6.7) и (5.66). Для определения сред­ невзвешенных значений 1\ и /г на графике зависимости / 1 ) 2 = = f{x) (рис. 46) проводятся линии, параллельные оси абсцисс, причем это сюущаствляетея таким образом, чтобы площади между ними и графиком l\,2='f (т) левее и правее точки пере­ сечения для каждого момента времени были бы одинаковы­ ми (площадь АБСА равна площади А'В'С'А' при т—6000 час). Для больших х (т > 20000 час) перемещаемость нейт­ ральной плоскости практически не оказывает влияния на точ­ ность расчета.

При возмущении начального температурного поля среды, на зависимость /г ,2 = /(t) существенное влияние оказывают на­ чальные разности температур на плоскостях .соприкосновения

156

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

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

вынимаемого

(второго) слоя

(от х = /х до х—1х

= 2,5

м). Поэто­

му рассмотрим более подробно динамику

температуры во

времени для

характерных

плоскостей этого

слоя

(характер­

ными плоскостями мы называем поверхности раздела и плос­ кость, разделяющую рассматриваемый слой на две равные части).

Температура на границе соприкосновения угольного массива с закладкой (при x = li), согласно решению (6.4), бу­ дет равна

ty(x, х) = tl-{tl - Q

К,

 

GO

(6.9)

т = 1

Анализ последнего уравнения показывает, что с увеличе­ нием времени функция ty(lv х) непрерывно повышается (так как

с повышением т член erf с (и) увеличивается). Минимальное зна­

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

момент времени (т =

+ 0):

'у&. о й » = я - w - й) т т к Г '

( 6 ' 1 0 )

Обозначим относительную

температуру любой

точки

угольного массива через &у , т. е.

 

 

7Т~р

" ^ У

(о - 1 1 )

Тогда решение (6.4) можно переписать следующим образом:

 

 

 

 

 

 

 

/н =

1

 

 

Л: — /, +

/,х.,

 

 

 

 

Х

2 7 ^ 7

=

 

-

(6.12)

где

 

 

 

 

 

 

 

 

 

3\. =

- г - гЛг - ег'с- —

 

-

'

 

 

 

1 + / С Х

'

2 / а у г

 

 

 

 

со

 

 

 

 

 

 

•9-2= 1

. у —

7 i(—hi)

 

erfc

 

 

о •

1

-г ATj.

 

 

 

 

 

2|/ аух

На рис. 47 представлены кривые изменении во времени: каждого члена выражения (6.12) в случае рассмотрения плос­ кости кровли вынимаемого* слоя —= 2,5 м). Из кривых В.ИДНО, что в начале процесса оба члена правой части урав­ нения (6.12) неуклонно растут. Однако, в начальный момент времени приращение второго члена меньше по сравнению с первым. При т=2500 час приращения становятся равными, а затем второй член обгоняет в приросте первый. Заштрихован­ ная область представляет разность между этими членами и дает представление о характере изменения функции •Э-у = f(x) •при х—7i =2,5 м. В данном случае рассматриваемая функция достигает максимального значения при т=2500 час. Когда

158

х-усо, ггу 'Становится равным нулю, т. е. температурное поле' восстанавливает свое первоначальное состояние.

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

0,!5

1

f

3

1 1

Р

0,05

 

Рис. 47. График зависимости

9 1 1 2 =

F(i)

уменьшается,

достигает

минимума, а

затем

начинает увели­

чиваться. Минимальная

температура на поверхности сопри­

косновения с

закладкой устанавливается сразу, в момент за­

ложения выработанного пространства, и величина ее 'Опреде­

ляется выражением

(6.10). По мере увеличения

координаты

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

растет

и при

х—Zi = 2,5

м составляет

2500

час.

 

 

 

Бели

члены с

erf с

выражения (6.12)

разделить и

умно­

жить

на величину

1\, получится

уравнение

в критериальной

форме

[91]:

х

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Кг

 

+К)

 

 

х

Ъу =

 

 

 

 

1 +

YFo1

 

 

 

 

 

 

 

 

 

 

-1 -\-2тжх

X erfc

' i

(6.13)

 

2V

Fox

Л59>

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