RU2498353C1 - Способ контроля изменения несплошности в массиве горных пород - Google Patents
Способ контроля изменения несплошности в массиве горных пород Download PDFInfo
- Publication number
- RU2498353C1 RU2498353C1 RU2012123762/28A RU2012123762A RU2498353C1 RU 2498353 C1 RU2498353 C1 RU 2498353C1 RU 2012123762/28 A RU2012123762/28 A RU 2012123762/28A RU 2012123762 A RU2012123762 A RU 2012123762A RU 2498353 C1 RU2498353 C1 RU 2498353C1
- Authority
- RU
- Russia
- Prior art keywords
- rock mass
- received pulse
- discontinuity
- changes
- change
- Prior art date
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относится к геофизике и может быть использовано в горной промышленности для контроля изменения состояния массива горных пород на более ранней стадии образования несплошностей, ведущих к динамическим проявлениям и разрушениям. Согласно заявленному способу контроля изменений несплошностей в массиве в качестве зондирующего сигнала используют серии отдельных одиночных прямоугольных акустических импульсов. Дополнительно определяют энергию каждого принятого импульса в выбранных частотных интервалах. Определяют отношение, а об изменениях несплошностей в массиве горных пород судят по невыполнению неравенства. Технический результат - повышение чувствительности контроля изменения напряженно-деформированного состояния массива горных пород. 1 ил.
Description
Изобретение относится к геофизике и может быть использовано в горной промышленности для ведения мероприятий, контролирующих зарождающиеся процессы разупрочнения горного массива, ведущих к динамическим проявлениям и, как следствие, человеческим жертвам.
Известен способ [1], в котором ведется прозвучивание массива шумовым стационарным сигналом, имеющим нормальное распределение со средним значение равным нулю, и при приеме измеряют дисперсию этого сигнала и период автокорреляционной функции. Считают, что положение максимума зоны опорного давления соответствует точке массива с максимальной дисперсией и минимальным периодом автокорреляции функции принятого сигнала.
К недостаткам следует отнести низкую чувствительность, т.к. этот способ интегральный и в нем оценивается дисперсия всего сигнала в целом. И, как следствие, он не способен различать мелкие разупрочнения, появляющиеся при различных нагрузках.
Более близким является способ контроля изменения напряженного состояния горного массива, в котором ведут прозвучивание участков массива сигналами, использующими серию из нескольких пачек равномерных импульсов, оценивают дисперсию составляющих спектра сигнала и по их изменению в последующих сериях судят об изменении напряженного состояния горного массива в целом [2].
К недостаткам следует отнести низкую чувствительность и сложность реализации, т.к. спектр таких сигналов весьма сложен и выделить составляющие гармоник является непростой задачей, что не позволяет вести контроль в реальном времени.
Заявленное решение направлено на повышение чувствительности контроля изменения несплошностей в массиве горных пород для выбора прогностических характеристик и ведения прогноза разупрочнения контролируемого участка горного массива в реальном времени.
Поставленная цель достигается тем, что в известном способе заключающимся в прозвучивании массива сигналами, приеме сигналов, прошедших контролируемой участок массива горных пород, в качестве зондирующего сигнала используют серии отдельных одиночных прямоугольных импульсов, дополнительно определяют энергию каждого принятого импульса в выбранных частотных интервалах и определяют отношение
где A - коэффициент, пропорциональный амплитуде принятого импульса;
k=tgβ, β - угол наклона фронтов принятого импульса;
τ - длительность принятого импульса;
а об изменениях несплошностей в массиве горных пород определяют по неравенству
Si,(i+p)<Fα(ri,ri+p)
где Fα - верхний α предел F - распределения со степенями свободы ri и ri+p;
α - уровень значности,
ri=bi(ni-1); ri+p=bi+p(ni+p-1);
γi=(ni-2)Ri-3; γi+p=(ni+p-2)Ri+p-3;
Сущность предложенного изобретения заключается в следующем.
На контролируемом участке горного массива устанавливают датчики, которые работают в режиме прозвучивания (излучающий и приемные). В качестве зондирующего сигнала используют отдельные одиночные прямоугольные акустические импульсы. Приемные датчики, подключенные к приемной аппаратуре, принимают сигнал, прошедший контролируемый участок горного массива. Принятый сигнал в своем спектре будет иметь особенности, т.е. характерные изменения параметров спектральной плотности. При каждом приеме импульса определяется отношение а
где A - коэффициент пропорциональный амплитуде принятого импульса;
k=tgβ, β - угол наклона фронтов принятого импульса;
τ - длительность принятого импульса;
Через промежуток времени t вновь излучают прямоугольный акустический импульс, измеряют его, и снова в спектре принятого сигнала будем иметь свои особенности. Особенности каждого такого спектра отражают физико-механическое состояние контролируемого участка массива горных пород, в том числе и изменение несплошности. Непосредственное сравнение спектров не дает полной уверенности в правильности выбора решения, ввиду недостаточной чувствительности и точности.
Для определения изменения несплошности контролируемого участка в массиве горных пород за время t, необходимо воспользоваться предложенным способом со следующей процедурой.
Во-первых, необходимо чтобы все излученные импульсы были приняты на приемной стороне. Во-вторых, длина трасс, по которым проходят импульсы, должна находиться в следующем соответствии:
Z1<Z2<Z3<…<Zl,
где Zi=V(τi+τu).
τi - время прохождения импульса от излучателя Is до i-го приемного датчика Di;
τu - длительность принятого импульса;
V - скорость распространения импульса в контролируемой среде.
Далее при первом излучении излучателем Is акустического импульса все приемное датчики принимают сигнал. Значит можно записать матрицу начального состояния А0:
По истечении времени t1 (t1>>τi+τu), снова излучают и принимают n импульсов по l трассами снова получают матрицу состояния
По истечении времени t2 получают матрицу состояния А2; По истечению времени t3 - А3 и т.д. В результате получаем матрицу Ach изменения состояния, т.к. в ней записаны все матрицы состояния А0, А1, А2 и т.д.
Для каждой строки матрицы состояния вычисляют
ri=bi(ni-1); ri+p=bi+p(ni+p-1);
γi=(ni-2)Ri-3; γi+p=(ni+p-2)Ri+p-3;
А для определения физических изменений, т.е. изменения трещиноватости контролируемого участка массива горных пород, необходимо сравнить в выбранных матрицах состояния соответствующие строки
где Fα - верхний α предел F-распределения со степенями свободы ri, и ri+P;
α - уровень значности.
Например, для определения изменения несплошности от t2 до t8, необходимо подставить в неравенство (1) соответствующие строки матриц состояния А2 и А8.
При выполнении неравенства (1) никаких существенных изменений несплошности на контролируемом участке массива горных пород не произошло. Если же указанное неравенство не выполняется, то за истекший период произошли изменения несплошности на контролируемом участке массива горных пород, которые могут быть зафиксированы неравенством (1) на выбранном уровне значности α.
Период повторения ti излучения импульсов должен выбираться исходя из скорости геомеханических процессов в контролируемом массиве горных пород. Чем выше скорость таких процессов, тем меньше период повторения ti и наоборот. Необходимость этого объясняется тем, что отслеживание изменения несплошности на ранней стадии существенно поможет выбрать прогностические характеристики для прогноза разрушения горного массива.
Длительность излучаемых импульсов выбирается одинаковой в сопоставляемых сериях и зависти от детальности контроля. Чем меньше контролируемый объект, тем короче выбирается длительность излучаемого импульса. Причем необходимо выбирать τu так, чтобы τu<<τi.
Положительный эффект заключается в использовании дифференциального подхода к определению изменения несплошности в массиве, т.к. контролируемый участок горного массива разбит на трассы, и изменение трещиноватости хотя бы в объеме одной трассы указывает на нарушение сплошности в контролируемом объеме массива горных пород.
В том числе выбор p для каждой трассы дает возможность менять чувствительность и определять значимо или незначимо произошли отличия между сравниваемыми временными интервалами.
Для сравнения предложенного способа и способа по [2] было проведено численное моделирование. Суть его в следующем импульсы, прошедшие контролирующий участок горного массива, имеют в своем спектре особенности, во временном понятии эти особенности выражаются в изменении переднего и заднего фронтов принятого импульса. Поэтому при численном моделировании учитывались изменения фронтов импульса с шагом один градус, т.е. прямоугольный импульс изменялся до трапеции, где боковые грани и «отвечали» за физико-механические изменения в контролируемом массиве. На фиг.1 приведены графики статистик по способу [2] линия 1, по предлагаемому способу - 2, а линия 3 это уровень выбранного классического F-критерия, по превышению которого считаются отличия значимы.
Из графиков видно, что предложенный способ чувствует различие практически с первого градуса (точнее со 2-го), в то время как способ [2] работает только с 12 градусов изменения фронтов и имеет провалы, которые отмечены на графике пересечением уровня критерия. Из приведенных вычислений можно судить о чувствительности предложенного способа.
ЛИТЕРАТУРА
1. Авт. св. СССР №1452984, 1989, БИ №3.
2. Патент РФ №2191411 от 20.10.02, БИ №29.
Claims (1)
- Способ контроля изменения несплошности в массиве горных пород, заключающийся в прозвучивании массива сигналами, приеме сигналов, прошедших контролируемой участок массива горных пород, отличающийся тем, что в качестве зондирующего сигнала используют серии отдельных одиночных прямоугольных акустических импульсов, дополнительно определяют энергию каждого принятого импульса в выбранных частотных интервалах и определяют отношение
где A - коэффициент, пропорциональный амплитуде принятого импульса;
k=tgβ, β - угол наклона фронтов принятого импульса;
τ - длительность принятого импульса;
q - интервал сравнения от i до i+p, где р=2,
а об изменениях несплошностей в массиве горных пород судят по невыполнению неравенства
Si,(i+p)<Fα(ri, ri+p),
где Fα - верхний α предел F-распределения со степенями свободы ri и ri+p;
α - уровень значности,
ri=bi(ni-1); ri+p=bi+p(ni+p-1);
γi=(ni-2)Ri-3; γi+p=(ni+p-2)Ri+p-3;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2012123762/28A RU2498353C1 (ru) | 2012-06-07 | 2012-06-07 | Способ контроля изменения несплошности в массиве горных пород |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2012123762/28A RU2498353C1 (ru) | 2012-06-07 | 2012-06-07 | Способ контроля изменения несплошности в массиве горных пород |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2498353C1 true RU2498353C1 (ru) | 2013-11-10 |
Family
ID=49683313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2012123762/28A RU2498353C1 (ru) | 2012-06-07 | 2012-06-07 | Способ контроля изменения несплошности в массиве горных пород |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2498353C1 (ru) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1645511A1 (ru) * | 1988-05-24 | 1991-04-30 | Предприятие П/Я Р-6542 | Устройство дл контрол изменени напр женного состо ни массива горных пород |
RU2046376C1 (ru) * | 1992-06-29 | 1995-10-20 | Малое предприятие "Геоакустик" | Способ оценки напряженного состояния горных пород |
US6098021A (en) * | 1999-01-15 | 2000-08-01 | Baker Hughes Incorporated | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method |
RU2191411C2 (ru) * | 2000-10-16 | 2002-10-20 | Вычислительный центр Дальневосточного отделения РАН | Способ контроля изменения напряженного состояния горного массива |
RU2192657C1 (ru) * | 2001-05-10 | 2002-11-10 | Вычислительный центр Дальневосточного отделения РАН | Способ контроля изменения напряженно-деформированного состояния горного массива |
RU2408035C2 (ru) * | 2007-12-11 | 2010-12-27 | Федеральное государственное унитарное предприятие "ВНИИОкеангеология" | Способ обнаружения вторичных нефтяных залежей |
-
2012
- 2012-06-07 RU RU2012123762/28A patent/RU2498353C1/ru not_active IP Right Cessation
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1645511A1 (ru) * | 1988-05-24 | 1991-04-30 | Предприятие П/Я Р-6542 | Устройство дл контрол изменени напр женного состо ни массива горных пород |
RU2046376C1 (ru) * | 1992-06-29 | 1995-10-20 | Малое предприятие "Геоакустик" | Способ оценки напряженного состояния горных пород |
US6098021A (en) * | 1999-01-15 | 2000-08-01 | Baker Hughes Incorporated | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method |
RU2191411C2 (ru) * | 2000-10-16 | 2002-10-20 | Вычислительный центр Дальневосточного отделения РАН | Способ контроля изменения напряженного состояния горного массива |
RU2192657C1 (ru) * | 2001-05-10 | 2002-11-10 | Вычислительный центр Дальневосточного отделения РАН | Способ контроля изменения напряженно-деформированного состояния горного массива |
RU2408035C2 (ru) * | 2007-12-11 | 2010-12-27 | Федеральное государственное унитарное предприятие "ВНИИОкеангеология" | Способ обнаружения вторичных нефтяных залежей |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2841561C (en) | Seismic geophysical surveying using a fibre optic distributed sensing apparatus | |
US9631972B2 (en) | Distributed fibre optic sensing | |
WO2008018061B1 (en) | A method and system for designating a target and generating target-related action | |
RU2014109546A (ru) | Способ определения уровня наполнения резервуара | |
US20180074019A1 (en) | Structure evaluation system, structure evaluation apparatus, and structure evaluation method | |
US11132542B2 (en) | Time-space de-noising for distributed sensors | |
RU2343502C2 (ru) | Способ и система определения положения наблюдаемого объекта по глубине в водной среде | |
EP3321673B1 (en) | Structure evaluation system and structure evaluation method | |
RU2510043C1 (ru) | Способ определения дальности до поверхности земли | |
CN114810213A (zh) | 用于煤与瓦斯突出的多源信息融合智能预警方法及装置 | |
US10458954B2 (en) | Structure evaluation system, structure evaluation apparatus, and structure evaluation method | |
RU2498353C1 (ru) | Способ контроля изменения несплошности в массиве горных пород | |
US7289388B2 (en) | Estimation of background noise and its effect on sonar range estimation | |
RU2631906C1 (ru) | Устройство для определения местоположения источника сигналов | |
US7239580B2 (en) | Noise adaptive sonar signal processor | |
RU90574U1 (ru) | Активный гидролокатор | |
Fialkowski et al. | Measurements and modeling of low-frequency near-surface scattering statistics | |
RU2480792C2 (ru) | Способ контроля изменения трещиноватости в массиве горных пород | |
RU2192657C1 (ru) | Способ контроля изменения напряженно-деформированного состояния горного массива | |
RU2191411C2 (ru) | Способ контроля изменения напряженного состояния горного массива | |
GB2463774A (en) | Radar system for detecting and analysing weather systems | |
US20230034759A1 (en) | System and method for acquiring marine seismic data using a vibrator | |
CN203551141U (zh) | 一种声波飞渡时间测量装置 | |
RU2288449C2 (ru) | Лазерный импульсный дальномер | |
RU2809469C1 (ru) | Способ и система сейсмоакустического контроля массива горных пород |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20140608 |