RU2661498C1 - Способ спектрально-акустического прогноза выбросоопасности угольных пластов - Google Patents
Способ спектрально-акустического прогноза выбросоопасности угольных пластов Download PDFInfo
- Publication number
- RU2661498C1 RU2661498C1 RU2017132120A RU2017132120A RU2661498C1 RU 2661498 C1 RU2661498 C1 RU 2661498C1 RU 2017132120 A RU2017132120 A RU 2017132120A RU 2017132120 A RU2017132120 A RU 2017132120A RU 2661498 C1 RU2661498 C1 RU 2661498C1
- Authority
- RU
- Russia
- Prior art keywords
- spectral
- amplitudes
- acoustic
- median
- coal
- Prior art date
Links
- 239000003245 coal Substances 0.000 title claims abstract description 41
- 238000013277 forecasting method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 58
- 230000003595 spectral effect Effects 0.000 claims abstract description 45
- 238000005065 mining Methods 0.000 claims abstract description 30
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 239000011435 rock Substances 0.000 claims abstract description 16
- 230000007246 mechanism Effects 0.000 claims abstract description 11
- 238000005553 drilling Methods 0.000 claims abstract description 10
- 231100001261 hazardous Toxicity 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 238000012935 Averaging Methods 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims abstract description 4
- 230000005534 acoustic noise Effects 0.000 claims description 5
- 238000009434 installation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 238000001228 spectrum Methods 0.000 description 11
- 238000012544 monitoring process Methods 0.000 description 8
- 238000013459 approach Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 241000566515 Nedra Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000000565 sealant Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F5/00—Means or methods for preventing, binding, depositing, or removing dust; Preventing explosions or fires
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21C—MINING OR QUARRYING
- E21C39/00—Devices for testing in situ the hardness or other properties of minerals, e.g. for giving information as to the selection of suitable mining tools
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Изобретение относится к горной промышленности и может быть использовано для прогноза динамических явлений типа внезапного выброса угля и газа, горного удара и им подобных. Техническим результатом является повышение достоверности текущего прогноза выбросоопасности при отработке угольных пластов. Предложен способ спектрально-акустического прогноза выбросоопасности угольных пластов, включающий непрерывное автоматическое измерение амплитуды генерируемого широкополосного акустического сигнала в горном массиве работающим механизмом, аналого-цифровое преобразование сигнала, формирование экспериментальных выборок из последовательности отсчетов аналого-цифрового преобразования, выполнение с ними процедуры быстрого преобразования Фурье и определение дискретного ряда амплитуд спектральных составляющих (гармоник), их усреднение во времени. При этом определяют текущее значение показателя выбросоопасности Rт по формуле Rт=(Smax,т-1,8)(gmax,т-a), где Smax,т и gmax,т – соответственно текущие максимальное значение выхода буровой мелочи и максимальное значение начальной скорости газовыделения при бурении контрольных шпуров, параметр a берут равным 5 для Воркутинского угольного месторождения, берут равным 4 для остальных бассейнов и месторождений восточных районов России. Затем определяют коэффициент относительной выбросоопасности инструментального метода прогноза Rо.в. контролируемой зоны пласта по формуле Rо.в.=Rт/b, где параметр b берут равным 21 для Воркутинского угольного месторождения, берут равным 6 для остальных бассейнов и месторождений восточных районов России. Определяют минимальное опорное значение медианы М0min дискретного ряда амплитуд спектральных гармоник шумов работающего механизма по формуле М0min=nmin⋅Δf, где Δf – интервал между соседними гармониками, находят критическое значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования по формуле Мкр=М0min/Rо.в. Измеряют текущее значение медианы Мт дискретного ряда амплитуд спектральных составляющих шумов и сравнивают его с критическим значением Мкр, при Мт≥Мкр зону пласта относят к выбросоопасной, при Мт<Мкр зону пласта относят к невыбросоопасной. 3 з.п. ф-лы, 1 ил., 1 табл.
Description
Изобретение относится к горной промышленности и может быть использовано для прогноза динамических явлений типа внезапного выброса угля и газа, горного удара и им подобных.
Анализ признаков выбросоопасности угольных пластов, зарегистрированных перед этим явлением в призабойном пространстве, а также известных моделей потери устойчивости горного массива при протекании этого явления, свидетельствует о том, что основными факторами, определяющими выбросоопасность, являются напряженное состояние призабойного пространства, внутрипластовое давление свободного газа и прочность угля. Поэтому известные способы прогноза динамических явлений, в том числе автоматизированного прогноза, основаны на контроле одного или нескольких этих факторов.
Известен способ текущего прогноза выбросоопасности угольных пластов по начальной скорости газовыделения и выходу буровой мелочи из скважин [Федеральные нормы и правила в области промышленной безопасности «Инструкция по прогнозу динамических явлений и мониторингу массива горных пород при отработке угольных месторождений». Утверждены приказом Федеральной службы по экологическому, технологическому и атомному надзору от 15.08.2016 г. №339. – 129 с.]. Метод основан на контроле основных факторов выбросоопасности: начальная скорость газовыделения контролирует газовый фактор, а выход буровой мелочи – напряженное состояние и прочность угля. Однако этот метод не является непрерывным, требует остановки ведения горных работ, достаточно трудоемок и продолжителен.
В качестве ближайшего аналога выбран способ акустического прогноза выбросоопасности угольных пластов, включающий генерирование акустических колебаний в горном массиве работающими механизмами, непрерывное измерение их амплитуд в области высоких и низких частот и оценку выбросоопасности пласта по отношению амплитуд высокочастотной и низкочастотной частей спектра акустических колебаний, которое является показателем выбросоопасности данного способа прогноза [Авторское свидетельство СССР №1222853, МПК Е21F5/00, Е21С39/00, опубл. 07.04.1986]. Преимуществом данного способа является непрерывность контроля выбросоопасности в процессе ведения горных работ. Но, поскольку на амплитуды спектральных компонент измеряемого акустического сигнала влияет напряженное состояние массива, а не давление газа в нем, данный способ контролирует преимущественно фактор напряженного состояния выбросоопасности. Именно поэтому он используется для прогноза динамических явлений типа внезапный выброс, горный удар и им подобных. По этой же причине ряд исследователей показатель опасности в форме отношения амплитуд из высокочастотной и низкочастотной областей спектра сигнала называют коэффициентом относительных напряжений [См, например, Копылов К.Н. Автоматизированная система контроля состояния массива горных пород и прогноза динамических явлений / К.Н. Копылов, О.В. Смирнов, А.И. Кулик, А.И. Пальцев // Безопасность труда в промышленности, 2015. - №8. - С. 32-37.].
Поскольку акустических методов контроля состояния горного массива много, а сущность рассматриваемого метода прогноза выбросоопасности основана на спектральном анализе «шумов» работающего горного оборудования, прошедших контролируемую область горного массива от источника до приемника, здесь и далее будем называть данный метод прогноза спектрально-акустическим.
При использовании этого метода для прогноза выбросоопасности идут по одному из двух путей:
По первому пути для того, чтобы скомпенсировать отсутствие учета давления газа и прочности перемятых угольных пачек в пласте, в качестве критического выбирают такое максимальное значение отношения амплитуд высокочастотной и низкочастотной частей спектра акустических колебаний Ккр, которое несколько ниже самого малого значения, из когда-либо замеренных перед внезапным выбросом угля и газа. Так, например, при использовании аппаратуры АК-1 нормативным документом регламентировалось устанавливать Ккр=3 [Предупреждение газодинамических явлений в угольных шахтах (Сборник документов) / Колл. авт. – М.: Государственное предприятие НТЦ по безопасности в промышленности Госгортехнадзора России, 2000. – С. 165-168]. Нечувствительность данного способа прогноза к газовому фактору выбросоопасности и наличию перемятых (с низкой прочностью) пачек угля определяют завышенный «запас надежности» этого способа прогноза и, как следствие, недостаточную точность (достоверность) и экономическую эффективность.
Согласно второму пути разрабатывают способы и реализующие их многофункциональные системы, в которых критическое значение показателя выбросоопасности при прогнозе спектрально-акустическим методом корректируют путем учета газового фактора и прочности угля (например, патент на полезную модель №34202, МПК Е21F5/00, 27.11.2003; патент на изобретение № 2231649, МПК Е21F5/00, опубл. 27.06.2004; патент на изобретение № 2250376, МПК Е21F5/00, опубл. 20.04.2005).
Таким образом, спектрально-акустический метод прогноза выбросоопасности и других типов динамических явлений используется как самостоятельно, так и в составе многофункциональных систем.
Однако реализация данного метода осложнена выбором рабочих частот, характеризующих низкочастотную и высокочастотную области спектра шумов работающего горного оборудования.
В настоящее время известны два подхода к определению рабочих частот. Первый реализуется, например, аппаратурой АК-1 (или ее модификацией АК-1М) и заключается в разбиении рабочего частотного диапазона (например, 20-1500 Гц) на поддиапазоны низких и высоких частот. Разбиение осуществляется с помощью фильтров нижних и верхних частот.
Так, например, нормативным документом при использовании аппаратуры АК-1 или АК-1М регламентируется провести оценочные разведочные наблюдения для выбора частот среза фильтров высоких частот (ФВЧ) на одно из трех значений 600, 800 или 1000 Гц, и частот среза фильтров низких частот (ФНЧ) на одно из трех значений 160, 200 или 300 Гц [Руководство по выполнению спектрально-акустического контроля (прогноза) выбросоопасности на шахтах Кузбасса. Утверждено Кузнецким управлением Гостехнадзора России 23.04.2002 г. // Кемерово, 2002, п. 2.3.2].
Недостатки этого подхода в следующем:
1. Не используются спектральные составляющие акустического сигнала, лежащие между частотами среза ФВЧ и ФНЧ. Поэтому, если изменение напряженного состояния горного массива при ведении горной выработки приведет к существенному изменению спектра акустического сигнала именно в этой области частот, оно окажется незамеченным.
2. Необходимость проведения специальных работ для выбора частот среза фильтров высоких и низких частот и периодической их корректировки в зависимости от изменения горно-геологических и горнотехнических условий.
3. Опыт применения этого подхода показал, что невозможно установить единое значение критической величины показателя выбросоопасности для всех шахт даже одного угольного бассейна. Поэтому критическое значение показателя выбросоопасности необходимо устанавливать экспериментально, однако методики выполнения этой процедуры нет.
Второй подход к выбору рабочих частот реализуется, например, системой акустического контроля состояния массива горных пород и прогноза динамических явлений (САКСМ) [Руководство по применению системы акустического контроля состояния массива горных пород и прогноза динамических явлений // М.: МНТЛ РИВАС. 2016. – 49 с.]. В соответствии с этим документом (п. 2.5 на с. 6) рабочий диапазон частот лежит в пределах 20-3500 Гц, а области низких и высоких частот устанавливаются следующим образом. Предполагается, что амплитудно-частотная характеристика (АЧХ) акустического сигнала в зоне размещения приемника акустических колебаний имеет максимум. Обработка сигнала осуществляется автоматизировано с помощью специально разработанного программного обеспечения. Определяется частота, на которой амплитуда сигнала имеет максимальное значение Амах. Левее этой частоты определяют частоты, сигналы на которых имеют соответственно амплитуды, равные 0,5Амах и 0,75Амах. Эти частоты являются граничными низкочастотной области спектра. Аналогично этому правее частоты, соответствующей Амах, определяют частоты, сигналы на которых имеют соответственно амплитуды, равные 0,75Амах и 0,5Амах. Эти частоты являются граничными высокочастотной области спектра.
Преимущество этого подхода, по сравнению с предыдущим, состоит в том, что границы высоких и низких частот не являются фиксированными, а автоматически корректируются в зависимости от АЧХ принимаемого сигнала.
Но данный подход имеет следующие недостатки:
1. Используется не весь спектр для определения показателя опасности (коэффициента относительных напряжений).
2. Подход предполагает, что форма АЧХ сигнала имеет явно выраженный максимум. Однако на практике регистрируются, в том числе, сигналы с убывающей формой АЧХ при изменении частоты от минимальной к максимальной в области рабочих частот. Для таких случаев алгоритм расчета предполагает принудительное искажение спектра сигнала, что влияет на показатель опасности.
3. Не существует единого критического значения показателя выбросоопасности для всех шахт даже одного угольного бассейна. Поэтому критическое значение показателя выбросоопасности необходимо устанавливать экспериментально, однако методики выполнения этой процедуры нет.
Задача изобретения – повышение достоверности текущего прогноза выбросоопасности (и возможно, других типов динамических явлений) при отработке угольных пластов.
Это достигается тем, что при осуществлении способа спектрально-акустического прогноза выбросоопасности угольных пластов, включающего непрерывное автоматическое измерение амплитуды генерируемого широкополосного акустического сигнала в горном массиве работающим механизмом, аналого-цифровое преобразование сигнала, формирование экспериментальных выборок из последовательности отсчетов аналого-цифрового преобразования, выполнение с ними процедуры быстрого преобразования Фурье и определение дискретного ряда амплитуд спектральных составляющих (гармоник), их усреднение во времени, при этом дополнительно инструментальным методом определяют текущее значение показателя выбросоопасности Rт по формуле Rт=(Smax,т-1,8)(gmax,т-a), где Smax,т и gmax,т – соответственно текущие максимальное значение выхода буровой мелочи и максимальное значение начальной скорости газовыделения при бурении контрольных скважин. Параметр a берут равным 5 для Воркутинского угольного месторождения, берут равным 4 для остальных бассейнов и месторождений восточных районов России.
Определяют коэффициент относительной выбросоопасности инструментального метода прогноза Rо.в. контролируемой зоны пласта по формуле Rо.в.=Rт/b, где параметр b берут равным 21 для Воркутинского угольного месторождения, берут равным 6 для остальных бассейнов и месторождений восточных районов России. Время цикла подвигания забоя разбивают на равные временные интервалы, определяют полусумму амплитуд всех усредненных спектральных составляющих сигнала 
в каждый j-й временной интервал одного цикла подвигания забоя, находят медианный интервал частот, в котором находится медиана путем определения номера nj спектральной гармоники, при которой выполняются условия 
, 
, для одного цикла подвигания забоя определяют номер nmin спектральной составляющей, соответствующей минимальному медианному значению частоты акустического сигнала работающего механизма. Затем определяют минимальное опорное значение медианы М0min дискретного ряда амплитуд спектральных гармоник шумов работающего механизма по формуле М0min=nmin⋅Δf, где Δf – интервал между соседними гармониками. Находят критическое значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования по формуле Мкр=М0min/Rо.в. Измеряют текущее значение медианы Мт дискретного ряда амплитуд спектральных составляющих шумов и сравнивают его с критическим значением Мкр: при Мт ≥ Мкр зону пласта относят к выбросоопасной, при Мт < Мкр зону пласта относят к невыбросоопасной.
Рабочий диапазон частот ограничивают максимальным значением, например 1000 Гц. Определение критического значения медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования проводят отдельно для каждого типа работающего оборудования. Установку приемника акустических колебаний в подготовительной выработке при определении опорного значения медианы дискретного ряда амплитуд спектральных составляющих шумов осуществляют на удалении от забоя выработки, близком к предельному, соответствующему отношению амплитуды акустического сигнала работающего оборудования к акустическим помехам, близкому к трем.
Описание иллюстрирующих материалов: на фиг. 1 показана амплитудно-частотная характеристика смоделированного сигнала при σт/σпр=0,2; 0,4; 0,6; 0,8 и 1,0, где σт и σпр – соответственно средние предельное (максимально возможное для контролируемого участка пласта, предшествующее разрушению массива из-за динамического явления) и текущее значения напряжений в массиве.
Обоснуем справедливость предлагаемого способа путем моделирования акустического сигнала с помощью суммы гармоник. Предположим, что аналоговый акустический сигнал, генерируемый режущим органом работающего оборудования, например, проходческого комбайна, принимается приемным преобразователем на расстоянии r от излучателя, оцифровывается, преобразуется с помощью БПФ в амплитудно-частотный спектр, который ограничивается в диапазоне частот 20-1000 Гц, причем интервал между соседними гармониками Δf равен, например, 20 Гц. Тогда частота i-ой гармоники будет равна
Ограничение частоты «сверху» обусловлено следующим. При контроле напряженного состояния впереди забоя подготовительной выработки источник и приемник акустических колебаний находятся позади контролируемой области горного массива. Однако, акустическая волна при своем движении не может быть сфокусирована в объеме среды (цилиндра, трубки), геометрические размеры поперечного сечения которой теоретически меньше половины длины волны [Савич А.И. О зоне «захвата» упругих волн // Труды Гидропроекта, 1971. - №21. - С. 29-40.]. Практически эта величина соизмерима с длиной волны. Следовательно, как минимум на эту глубину угольный пласт «прозвучивается» звуком, излучаемым режущим органом комбайна. Скорость звука в каменном угле для различных угольных бассейнов заметно отличается. Так, для Карагандинского угольного бассейна скорость продольного и поперечного звука минимальна (из основных угольных бассейнов на территории бывшего СССР) и соответственно равна 1,2-1,5 км/с и 0,8-1,0 км/с [Азаров Н.Я. Сейсмоакустический метод прогноза горно-геологических условий эксплуатации угольных месторождений / Н.Я. Азаров, Д.В. Яковлев // М.: Недра, 1988. – 199 с.]. Следовательно, длины волн продольного и поперечного звука частотой 1000 Гц (соответствует минимальной длине волны в выбранном рабочем диапазоне частот) здесь будут соответственно равны 1,2-1,5 м и 0,8-1,0 м. Для Печерского угольного бассейна скорость продольного и поперечного звука максимальна и соответственно равна 2,3-3,6 км/с и 1,0-1,3 км/с. Для этих скоростей длины волн продольного и поперечного звука частотой 1000 Гц будут находиться соответственно в пределах 2,3-3,6 м и 1,0-1,3 м. Таким образом, чтобы глубина контроля массива впереди забоя была не менее одного 1 м, использовать частоты, например, свыше 1000 Гц нецелесообразно.
Известно, что амплитуду i-ой гармоники акустического сигнала на расстоянии r от источника можно описать следующим образом:
где Ai0 – амплитуда i-ой гармоники у источника; F(r) – функция, учитывающая диаграмму направленности источника сигнала (например, F(r)=1 – для плоской волны; F(r)=1/r – для сферической волны); α – коэффициент затухания звука.
Известно также, что для твердого тела коэффициент затухания звука в первом приближении прямо пропорционален частоте и обратно пропорционален действующим на тело средним текущим напряжениям α ~ f/σт. Поэтому можно записать следующее равенство [Шадрин А.В. Основы автоматизированного непрерывного ГДЯ-мониторинга на угольных шахтах Кузбасса / А.В. Шадрин, В.А. Коноваленко // Вестник КузГТУ, 2001. - №3. - С.28-31]:
где α0 – затухание на некоторой частоте f0 при отсутствии напряжений (в разгруженном состоянии), м-1; β – безразмерный коэффициент пропорциональности, определяемый акустическими свойствами массива; fi – частота i-ой гармоники, Гц; σпр и σт – средние соответственно предельное (максимально возможное для контролируемого участка пласта, предшествующее разрушению массива из-за динамического явления) и текущее значения напряжений в массиве, Па; 
, м-1Гц-1.
Тогда для выбранного нами диапазона частот сумма дискретного ряда амплитуд спектральных составляющих акустического сигнала 
с учетом (1-3) будет равна:
Предположим, что АЧХ акустического сигнала имеет максимум и ее компоненты у излучателя могут быть описаны двумя экспоненциальными функциями (возрастающей и убывающей) в следующем виде:
где А0 – амплитуда гипотетической «нулевой» гармоники при i=0, В; параметры ξ и η – определяют скорости изменения экспонент.
Определим параметры ξ и η из условия, что А0=1 В и амплитуды следующих гармоник удовлетворяют условию А14,0=А15,0 ≈ 1 В (условие «сшивания» двух экспонент). Тогда (5) примет вид:
Подставим (6) в (4). Результаты расчета амплитуд отдельных гармоник представлены на фиг. 1 при следующих значениях входящих величин: α0=1,3 м-1; β=0,07; f0=500 Гц; F(r)=1; r =10 м.
Из рисунка видно, что с ростом напряжений амплитуды высокочастотных гармоник возрастают сильнее низкочастотных.
Для приведенных параметров акустического сигнала определим зависимость текущего значения медианы дискретного ряда амплитуд спектральных составляющих Мт от отношения 
При этом в соответствии с определением под медианой будем понимать корень уравнения A(f,r)=0,5, иначе говоря, в нашем случае Мт – это значение частоты гармоники fn, при которой выполняется условие:
Введем понятие показателя выбросоопасности спектрально-акустического метода (коэффициента относительных напряжений) К, равного отношению текущего значения медианы Мт к опорному М0, определенному аналогично Мт на участке проводимой горной выработки, на котором предварительно был выполнен прогноз инструментальным методом по начальной скорости газовыделения и выходу буровой мелочи при бурении контрольных шпуров:
К=Мт/М0. (8)
В таблице в качестве примера представлены значения показателя выбросоопасности для отношений текущего и предельного значений средних напряжений при условии, что в качестве опорного М0 взято значение медианы М0,2 при 
.
Таблица
Зависимость показателя выбросоопасности спектрально-акустического метода прогноза от напряженного состояния
| σ т /σ пр | К=М т /М 0,2 |
| 0,2 | 1,0 |
| 0,4 | 1,8 |
| 0,6 | 2,2 |
| 0,8 | 2,4 |
| 1,0 | 2,6 |
Из таблицы видно, что при рассмотренной модели акустического сигнала, имитирующего шум работающего комбайна, увеличение отношения текущих напряжений к предельным в 5 раз привело к росту показателя выбросоопасности (коэффициента относительных напряжений) К в 2,6 раза. При этом для определения показателя выбросоопасности использовался весь спектр акустического сигнала, что исключает появление ошибки в прогнозе из-за того, что изменение напряженного состояния призабойного пространства привело к изменению амплитуд спектральных составляющих, находящихся за пределами рабочего диапазона частот.
Предлагаемый способ реализуется непосредственно следующим образом. Спектрально-акустический прогноз выбросоопасности спектрально-акустическим методом начинают с определения критерия выбросоопасности этим методом в контролируемой выработке. Для этого выбирают экспериментальный участок, на котором выполняют прогноз выбросоопасности вначале инструментальным методом, характеризующимся высокой степенью достоверности прогноза, а затем спектрально-акустическим. Из инструментальных методов в настоящее время на угольных шахтах восточных районов России наибольшей достоверностью обладает метод текущего прогноза по начальной скорости газовыделения и выходу буровой мелочи при бурении контрольных скважин (шпуров). В соответствии с этим методом критерий выбросоопасности для угольных шахт России описывается следующим выражением [Федеральные нормы и правила в области промышленной безопасности «Инструкция по прогнозу динамических явлений и мониторингу массива горных пород при отработке угольных месторождений». Утверждены приказом Федеральной службы по экологическому, технологическому и атомному надзору от 15.08.2016 г. №339. – 129 с.]:
R=(Smax-1,8)(gmax-a) – b=0, (9)
где R – безразмерный показатель выбросоопасности; Smax – максимальное значение выхода буровой мелочи, л/м; gmax – максимальное значение начальной скорости газовыделения, л/мин⋅м; а=5 – для Воркутинского месторождения; а=4 – для остальных бассейнов и месторождений восточных районов России; b=21 – для Воркутинского месторождения; b=6 – для остальных бассейнов и месторождений восточных районов России.
При R ≥ 0 зона пласта относится к выбросоопасной.
Из (9) определяют следующим образом критическое значение показателя выбросоопасности Rкр инструментального метода прогноза:
Rкр=(Smax-1,8)(gmax-a)=b, (10)
Далее определяют текущее значение показателя выбросоопасности Rт по следующей формуле:
Rт=(Smax,т-1,8)(gmax,т-a), (11)
где Smax,т и gmax,т – текущие соответственно максимальное значение выхода буровой мелочи и максимальное значение начальной скорости газовыделения.
Затем определяют коэффициент относительной выбросоопасности инструментального метода прогноза Rо.в. контролируемой зоны пласта как отношение текущего и критического значений показателя выбросоопасности:
После этого на этом же участке пласта определяют опорное значение медианы М0 дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования, например, проходческого комбайна. Для этого на участке проводимой горной выработки, оснащенной аппаратурой для спектрально-акустического прогноза, например, САКСМ, в соответствии с нормативным документом, например, инструкцией по эксплуатации, во время одного цикла подвигания забоя измеряют амплитуду акустического сигнала, генерируемого работающим механизмом, проводят аналого-цифровое преобразование сигнала, формируют экспериментальные выборки из последовательности отсчетов аналого-цифрового преобразования, выполняют с ними процедуру БПФ, проводят их усреднение во времени, например, за 15-секундные интервалы времени, в результате чего получают дискретный ряд усредненных амплитуд спектральных составляющих Aij (i∈[1; N], j∈[1; P],), где Aij – амплитуда i-ой гармоники в j-й 15-и секундный интервал; P – количество 15-и секундных интервалов за время цикла подвигания забоя. Далее определяют полусумму амплитуд всех усредненных спектральных составляющих сигнала 
в каждый j-й 15-и секундный интервал и находят медианный интервал частот, в котором находится медиана путем определения номера nj спектральной составляющей (гармоники), при которой выполняются условия:
Для простоты принимают опорное значение медианы М0j, равным меньшему значению частоты медианного интервала:
М0j=nj⋅Δf, (14)
где Δf – интервал между соседними гармониками.
Измерения проводят в течение одного цикла подвигания забоя, по результатам которого определяют минимальное опорное значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования М0min. Причем процедуру определения минимального опорного значения медианы М0min проводят отдельно для каждого типа работающего горного оборудования (комбайна, струга, буровой установки, электросверла, отбойного молотка), и из обработки исключают результаты нескольких, например, трех первых 15-и секундных интервалов.
Далее определяют (вводят понятие) коэффициента относительной выбросоопасности спектрально-акустического метода Ко.в. как отношение текущего значения медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования к ее критическому значению:
Предполагая, что коэффициенты относительной выбросоопасности инструментального и спектрально-акустического метода примерно равны, находят критическое значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования:
Мкр=М0min/Rо.в. (16)
Затем осуществляют текущий прогноз выбросоопасности спектрально-акустическим методом путем измерения текущего значения медианы дискретного ряда амплитуд спектральных составляющих шумов и сравнения его с критическим значением. При условии Мт≥Мкр зону пласта относят к выбросоопасной, а если Мт<Мкр зону пласта относят к невыбросоопасной.
По мере подвигания забоя подготовительной или очистной выработки критическое значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования корректируют в случае изменения горно-геологических условий. При кажущемся постоянстве горно-геологических условий коррекцию критического значения медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования проводят не более чем через, например, 300 м подвигания забоя.
При выборе места установки приемника акустических колебаний в подготовительной выработке при определении опорного значения медианы дискретного ряда амплитуд спектральных составляющих шумов следует учесть следующее. Он должен быть установлен на удалении от забоя выработки, близком к предельному, например, на расстоянии 25 м. Это делается с тем, чтобы получить минимальную величину критического значения медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования, что позволит избежать при прогнозе ошибки первого рода. Предельное расстояние определяется коэффициентом затухания звука в конкретной выработке и уровнем акустических шумов, замеренным до начала работы горного оборудования, акустическое излучение которого используются для прогноза выбросоопасности. Это расстояние соответствует отношению амплитуды акустического сигнала работающего оборудования к акустическим помехам, близкому к трем.
Для осуществления предлагаемого способа применяют серийно выпускаемое оборудование: для измерения начальной скорости газовыделения при бурении контрольных шпуров – ручное электросверло, комплект витых составных штанг длиной до 6,5 м, герметизатор скважины и измеритель начальной скорости газовыделения, например, типа ИГ-1; для измерения выхода буровой мелочи – мерный сосуд; для спектрально-акустического прогноза (контроля) выбросоопасности – например, систему акустического контроля состояния массива горных пород и прогноза динамических явлений (САКСМ), состоящую из аппаратуры регистрации и передачи на поверхность акустического сигнала (АРАС), программного обеспечения акустического контроля массива горных пород и прогноза динамических явлений (программа АКМП-РИВАС), персонального компьютера, принтера и устройства бесперебойного питания. Связь подземной части аппаратуры с наземной осуществляется линией связи, в качестве которой может использоваться, например, свободная пара проводов телефонного кабеля шахты.
Алгоритм цифровой обработки дополнительно к выполняемым программой АКМП-РИВАС операциям включает определение критического значения медианы дискретного ряда усредненных во времени амплитуд спектральных составляющих Мкр по формулам (13)-(16), а затем текущего значения медианы дискретного ряда усредненных во времени амплитуд спектральных составляющих Мт по формулам (13) и (14), в которых индекс j заменен на индекс т (текущий), и отнесение зоны пласта к выбросоопасной, если Мт≥Мкр, а если Мт<Мкр, то зону пласта относят к невыбросоопасной.
Пример. По результатам бурения контрольной скважины определили инструментальным методом текущего прогноза выбросоопасности по начальной скорости газовыделения и выходу буровой мелочи максимальное значение начальной скорости газовыделения gmax и максимальное значение выхода штыба Smax. По ним по формуле (11) определили текущее значение показателя выбросоопасности инструментального метода Rт=4,8. По формуле (12) определили коэффициент относительной выбросоопасности инструментального метода Rо.в=Rт/b=4,8/6=0,8. Далее приступили к выполнению цикла подвигания забоя при ведении, например, подготовительной выработки проходческим комбайном. Длительность цикла разбили на 15-секундные интервалы времени, по окончании каждого из которых определяли из условия (13) номер nj спектральной составляющей, соответствующей медианному значению частоты АЧХ акустического «шума» комбайна. Проанализировав весь цикл подвигания забоя, нашли, что номер nmin спектральной составляющей, соответствующей минимальному медианному значению частоты АЧХ акустического «шума» комбайна, nmin=12. По формуле (14) определили минимальное опорное значение медианы М0min=nmin⋅Δf=12⋅20=240 Гц. Затем по формуле (16) нашли критическое значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего комбайна Мкр=240/0,8=300 Гц.
После этого перешли к прогнозу выбросоопасности с помощью спектрально-акустического метода в автоматическом режиме. Во время работы проходческого комбайна, для которого определили значение Мкр=300 Гц, измерили в 15-и секундные интервалы времени следующие текущие значения медианы дискретного ряда амплитуд спектральных составляющих шумов работающего комбайна Мт=180; 220; 240; 260; 200; 280 Гц, сравнили их с критическим значением медианы Мкр=300 Гц и пришли к выводу, что призабойное пространство впереди выработки выбрособезопасно до выполнения следующего 15-секундного прогноза.
Claims (15)
1. Способ спектрально-акустического прогноза выбросоопасности угольных пластов, включающий непрерывное автоматическое измерение амплитуды генерируемого широкополосного акустического сигнала в горном массиве работающим механизмом, аналого-цифровое преобразование сигнала, формирование экспериментальных выборок из последовательности отсчетов аналого-цифрового преобразования, выполнение с ними процедуры быстрого преобразования Фурье и определение дискретного ряда амплитуд спектральных составляющих (гармоник), их усреднение во времени, отличающийся тем, что определяют текущее значение показателя выбросоопасности R т по формуле
R т=(S max,т-1,8)(g max,т-a),
где S max,т и g max,т - соответственно текущие максимальное значение выхода буровой мелочи и максимальное значение начальной скорости газовыделения при бурении контрольных шпуров, параметр a берут равным 5 для Воркутинского угольного месторождения, берут равным 4 для остальных бассейнов и месторождений восточных районов России,
определяют коэффициент относительной выбросоопасности инструментального метода прогноза R o.в. контролируемой зоны пласта по формуле
R o.в.=R т/b,
где параметр b берут равным 21 для Воркутинского угольного месторождения, берут равным 6 для остальных бассейнов и месторождений восточных районов России,
время цикла подвигания забоя разбивают на равные временные интервалы, определяют полусумму амплитуд всех усредненных спектральных составляющих акустического сигнала работающего механизма 
в каждый j-й временной интервал одного цикла подвигания забоя, находят медианный интервал частот, в котором находится медиана, путем определения номера n j спектральной гармоники, для которой выполняются условия 
, 
, для одного цикла подвигания забоя определяют номер n min спектральной составляющей, соответствующей минимальному медианному значению частоты акустического сигнала работающего механизма, определяют минимальное опорное значение медианы M 0min дискретного ряда амплитуд спектральных гармоник шумов работающего механизма по формуле
M 0min=n min⋅Δf,
где Δf – интервал между соседними гармониками,
находят критическое значение медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования по формуле
М кр=M 0min /R o.в,
измеряют текущее значение медианы М т дискретного ряда амплитуд спектральных составляющих шумов и сравнивают его с критическим значением М кр, при М т≥М кр зону пласта относят к выбросоопасной, при М т<М кр зону пласта относят к невыбросоопасной.
2. Способ спектрально-акустического прогноза выбросоопасности угольных пластов по п. 1, отличающийся тем, что рабочий диапазон частот ограничивают максимальным значением, например 1000 Гц.
3. Способ спектрально-акустического прогноза выбросоопасности угольных пластов по п. 1, отличающийся тем, что определение критического значения медианы дискретного ряда амплитуд спектральных составляющих шумов работающего горного оборудования проводят отдельно для каждого типа работающего оборудования.
4. Способ спектрально-акустического прогноза выбросоопасности угольных пластов по п. 1, отличающийся тем, что установку приемника акустических колебаний в подготовительной выработке при определении опорного значения медианы дискретного ряда амплитуд спектральных составляющих шумов осуществляют на удалении от забоя выработки, близком к предельному, соответствующему отношению амплитуды акустического сигнала работающего оборудования к акустическим помехам, близкому к трем.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2017132120A RU2661498C1 (ru) | 2017-09-13 | 2017-09-13 | Способ спектрально-акустического прогноза выбросоопасности угольных пластов |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2017132120A RU2661498C1 (ru) | 2017-09-13 | 2017-09-13 | Способ спектрально-акустического прогноза выбросоопасности угольных пластов |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2661498C1 true RU2661498C1 (ru) | 2018-07-17 |
Family
ID=62917311
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2017132120A RU2661498C1 (ru) | 2017-09-13 | 2017-09-13 | Способ спектрально-акустического прогноза выбросоопасности угольных пластов |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2661498C1 (ru) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113187559A (zh) * | 2021-05-31 | 2021-07-30 | 煤炭科学研究总院 | 矿井安全预警方法及装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU861648A1 (ru) * | 1980-01-10 | 1981-09-07 | Всесоюзный научно-исследовательский институт горной геомеханики и маркшейдерского дела | Способ определени выбросоопасных участков угольного пласта при ведении горных работ в очистных и подготовительных выработках |
| SU1222853A1 (ru) * | 1984-10-23 | 1986-04-07 | Днепропетровский Ордена Трудового Красного Знамени Горный Институт Им.Артема | Способ акустического прогноза выбросоопасности угольных пластов и устройство дл его осуществлени |
| SU1657659A1 (ru) * | 1989-07-03 | 1991-06-23 | А.М.Морев, В.С.Маевский Ю.В.Деев и С.И.Никишин | Способ дегазации угольного пласта |
| RU2231649C1 (ru) * | 2002-11-10 | 2004-06-27 | Кемеровский государственный университет | Способ текущего прогноза внезапных выбросов угля и газа |
| RU2250376C2 (ru) * | 2003-06-16 | 2005-04-20 | Кемеровский государственный университет | Способ автоматизированного текущего прогноза внезапных выбросов угля и газа |
| RU2467171C1 (ru) * | 2011-06-01 | 2012-11-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Новосибирский национальный исследовательский государственный университет" (Новосибирский государственный университет, НГУ) | Способ диагностики опасных ситуаций при подземной добыче каменного угля и методика прогноза параметров зон трещиноватости, образованной гидроразрывом пласта |
| CN103089305A (zh) * | 2013-01-14 | 2013-05-08 | 中国矿业大学(北京) | 一种煤层煤与瓦斯突出危险性评估方法 |
-
2017
- 2017-09-13 RU RU2017132120A patent/RU2661498C1/ru not_active IP Right Cessation
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU861648A1 (ru) * | 1980-01-10 | 1981-09-07 | Всесоюзный научно-исследовательский институт горной геомеханики и маркшейдерского дела | Способ определени выбросоопасных участков угольного пласта при ведении горных работ в очистных и подготовительных выработках |
| SU1222853A1 (ru) * | 1984-10-23 | 1986-04-07 | Днепропетровский Ордена Трудового Красного Знамени Горный Институт Им.Артема | Способ акустического прогноза выбросоопасности угольных пластов и устройство дл его осуществлени |
| SU1657659A1 (ru) * | 1989-07-03 | 1991-06-23 | А.М.Морев, В.С.Маевский Ю.В.Деев и С.И.Никишин | Способ дегазации угольного пласта |
| RU2231649C1 (ru) * | 2002-11-10 | 2004-06-27 | Кемеровский государственный университет | Способ текущего прогноза внезапных выбросов угля и газа |
| RU2250376C2 (ru) * | 2003-06-16 | 2005-04-20 | Кемеровский государственный университет | Способ автоматизированного текущего прогноза внезапных выбросов угля и газа |
| RU2467171C1 (ru) * | 2011-06-01 | 2012-11-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Новосибирский национальный исследовательский государственный университет" (Новосибирский государственный университет, НГУ) | Способ диагностики опасных ситуаций при подземной добыче каменного угля и методика прогноза параметров зон трещиноватости, образованной гидроразрывом пласта |
| CN103089305A (zh) * | 2013-01-14 | 2013-05-08 | 中国矿业大学(北京) | 一种煤层煤与瓦斯突出危险性评估方法 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113187559A (zh) * | 2021-05-31 | 2021-07-30 | 煤炭科学研究总院 | 矿井安全预警方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10669846B2 (en) | Apparatus, computer readable medium, and program code for evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
| US10698427B2 (en) | System and method for assessing sand flow rate | |
| US10415373B2 (en) | Submersible pump monitoring | |
| CN102292518B (zh) | 利用分布式感测装置进行井监测 | |
| US10180061B2 (en) | Methods of evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
| EP3433643B1 (en) | Method and device for estimating sonic slowness in a subterranean formation | |
| US11231512B2 (en) | Apparatus and methods of evaluating rock properties while drilling using acoustic sensors installed in the drilling fluid circulation system of a drilling rig | |
| CN112925021B (zh) | 用面波探测地层横波速度径向分布的测井方法 | |
| Young et al. | Seismic spectroscopy in fracture characterization | |
| Khoshouei et al. | Application of Acoustic Emission (AE) in mining and earth sciences: a review | |
| WO2019040639A1 (en) | SYSTEM AND METHOD FOR EVALUATING SAND FLOW | |
| US9260962B2 (en) | Reducing distortion in amplified signals in well logging tools | |
| RU2661498C1 (ru) | Способ спектрально-акустического прогноза выбросоопасности угольных пластов | |
| CN105008663B (zh) | 揭示未钻探的地质地层中孔隙压力的异常不连续性界面的方法以及实现其的系统 | |
| WO2020097090A1 (en) | Apparatus and methods of evaluating rock properties while drilling using acoustic sensors installed in the drilling fluid circulation system of a drilling rig | |
| CN106801602A (zh) | 利用随钻测量工具的压力波信号实时监测气侵的方法 | |
| Deb et al. | Estimation of blast induced peak particle velocity at underground mine structures originating from neighbouring surface mine | |
| Zastrelov et al. | Determination of geological conditions of gassy coal seams on the basis of seismic acoustic profiling in underground mine workings | |
| RU2310758C1 (ru) | Способ контроля состояния призабойной части массива | |
| Shadrin et al. | The methodology for working out the techniques for geophysical monitoring of the process of directional hydraulic fracturing of a hard roof and the process of local hydrotreating of a coal seam face working space, its implementation within the framework of RSF project | |
| US10253620B1 (en) | System for kick detection during a drilling operation | |
| Zberovskyi | Control of the mud pulse method the loosening of coal layers by amplitude-frequency recommendation of acoustic signal by the APSS-1 system | |
| CN117287261A (zh) | 一种煤层突出危险性频谱声学预测方法 | |
| Reuter et al. | Acoustic Noise in Coal Seam as an Indicator of Geodynamic Behavior in Longwall | |
| Shadrin et al. | Researching acoustic radiation spectrum of the mining equipment used for controlling stress state of a face space |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20190914 |