книги / Механика деформирования и разрушения горных пород
..pdfгде a = 2sin р/(1 —sin р);
К0 = 1/ [ос(1 —m pJ] [2mp\ipC1/Kp — ainppa — 2Kcos p/sin р].
Упругая область. В упругой области материал скелета деформируется по закону Гука. Запишем уравнения теории упругости для пористой среды, используя значения средних r fJ напряжений
еу- е 8У= 1/(2С)(Гу-Г8у), |
(4.66) |
|
где Г = (1 — т)а — тр\ а=\-/3сц; е= АТ; |
K = ( \ —2v)/E-, |
v; Е — |
соответственно коэффициент Пуассона |
и модуль упругости. |
Сучетом зависимостей (4.59), (4.66) компоненты напряжений
вскелете равны
|
аг c = {Pi ± С 4/г2)—\аетс1(\ — те)С\пг, |
(4.67). |
|
tfe J |
|
где |
С4— произвольная постоянная. |
|
Радиус зоны пластических деформаций определяется из |
||
условий непрерывности напряжений на границе |
раздела |
|
двух |
сред: |
|
|
(1 -п }р)(ав + а г) |г=Г| =(1 - /и е)р, - \х етеС\пг0. |
(4.68) |
Учитывая напряжения (4.65), уравнение (4.68) запишем в виде
( \ - т е)(2 + а.)[ро + К0)(г01а)°-Ко] -
—amppa + 2Kcos р/( 1—sin р) +
+ трцрКе/{Крце){рь- р а)( 1 + 2/а + In ( r j a ) ) / / ( ^ Кр)In {rja) -
- In (гJ b)] =(1 - m e)pl - \ i eme(pb- p a)\nr0j [\ieKe/(\ieKp)ln {rja)~
-ln (r0/6)]. (4.69)
Упругопластическое распределение напряжений при нерав нокомпонентном поле напряжений в нетронутом массиве. В общем случае распределение напряжений в нетронутом массиве будет отличаться от гидростатического. Горизонтальная составля ющая напряжений при действии гравитационных сил обычно меньше вертикальной. Однако в зоне действия тектонических сил горизонтальная составляющая может превосходить вер тикальную.
Обобщим приведенное выше решение задачи на случай неравнокомпонентных напряжений в нетронутом массиве. Поста новка задачи состоит в следующем. В массиве пород сооружается выработка или производится бурение скважины круглого сечения радиуса г = а. Массив пород предполагается упругим и изотроп ным, в котором под давлением находится газ. Под действием напряжения и давления газа около контура выработки или скважины образуется область предельных деформаций.
Напряженное состояние рассматривается в плоскости, пер пендикулярной к оси выработки или скважины. Тогда задача сводится к исследованию напряженного состояния в газона сыщенной бесконечной плоскости с отверстием. При этом граничные условия на контуре скважины имеют вид (4.55) и (4.57).
На бесконечности для газонасыщенной плоскости они
записываются |
|
|
<*?'=PÙ |
o ^ = XPl-, т£> = 0, |
(4.70) |
где р и X— соответственно |
величина гравитационной |
состав |
ляющей напряжений и коэффициента бокового распора в не тронутом массиве.
Как и выше, в первом приближении характеристики пласта т, К, |1 в упругой и пластической областях принимаются постоянными.
Газовая фаза. Изменение давления вокруг выработки или скважины, как и выше, характеризуется зависимостью (4.59).
Компоненты напряжений в скелете в пластической области
определяются |
формулами (4.65). |
|
Условия на границе L раздела упругой и пластической |
||
областей записываются в виде |
|
|
|
4Re<l>*(z) = A l (zz)a/2 + Bl + Dl In (zz/a); |
|
2 |
[гФ* (z)+ 'P* (z)] = A 2(z/z) (zz)“/2 + B2z/z, |
(4.71) |
где A i =(2 + a)(p0 +A’oJ/a®;
=[2ATcosp —2mppasin p)/[(l —mp)(l —sin p)] —/£0(a+2) —
-m y C i/[ ( 1— mp)Kp~\(2+ a)/a;
Dx = 2mp\xpCxl\(\ —Я1Р)ЛТР];
A 2 = v.{p0 + K0)la«;
B2 = — K0<x — ctmp/(\ —mp)pa + 2Kcos p/(l —mp){\ —sin p)— —w p|ipC1/(l —mp)Kp\
Ф*(г), T*(z) — функции комплексного переменного, .через ко торые выражаются компоненты напряжений в упругой области.
Для их определения имеем граничные условия (4.71) на неизвестном контуре L и условия на бесконечности (4.70). Отобразим контур L области предельных деформаций на единичный круг у, с помощью функции
X |
|
z= co(p) = c0<p + X о/ф"’ |
(4.72) |
J=I |
|
где с0, cj — неопределенные коэффициенты; |
ср — комплексная |
переменная. |
|
Тогда получим
4Re<D(ф)= Bi + A t[to(ф)со (ф) ] 0 / 2 + ZMn (ю(ф)ш(ф)/а,); (4.73) 2[Ф'(ф)/оУ(ф)юф + \1>(ф)] = Л2 [со (ф)]в / 2 “ 1 [со (ф)]а / 2 + ,+
+ 52 со(ф)/со(ф).
Решение граничной задачи (4.73) при условиях на бес конечности (4.70) построено [32]. Область пластических дефор маций имеет формулу эллипса:
со(ф)= с0 ф-(-с1 /ф, |
(4.74) |
где коэффициенты с0 и сх находятся из соотношений:
(1 +A.)/?i = A I C Q -\T B i +2Z)ilnc0;
л |
(l - ^)P i = A2c1ce0- 1+ 2B2clc ô l. |
(4.75) |
Приведенное решение имеет смысл, если область предельных деформаций полностью охватывает контур выработки, для чего необходимо выполнение условия с0—с, ^ 1.
Используя зависимости (4.59), (4.62) и (4.63), распределение давления газа р в пластической и упругой зонах можно записать в следующем виде:
Р=Ра + |
\ipK‘(pb-p„)\n(rla) |
(4.76) |
|
(цр* 7 { \ i ‘ K p ) In (rja)- In (rjb))’ |
|||
|
|
||
Р=Рь + |
iPb-Pa) |
(4.77) |
|
PpKel(\il'Kp) ln(rja)— \n(rjb) |
Запишем зависимости для определения дебита скважины. При упругопластическом деформировании массива пород значе
ние дебита Q при |
известной скорости W, может быть |
легко |
найдено: |
|
|
|
2пК‘(рь- Р') |
(4.78) |
|
Qp= n'ln (rjb)—\ipKeIKp\n (гJа) |
|
При ц.р = |1с и |
К е = К р из выражения (4.78) следует |
зави |
симость |
|
|
|
2пК‘(рь-ра) |
(4.79) |
|
ц'1п alb |
|
|
|
совпадающая с решением для дебита при упругом дефор мировании пород.
Пример. Найдем размер зоны предельных деформаций вокруг вертикальной
скважины, пересекающей газоносный пласт на глубине |
Н = |
1 км. |
Исходные |
|||||
данные для расчета: уЯ=25 МПа; Х=1; Ь / а = |
20; /?«,= 0; |
К |
с = |
2 |
МПа; т р = 0,2; |
|||
ш*=0,02; |
ра= 10 МПа; |
/>* = 20 МПа; \ip=\xe; |
K e =e ' i 0 ; |
K |
p |
= e 2 S . |
Сцепление |
|
пористого |
массива К = |
( \ — т р ) К с . |
|
|
|
|
|
|
Расчет по формуле (4.69) показывает, что радиус зоны предельных деформаццй вокруг скважины без учета и с учетом ослабляющего влияния давления газа составляет го = 2,04я и /*о = 4,0д. Таким образом, в результате ослабляющего давления газа радиус зоны предельных деформаций увеличива ется в 1,96 раза.
Оценим изменение дебита газа в скважине для проведенных выше горно-геологических и горно-технических условий заложения скважины. Расчет по формулам (4.78) и (4.79) показывает, что наличие предельных деформаций радиуса зоны г0—\а увеличивает дебит газа в Qp= \£6Q e.
Приведенные методы расчета могут быть использованы для прогнозирования напряженно-деформированного состояния массива вокруг выработок и скважин, а также обоснования технологических решений при добыче газа и нефти и решения проблем гидрогеологии.
5. УПРУГОПЛАСТИЧЕСКОЕ РАСПРЕДЕЛЕНИЕ НАПРЯЖЕНИЙ ВОКРУГ ВЫРАБОТКИ С УЧЕТОМ АНИЗОТРОПИИ ПРОЧНОСТНЫХ
СВОЙСТВ ГОРНЫХ ПОРОД |
' |
5.1.ПРОЧНОСТЬ АНИЗОТРОПНЫХ ГОРНЫХ ПОРОД
Горные породы в зависимости от горно-геол<ргических условий залегания, характера расположения и числа трещин, слоистости, наличия различного рода включений и ряда других факторов проявляют свойства упругой и пластической анизот ропии.
Наиболее существенным фактором, влияющим на прочност-
ьые и деформационные характеристики массива и процесс его деформирования, является трещиноватость горных пород.
По характеру происхождения различают естественную и ис кусственную трещиноватость. Естественная трещиноватость связана с особенностями образования пород, а искусственная формируется вследствие принятой технологии добычи полезных ископаемых и сооружения выработок, процесса деформирования массива и влияния внешней среды.
Естественная трещиноватость, в зависимости от ориентации трещин, может быть упорядоченной или хаотичной. Наличие трещин в массиве приводит к его ослаблению й снижению прочностных и деформационных характеристик. Прочностные характеристики породы в массиве будут значительно меньше, чем в куске.
Другим фактором, существенно влияющим на прочностные и деформационные свойства массива, является слоистость. Обычно различают микрослоистость, т. е. слоистость в преде лах одной литологической разности, и макрослоистость,
образующуюся вследствие наличия в массиве различных ли тологических разностей.
Микрослоистость и упорядоченная естественная трещино ватость приводят к появлению у породы анизотропии меха нических свойств.
Испытания образцов горных пород при одноосном сжатии на нагрузку, направленную перпендикулярно и параллельно слоисто сти, показывают, что прочность их в этих направлениях различна. Так, в табл. 5.1 приведены [16] экспериментальные данные по определению прочностных и деформационных характеристик пород при одноосном сжатии и растяжении их перпендикулярно и параллельно слоистости. Для испытаний использовались породы рудных месторождений различной степени устойчивости, которые охватывали все основные типы пород. Испытания выполнены для устойчивых пород, средней устойчивости, неустойчивых и группы весьма неустойчивых пород.
Анализ этих данных показывает, что пределы прочности на одноосное сжатие перпендикулярно и параллельно сло истости зависят от типа пород. Коэффициент анизотропии, равный отношению предела прочности на сжатие перпен дикулярно слоистости к пределу прочности параллельно сло
истости, |
для |
пород средней устойчивости изменяется |
от 2,5 |
до 1,77, |
а для неустойчивых от 2,66 до 2,57. Значение |
этого |
|
коэффициента |
для весьма неустойчивых пород равно |
3. |
Анизотропия горных пород сказывается также на величине статического модуля упругости, полученного при испытаниях пород параллельно и перпендикулярно слоистости. На дина мический модуль упругости влияние слоистости сказалось только для группы устойчивых пород. Значения коэффициента Пуассона оказались независимыми от слоистости пород.
Предельные свойства пластически анизотропных пород и материалов в настоящее время изучены недостаточно, имеются только отдельные предложения по выбору тех или других условий пластичности.
Условие пластичности анизотропной среды И. И. Гольден-
блат и В. А. Копнев (1968 г.) предлагают принимать |
в виде |
h + y /Ï 2 = U |
(5.1) |
где Il = Cijo ij; /2 = Сд аа |-/-ал/— линейная и квадратичная формы
напряжений, записанные в |
тензорной форме; С0, Cijki— неко |
||
торые постоянные. |
|
|
|
Попытки обобщения условий пластичности Мизеса и Кулона |
|||
на анизотропные среды |
приведены |
в работе |
Д. Д. Ивлева |
(1966 г.). |
анализе |
поведения |
слоистых сред |
В. В. Соколовским при |
предложено наряду с обычным условием предельного состояния вводить специальное условие, учитывающее их слоистость.
На и м е н о в а н и е
ка т его р и и
усто й ч и в о с ти
Устойчивые, 1
Средней устой чивости, 2
Неустойчивые,
3
Весьма неус тойчивые, 4
На п
ра в л е
ние |
П р е д ел п р о ч |
|||
уси л и й |
|
н о сти , |
М П а |
|
1 |
160—330 |
|||
il |
250 |
|
|
|
1 |
100— 160 |
|||
il |
40—90 |
|
||
1 |
40—90 |
* |
||
il |
15—35 |
|
||
1 |
30 |
и |
ниже |
|
il |
1 0 |
и |
ниже |
О д н о о сн о е |
сж ати е |
|
|
М о д у л ь у п р у го |
К о эф ф и ц и ен т |
||
сти |
£ 1 0 4, М П а |
П у а сс о н а v |
|
7,5— 11 |
|
0,18—0,29 |
|
|
— |
|
0 , 2 2 |
(N 40 |
г*^ 00 |
|
0,24—0,26 |
1 |
|
|
0,28—0,29 |
|
|
|
|
5,4—6,1 |
|
0,27—0,28 |
|
|
— |
|
0,26—0,34 |
Не |
более |
2,8 |
0,3 |
Не |
более |
0,6 |
0,35 |
|
О д н о о сн о е р астя ж ен и е |
|
||
П р е д ел п р о ч |
М о д у л ь у п р у го |
Д и н ам и ч еск и й |
5 |
|
м о д у л ь у п р у го |
|
|||
н о сти , М П а |
сти |
£ 1 0 4, М П а |
сти Е -104, М П а |
|
1 0 — 2 0 |
4,9—6,5 |
0,84—0,97 |
0,14 |
|
7,5 |
8 ,6 |
|
1,4 |
|
2,5—3,5 |
1,4—2,4 |
0,75 |
0,43 |
|
5— 10 |
3,8—4,1 |
0,75—0,83 |
|
|
_ |
0 ,8 |
—2 , 2 |
0,62—0,67 |
0,45 |
— |
2 ,1 |
—4,1 |
0,55—0,73 |
|
1,5—5 |
0,7 |
_ |
0,24 |
0,5 |
1,5 |
|
0,23 |
|
Основываясь на приведенных выше экспериментальных данных, а также анализе исследований, перейдем к выбору условий пластичности при пластическом деформировании ани зотропных горных пород. Анализ экспериментальных данных показывает, что предел прочности на сжатие а сж зависит от направления слоистости пород по отношению к действующей нагрузке
|
|
|
Осж = стс*(0), |
(5.2) |
||
где |
0— угол |
между |
слоистостью |
образца и |
направлением |
|
действия нагрузки. |
|
|
|
|
||
|
Тогда сцепление К горных пород также зависит от коор |
|||||
динаты 0, поскольку |
оно |
связано |
с а сж зависимостью |
|||
|
|
K=K(Q)=(\ - |
sin р)/[(2cos р) стсж (0)], |
(5.3) |
||
где |
р — угол |
внутреннего |
трения |
горной породы. |
||
|
Обобщим наиболее распространенные условия пластичности |
горных пород на анизотропные среды, учитывая, что предел
прочности их ' |
зависит от направления приложения нагрузки |
по отношению |
к слоистости пород. |
На основе большого числа экспериментальных исследований для описания предельного состояния горных пород в работе [32] предложено экспоненциальное условие пластичности, ко торое можно обобщить на анизотропные среды следующим
образом: |
|
х = тЦ 0)еЛ |
(5.4) |
где т= 1/2 (cfj- а з ); |
c = a 3/<Si; |
х°п (0) = стсж (0)/2; |
Л — параметр, |
значения которого |
приведены |
в вышеуказанной |
работе. |
При практических расчетах получило распространение усло вие пластичности Кулона в форме наклонной линейной оги бающей кругов главных .наибольших напряжений. С учетом связи (5.3) обобщим это условие на анизотропные среды,
которое для плоской деформации можно записать так |
|
(ae-CTr)2 + 4T?e= sin2 р [ств + <тг + 2ЛТ (0) ctg р] 2. |
(5.5) |
Для некоторых горных пород может быть использована также горизонтальная огибающая кругов напряжений, уравне ние которой для анизотропных сред имеет вид
(ств —а г)2 + 4т2в = 4 £ 2(0). |
(5.6) |
При расчете пластического деформирования пород вокруг выработок, сооружаемых в массивах с различными пределами на сжатие параллельно и перпендикулярно слоистости, не обходимо учитывать угол между действующими главными
напряжениями и направлением слоистости. Для учета этого фактора в работе [16] предлагается использовать паспорт прочности выработки, с помощью которого учитывается на правление действия главных напряжений и расположение пределов прочности пород в зависимости от их слоистости.
Для аппроксимации сцепления пород около контура вы работки, сооружаемой в слоистом массиве, используем зави симость
|
К(в) = К |
|
т |
(а„ cosnQ + b„ sin«0) |
} |
|
(5.7) |
|
|
1+5 £ |
|
||||||
|
|
п = 1 |
|
|
|
|
||
где а„, Ьп— коэффициенты; |
0— угловая |
координата; 5 — коэф |
||||||
фициент, характеризующий |
степень анизотропии |
пород. |
|
|||||
В первом приближении в зависимости (5.7) можно огра |
||||||||
ничиться двумя слагаемыми |
|
|
|
|
|
|||
|
|
K{Q) = K(\+bcos2Q). |
|
|
(5.8) |
|||
Обозначим сцепление при сжатии пород перпендикулярно |
||||||||
слоистости через К±, а параллельно К^. |
|
|
|
|||||
Тогда из зависимости (5.8) имеем |
|
|
|
|
||||
|
|
Л:(0) = /^(Ц-0) = /<Гх; |
|
|
(5.9) |
|||
|
|
АГ(л:/2) = АГ(1 —5)= AT,, |
|
|
||||
|
|
|
|
|
||||
Отсюда |
следует |
|
|
|
|
|
|
|
|
K = (*i + *ll)/2; |
6 = (* х -* ц )/(*! + *„). |
|
(5.10) |
||||
Если ввести в рассмотрение коэффициент анизотропии |
||||||||
прочностных свойств |
|
|
|
|
|
|
||
|
|
* . = *±/*11, |
|
|
|
(5.11) |
||
то из соотношения |
(5.10) получим |
|
|
|
|
|||
|
*=(*± + *ц)/2; |
8=(Кв—1)/(*0+1). |
|
|
(5.12) |
|||
Таким образом, введенные выше коэффициенты К и 5 за |
||||||||
висимости |
(5.10) имеют |
физический |
смысл. |
|
Коэффициент |
|||
К характеризует значение среднего сцепления породы, |
а 5 — |
|||||||
разброса, вызванного анизотропией прочностных свойств. |
||||||||
Значения |
коэффициента |
5, |
вычисленные по формуле |
(5.10), |
и прочности (по Н. Ф. Ренжиглову) некоторых горных пород Восточного Донбасса приведены в табл. 5.2.
Анализ данных (см. табл. 5.2) показывает, что закономер ности распределения пределов прочности для угольных ме сторождений являются такими же, как для рудных.
|
|
|
Предел прочности |
на одноосное |
|
№ |
Литологическая |
сжатие, |
МПа |
Значения |
|
|
|
||||
пробы |
разность |
перпендикулярно |
параллельно |
Ô |
|
|
|
|
|
||
|
|
|
слоям |
слоям |
|
1 2 |
Песчанистый сланец |
662 |
527 |
0,13 |
|
13 |
Глинистый |
сланец |
470 |
284 |
0,25 |
14 |
|
» |
615 |
537 |
0,07 |
15 |
|
» |
322 |
250 |
0,09 |
16 |
Песчано-глинистый сла |
6 6 6 |
394 |
0,26 |
|
17 |
нец |
сланец |
279 |
170 |
0,24 |
Глинистый |
|||||
18 |
|
» |
123 |
59 |
0,35 |
19 |
Песчанистый сланец |
390 |
360 |
0 , 2 |
|
2 0 |
Глинистый |
сланец |
264 |
225 |
0,08 |
2 1 |
|
» |
566 |
460 |
0 ,1 |
2 2 |
|
» |
460 |
330 |
0,16 |
5.2. УРАВНЕНИЯ ПЛОСКОЙ ДЕФОРМАЦИИ ПЛАСТИЧЕСКИ АНИЗОТРОПНЫХ ГОРНЫХ ПОРОД
При исследовании процессов деформирования горных пород и сплошных тел большое распространение получила задача о плоской деформации. Теория плоской деформации является одним из наиболее разработанных разделов теории пластич ности горных пород и сплошных тел.
Рассмотрим уравнения пластического деформирования среды при плоской деформации с учетом пластической анизотропии горных пород. Для плоской деформации условие предельного , состояния и два уравнения равновесия образуют замкнутую
систему, содержащую три неизвестных.
Рассмотрим вначале систему при экспоненциальном условии (5.4) предельного состояния. При плоском деформированном
состоянии |
условие (5.4)- отождествляется |
подстановкой |
вида |
а * > = ЛГ(ф) [Z)(c)±fi(c)cos2(p]; zxy= K((fi) £(c)sin2cp, |
(5.13) |
||
а,, J |
|
|
|
где В (с), |
D(c) — функции параметра |
с; К (ф) — параметр, |
характеризующий анизотропию в условии предельного состо яния; (р — угол между первым главным напряжением и осью Ох.
Для |
условия пластичности |
|
|
|
|
|
'г = х" (ф) еХс |
(5.14)' |
|
значения переменных В, |
D, К, |
с равны |
|
|
|
В (c) = exp(/lc); |
D (с) = ехр(Лс) (1 + с)/(1 —с); |
|
|
|
AT(<p)= xS (ф); |
с = а 31а1. |
209 |
|
14 Заказ |
3356 |
|
|
Подставив соотношения (5.13) в уравнения равновесия, получим ^
K(D' + B' cos2(pdcjôx + KB' sin2фдс/ду+[К' (D+
+В cos 2ф) — 2KB sin 2ф] 8ф/дх + В(2К cos 2ф +
+К' sin 2ф5ф /ôjn= Л';
КВ' sin2<pdc/dx+K(D' — B' cos 2ф) дс/ду + В (2К cos 2ф +
+ К' $1п2фЗф/5л' + [5 (2ATsin 2ф —К' со$2ф +
|
+ DK']d(f>/Ôy=Y, |
|
(5.16) |
||
где Z)' = dZ>/dc; B' — dB/dc, |
X, |
Y — массовые |
силы. |
вид |
|
Уравнения |
характеристик |
этой системы имеют |
|||
dу |
К' (BD' - B'D) cos 2ф - 2 KBD' sin 2ф ± |
|
|||
d* |
2KB{B' + D' соз2ф)+ |
* |
|||
i ± V 4 К 2В2 [{D')2- ( B ')2] + {K')2 { B ' D - B P f |
|||||
Очевидно, |
+K'(BD' + B'D) ып2ф |
|
‘ 1 ' |
||
если выражение, |
стоящее под |
знаком |
радикала |
в связи (5.17) больше нуля, то система уравнений (5.16) является гиперболической. В дальнейшем будем предполагать,
что это условие выполняется. |
урав |
|||
|
При |
отсутствии |
анизотропии ЛТ(ф) = К = const система |
|
нений |
упрощается. |
Этот случай системы исследован |
ранее |
|
и |
приведен в работе [32]. |
|
||
с |
По |
аналогии с |
плоской в случае осесимметричной задачи |
|
использованием |
условия Хаара — Кармана может |
быть |
проведено исследование системы уравнений теории предельного равновесия.
Перейдем к рассмотрению системы уравнений теории преде льного равновесия анизотропной среды при условии пластич
ности |
Кулона |
|
|
{cre- o 'f)2 + 4T?e = sin2 р [аг + ств + 2/:(ф) ctgр]2. |
(5.18) |
При плоском деформированном состоянии условие преде |
||
льного |
состояния (5.18) отождествляется подстановкой |
вида |
* > = с т ± [ а sin р + X (ф )C O S р] сов2ф;
СУу J
(5.19)
Тду = [о sin р+ К (ф) cos р] sin 2ф,
где ст= 1/2 (а* + а у).
Внося выражение (5.19) в уравнения равновесия, получим систему
Ü\ = da/dx+a2dcyldy + a3d(p/dx + a4.ô(p/dy — 0;
=dG/ôx + b28cy/dy + b2d(pldx+b^d(pldy = 0,