CC BY
36
ФИЗИЧЕСКИЕ ПРОЦЕССЫ ГОРНОГО ПРОИЗВОДСТВА
УДК 622.822.2: 622.271: 519.6
Н.В. Трушникова
ОБРАТНАЯ ЗАДАЧА ЛОКАЦИИ ОЧАГА САМОВОЗГОРАНИЯ УГЛЯ ПО ИЗМЕРЕНИЯМ ПОТЕНЦИАЛА ТЕРМОЭЛЕКТРИЧЕСКОГО ПОЛЯ НА ПОВЕРХНОСТИ УГОЛЬНОГО РАЗРЕЗА
Одним из эффективных технологических процессов добычи углей является открытый способ. Преимущества открытых горных работ по сравнению с подземными разработками заключается в том, что на карьерах обеспечивается достаточно высокая безопасность труда, производительность работ выше, а себестоимость добываемых углей в 2-4 раза ниже.
Однако открытые способы добычи имеют свои недостатки: необходимость отчуждения
больших земельных участков; понижение уровня грунтовых вод на больших площадях; зависимость от климатических условий, что наиболее сильно проявляется в районах сурового климата.
Наряду с этим следует отметить, что при проведении вскрышных работ нарушается структура горной среды (увеличивается трещиноватость и кусковатость породы). Это может приводить к процессам самонагревания и самовозгорания угля в уступах рабочих площадок [1-4]. Для Кузбасса и других угольных регионов, актуальной задачей является определение расположения очагов самонагревания и самовозгорания углей в угольных пластах и рабочих уступах карьерного поля. Из-за отсутствия информации о точном расположении очагов тепловыделения и их геометрической формы затруднительно построить правильную стратегию ликвидации эндогенных пожаров, которые могут продолжаться по нескольку лет, что приводит к экономическим потерям и экологическим проблемам.
Поэтому определение основных параметров очагов самонагревания и самовозгорания в угольных пластах и уступах при отрытых горных разработках является актуальной как с научной, так и с практической точки зрения. Одним из методов решения подобных проблем является теория обратных задач и применение численных алгоритмов нахождения искомых параметров. Решение этих обратных задач следует осуществлять последовательно: первоначально определяется «геометрия» (форма) области самонагревания, затем восстанавливается распределение температуры (или градиента электрического потенциала) на поверхности очага.
Математическая модель, численные алгорит-
мы и тестовые расчеты по первой задаче рассмотрены работах [5,6], для решения второй проблемы может быть использован математический аппарат, предложенный в [7,8].
В настоящей работе рассматривается обратная задача восстановления поверхности объемного теплового источника в массиве по измеренному экспериментально градиенту потенциала термоэлектрического поля на дневной поверхности угольного разреза. Проводится сравнение полученных результатов с данными натурного эксперимента из работы [3]. При изложении результатов используются данные экспериментов по измерениям, как температуры угольной среды, так и потенциала электрического поля [3]. Следует отметить, что в работах [2-4] из теоретического анализа и обработки данных экспериментов получены аналитические зависимости между вариациями температуры и потенциала электрического поля. Это позволяет определять связи между градиентами потенциала электрического поля и распределением температуры на дневной поверхности угольного массива.
Рис. 1. Схематичное изображение очага самонагревания в угольном уступе
На рис. 1. приведено схематичное изображение очага самонагревания на уступе угольного пласта, расположенного на глубине Н от дневной поверхности. Предполагается, что область экспериментальных измерений расположена в квадрате
[-М,М] х [-М,М] , а распределение температуры внутри очага варьируется в пределах от 30° до 95°. Такое предположение о распределении температуры внутри очага является достаточно условным, основано на экспериментальных данных и требует дополнительных исследований (например, решением задач, рассмотренных в [7,8]).
Ь(р(х,у)
Рис. 2. Геометрическая интерпретация задачи в пространстве Я3
На рис. 2. представлена геометрическая иллюстрация решаемой пространственной задачи: источник тепловыделения расположен в области Б с границей ¡(р, 0) , где р 0 - значения соответствующих углов в сферической системе координат (р, 0, ¡) , привязанной к декартовой
системе (£, Т, п); Ар(х,у) - градиент потенциала р, измеряемый на дневной поверхности
(г=0) в области [-М,М] х [-М,М].
Решение обратной задачи сводится к восстановлению поверхности области Б по измеряемому градиенту потенциала электрического поля Ар на дневной поверхности:
др( сo, Уo, Ч)
Лр( х0, у0,0) =
ді
го = о
где Хо, уо - координаты точек измерений. Искомая функция поверхности очага тепловыделения ¡=1(р, 0) определяется из нелинейного
интегрального уравнения первого рода (см. [5,6]): 2п п
У\( х0 ’Уо) = / Р к( х0 'У0 Р’вМ (р,0))йв,
0 0
(1)
где
(Хо, Уо) = -АР( Хо, Уо,0) / к,
к - эффективная плотность электрических зарядов на поверхности Б;
х~>-
К(Х0 ,у0 ,р’в’1(р, в)) = 5ІП в Х
(а2 -Ь2)уІ£2(Р’в)-2Ы(Р’в)+а2
Х
(а 2[НЬ + с(2а 2 - 3Ь 2 )]
+ [Н(а2 -2Ь2)-сЬ(5а2 -6Ь2)]х
х 1(Р’в)+с(а2 -Ь2)12(р’в))+(Н-3Ьс) д/12 (Р’в) - 2Ь1( Р’в)+а2 +1( Р’в) - Ь
Лп
а-Ь
+
а[НЬ + с(2а2 -3Ь2)]} а2 - Ь2 '
а2 = х02 + у02 + Н2 ’
Ь = Х0 ъ1п в свъ р + У0 8ЇП в ъ1п р + Н COS в’ с = ео8в, а >1
Так как входные данные задачи: Ар(хо,уо,0) определяются из эксперимента с погрешностью 8 > 0 и задача некорректна, то для нахождения приближенного решения уравнения (1) целесообразно применение регуляризирующих по Тихонову алгоритмов [5].
Численный алгоритм решения нелинейного уравнения (1) и модельные расчеты, иллюстрирующие его эффективность, приведены в [6].
Результаты обработки реального эксперимента измерения температуры стенок скважин
Глубина скважины, М Номера скважин
1 2 3 4 5 6 7
Температу] ра Т, °С
0.25 15 20 23 40 25 15 15
0.50 20 20 20 40 28 13 13
1.00 19 18 18 93 32 10 8
1.50 15 30 21 48 60 9 6
2.00 8 21 20 30 70 9.5 6
2.50 18 20 18 18 44 10 6
3.00 19 18 18 18 23 10.5 6
Важным вопросом является проверка адекватности построенных математических моделей, а также их анализ и сравнение получаемых результатов при фиксированной погрешности измерений 8 Основным критерием состоятельности
1
Х
Рис. 3. Аппроксимация границы очага тепловыделения по данным
эксперимента (скважины с №№ 1 — 7)
Др(лг,0), В!М
о
-0.5
-1.5
-2
-2.5
-3
-6
-4
-2
Рис. 4. Распределение градиента потенциала Ар(Х,0) и результаты восстановления границы теплового источника 2(£) для пласта "Мощный " (профиль 35) .
моделей может быть проведение натурного эксперимента, например, бурение контрольных скважин, а также согласованность численных результатов с данными эксперимента.
На стенках скважин производится измерение температуры горного массива по глубине скважины, а также регистрация основных компонент газовой фазы: 02, С02 и СО, выделяемых при окислительно-восстановительных процессах в угольной массе [1-3].
Результаты обработки реального эксперимента измерения температуры стенок скважин, проведенного на пласте «Мощный» (профиль 35), при наличии теплового источника в угольном массиве [3] приведены в таблице.
Анализ экспериментальных данных показывает, что сечение очага тепловыделения (Уо=0) расположено между скважинами № 3 и № 6, на
глубине от о.5 М до 2.5 М . Температура атмосферы (воздуха) То над дневной поверхностью считалась постоянной и, приблизительно, равнялась 15°С .
Построим приближенно границу источника тепловыделения, исходя из следующих соображений. В точки границы будем включать значения координат (^, 2), если температура в этих узлах более 30°С.
На рис. 3 представлена "аппроксимация" границы источника тепло- выделения - четырехугольник, построенный на основе данных из таблицы.
Используя формулу связи между вариациями потенциала электрического поля 8р(Х) и изменениями поля температуры 8Т (Х) в массиве [4], были найдены приближенные значения потенциа-
ла ~(хі ,0) и градиента потенциала Лр( Х,0).
Приведем результаты первичной обработки этого эксперимента, методом наименьших квадратов, используя аппроксимацию данных многочленом четвертой степени.
Пусть массивы х1, у1 - координаты скважин и значения потенциала р(хі ,0), соответственно, У- аналитическая зависимость при аппроксимации градиента потенциала Л р ( Хі, 0).
Имеем
У = 0.201515Х4 + 0.066667Х3 -
-10.648485Х2 -2.780952Х +
+ 123.982664 мВ/м;
Х1 = [-6 - 4 - 2 0 2 4 6] м;
у1 = [0 32 51.2 160 0 0] мВ/м
При этом относительная среднеквадратическая погрешность аппроксимации экспериментальных данных Я1 составляла приблизительно 32% .
Результаты численных расчетов решения нелинейного интегрального уравнения (1) представлены на рис. 4.
Сравнение полученных результатов (рис. 3 и 4) показывают удовлетворительное согласование "восстановленной" границы теплового источника х(^) с данными эксперимента.
Это может служить подтверждением адекватности математической модели, а также эффективности методики для определения границы очага тепловыделения в угольном массиве, разработанной в работах [5,6].
Оценка границы очагов тепловых источников может быть использована на практике: для оперативного прогнозирования развития (и контроля) процессов самонагревания и самовозгорания, а также их ликвидации (тушения) при минимальных экономических затратах (например, бурение минимального количества скважин для подачи антипирогенов).
СПИСОК ЛИТЕРАТУРЫ
1. Линденау, Н.И. Происхождение, профилактика и тушение эндогенных пожаров в угольных шахтах / Н.И. Линденау, В.М. Маевская, В.Ф. Крылов.- М.: Недра, 1977. - 319 с.
2. Физический контроль массивов горных пород / Б.Г. Тарасов, В.В. Иванов, В.В.Дырдин, А.Н. Фокин. - М.: Недра, 1994. - 240 с.
3. Кроль, Г.В. Разработка электрометрического способа контроля самонагревания и самовозгорания каменного угля на разрезах: дис. ... канд. техн. наук: 05.26.01: защищена 16.12.83: утв. 16.05.84 / Кроль Георгий Васильевич; ВостНИИ. - Кемерово, 1983. -195 с. - Библиогр.102 назв.
4. Алексеев, Д. В. Механизм формирования квазистационарного электрического поля в массиве горных пород при наличии естественных и техногенных тепловых источников / Д.В. Алексеев, П.В. Егоров // ФТПРПИ. -1994. - № 5. - С. 3 - 7.
5. Трушникова, Н.В. Об определении формы объемного теплового источника в массиве горных пород по измерениям электрического потенциала на поверхности / Н.В. Трушникова // Геодинамика и напряженное состояние недр Земли. Тр. междунар. конф. - Новосибирск. - ИГД СО РАН. - 2001. - С. 165-167.
6. Трушникова, Н.В. Определение некоторых параметров объемных тепловых источников в массиве горных пород по распределению градиента электрического потенциала на поверхности / Н.В.
Трушникова // Вестн. Кузбасского гос. тех. унив., 2006, № 6. - С. 3 - 7.
7. Berntsson, F. Numerical solution of a Cauchy problem for the Laplace equation / F. Berntsson, L. Elden // Inverse Problems. -2001. - № 17. - Р. 839 - 853.
8. Kress, R. Nonlinear integral equations and the iterative solution for an inverse boundary value problem / R. Kress, W. Rundell // Inverse Problems. -2005. - № 21. - Р. 1207 - 1223.
□ Автор статьи:
Трушникова Надежда Васильевна - ст. преп. каф. математики КузГТУ.
Тел 8-3842-39-63-19.
8-3842-53-05-63
