ГОРНЫЕ МАШИНЫ И ОБОРУДОВАНИЕ
УДК 534.6.08
Е.В. Денисова, И.В. Тищенко, А.П. Хмелинин, Я.З. Бадмаева
О ВОЗМОЖНОСТИ КОНТРОЛЯ ПРОЦЕССА РАЗРУШЕНИЯ ПОРОДНОГО МАССИВА БУРОВОЙ МАШИНОЙ ПО ВРЕМЕНАМ ЗАДЕРЖКИ СОЗДАВАЕМОГО ЕЮ АКУСТИЧЕСКОГО СИГНАЛА
Акустические методы, основанные на исследовании параметров создаваемого в процессе движения буровой машины в породном массиве сигнала, представлены в [1-5] V Они подразумевают амплитудный и спектральный анализ принятой от излучателя, расположенного на самой машине, информации, по которой судят о местоположении машины в породном массиве. Ниже предложено использовать известный способ измерения трехмерных координат источника динамического воздействия на породный массив по временам задержки создаваемых при этом на поверхности массива акустических сигналов примени-
Для подтверждения возможности реализации такого способа вычисления координат буровой машины были выполнены следующие экспериментальные исследования.
Экспериментальные исследования проводились на полигоне ИГД СО РАН "Зеленая горка". Вблизи места проведения исследований отсутствуют высоковольтные сети, а так же другие подземные коммуникации. Массив можно принять однородным. Во входной приямок запускалась машина "Тайфун-2" (габаритные размеры: длина 550 мм, диаметр 40 мм, энергия удара 22 Дж) на глубину 0,5 м. Акустоэлектрические преобразова-
ноутбук
т
Рис. 1. Схема эксперимента по измерению времен задержки акустического сигнала, создаваемого движением буровой машины в породном массиве
тельно к решению задачи определения местоположения буровой машины. Этот способ подразумевает использование сигнала, создаваемого машиной в качестве задающего для вычисления относительно него значений времен задержки прихода акустического сигнала на приемные каналы, расположенные на поверхности земли.
1 Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России на 2009-2013 гг.»
тели (1...4), в качестве которых использовались акселерометры КД29, размещались на поверхности земли в вершинах воображаемого квадрата со стороной 1 м, крепление акселерометров осуществлялось при помощи шпильки. Акустоэлектриче-ский преобразователь 5 (акселерометр КД29) с помощью магнита крепился к металлическому штырю, соединяемому с задней гайкой машины. Сигналы, принятые с пяти акселерометров поступали на входы многоканального измерительного комплекса PULSE и далее на ноутбук. Для обра-
ботки экспериментальных данных применялось программное обеспечение Pulse Lab Shop Version 10.1. Из серии записанных в процессе движения машины сигналов для каждого канала выбирают один импульс с одинаковым порядковым номером. Это позволит исключить из анализа те импульсы, которые соответствуют началу и концу процесса бурения, когда давление подаваемого воздуха не будет постоянно и иметь максимальное значение. Для выбранных по каждому из каналов импульсов, высчитывалось время задержки относительно сигнала, пришедшего с самой машины (5-ый акселерометр), которое заносилось в таблицу с данными для соответствующей координаты. Далее с использованием математического аппарата вычисляют три координаты машины (x, y, z).
За исходные данные для вычисления координат машины используют расстояния между аку-
стоэлектрическими преобразователями 1____4:
а - расстояние между акустоэлектрическими преобразователями: АВ=СВ=ЛВ=БС=1 м;
h - глубина запуска равна расстоянию AK=0,5
м;
l - длина машины равна 0,55 м;
Начало координат принято в точке А.
^ - время прихода опорного сигнала на 5-ый акселерометр (для вычислений принимают равным 0);
tА, tB, tc, tD - времена задержек распространения импульсного акустического сигнала в точках Л, Б, С, D соответственно относительно сигнала с машины.
При выполнении расчетов учитывалась дополнительная задержка сигнала с самой машины при прохождении им пути равного сумме длин рабочего органа и металлического стержня, на который крепился акселерометр. Расчетное значение времени прохождения такого пути в металле равно 0,14 мс. Схема эксперимента представлена на рис. 1.
При заглублении буровой машины в грунтовый массив на всю длину l и установке акселерометров 1...4 согласно схеме на рис.1, координаты точки О, где формируется импульсный акустический сигнал, относительно начала координат бу-а
дут равны: x = —, z = h, y = 0. Зная эти координаты, находят значения скоростей Va, Vb, Vc, Vd распространения импульсного акустического сигнала до каждого акселерометра 1...4 из системы уравнений (в общем виде):
\(OA)2 = x2 + y2 + z2 (OB)2 = x2 + (a - y)2 + z2 (OC)2 = (a - x)2 + (a - y)2 + z2 (OD)2 = (a - x)2 + y2 + z2 откуда вычисляют:
^ = %
- X
- °% г‘- °%
Далее при движении работающей буровой машины, когда ее местоположение неизвестно, по известным скоростям Уа, Уб, Ус, Уб распространения упругих волн в грунте и временам Iа, IБ, Iс, ^D задержки их распространения в нем, после ударного воздействия машины на грунтовый массив, происходит перемещение точки О, где формируется импульсный акустический сигнал, в точку О* удара. Далее вычисляют расстояния от точки О* удара до акселерометров 1...4:
° * А - Уа х хА ° * Б - Уб х 1Б ° * С - Ус х ^С ° * Б - уб X гв
и после этого определяют координаты машины из системы уравнений:
\° * А)2 - х2 + у2 + 22 (1)
(° * Б)2 - х2 + (а - у)2 + 22 (2)
(° * С)2 - (а - х)2 + (а - у)2 + 22 (3)
(° * Б)2 - (а - х)2 + у2 + 22 (4)
Из каждых трех уравнений последовательно вычисляют три координаты. Например, из уравнений (1), (2) и (3) можно получить:
x=
У = z
(O * A)2 - (OD)2 + а2 2b
(O * A)2 - (OB)2 + a2 . 2a
= V(O * A)2 - У2 - x2
Аналогичным образом проводят вычисления и для систем уравнений (1,2,4), (2,3,4), (1,3,4). В итоге получают 4 набора координат (х1, у1, ^1)...( х4, у4, г 4), по которым определяют средние арифметические значения трех координат машины:
х, +... + х4
Уі +... + У 4
У = ;
4
Вычисление времен задержек распространения импульсного акустического сигнала до аку-стоэлектрических преобразователей 1... 4 может производиться несколькими методами.
1. Вычисление времен по началу прихода импульса. В этом случае один из импульсных аку-
X =
z, +... + z 4
Таблица 2.Расчетные значения трехмерных координат машины
Величина и направление отклонения машины, м Реальные значения координат, м Расчетные значения скорости распространения акустического сигнала, при обработке по уровню к ■ итах, м/с
к=0,5 к=0,1 к=0,7
х у г х у г х у г х у г
0,1 вправо 0,40 0,00 0,50 0,30 0,02 0,48 0,12 -0,07 - 0,32 0,02 0,48
0,1 влево 0,60 0,02 0,50 0,70 0,05 0,53 - - - 0,68 0,05 0,50
стических сигналов используют как опорный (сигнал с машины) и вычисляют значения времен задержек между моментом времени t°, когда уровень абсолютного значения опорного сигнала превысил пороговое значение, и моментом времени, когда уровень абсолютного значения сигнала соответствующего акустоэлектрического преобразователя 1...4 превысил пороговое значение. Недостатки такого метода заключаются в низкой точности определения времен задержки из-за высокого уровня шума в полосе полезного сигнала и, как следствие, требуется достаточная величина порогового значения. Но при этом данный метод работает даже при сильных искажениях сигнала вследствие дисперсии среды и возможных ее неоднородностей.
2. Вычисление времен задержек распространения импульсного акустического сигнала по уровню 0,5 или 0,7 от своего максимального значения.
Таблица 1.Рассчитанные скорости распростране-
ния акустического сигнала
Обозна чение скорости Расчетные значения скорости распространения акустического сигнала, при обработке по уровню к ■ Umaх:, м/с
к=0,5 к=0,1 к=0,7
Уа 476 511 450
УЬ 324 359 312
Ус 350 408 339
Уё 453 544 436
В этом случае один из сигналов также используют как опорный и вычисляют значения времен задержек между моментом времени, когда уровень абсолютного значения опорного сигнала принимает свое значение 0,5■Umaх или 0,7■Umaх, и моментом времени, когда уровень абсолютного значения сигнала соответствующего акустоэлектри-ческого преобразователя 1...4 принимает то же значение 0,5■Umaх или 0,7итах. На точность данного метода оказывает влияние искажение формы сигнала вследствие дисперсии среды, либо наличия в ней различного типа неоднородностей. Точность по сравнению с первым методом повышается за счет снижения вероятности ложного срабатывания устройства от воздействия шумовой по-
мехи большой амплитуды.
В табл. 1 представлены результаты рассчитанных скоростей распространения акустического
сигнала в массиве, а в табл. 2 — координат машины, на основе практически измеренных времен задержки распространения акустического сигнала в породном массиве.
Экспериментально установлено, что скорость распространения акустического сигнала в породном массиве зависит от расстояния между акусто-электрическим преобразователем и местом приложения ударной нагрузки. Так, при измерениях в точках В и С скорость распространения акустического сигнала, создаваемого машиной ниже, чем при ее измерении в точках А и Б. Это можно объяснить тем, что с увеличением расстояния от места удара до точки, в которой происходят измерения возникает большая вероятность нахождения на пути распространения акустической волны различного типа неоднородностей (например, твердых породных включений). В связи с этим измерения координат машины целесообразно производить, когда машина находится примерно в центре между акустоэлектрическими преобразователями, это позволит избежать погрешностей при последующем вычислении ее координат.
Анализ рассчитанных координат машины по измеренным временам задержки акустического сигнала в массиве (см. данные табл. 2) показал, что наибольшая точность измерений достигается для координат у и 2, она составляет единицы сантиметров или 4-6%. В то время как для координаты х точность измерений составляет 15-20%, это может быть связано с погрешностями при измерении скоростей распространения акустического сигнала в породном массиве из-за наличия в нем неоднородностей. Установлено, что на точность вычисления координат машины влияет метод обработки экспериментальных данных. При обработке данных по уровню 0,7итах расчетные значения координат наиболее близки к реальным, а при обработке данных по уровню 0,1- итах полученные в результате математических расчетов координаты машины резко отличаются от реальных значений. Это объясняется высоким уровнем шума в полосе полезного сигнала.
Выводы. Экспериментальные исследования и их математическая обработка показали, что по времени задержки распространения акустического сигнала, создаваемого движением буровой маши-
ны в породном массиве, относительно времени можно определять ее трехмерные координаты с
прихода опорного сигнала (с самой машины) точностью до единиц сантиметров.
СПИСОК ЛИТЕРАТУРЫ
1. Денисова Е.В., Тищенко И.В., Конурин А.И. О результатах испытаний многоканального комплекса для контроля геомеханических процессов // Вестник КузГТУ, Ю12. - №6. - С. 16-19.
2. Патент на полезную модель №118765 Измеритель координат импульсного акустического сигнала в среде. Патент РФ на полезную модель / В.Н. Опарин, Е.В. Денисова, С.Ю. Гаврилов, А.И. Конурин, Н.С. Полотнянко // Опубл. Бюл. - Ю12. - №21.
3. USA, United States Patent № 6,886,644. Apparatus and method for horizontal drilling / Vermeer Manufacturing Company. USA. Appl. No.: 1G/224,2G5., Publ. May 3, 2GG5.
4. USA, United States Patent №8,213,264. Method and device of measuring location, and moving object / Samsung Electronics Co., Ltd., USA. Appl. No.: 12/656,G24. Publ. July 3, 2G12.
5. USA, United States Patent № 8,264,9G9. System and method for depth determination of an impulse acoustic source by cepstral analysis / The United States of America as represented by the Secretary of the Navy, USA. Appl. No.: 12/698,679. Publ. September 11, 2G12.
□ Авторы статьи
Денисова Екатерина Вячеславовна, канд.техн.наук старший научный сотрудник Института горного дела СО РАН, Тел.(383) 217-09-52, е-шаД: slimthing@mail.ru
Тищенко Игорь Владимирович, канд.техн.наук, научный сотрудник Института горного дела СО РАН, телефон (383)217-05-13
Хмелинин Алексей Павлович, аспирант Института горного дела СО РАН,, телефон (383) 217-09-52
Бадмаева Янжима Зандановна, магистрант НГТУ
УДК 534.61
Е.В. Денисова, С.Ю. Гаврилов, А.И.Конурин
ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ И ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ АКУСТИЧЕСКОГО ПОЛЯ, СОЗДАВАЕМОГО НА ПОВЕРХНОСТИ ГРУНТОВОГО МАССИВА ПРИ ДВИЖЕНИИ В НЕМ
ПНЕВМОУДАРНОЙ МАШИНЫ
Использование акустического сигнала, создаваемого движением пневмоударной машины в грунтовом массиве, для определения ее местоположения 2 позволило бы существенно повысить надежность и экономичность современных систем навигации для бестраншейных технологий сооружения подземных переходов. Применение машины в качестве излучателя полезного сигнала исключает необходимость размещения в ее рабочем органе дополнительных устройств, что повышает надежность системы навигации к ударным нагрузкам, и исключает необходимость изготовления дополнительного передающего устройства, что в целом снижает затраты на реализацию системы [1].
Экспериментальные исследования акустического поля, создаваемого движением в грунтовом массиве пневмоударной машины, проводились на полигоне ИГД СО РАН "Зеленая горка". Вблизи места проведения исследований отсутствуют вы-
2 Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России на 2009-2013 гг.»
соковольтные сети, а так же другие подземные коммуникации. Массив, в который запускалась машина можно принять однородным. Во входной приямок помещается пневмоударная машина «Тайфун-2» (габаритные размеры: длина 55G мм, диаметр 4G мм, энергия удара 22 Дж) на глубину 5GG мм строго под измерительным акселерометром А6 (схема эксперимента на рис. 1). Многоканальный измерительный комплекс PULSE располагается на поверхности земли, акселерометры А1...Ап подключены к 11 каналам блока 3G53 и устанавливаются на выровненной поверхности земли с шагом iGG мм друг относительно друга вдоль оси х (фотография эксперимента на рис. 2). Для измерения поля, создаваемого машиной вдоль оси ее движения, акселерометры А1...А11 перемещают вдоль оси у (точки с координатами -3GG, -2GG, -1GG, G, 1GG, 2GG, 3GG, 4GG, 5GG, 6GG, 7GG и 8GG мм по схеме на рис. 1). Это позволит оценить степень изменения амплитуды акустического сигнала при произвольном расположении приемных каналов относительно пневмопробойника, характерном для проходки реальной скважины.