RU2106493C1 - Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления - Google Patents

Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления Download PDF

Info

Publication number
RU2106493C1
RU2106493C1 RU95116415A RU95116415A RU2106493C1 RU 2106493 C1 RU2106493 C1 RU 2106493C1 RU 95116415 A RU95116415 A RU 95116415A RU 95116415 A RU95116415 A RU 95116415A RU 2106493 C1 RU2106493 C1 RU 2106493C1
Authority
RU
Russia
Prior art keywords
rock
invariant
rocks
well
bore
Prior art date
Application number
RU95116415A
Other languages
English (en)
Other versions
RU95116415A (ru
Inventor
А.А. Репко
В.Н. Рева
Ю.М. Карташов
А.В. Синицын
Original Assignee
Научно-исследовательский институт горной геомеханики и маркшейдерского дела
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Научно-исследовательский институт горной геомеханики и маркшейдерского дела filed Critical Научно-исследовательский институт горной геомеханики и маркшейдерского дела
Priority to RU95116415A priority Critical patent/RU2106493C1/ru
Publication of RU95116415A publication Critical patent/RU95116415A/ru
Application granted granted Critical
Publication of RU2106493C1 publication Critical patent/RU2106493C1/ru

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

Использование: для обеспечения структурной устойчивости горных выработок и безопасного выполнения подземных работ в горной и других отраслях промышленности. Сущность изобретения: способ оценки предельного напряженного состояния пород включает бурение скважины из горной выработки в зоне опорного давления, измерения радиальных и осевых смещений пород и интерпретацию замеров. Скважину бурят в зоне опорного давления относительно элементов залегания пласта, измеряют радиальные и осевые смещения пород, в качестве характеристик напряженного состояния выбирают первый инвариант тензора деформаций и второй инвариант девиатор деформаций после чего строят график зависимости первого инварианта тензора деформаций от времени и по точке перехода первого инварианта через максимум определяют момент времени наступления предельно о состояния горных пород, а соответствующее ему значение относительной энергетической прочности материала рассчитывают по определенной формуле. Устройство для измерения деформаций горных пород снабжено передвижным репером, дополнительным преобразователем для измерения смещений пород вдоль оси скважины, сердечник которого связан с торцом передвижного реперп через подвижную головку устройства и направляющим устройством относительно оси скважин, при этом передвижной репер размещен на оси скважины за корпусом устройства, а дополнительный преобразователь встроен в торце устройства, причем передвижной репер и направляющее устройство образуют базу устройства при определении деформации пород вдоль оси скважины. 3 ил.

Description

Изобретение относится к годному делу и может быть использовано для решения различных геомеханических задач, в частности, прогнозирования статической и динамической устойчивости горных выработок, пройденных, главным образом, в горных породах осадочного происхождения.
В практике исследований напряженно деформированного состояния горных массивов известен способ Лимана (см. например, Крупенников Г.А., Филатов Н.А., Амусин Б. 3, и др. Распределение напряжений в горных массивах. М.: Недра, 1972). Сущность этого способа состоит в определении деформаций на поверхности центральной скважины, возникших в результате полной разгрузки пород после обуривания центральной скважины коаксиально расположенной скважиной, и пересчете деформаций в напряжения. Кроме того, известен способ определения напряжений при помощи полной разгрузки пород по схеме Хаста (см. например, Крупенников Г.А., Филатов Н.А., Амусин Б.З. и др. Распределение напряжений в горных массивах. М., Недра, 1972). В отличие от способа Лимана в центральной скважине изменяют радиальные смещения стенок.
Недостатком указанных способов является сложная технология разгрузки керна с центральной скважиной. Оба способа применяются, главным образом, в скальных породах.
Наиболее близким к предлагаемому по технической сущности является способ определения напряжений в осадочных горных породах при помощи буровых скважин (см. Руководство по применению метода буровых скважин для определения напряжении в осадочных горных породах, АН СССР, Сибирское отделение института горного дела, Новосибирск, 1969). Сущность метода состоит в следующем. Из горной выработки бурят скважину, в которой устанавливают деформометр, ориентированный относительно элементов залегания пласта. Деформометр фиксирует радиальные перемещения нескольких точек контура скважин. Отсчеты снимают через определенные интервалы времени. Затем измеренные величины смещений пересчитывают в напряжения с помощью формул решения плоской задачи теории упругости с учетом свойств ползучести горных пород.
Трудность применения теоретических решений для определения компонент напряжений состоит в том, что для этого необходимо знать уравнение состояния пород, отражающее их фундаментальные свойства - упругость, пластичность и вязкость в реальных условиях деформирования пород. В то же время учет указанных особенностей предопределяет усложнение модели среды, и в математическом плане представляет существенно более сложную задачу по сравнению с плоской задачей теории упругости. Способ буровых скважин не позволяет определить предельные значения компонент тензора напряжений.
Для измерения деформаций пород известны различные устройства.
Наиболее близким к предлагаемому по технической сути является устройство для измерения деформаций горных пород (см. патент России N 1382955, E 21 C 39/00, 1986).
Устройство дает возможность измерять деформации пород только в двух ортогональных направлениях. Устройство содержит корпус, закрепленный в нем полый установочный стержень, штоки, установленные внутри стержня с возможностью осевого перемещения, упругие элементы, размещенные в корпусе, ползуны, связанные со штоками, и элементы передачи информации.
Недостатком устройства является отсутствие возможности измерять третью компоненту деформации пород (вдоль оси скважины).
Основная задача, которую решает предлагаемый способ - повышение точности определения напряженного состояния пород в горном массиве путем установления инвариантных (деформационных) характеристик пород и момента перехода объемных деформаций пород через максимум.
Устройство для реализации этого способа решает задачу расширения функциональных возможностей устройств для измерения деформации.
Согласно способу определения предельного состояния пород устанавливают инвариантные деформационные характеристики пород - первый инвариант тензора деформаций I1(Tε) и второй инвариант девиатора деформаций I2(Dε) , для чего из горной выработки в зоне опорного давления бурят измерительную скважину, в очищенную скважину помещают 3-х компонентный деформометр, ориентированный относительно главных компонент напряжений, и измеряют в определенный момент времени радиальные и осевые смещения пород, затем через установленный промежуток времени измерения повторяют, после чего строят графики зависимости первого инварианта тензора деформаций от времени и по точке перехода первого инварианта через максимум определяют момент времени наступления предельного состояния пород, а соответствующее ему значение относительной энергетической прочности материала рассчитывают по формуле:
Figure 00000002

где Cμ - параметр, зависящий от структуры горных пород, после этого производят оценку предельного напряженного состояния массива.
Сущность предлагаемого способа состоит в следующем.
Для решения поставленной задачи в породном массиве измеряют компоненты тензора активных деформаций некоторого элементарного объема среды, для которого выполняется условие λ ≪ L ≪ Z , где λ - масштаб естественной неоднородности массива (трещины, кристаллы и т.д.), Z - макроскопический масштаб задачи (размер зоны влияния выработки) и L - размер элементарного объема, принимаемый равным базе измерения смещения пород. По компонентам деформаций определяют I-ый инвариант тензора деформаций I1(Tε) и корень квадратный из II-го инварианта девиатора деформаций
Figure 00000003
. Указанные инварианты широко применяются в теории пластичности (Качанов Л.М. Основы теории пластичности., М.; Недра, 1969, с.418).
По физическому смыслу первый инвариант I1(Tε) выражает относительное изменение объема, а корень квадратный из второго
Figure 00000004
- относительное суммарное изменение формы малого элемента среды. По указанным инвариантам рассчитывают характеристики с более высокой степенью обобщения, а именно: относительное гидростатическое давление:
Figure 00000005

и относительную энергетическую прочность материала
P*= CμK 2 Б , (2)
где Cμ - параметр, зависящий от структуры горных пород, а именно от коэффициента поперечных деформаций материала. На основании теории размерностей можно показать, что при равенстве третьего (кубического) инварианта тензора деформаций I3(Dε) = 0 , представление KБ в формуле (I) единственно. Структура формулы (I) не изменяется при расчете KБ через компоненты напряжений (при пропорциональном нагружении). Относительное гидростатическое давление KБ изменяется oт
Figure 00000006
. При одноосном сжатии KБ = +1, при одноосном растяжении KБ = -1, а при чистом сдвиге KБ = 0. Чем больше условия нагружения элемента породы приближаются к равномерному всестороннему сжатию, тем "мягче" схема нагружения этого элемента, и, наоборот, чем больше условия нагружения приближаются к всестороннему растяжению, тем жестче схема нагружения. Указанный показатель не совпадает с известным в геомеханике понятием о "жестком" и "мягком" нагружении образца породы испытательной машиной. Показатель P*, называемый относительной энергетической прочностью материала в пределах малого объема породы, выражает соотношение между удельной потенциальной энергией изменения объема и удельной потенциальной энергией изменения формы рассматриваемого малого объема породы.
Определив показатель KБ, строят график зависимости первого инварианта тензора деформаций I1(Tε) от времени. Как показывают лабораторные испытания образцов пород, в условиях одноосного и многоосного нагружения функция I1(Tε) представляет собой немонотонную кривую с ярко выраженным перегибом, отражающим момент начала увеличения объема материала. По этой точке определяют момент времени t* наступления предельного состояния материала, а соответствующую величину относительной энергетической прочности рассчитывают по формуле (2). При этом под предельным состоянием материалов подразумевается возникновение в них заметных пластических деформаций и разрушения (см., например. Филин А.П. Прикладная механика твердого деформируемого тела, том. 1, М.; Недра, 1975, с.832).
Физическая сущность способа заключается в использовании свойств горных пород изменять свой объем соответственно характеру нагружения.
На фиг. 1 показана схема расположения измерительной скважины с деформометром относительно выработки; на фиг. 2 показан график зависимости первого инварианта тензора деформаций I1(Tε) элементарного объема породы от времени t; на фиг. 3 - устройство для осуществления способа (трехкомпонентный скважинный деформометр). Способ осуществляют следующим образом. На изучаемом участке (см. фиг.1) горного массива выбирают выработку 1 подверженную влиянию опорного давления. Из стенки выработки бурят скважину 2, ориентированную относительно элементов залегания пласта 3. Скважину тщательно очищают от буровой мелочи и пыли. В скважину досылают 3-х компонентный скважинный деформометр 4 и в намеченной точке массива измеряют смещения стенок скважины в трех ортогональных направлениях X, Y, Z, обусловленных влиянием прямого упругого последействия и других факторов. Интервал времени между отдельными измерениями определяют в зависимости от скорости смещения стенок. Измерения прекращают, когда соотношение между наименьшей и наибольшей компонентами смещений превышает единицу. Затем по измеренным смещениям рассчитывают компоненты деформаций пород, далее - первый инвариант тензора деформаций I1(Tε) и корень квадратный из второго инварианта девиатора деформаций
Figure 00000007
. После чего строят график (см. фиг.2) зависимости первого инварианта деформаций от времени I1(Tε) = φ(t) . Установив по указанному графику точку "А"^ перехода первого инварианта деформаций через максимум I1(Tε)max определяют момент t* времени наступления предельного состояния пород. Соответствующее значение относительной энергетической прочности материала определяют по формуле (2).
Примером реализации способа может служить оценка эффективности мероприятий по предупреждению пучения пород в окрестности горной выработки. Применяемый для этой цели традиционный способ измерения смещений пород при помощи глубинных реперов не позволяет получить информацию о напряженно-деформированном состоянии пород в глубине массива и требует больших затрат времени. В соответствии с предлагаемым способом в пучащем слое пород из выработки проходят при помощи станка для кернового бурения (типа НКР) скважину диаметром 59 мм. Длину скважины выбирают такую, чтобы скважина пересекла зону концентрации напряжений ("пик" опорного давления). Для получения сравнительно быстро, примерно, в течение одного часа, необходимой информации деформометр располагают в районе "пика" опорного давления, ориентировочно на расстоянии 1 м от стенки выработки. Затем фиксируют время и производят первоначальные отсчеты по трем датчикам скважинного деформометра. Относительно этих первоначальных отсчетов в дальнейшем определяют смещения пород в трех направлениях. При измерениях в относительно слабых породах (глинистые, углистые сланцы, соли и др.) интервал времени между измерениями принимают равными 5 мин. Полученные графики смещения пород во времени являются базой для расчетов деформаций и инвариантов
Figure 00000008
. Далее строят графики зависимости I1(Tε) = φ(t) . По этому графику фиксируют время t* перехода функции I1(Tε) через максимум I1(Tε)max . Соответствующее этому моменту величину показателя P* рассчитывают по формуле (2). Если в течение отведенного времени для измерений предельное состояние пород не наступило, т.е. график зависимости I1(Tε) = φ(t) оказался монотонным, то максимальную величину относительной энергетической прочности определяют по последнему измерению. После проведения мероприятий проходят по этому же пласту вторую скважину и повторяют циклы измерений смещений аналогично описанному. В заключение сравнивают максимальные величины относительной энергетической прочности пород и по ним оценивают эффективность мероприятий.
На фиг. 3 показано устройство для измерения деформаций горных пород.
Устройство содержит корпус 1, полый установочный стержень 2, закрепленный в корпусе 1 между перегородками 3 при помощи гаек 4, ползуны 5 и 6, размещенные на стержне 2 и выполненные в виде втулок с крестообразными пазами, упругие элементы 7 с гайками 8 регулировки, установленными на стержне 2 между ползунами 5 и 6 и перегородками 3, штоки 9, размещенные внутри стержня 2 и связанные посредством штифтов 10 с ползунами 5 и 6, стойки 11 и 12, закрепленные на стержне 2 параллельно соответствующим пазам ползунов 5 и 6 и повернутые относительно друг друга на 90o, направляющее устройство 13 передвижной репер 14, торец которого сопряжен с сердечником одного из трех преобразователей 15 линейных перемещений через подвижную головку 16 устройства при измерении продольных деформаций, два других преобразователя 15 взаимодействуют с концами штоков 9 при измерении поперечных деформаций, и с подвижной головкой 26, элементы передачи деформаций, выполненные в виде четырех пар рычагов 17, 18 (две пары рычагов, расположенные перпендикулярно плоскости чертежа, не показаны). Одни концы рычагов 17 и 18 соединены при помощи подшипников качения 19, предназначенных для контакта со стенками скважины, а вторые шарнирно прикреплены к концам стоек 11 или 12 и ползунам 6 или 5 соответственно. Один из пазов ползунов 5 и 6 образован щечками 20 с отверстиями для закрепления концов рычагов 17 и 18, во втором пазу размещаются стойки 11 и 12. Концы длинных 17 и коротких 18 рычагов двух пар, расположенных в одной плоскости, присоединены соответственно к стойке 11 и щечкам 20 ползуна 6, а концы рычагов аналогичных двух пар, расположенных в перпендикулярной плоскости, прикреплены соответственно к стойке 12 и щечкам (не показаны) ползуна 5.
Направляющее устройство 13 прикреплено к корпусу 1 и представляет собой систему, включающую восемь пар рычагов 17 и 18, соединенных со стороны стенок скважины подшипниками 19. Другие концы рычагов присоединены к стойкам 11 и подвижным ползунам 21. Подвижные ползуны соединены упругим элементом 22. Аналогичную конструкцию имеет передвижной репер 14. С последним контактирует головка 16, включающая упругий элемент 23, перегородку 24 и гайку 25 для регулировки. Головка 16 соединена с преобразователем 15 линейных перемещений.
Каждый преобразователь 15 линейных перемещений представляет собой индуктивный датчик с частотным выходом сигнала, состоящий из корпуса с установленной в нем катушкой индуктивности с выступающим сердечником, пружины и держатели (на черт. не показаны). Катушка с навесными электронными элементами залита компаундом.
Устройство работает следующим образом.
В скважине в начале при помощи специального досыльника устанавливается передвижной репер 14 и основное устройство до контакта с передвижным репером. Во время движения репера и устройства в скважине подшипники 19 рычагов 17 и 18 автоматически прижимаются к стенкам за счет усилия упругих элементов. При достижении заданной глубины в скважине передвижной репер и устройство фиксируются в определенном положении. Деформации скважины в поперечном направлении воспринимаются через подшипники 19 двумя парами рычагов. Короткие рычаги 18, укрепленные на стойке 11, поворачиваются относительно этой стойки, передают движение длинным рычагом 17, а последние толкают ползун 6. В свою очередь, ползун 6 через штифт 10 сообщает движение штоку 9, а шток воздействует на сердечник преобразователя 15 линейных перемещений. Одновременно аналогичное движение совершает другая пара рычагов (не показана), перемещая ползун 5 по установочному стержню 2.
При этом движение штока 9 передается другому преобразователю 15 линейных перемещений, установленному с другой стороны стержня 2.
Одновременно с перемещениями рычагов в поперечном направлении перемещается головка 16, передавая осевые смещения пород сердечнику преобразователя линейных перемещений 15. Смещения пород в осевом направлении происходит в результате деформации пород в интервале базы измерений, определяемой расстоянием между передвижным репером 14 и направляющим устройством 13. Движение ползунов 5 и 6 и соответственно пар подшипников 19 в противоположных направлениях вдоль оси скважины, а также движение головки 16, позволяет измерять деформации в скважине в трех взаимно ортогональных направлениях. Благодаря указанному увеличению информативности достигается расширение функциональных возможностей устройства.
Использование в качестве преобразователя линейных перемещений индуктивного датчика с частотным выходом сигнала позволяет производить измерения относительных деформаций с высокой точностью порядка 10-5 и передавать сигнал на расстояние до 2000 м.
Предлагаемый способ позволяет оценить напряженное состояние массива горных пород в естественных условиях без привлечения математического аппарата для решения объемных задач. Исходными данными для его реализации являются зависимости смещений пород в пространстве и во времени. Способ может быть использован для решения сложных геомеханических задач, как-то: определение структурной устойчивости горных выработок; расчет ресурса работоспособности предохранительных целиков; прогнозирование динамических явлений; оценка эффективности технологических приемов при сонации предприятий добывающей промышленности.
Для осуществления способа предлагается устройство, обладающее следующими достоинствами: возможность измерения смещений пород в трех направлениях; высокая точность измерения (1•10-5) и передача сигнала на расстояние до 2000 м; компактность.

Claims (2)

1. Способ оценки предельного напряженного состояния горных пород в объемном поле сжимающих напряжений горного массива, включающий бурение скважин из горной выработки, измерение смещений пород через определенный момент времени с помощью деформометра и интерпретацию замеров, отличающийся тем, что скважину бурят в зоне опорного давления относительно элементов залегания пласта, измеряют радиальные и осевые смещения пород, в качестве характеристик напряженного состояния выбирают первый инвариант тензора деформаций J1(Tε) и второй инвариант девиатора деформаций J2(Dε) после чего строят график зависимости первого инварианта тензора деформаций от времени и по точке перехода первого инварианта через максимум определяют момент времени наступления предельного состояния горных пород, а соответствующее ему значение относительной энергетической прочности материала рассчитывают по формуле
Figure 00000009

где Cμ - параметр, зависящий от структуры горных пород,
после этого производят оценку предельного напряженного состояния пород.
2. Устройство для измерения деформаций горных пород, содержащее корпус, закрепленный в нем установочный стержень, стойки, закрепленные на противоположных концах стержня и повернутые одна относительно другой на 90o, ползуны в виде размещенных на стержне втулок с двумя крестообразными пазами каждая, штоки, установленные внутри стержня и связанные с ползунами, преобразователи линейных перемещений, соединенные с торцами штоков, и элементы передачи деформации, связанные с ползунами, отличающееся тем, что оно снабжено передвижным репером, дополнительным преобразователем для измерения смещений пород вдоль оси скважины, сердечник которого связан с торцом передвижного репера через подвижную головку устройства, и направляющим устройством относительно оси скважины, при этом передвижной репер размещен на оси скважины за корпусом устройства, а дополнительный преобразователь встроен в торце устройства, причем передвижной репер и направляющее устройство образуют базу устройства при определении деформаций пород вдоль оси скважины.
RU95116415A 1995-09-26 1995-09-26 Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления RU2106493C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU95116415A RU2106493C1 (ru) 1995-09-26 1995-09-26 Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU95116415A RU2106493C1 (ru) 1995-09-26 1995-09-26 Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления

Publications (2)

Publication Number Publication Date
RU95116415A RU95116415A (ru) 1997-09-27
RU2106493C1 true RU2106493C1 (ru) 1998-03-10

Family

ID=20172295

Family Applications (1)

Application Number Title Priority Date Filing Date
RU95116415A RU2106493C1 (ru) 1995-09-26 1995-09-26 Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления

Country Status (1)

Country Link
RU (1) RU2106493C1 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2485313C1 (ru) * 2012-01-10 2013-06-20 Учреждение Российской академии наук Институт горного дела Сибирского отделения РАН Способ оценки напряженного состояния горных пород
RU2763565C1 (ru) * 2021-07-07 2021-12-30 федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский горный университет» Устройство для оценки напряжённо-деформированного состояния горного массива

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Руководство по применению метода буровых скважин для определения напряжений в осадочных горных породах", АН СССР, Сибирское отделение института горного дела, Новосибирск, 1969. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2485313C1 (ru) * 2012-01-10 2013-06-20 Учреждение Российской академии наук Институт горного дела Сибирского отделения РАН Способ оценки напряженного состояния горных пород
RU2763565C1 (ru) * 2021-07-07 2021-12-30 федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский горный университет» Устройство для оценки напряжённо-деформированного состояния горного массива

Similar Documents

Publication Publication Date Title
Kuwano et al. On the applicability of cross-anisotropic elasticity to granular materials at very small strains
CN103267678B (zh) 一种油气藏岩体力学地下原位模型恢复的方法及装置
CN103278389B (zh) 一种岩石动、静态弹性参数同步测量的方法
Yilmaz A new testing method for indirect determination of the unconfined compressive strength of rocks
CN103278614B (zh) 一种动静态岩石力学参数校正方法及装置
Viggiani Small strain stiffness of fine grained soils
Zhang et al. Deformation and failure characteristics of sandstone under uniaxial compression using distributed fiber optic strain sensing
Gstalder et al. Measurement of some mechanical properties of rocks and their relationship to rock drillability
US5345819A (en) Method and apparatus for wellbore stability analysis
Aziz et al. Mechanical rock properties estimation for carbonate reservoir using laboratory measurement: A case study from Jeribe, Khasib and Mishrif Formations in Fauqi Oil Field
Zhou et al. Empirical ratio of dynamic to static stiffness for propped artificial fractures under variable normal stress
RU2106493C1 (ru) Способ оценки предельного напряженного состояния горных пород и устройство для его осуществления
Voznesenskii et al. On the evaluation of rock integrity around mine workings with anchorage by the shock-spectral method
Benson et al. Modulus testing of rock at the Churchill Falls underground powerhouse, Labrador
Yu James K. Mitchell Lecture. In situ soil testing: from mechanics to interpretation
Davies et al. Determination of geomechanical properties of a typical Niger Delta reservoir rock using geophysical well logs
Teufel Acoustic emissions during anelastic strain recovery of cores from deep boreholes
Heuzé Sources of errors in rock mechanics field measurements, and related solutions
Nikolenko et al. Effect of discontinuities on elastic wave velocities in high-stress rock samples: Experimental research using ultrasonic interferometry
Haque et al. Development of Axial Load Transfer (TZ) Analytical Model for the PSC Piles
Koopmans et al. The effect of stress on the determination of deformation modulus
Amadei et al. Methods of in situ stress measurement
Muralha et al. Laboratory and Field Testing of Rock Masses for Civil Engineering Infrastructures
Naeimipour et al. Applications of rock strength borehole probe (RSBP) in underground openings
Asslan et al. Factors Influencing Small-Strain Stiffness of soils and its Determination