Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика деформирования и разрушения горных пород

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
8.26 Mб
Скачать

где 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,