Федеральное агентство по образованию

Государственное образовательное учреждение

высшего профессионального образования

“Якутский государственный университет им. М.К. Аммосова”

Технический институт (филиал) в г. Нерюнгри

Контрольная работа № 1

на тему: «Определение геотермии горного массива»

Вариант 5

Выполнил: ст. гр. ТиТР - 06

Денисов Д.С.

Проверил: преподаватель

Скоморошко Ю.Н.

Нерюнгри 2008г.

Определение геотермии горного массива

Цель занятия – построить температурный профиль горного массива по глубине (в гелиотермозоне, криолитозоне) и оценить мощность СТС, а также мощность распространения вечномерзлых горных пород.

Теоретические положения:

Температурное поле верхней части земной коры определяется взаимодействием внутренних и внешних источников тепла. Внутренние источники тепла относительно стабильны, т.к. связаны с постоянно действующими факторами (радиоактивный распад, гравитационная дифференциация вещества и т.д.). Эти источники вызывают повышение температуры пород с глубиной. Внешние источники (основным из которых является переменная во времени солнечная радиация) вызывают периодические температурные колебания горного массива, затухающие на определенной глубине от поверхности Н₀, называемой глубиной гелиотермозоны или глубиной нейтрального слоя.

Температурный режим поверхности Земли в конкретном районе определяется как:

, (1.1)

где т_ср – среднегодовая температура почвы, С;

т_ср=t_ср+2;

t_ср – среднегодовая температура воздуха, С;

А_т – амплитуда колебаний температуры почвы,С;

А_т А_t – 2,5;

А_т – амплитуда колебаний воздуха, С;

 - время, изменяется от 0 до 8760 – продолжительность года в часах,

Для полуограниченного массива амплитуда годовых колебаний температуры пород на глубине Н определяется по известной формуле:

А_т (Н)=А_техр С, (1.2)

где а - коэффициент температуропроводности, м²/ч;

а=3600 ;

_ коэффициент теплопроводности пород, Вт/(мК);

с – удельная теплоемкость пород, Дж/(кгК);

 – плотность пород, кг/м³.

Запаздывание колебаний температуры пород по отношению к изменениям температуры воздуха для полуограниченного массива имеет вид:

_тч

Тогда изменение температуры пород в пределах гелиотермозоны с учетом зависимостей (1.2), (1.3) приблизительно описывается уравнением:

т_нт_срА_техр

Глубину гелиотермозоны можно определить из выражения (1.2)

Н_о=. (1.5)

где А_т(Н_о) – амплитуда пород на глубине Н_о, для расчетов можно принять А_тН_оС.

Изменение температуры пород при углублении на 1м называется геотермическим градиентом q_г (G). Тепловой поток в недрах Земли q связан с геотермическим градиентом законом Фурье

_п_г

Знак минус в формуле говорит о том, что вектор геотермического градиента направлен сверху вниз (в сторону увеличения температуры), а тепловой поток – снизу вверх (направление теплопередачи).

Поэтому, геотермический градиент можно определить следующим образом

_г_п

Средний удельный тепловой поток из недр Земли к ее поверхности составляет ^Втм^

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

ТТ_о_гНН_о

где Т_о- температура пород на глубине нейтрального слоя Н_о и вычисляется по формуле (1.4).

Исходные данные

Теплофизические свойства пород

Таблица 1.1

Наименование	λπ ,Вт/(мК)	С*103, Дж/(кгК)	ρ,кг/м3
Алевролит	1,9	0,83	2540
Гранит	3,5	0,67	2600
Гипс	1,1	1,05	2320
Глина	1,4	0,78	1900
Кварц	2,7	0,96	2500
Песчаник	2,9	0,82	2300
Сланец глинистый	1,75	0,75	2000

Данные по варианту

Показатели	Вариант 5
tср, 0С	-8
At, 0C	19
Мощность пород, м:
Глина	20
Алевролит
Глин. Сланец	20
Песчаник	180
Кварцит
Гипс	220

Порядок проведения работы

1. Рассчитать по формуле (1.1) и построить график изменения текущей температуры поверхности по заданным t_ср и А_t в функции времени на период один год.

время	0	730	1460	2190	2920	3650	4380	5110	5840	6570	7300	8030
T(τ)	-6	2,3	8,3	10,5	8,3	2,2	-6,0	-14,3	-20,3	-22,5	-20,3	-14,3

2. Вычислить годовое изменение температуры пород на разных глубинах (2, 5, 8,10,13,15,17 метров и т.д.) в пределах гелиотермозоны по формуле (1.4).

	Месяц	Время, час	Глубина, м
			0	2	5	8	10	13	15	17
3	март	0	-6,0	-11,7	-9,2	-6,3	-5,7	-5,8	-5,9	-6,0
4	апрель	730	2,3	-7,6	-9,1	-6,9	-6,1	-5,8	-5,9	-6,0
5	май	1460	8,3	-3,2	-8,2	-7,2	-6,4	-5,9	-5,9	-5,9
6	июнь	2190	10,5	0,5	-6,7	-7,2	-6,6	-6,0	-5,9	-5,9
7

июль	2920	8,3	2,5	-5,0	-6,9	-6,6	-6,1	-6,0	-6,0
8	август	3658	2,2	2,1	-3,6	-6,3	-6,5	-6,2	-6,0	-6,0
9	сентябрь	4380	-6,0	-0,3	-2,8	-5,7	-6,3	-6,2	-6,1	-6,0
10	октябрь	5110	-14,3	-4,4	-2,9	-5,1	-5,9	-6,2	-6,1	-6,0
11	ноябрь	5840	-20,3	-8,8	-3,8	-4,8	-5,6	-6,1	-6,1	-6,1
12	декабрь	6570	-22,5	-12,5	-5,3	-4,8	-5,4	-6,0	-6,1	-6,1
1	январь	7300	-20,3	-14,5	-7,0	-5,1	-5,4	-5,9	-6,0	-6,0
2	февраль	8030	-14,3	-14,2	-8,4	-5,7	-5,5	-5,8	-6,0	-6,0
3	март	8760	-6,0	-11,7	-9,2	-6,3	-5,7	-5,8	-5,9	-6,0

а глины	0,003
Аср	16,5
Тср	-6
H0=	16

	dh	qг	a
глина	10	0,050	0,0034
глин. сл.	30	0,040	0,0042
песчаник	20	0,024	0,0055
гипс	0	0,063	0,0016
H0=	16