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

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

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

Минеральные

трещины

 

Открытые

трещины

 

ст2, МПа

 

 

 

 

 

b, мкм

Гм, 1/м

h, мкм

7'0,

1/м

мД

0

5

5

10

5

25

5 —7

50

5— 10

0

5

5

5

10-

5

25

5—7

50

5 —10

Невыбросоопасные песчаники

_

2—5

20—30

15—

18

30—40

5—

10

4050

10

50 - 60

1-Г 7 о

Выбросоопасиые песчаники

5— 10

5—7

5— 10

18—20

10— 15

5—10

15 20

10—15

2 0 - 25

15 18

 

1

0,01

23

2,1

1

2

1,8

2 —4

3,4

 

5

4,2

5— 10

1

 

30

3,5

5—7

2,1

 

10

4,3

10— 12

6,4

Таким образом, в условиях неравномерного объемно-на­ пряженного состояния в ВО и НВО песчаниках отмечается развитие остаточных деформаций разного механизма. Остаточ­ ная деформация ВО песчаников сопровождается переориен­ тировкой глинистых минералов цемента и обломочного матери­ ала за счет поворота, а иногда — изгиба отдельных пластинок слюд. Переориентировка способствует образованию и росту межзерновых трещин, разобщающих зерна и цементирующую массу, их раскрытию и разрыву в направлении продольного напряжения.

Остаточная деформация НВО песчаников начинается с двойникования кристаллов кальцита, слагающих цемент этих пород, и развития внутрикристаллических трещин в породообразую­ щих зернах. В результате двойникования зерен цемента образу­ ются каналы; nopbf и разрывы сплошности зерен. Внутрикристаллические трещинки развиваются в основном вдоль двойни­ ковых швов под углом к ним и иногда хаотически. Как правило, перечисленные микротрещины короткие, волосяные, заполненные в ходе деформации тонкодисперсным веществом.

Дальнейший рост и раскрытие внутрикристаллических тре­ щин приводит к развитию межзерновых и секущих зерна

трещин,

которые

и увеличивают трещинную проницаемость

на два

порядка

по сравнению с исходной породой.

В случае ВО

песчаника, проницаемость которого высока,

в зону опорного давления поступает большое количество газа. Необратимое деформирование массива ВО песчаника несущест­ венно увеличивает его проницаемость, дегазации зоны опорного давления не происходит* и газ может совершать работу по выбросу нарушенной деформацией части массива.

Результаты по фильтрации песчаников помогают в какой-то мере разобраться в механизме явлений выброса. ВО песчаники имеют более высокую проницаемость и соответствующую газоотдачу. При механическом его деформировании силами горного давления проницаемость не очень изменяется и не приводит к сильной дегазации деформированной области пласта, а приток газа в эту область по-прежнему остается большим. Деформированный объем пласта, потеряв свою прочность, особенно прочность на отрыв, и постоянно вы­ водимый из массива газ могут существенно способствовать возникновению разрушения пласта.

В случае НВО песчаника процесс может быть представлен несколько иначе. Начальная проницаемость низкая, газоотдача также низкая. Деформированный объем пласта резко увеличи­ вает свою проницаемость, пласт интенсивно дегазируется, подводимый из массива газ (а его подводится из-за более низкой начальной проницаемости НВО песчаника значительно меньше по сранению с ВО песчаником) успевает уходить и в деформационной зоне пласта не способен создавать давлений, необходимых для разрушения.

4.3. УСЛОВИЯ ПРЕДЕЛЬНЫХ СОСТОЯНИЙ НАСЫЩЕННЫХ ПОРИСТЫХ ГОРНЫХ ПОРОД

При горных ударах, выбросах угля, соли и газа происходит разрушение газонасыщенных массивов горных пород. Ис­ следование процессов деформирования и разрушения пористых насыщенных пород имеет большое значение для выявления закономерностей формирования и протекания динамических проявлений горного давления, выбора рациональных способов добычи нефти и газа, изучения механизма разрушений в земной коре при землетрясениях и решения других геомеханических задач.

Перейдем к рассмотрению предельного деформирования насыщенной пористой среды.

Двухфазная среда. Насыщенную газом (жидкостью) по­ ристую среду с точки зрения механики сплошной среды можно рассматривать как двухфазную. Одной из фаз явля­ ется газ (жидкость), а другой — твердые частицы скелета. При этом все пространство элементарного макрообъема запол­

нено двумя сплошными

средами, которые проникают

друг

в друга.

микрочастицы

уложены случайным

В пористых средах

образом, поэтому в

элементарном

макрообъеме

AV=

= Дх1Д.т2Дл'3 имеют место случайные тензорные поля. Струк­ туру двухфазной горной породы будем характеризовать от­ носительным объемным содержанием фаз. Пористость среды

с хаотической структурой может быть определена на основе усреднения по объему AV

W

где Х ( М , со)= 1, если точка М принадлежит поровому простран­ ству; л(М , со)= 0, если точка М принадлежит скелету; со харак­ теризует случайный характер принадлежности точки М одной или другой фазе.

Усреднение может быть выполнено по объему AV или плоскости AS. Будем рассматривать в дальнейшем массивы насыщенных пород в виде однородных изотропных сред.

Построение уравнений состояния твердой фазы для двухфаз­ ной среды осложняется наличием двух систем напряжений, одна из которых представляет собой гидростатическое сжатие

сплошного материала

под

действием

порового

давления,

а вторая — деформацию

скелета

имеющимися в нем

напряже­

ниями. Средние напряжения

в

скелете

обозначим

через а,;,

а полные средние напряжения в двухфазной среде через Ги. Связь между напряжениями в этом случае можно записать так:

Г0= (1—«) CT,V—«р50,

(4.7)

где п — средняя пористость среды; р — давление газа (жидкости). Это соотношение является фундаментальным и выполняется для любой двухфазной среды. Нужно отметить, что в связи (4.7) учитывается объемное соотношение фаз, в котором

твердая фаза

составляет 1—и и

газовая (жидкая)

п частей.

Для описания предельного состояния однородных однофаз­

ных горных

пород

наибольшее

распространение

получила

в настоящее время

теория прочности Мора

 

 

 

т=/(ст),

(4.S)

где'

 

 

 

 

 

T= 0,5(G !—ст3); а = 0,5(ст!+ ст3).

(4.9)

Заметим, что в случае прямолинейной огибающей условие

(4.8) в главных напряжениях записывается

 

 

Gi —а 3 = (а1+ a 3)sin p-h2/Ccos р,

(4.10)

где К, р— соответственно сцепление и угол внутреннего трения изотропной однофазной среды.

В механике грунтов при рассмотрении двухфазных сред получили распространение так называемые эффективные (фик­ тивные) аФ напряжения:

 

a«f = ( l - « ) ( a lV+p5lV).

 

(4.11)

Среднее напряжение для двухфазной среды в этом случае

имеет вид

 

 

 

 

 

 

Г

 

(4.12)

Из вышеизложенного следует, что напряженное состояние

скелета в

двухфазной среде может быть описано

[32, 43 ]

с помощью

средних

или эффективных аФ

напряжений.

К. Терцаги, анализируя результаты испытаний водонасыщен­

ных песков

и других

грунтов, пришел (1961

г.) к

выводу,

что условия разрушения этих тел зависят только от эффек­ тивных напряжений.

В связи с этим для описания предельного состояния насыщенных пористых грунтов К. Терцаги предложено ис­ пользовать предельное условие, которое справедливо для рассматриваемой среды при отсутствии жидкости, и заменить в нем напряжения Г,, на эффективные аФ. Сформулированное выше положение о предельном условии называется принципом К. Терцаги.

Для обоснования выбора вида напряжений в условии пластичности двухфазной среды рассмотрим вначале прямо­ линейную огибающую кругов главных наибольших напряжений (4.10).

Запишем условие (4.10) предельного состояния применитель­ но к изотропной насыщенной пористой среде, используя

принцип К. Терцаги:

 

Т(Л = а (Л sin р+ К cos р,

(4.13)

где т(/, = (1 —«)т = 0,5( 1- и ) ( а Г - а 3);

 

ст(/) = (1 —и) [(<?! + а 3)/2+/?].

(4.14)

Для формулировки условия предельного состояния двухфаз­ ной среды можно воспользоваться также средними напряже:

ниями Tij в двухфазной среде.

(4.10) предель­

Используя связь (4.7), запишем [32] условие

ного состояния относительно напряжений а 0-

в скелете

а! —а 3 = (а! + a 3)sin р+ 2АГ/(1 —/7)cos р —2л?/( 1—п)рsin р

 

(4.15)

и эффективных напряжений

 

т(Л = а (Л sin р+ К cos р—р sin р.

(4.16)

Сопоставив условия пластичности (4.13) и (4.16) с эффек­ тивными напряжениями, заметим, что они отличаются слага­ емым в правой части. Оценка записанных выше условий пластичности выполняется путем сопоставления эксперимен-

 

 

 

 

 

 

т</>

т(/>,

т(Л,

 

I,

о3,

А

<Ji,

а = 0,5(а, - а 3),

ст(л

эксперимен­

рассчитан­

рассчитан­

тальные,

ные по фор­

ные по

пред­

МПа

МПа

МПа

МПа

МПа

МПа

МПа

муле (4.13)

лагаемой

фор­

 

 

 

 

 

 

 

Терцаги,

муле (4.16),

 

 

 

 

 

 

 

МПа

МПа

5

2,5

0

12,5

7,5

4,2

2,8

2,5

2,5

 

2.5

2,5

2

7,5

5

3,9

1,4

2,4

1,4

 

6,9

5

0

18,8

11,9

6,7

3,9

3,8

3,8

 

5

5

2

15

10

6,7

2,8

3,8

2,8

 

5

10

7,8

20

15

12,8

2,8

6,8

2,9

 

7,9

10

4

25,8

17,9

12,3

4,4

6,6

4,6

 

тальных и расчетных значений касательных напряжений для шурабского бурового угля.

Значения показателей п, К, р, найденных на основе ап­

проксимации

экспериментальных

данных (рис. 4.1, я)

условием

предельного

состояния (4.15),

составили л = 0,438;

#=0,505

МПа; р = 30°

Сопоставление экспериментальных и

расчетных

значений т свидетельствует об их удовлетворительной сходи­ мости.

Перейдем теперь к анализу условий (4.13) и (4.16) преде­ льного состояния, записанных относительно напряжений стФ. Экспериментальные и расчетные значения наибольших касатель­ ных эффективных напряжений для шурабского угля с исполь­ зованием условий предельных состояний (4.13) и (4.16) при­ ведены в табл. 4.3.

Анализ экспериментальных и расчетных значений напряже­ ний т(Л показывает, что с экспериментальными данными согласуются только расчетные значения т(Л, найденные по формуле (4.16). Значения расчетных напряжений т(Л, полученные

по формуле

(4.13), основанной на

использовании принципа

К. Терцаги,

в большинстве случаев

при р ф 0 существенно

превышают экспериментальные данные. Таким образом, в даль­ нейшем для описания предельного состояния насыщенных пористых двухфазных изотропных гррных пород будем пользо­ ваться условием предельного состояния (4.10) относительно средних напряжений Г,,, записанных через напряжения в скелете Gij (4.15) или эффективные напряжения а |л (4.16).

Предельное условие Мора (4.8) для полных напряжений (4.7) в двухфазной насыщенной среде можно обобщить в виде:

т = -р 1 -/[(1 -и )ст -й р ].

(4.17)

Анализ экспериментальных данных (рис. 4.1), а также усло­ вий (4.15), (4.16) предельного состояния показывает, что для двухфазной среды при наличии давлении р газа (жидкости)

наибольшее касательное напряжение, при котором происходит разрушение материала, оказывается меньшим, чем соответст­ вующее касательное напряжение для того же материала при отсутствии давления. С точки зрения теории прочности Мора наличие давления р газа (жидкости) приводит уменьшению среднего давления и, следовательно, смещению точки по огибающей главных наибольших напряжений к началу коор­

динат,

т. е.

к снижению т.

В

работе

[32 ] приведены экспериментальные данные по

предельному деформированию ряда горных пород. Показано, что угол наклона экспериментальных линий (т—р) зависит от относительного содержания п газовой (жидкой) фазы.

На пористость в предельном состоянии оказывает влияние развитие процесса деформирования горных пород. Под дей­ ствием нормальных напряжений происходят объемное упругое сжатие среды и уменьшение ее пористости. В допредельном, предельном и запредельном состояниях на этот процесс накладывается второй конкурирующий процесс накопления повреждений под действием сдвиговых и нормальных напряже­ ний. В результате этого среда пластически разрыхляется, объем пор в ней увеличивается.

На основании статистической модели [31, 32] прочности, горных пород можно полагать, что разрыхление пород приводит к проявлению микротрещин и, следовательно, до­ полнительной пористости. Тогда можно считать, что по­ ристость в предельном состоянии горных пород приближенно равна истинному ее значению п0 в начальный момент времени и дополнительной пористости, вызванной процессом

деформирования. Значение ее можно

приближенно найти

по уравнению

 

л = /7о + 0р,

(4.18)

где 0Р — пластическое разрыхление горных пород в предельном состоянии.

Для наиболее простого случая линейной аппроксимации огибающей кругов главных наибольших напряжений могут быть написаны соответствующие условия предельного состо­ яния. При плоской деформации они выглядят таким образом:

(стя- а ,,)2Н4т2,.= sin2 р[стЛ.+ ст},+2/(1 —«)(AT—«/?tgр)ctgр]2,

(4.19)

где К, р — соответственно средние значения сцепления и угла внутреннего трения двухфазной среды.

Значения показателей К, р, п для ряда пород, полученные в результате аппроксимации экспериментальных данных зави­ симостью (4.15), приведены в работе [32]. Анализ этих данных

показал, что значения показателя п значительно превосходят значение прочности. Так, установлено, что значение параметра п, найденное в результате аппроксимации экспериментальных данных, для мрамора равно 0,66. Значение пористости мрамора

в

недеформируемом состоянии [31, 32] — п0 = 0,001;

0,009,

а

с учетом пластического разрыхления л = 0,001; 0,084.

Таким

образом, имеется существенное различие между значением параметра л, полученного в результате аппроксимации, и значе­ нием пористости мрамора.

Установленное различие может быть объяснено на основе статистической модели [31] прочности. В результате пластичес­ ких деформаций по элементарной площадке происходит сдвиг с отрывом, т. е. имеет место раскрытие микротрещины. При пластических деформациях контакт между частицами в этом случае происходит не по всему скелету, а по его части.

Следовательно, соотношение (4.7) может быть

обобщено [32]

таким образом:

 

 

 

Гу = 5(1 -л)<Тц-[1 -5 (1 -п)]рЬф

(4.20)

где S — площадь контактов

в скелете породы

(0<5<1).

Введя обозначение

 

 

 

т= 1-5(1 -л ),

(4.21)

соотношение (4.20)

запишем

так:

 

 

Гу = (1

t n ) o (j — m p b j j .

(4.22)

Если пористость

породы

оценивается значением (4.18), то

т можно представить в виде

 

 

ш = шо + 50р,

(4.23)

где m0= l - S ( l - n 0).

Таким образом, соотношение (4.22) по виду совпадает со связью (4.7). Однако значение показателя ш учитывает пори­ стость среды и ее микроструктуру, которая формируется на фоне пластических деформаций и развития процесса разруше­ ния. Учитывая это в дальнейшем, под показателем т будем понимать условное относительное объемное содержание газо­

вой (жидкой) фазы.

мрамора //? = 0,66.

В рассмотренном выше примере для

При пористости,

оцениваемой значением /2

= 0,05,

5 = 0,36.

В дальнейшем

будем предполагать, что

между

напряжени­

ями в двухфазной среде имеется связь (4.22). Тогда условия

пластичности (4.15), (4.16), (4.17)

сохраняют свой вид,

однако

в них

вместо п следует использовать т.

данных

На

основе многочисленных

экспериментальных

в работе [31] предложено условие (1.2) пластичности, наиболее полно описывающее поведение горных пород. Это условие

может быть обобщено на насыщенную пористую среду сле­ дующим образом:

x - t |,/ ( l - m)ex p L {;-'";°--"'’l .

(4.24)

[_ (1 —m)Gi —nipJ

 

Уравнения предельного состояния двухфазной газонасыщенной среды с учетом адсорбции угля. При анализе напряженнодеформированного состояния газонасыщенного трещиноватого пористого массива угля ^может быть опцеана модель двухфаз­ ной среды, с помощью которой удается учесть влияние свободного и сорбированного газа.

Будем рассматривать двухфазную среду, которая харак­ теризуется наличием трещин, разделяющих ее на элементарные блоки. В трещинах и макропорах этой среды происходит движение газа, которое, в основном, описывается законом Дарси. В микропорах основной формой переноса является диффузия газа.

Для описания предельного состояния такой среды будем считать, что она состоит из скелета, трещин и макропор, удельное содержание которых характеризуется т ь и макропор с содержанием т2. При этом предполагается, что трещины и макропоры не связаны с микропорами. Поэтому давление газа в них будет различным и характеризоваться значением /?ь а в микропорах^—р2. Тогда среднее напряжение в такой двухфазной среде будет

Гу = { l - m 1- m 2)aij- m lp l - m 2p2.

(4.25)

Будем полагать, что прочность двухфазной среды определя­ ется физико-механическим взаимодействием газа с пористой трещиноватой структурой пород. Тогда с учетом связи (4.25) могут быть обобщены условия (4.15), (4.16), (4.17) предельного состояния. Например, условие предельного состояния (4.19) запишется так

(а* - а.у ) + 4т^,. = sin2 р {а* + сту + 2 ctg р [A —

+ ю2р2) tg р]/

( 1 - т , —ш2)}2.

(4.26)

Уравнения предельного состояния многокомпонентных сред и горных пород. При анализе насыщенных пористых трещинова­ тых горных пород могут быть использованы основные положе­ ния механики гетерогенных многокомпонентных сред.

С точки зрения механики гетерогенных сред насыщенный массив горных пород будем моделировать многофазной мно­ гокомпонентной сплошной средой. Эта среда содержит скелет, поры и трещины которого насыщены жидкостью и газом. Материал скелета имеет различные размеры зерен, которые

с помощью петрографического анализа могут быть разделены на ряд компонент. Эти компоненты имеют различные прочност­ ные и деформационные характеристики. Заметим, что резуль­ таты [32, 43] экспериментальных исследований по предельному деформированию не насыщенных газом (жидкостью) одноком­ понентных и двухкомпонентных горных пород показывают их различие. Удельное относительное содержание г'-й компоненты будем характеризовать величиной /и;. Количество компонент в твердой фазе обозначим через г, а в жидкой и газовой — S. Тогда удельное содержание различных фаз в гетерогенной среде характеризуется

 

 

1 > £ = 1 .

 

 

 

(4-27)

 

1= 1

 

 

 

 

 

Среднее напряжение Ги в многофазной многокомпонентной

среде имеет

вид

 

 

 

 

 

 

 

Гij= £

w ,ay-

£

miPi,

 

(4.28)

 

i = 1

i =

1

 

 

 

где Gij— напряжение в скелете /-й

компоненты

твердой

фазы;

Pi — давление

жидкости или газа

в

/-й

компоненте жидкости

и газа.

 

скелета

в

общем

случае

имеют

Компоненты материала

различные законы деформирования. Например, наиболее проч­ ная часть скелета может деформироваться по закону Гука, а другие составляющие его находятся в пластическом состоянии или испытывают деформационное разупрочнение. Математичес­ кая модель подобной среды при отсутствии жидкости и газа рассмотрена в разд. 3. Давление в жидкой и газовой фазах Pi в общем случае предполагается различным.

Построение условий предельного состояния многофазных многокомпонентных горных пород выполним с помощью средних Г., напряжений. Для этого используем условие пре­

дельного состояния однородной однофазной среды, в котором напряжения заменяются средними Г0напряжениями (4.28). Запишем, например, условие предельного состояния Кулона

(4.19) для

многофазной

многокомпонентной

среды

 

 

2

/

г

\2

 

 

_1=1

ч

=

Z Wit*,

= sin2p^]T

+

,

 

 

 

 

- 2

£

miPi+ 2К ctg pI

 

(4.29)

 

 

i =

1

J

 

 

где

a*, o yl, i xyl — компоненты напряжений

в

скелете;

К,

р— соответственно средние для скелета

среды

сцепле­

ние

и

угол

внутреннего трения, которые могут

быть

записаны

в виде

 

 

 

 

 

К= £ т{Ка

р= X mi Рь

(4.30)

 

 

 

1 = 1

1 = 1

 

при

этом

р, являются

значением сцепления и

угла

внутреннего трения в /-й компоненте скелета.

Предельное условие Мора (4.17) многофазной многоком­ понентной среды может быть обобщено следующим образом:

1/2 £ m i ( o \ - a i ) = f

1/2 £ /и,(а‘ + а '3) - £

т,/>,

i=i

L i=i

î=i

 

где а \, а'3 — компоненты

главных

напряжений

в скелете.

Рассмотрим частный случай многокомпонентной среды, имеющий прикладное значение при анализе процессов фильтра­

ции

нефти

и газа. Будем считать, что среда состоит

из

скелета,

содержит газовую фазу, а также

жидкостную

в виде нефти и

воды. Удельное

соотношение различных

фаз

этой

среды

характеризуется

величинами

mH, т г, т в,

означающими соответственно удельные относительные соде­ ржания нефти, газа и воды. Предполагается, что жидкая и газовая фазы образуют многокомпонентную смесь, давление

в

которой равно р. Тогда средние напряжения

Г,7

в многокомпонентной среде

 

 

Гij= (1—ти—тв — тг)оц —(ти + т в+ т г)р.

(4.32)

Используя предложенную идею построения условий преде­ льного состояния многофазных многокомпонентных сред, усло­ вие Кулона (4.29) можно записать в форме (4.19), где параметр п = тн + тв + т г характеризует суммарное влияние жидкой и га­ зовой фаз.

4.4. ИССЛЕДОВАНИЕ СИСТЕМЫ УРАВНЕНИЙ ПРЕДЕЛЬНОГО ДЕФОРМИРОВАНИЯ ДВУХФАЗНЫХ СРЕД И ГОРНЫХ ПОРОД ПРИ НЕЛИНЕЙНОЙ ФИЛЬТРАЦИИ

При изучении процессов пластического дефор­ мирования насыщенных горных пород необходимо рассматри­ вать уравнения движения для твердой и газовой фаз.

Образование областей предельного состояния в массивах пород с пластическим разрыхлением или уплотнением среды меняет взаимосвязь между твердой и газовой (жидкой) фазами. Приобретенная пористость среды в результате ее пластического деформирования зависит от распределения напряжений в ске­ лете и давления газа (жидкости) в порах, снижает прочностные