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

книги / Механика горных пород при разработке месторождений углеводородного сырья

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

ГЛАВА 2

МОДЕЛИ МЕХАНИКИ ГОРНЫХ ПОРОД, ИСПОЛЬЗУЕМЫЕ ДЛЯ РАСЧЕТОВ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ГОРНЫХ МАССИВОВ МЕСТОРОЖДЕНИЙ УГЛЕВОДОРОДОВ

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

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

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

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

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

2.1. ОСНОВНЫЕ МЕХАНИЧЕСКИЕ МОДЕЛИ ГОРНЫХ ПОРОД

Различные механические модели горных пород в той или иной степени идеализируют основные свойства поведения пород под нагрузкой, такие как упругость, пластичность и ползучесть. Эти модели весьма широко представлены в научно-технической литературе [1, 3, 4, 6, 7, 8], поэтому ниже они освещаются весьма кратко. При этом авторы данной работы придерживаются обо­ значений, принятых в европейской научно-технической литера­ туре, в частности, в работах В. Виттке [34, 35] и его учеников.

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

Рис. 2.1.1. Деформируемость изотропной горной породы:

1 - упругие деформации; 2 - пластические деформации; / - высота образца

никают исключительно упругие деформации, которые не зависят от времени и прямо пропорциональны нагрузке в соответствии с законом Гука (рис. 2.1.1).

При описании упругих деформаций пород исходят из моделей их структуры. Для пород с неориентированной зернистой струк­ турой возникающие упругие деформации не зависят от направ­ ления действия нагрузки. В этом случае для описания напря­ женно-деформированного состояния достаточно двух констант материала: модуля упругости Е и коэффициента Пуассона v, так как модуль сдвига можно выразить через Е и v: G = £/[2(1 + v)]. В соответствии с законом Гука напряжения и деформации в произвольной декартовой системе координат (х, у} z) можно свя­ зать между собой следующими соотношениями:

=

Е

[ ( 1 " VK + vey + v e 2];

1 - V -

 

2v2

<т = °У

1

Е >1

МС > МС

[ v e , +(1

v)ey + v e ,] ;

(Т =

 

Е

1

К]

< 1 to<

тХУ

[ve* + vey + (1 - v )6 j];

Е

( 2.1.1)

ю Т 7 )Уху’

 

Уyz >

 

2(1 + V)

Tzr

2(l + v )Yzr*

В матричной форме записи выражения (2.1.1) будут выгля­

деть следующим образом:

 

н =

(2.1.2)

где

 

~~ х ' ® у ' ®

ху' ^уг' ^гу )

{в} =

 

£ ( 1

-

1

<

1

 

E v

1

1 >

 

E v

1 -

v

-

Н =

0

0

0

v ) ГО<

гм>СМ

2 v 2

 

E v

 

1 -

v

-

2 v 2

£ ( 1

-

v )

1

1 >

МГ > СМ

 

E v

 

> см 1 > 1

 

 

0

 

 

 

0

 

 

 

0

 

 

 

£ v

 

 

1

1 >

СМ>см

 

 

E v

 

1

-

v -

2 v 2

 

£ ( 1 -

v )

1

-

v -

2 v 2

 

 

0

 

 

 

0

 

 

 

0

 

0

0

0

0

0

0

0

0

0

£

0

0

2(1 +

V )

 

£

0

0

2(1 + V )

£

0 0

2(1 + v)

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

Для таких пород напряженно-деформированное состояние описывается с помощью модели трансверсально-изотропного уп­ ругого тела [34]. Для этого необходимо пять констант материала:

 

 

ст

 

Рис.

2.1.2. Анизотропная

 

/ I

*-----------71

 

деформируемость породы

 

 

со слоистой структурой:

/

I

/

I

1 -

плоскость изотропии;

 

2 - слоистость

модули упругости Ei и Е2, характеризующие сжимаемость соот­ ветственно перпендикулярно к плоскости изотропии и в плоско­ сти изотропии, два коэффициента Пуассона vj и v2, а также мо­

дуль сдвига

G2 для напряжений сдвига в

плоскости изотропии.

Если принять, что в декартовой системе координат (,х', у \ / )

ось

/ совпадает

с направлением наибольшей

деформируемости,

то

можно записать следующее соотношение между напряжениями и деформациями:

 

И

= 1

 

 

 

 

(2.1.3)

 

 

 

 

 

 

при

 

 

 

 

 

 

 

1 -nvl

V ,

+ nvl

E ^

0

0

0

 

(1 + vi)m

1(1 + v{)m

1m

 

 

 

 

v1+nvi

1 -nvl

E ^

0

0

0

 

(1 + Vj)m

1(1 +v,)m

1m

 

 

 

 

[£>1 =

 

m

£21 ~mV|

0

0

0

(2.1.4)

m

 

 

0

 

0

0

2(1+v,)

0

0

 

0

 

0

0

0

G2

0

 

0

 

0

0

0

0

G2

 

где п = £,£2; от = 1 -

v, -

2п\\.

 

 

 

 

 

Также можно записать обратное соотношение

И - И" К).

где

1

J4 1\<

_ V2

е2

0

0

0

1 \<

_ v2

 

е2

1

_ v2

же2

_ V2

1

Т2

е2

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

2(1+V,)

0

0

 

 

 

0

1

0

 

0

0

1

G2

 

 

Инженерные задачи обычно решаются в некоторой неподвиж­ ной в пространстве системе координат (х, у , z), которая, как пра­ вило, не совпадает с системой координат (х1, у , z/). Связь между ними может быть установлена с помощью двух углов, которые задают направление горизонтали (угол простирания ос и угол падения Р) плоскости изотропии относительно системы коорди­ нат (xt у уг) (рис. 2.1.3). Связь между напряжениями и деформа­ циями в неподвижной системе координат (xt */, z) и системе ко­ ординат, связанной со слоистостью (х?у у'у z'), устанавливается с помощью матриц преобразования [Т\ и [Т*]:

{&} = [ГИа};

{£'} = [Г Н в}.

Матрицы [7] и [Т*] имеют следующий вид:

 

2

0

2 1 ,7 7 * ,

0

1 ?

щ

 

 

 

 

 

 

 

2

„ 2

2

1 2 77*2

277*27*2

1 ’

77*2

п 2

 

 

 

 

 

 

 

 

2

 

2 1 3 т

3

277*3 7*з

13

щ

« 3

 

 

 

 

1 , 1 2

m i m 2

0

1,77*2

+

1 2 т ,

277*, 7*2

^ 2 ^ 3

m 2 m 3

Я Л

1 2 77*3

+

1 3 Т П 2

77*27*3 + Ш 3 П 2

M l

щ

щ

0

1 , 7 7 * 3 + 1 3 77*,

77*, 7*3

 

 

 

 

 

 

 

(2.1.5)

(2.1.6)

0

2 ^ 2 1 г

2 ” з 1 з

” 2 ! .

” 2 l 3 + Я з 1 :

« з ! .

(2.1.7)

Рис. 2.1.3. Связь глобальной (х, у, г) и локальной (х\ у', г*) систем координат:

1 - линия падения; 2 - горизонталь; 3 - плоскость изотропии; а -

угол простира­

 

 

 

 

ния; р - угол падения

 

 

 

1 ?

 

щ

0

1 , 777,

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

„ _ 2

„ 2

 

 

 

 

 

 

 

 

ТП>2

2

 

 

щ

щ

^ * 2 1 г

 

 

 

 

„ 2

 

 

щ щ

 

 

1 з

 

Щ

 

1 3 ™ з

” з 1 з

2 1 , 1 а

 

2

0

1 ,7 7 1 2

+

1 2 777,

Щ

П 2

« Л

 

 

 

 

 

 

 

 

2

1 2 1 3

Ъ п ^ Щ

 

1 2 Т ^ з

+

1 3 7772

T n j n j

+ Ш з П 2

^ М з + ” 3 1 2

 

 

 

 

 

 

 

 

 

. 2 M

I

2 щ щ

0

1 ,7 7 7 3

+

1 3 777,

Щ П г

” a l i

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

(2.1.8)

 

 

 

 

 

 

 

 

 

 

li

= sin ос;

 

т\ =

cos а;

 

 

 

Ь

= cos Р cos а;

 

т2 = -

cos Р sin а;

п2 = -

sin Р;

13 = -

 

sin р cos а;

 

га3 = sin Р sin а;

п3 = -

cos р.

При этом

 

 

 

 

 

 

 

 

 

 

 

 

[Г!’1= [Г]т

 

 

(2.1.9)

Если выражения (2.1.5) и (2.1.6) подставить в (2.1.3), то по­ лучим

[Л М = [D'][f]{e}.

При умножении на [Г]"1находим

{а} = П Г '-т [Г ].{г).

Окончательно с учетом (2.1.9) получим

 

{о} = [ D m ,

( 2. 1. 10)

где [D] = [ f] T \D ']\T ).

 

Аналогично можно получить зависимость

 

М = [Т\т[07'[Т \{°}-

(2.1.11)

Выражения (2.1.10) и (2.1.11) описывают трансверсально­ изотропную зависимость между напряжениями и деформациями в неподвижной в пространстве системе координат (х, у} z) при наличии плоскостей изотропии, ориентированных относительно глобальной системы координат с помощью углов а и р .

Предположение о том, что породы со слоистой или линейной структурой проявляют трансверсально-изотропные упругие свой­ ства, не всегда соответствует встречающимся в природе услови­ ям. Например, если скальная порода расчленена тремя взаимно ортогональными системами трещин, которые частично раскрыты и заполнены, то будет проявляться зависимость деформируемо­ сти от направления действия нагрузки. В этом случае имеет ме­ сто более высокая степень анизотропии - ортотропия. Для ее описания необходимо девять независимых констант упругости - это модули упругости Ей Е2 и Ез для трех взаимно ортогональных направлений, три коэффициента Пуассона vb v2 и v3, а также три модуля сдвига Gj, G2 и G3. Из теории упругости известно, что может понадобиться до 21 постоянной упругости для точного описания анизотропного поведения породы. Современные чис­ ленные методы расчета позволяют без труда учесть в расчетах любую степень анизотропии породы. Однако на практике экспе­ риментально определить такое число постоянных упругости вряд ли возможно из-за высокой сложности и стоимости работ. По­ этому на практике чаще всего применяется упругая изотропная или трансверсально-изотропная модель горного массива.

2.2. КРИТЕРИИ РАЗРУШЕНИЯ И НЕУПРУГИЕ МОДЕЛИ ГОРНЫХ ПОРОД

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

Q*

Рис. 2.2.1. Испытания образцов на трехосное сжатие и диаграмма напряжения Мора

Прочность пород изотропной структуры в большинстве случа­ ев определяется в испытаниях на одноосное или трехосное сжа­ тие на цилиндрических или прямоугольных образцах. Трехосные испытания проводятся таким образом, что минимальное и сред­ нее главные напряжения равны: а3 = а2 < Сть Значения напряже­ ний CTI и а3, замеренные в момент разрушения, можно изобразить на круговой диаграмме Мора (рис. 2.2.1). Огибающая кругов Мора, построенная по результатам отдельных испытаний, будет представлять собой паспорт прочности горной породы, т.е. точки, находящиеся выше огибающей, будут соответствовать таким па­ раметрам напряженного состояния, при которых предел прочно­ сти породы превышен.

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

т = atg фС + cG.

(2.2.1)

Это простейшая форма критерия разрушения Кулона-Мора. Его параметры cG и фс называются соответственно коэффициент сцепления и угол внутреннего трения. Из диаграммы Мора сле­ дует, что критерий разрушения реализуется на площадках, кото­ рые наклонены к наибольшему главному напряжению под угла­ ми р = ±(45° - фс/2) и называются площадками скольжения. При испытаниях на трехосное сжатие площадка скольжения на­ клонена к горизонту по углом a = 45° + фС/2 (см. рис. 2.2.1).

Если уравнение (2.2.1) записать через главные напряжения, то получим выражение

 

e l + sin(pG

2cos фс

(2.2.2)

 

z--- :------------------ :-----

 

1 - sin <pG

1 - sin <pG

 

В координатах

- a3 условие разрушения Кулона-Мора так­

же изображается прямой линией (рис. 2.2.2). Угол наклона этой прямой является функцией угла внутреннего трения фС, а отре­ зок на оси ai представляет собой прочность на одноосное сжатие

®dG-

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

образом:

 

(а* - Gy)2+ 4 т2ху = (a* + Gy +2 cGsin фс)2 sin2 фс.

(2.2.3)