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

книги / Физические свойства осадочных пород при высоких давлениях и температурах

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

Физические свойства песчаника Бнрп (США), насыщенного 2 н. раствором КаС1 при различных

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

г эф' к г *см>

 

 

 

 

 

 

Источник данных

Параметры

 

50

70

150

280

300

420

450

600

700

750

900

1 2 0 0

1500

 

0

 

 

3000

3600

3840

4200

4360

4430

4500

4580

4700

4800

ВНИИГеофп-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зпка

ур , м/сек

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3200

3600

3750

4100

4150,

4180

[102]

 

 

 

 

 

 

 

4000

 

 

 

 

 

 

 

 

 

1,0

1,14

2,40

4,44

5,30

6,71

ВНИИГеофп-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зика

А^п, %

1.0

2,5

3,5

5,9

[43, 101,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

106]

 

1.0

1,5

3,0

 

— —

р„, 100 смЪ/кГ

— 200

рП1 ом-м

1,70 1,85

 

114,0

84,0

51,0

41,0

27,8

ВИИИГеофн-

 

 

 

 

 

 

 

 

 

 

 

 

знка

1-Т

О

 

93

6 0 -

64

 

 

40—42

 

[43,

101, 106

О1

 

 

 

 

 

 

 

и

др.]

 

С"

 

 

70

 

 

 

 

 

 

 

 

 

 

1,90

1,96

2.03

1 2,07

2,11 2,14

2,16

2,17

ВНИИГеофи-

 

 

 

 

 

 

 

 

 

 

 

|

злка

переходе от толщ-покрышек к скоплениям нефти и газа в коллек­ торах, исследования должны проводиться при одновременном повышении всех параметров (рс, р 1 и I).

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

Результаты измерения некоторых физических параметров образца песчаника Бири на установке УФС ВЫИИГеофизика даны в табл. 4 и 5. Сравнение приведенных значений с опубликован­ ными данными показывает хорошую сходимость как абсолютных величин параметров, так и их изменений под действием давления.

Глава III

ВЛИЯНИЕ ВСЕСТОРОННЕГО ДАВЛЕНИЯ НА СЖИМАЕМОСТЬ, ПОРИСТОСТЬ И ПЛОТНОСТЬ ОСАДОЧНЫХ ГОРНЫХ ПОРОД

1. СЖИМАЕМОСТЬ ОСАДОЧНЫХ ПОРОД ПРИ ВСЕСТОРОННИХ ДАВЛЕНИЯХ

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

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

1 Образцы песчаника любезно были псредапы паучно-псслодовательским центрам фирмы «Шелл» профессору |С. Г. Комарову!» Оппсаппе шлифа про­ ведено в МИНХ п ГП*

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

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

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

В. М. Добрынин [35, 36], анализируя некоторые из указан­ ных теоретических расчетов и сопоставляя их с имеющимися экспериментальными результатами, пришел к важному выводу о том, что для зернистых горных пород в отличие от сплошных сред сжимаемость зависит от величины давления на скелет, умень­ шаясь при увеличении давления по закону

|3П= сопзЪ •/г 1/*.

(12)

Однако следует отметить, что уравнение (12) было получено на основании обобщения решений, где использовались в основном результаты контактной задачи Герца для случая упругого сжатия двух шаров под действием сосредоточенной силы Р, т. е. решений, где не учитывались возникающие на контактах тангенциальные напряжения. Согласно уравнению (12), при р -> 0 рп -> °о, при р оо или достаточно больших давлениях (5П= 0. Это противо­ речит физической сущности понятия сжимаемости, вследствие чего можно предположить, что уравнение (12) применимо в огра­ ниченном диапазоне давлений.

Полное изменение объема пористой однородной и изотропной среды является функцией напряжения в скелете породы и жид­ кости. В дифференциальной форме относительная объемная де­ формация представляется в виде следующей системы уравнений [35, 94, 108, 1111:

где Ус'к, К и К — объемы скелета породы,-се твердой фазы и пе­ рового пространства.

В. М. Добрынин. [361 рекомендует рассматривать напряжение в скелете породы как функцию разности давлений ре р; = р,,. В этом случае система уравнений (13) будет иметь вид

Ш

, - у Ь ( - Г О Й И - ) * й (л - й ) +

 

(14)

(Щ - Ь й т а Э а ) * . * * * - * (15)

т

( т

&

г к

 

, - < * - » > + . < * •

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 16)

 

 

 

 

 

 

 

 

( д Г с к \

 

, _ } _ ( ™

к \

 

(17)

 

_ ^ _ / З П к \

 

 

■ 1

 

 

.

 

Уск

\

дрс

/ ( Р е- Р,-).

<

Уск

\

аРс

) Рр

I

Уск. V аР1

/ рс, I ’

 

 

Уп

\

дРе

){Ре-Р[). *

 

У'п

\

дРо

/

Р/1 '

'

Уп \

ЭР1

/ Р с, Г

 

^ '

 

 

 

Ут

V

дре ) ( р е- Р1),

=

± . ( * Б \

1

(

дУт

\

(19)

Обозначим

I

V* \

дре ) Р(,

Ут \

дР1

) Р(

 

 

 

1

 

дУ-ы__ \

 

= —3 ' -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(20)

 

 

 

 

 

ск \ П Р е ~ Р д

) Р1, (

 

1ск’

 

 

 

 

 

 

 

 

 

Уы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(

дРе

1

 

 

_

_ й .

 

 

 

 

(21)

 

 

 

 

 

У1

V

) {РгР1), I

 

1 г ’

 

 

 

 

 

 

 

 

 

 

I (

еЪ

 

\

 

= - _ з ' -

 

 

 

 

(22)

 

 

 

 

 

У'„

\Э{ре—Р1) )р[, <1

 

 

 

 

 

 

 

 

 

 

 

1

(

дУт

\

 

=

_3".

 

 

 

 

 

 

 

 

 

 

у ;

\а (Р е ~ Р 1 ) ) р, , (

 

1 т ’

 

 

 

 

 

У'ск

V дре

=

_ в " .

 

)(»'->,). I

Рск’

 

- А - № )

=

_Р /

(25)

К

V «Ре

*

' п'

 

Если принять сжимаемость скелета рек и пор Рп как относи­ тельное изменение соответствующих объемов на единицу измене­ ния всестороннего давления при постоянном норовом давлении — уравнения (20) и (22), — а сжимаемость твердой фазы рт — при постоянном дифференциальном давлении (уравнение (21)), то для решения системы уравнений (13) необходимо определить еще зна­

чения производных рп,

Рек и рт. Величину Рт можно найти, за­

менив в уравнении (23)

значение Ут на Ус'к — У„ с учетом равен­

ств (20) и (22).

'

а-

1

( ЭК

\

1

( д(У’ы -Уп\

=

1 т

У х

V д ( Р е - Р . )

) Р{, I

У ^ к - У п

V а ( Р е - Р 1 ) ) Р(, ,

 

 

 

= “

Г = Ц - (Р с к - а д .

(26)

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

( * У т’ \

 

= (

)

= 3 —

(27)

\ Ут )(Рй- Н ) ' г

^

^

7 ск

Г

 

Условие однородности скелета и вытекающее из этого равен­

ство (27) использовано щля

оценки

деформаций

пористой

среды

в работах [36,

94].

 

 

 

 

С учетом уравнений (20) и (27) значения^ Рек и Рп можно опреде­

лить, рассматривая равенство (24)

при

условии ре р1 =

соп§1.

\ У ск’

1 '

'{ Р е ~ р1 ) ' г

'

= -Рт*Л -

(28)

' [Ре~рС)'

 

Следовательно, на основании выражения (28) можно записать

Рт =•-• Рек =• Рп :

__ 1

/ ^ ; \

 

=

(

^ск \

=

У'т\ аре /(р^-р.), I

VI*

\

дРе / ( Р с-Р ,),

I

 

_____ 1

( агп \

 

 

 

(29)

 

~

Ут,

VРе ) СРе~р1)<

 

 

 

 

 

 

В окончательном виде связь напряжение — деформация в слу­ чае объемного изменения однородной пористой среды с учетом

равенств (20), (26) и (28) можно выразить следующей системой уравнений:

-

( 4 § г ) , = Р« л (рс -

Р д +

Рт

(30)

-

( ^ • ) / = м

(Ре -

Р1) рт лрр

(31)

- ( т г ) ,

-

/спРп)

(/>« -

Л) + Рт

(32)

Таким образом, объемная деформация пористой среды выра­ жается четырьмя независимыми величинами: Рек* (Зп, Рт и &п«

При выводе равенства (27) предполагается постоянство коэф­ фициента пористости породы при ре р 1 = соп81. Однако для

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

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

 

 

 

о

ЛЯ

1

3

 

 

 

 

(33)

 

 

 

3 -тг* 1+ с

г

1+

с '

 

тогда

как относительная

деформация

твердой

фазы

 

 

/ ау1 \

ч ап.

1

от

з

 

 

(34)

 

\ у ;

Л

° я

1 +

ЗеЗ г а

г •

1+ 3 е з

 

 

Здесь

с — доля

зерен

радиусом Я.

 

 

 

 

 

 

Из сравнения выражений (33) и (34) следует, что для рассмат­ риваемой модели равенство относительных деформаций объема скелета и твердой фазы не соблюдается, а уравнение (27) является частным случаем при условии йЯ/Я = йг/г, т. е., когда скелет однородный.

Поскольку

Уп

У?+Уп’

4 Г»

то после логарифмирования и дифференцирования уравнения

(35)

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

№ ) (* ,- * > . / =

 

 

( ^ ) ( р е-рЛ. <+

( - Т г ) ( Ре-Р1), г

<36>

 

 

 

 

 

к'п

( а <п \

I

( ду; \

 

 

/07ч

 

 

\ ^’ск

/ {Ре~Р() . 1

1 ~ * п \ 1‘п

' ( р е~р1)’ 1

1

^

у т ' ( Р е- Р { ) - 1 ‘

 

 

Учитывая равенства (14), (15) и (28)

при (ре — р{) =

сопзЬ,

можно записать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ж п \

 

 

 

 

1

 

/

йк'г,

 

 

 

 

 

(38)

 

 

(

^п" ) ( Р е- Р 1), I

 

 

1 - * п

(

к Ц )(р е-р,),

г

Рт * Р 1 '

 

 

 

 

 

 

 

 

 

/

^ с к

\

 

^

 

А-п

/

^

 

 

__ О

л

,

 

 

 

 

 

__

 

 

 

 

 

 

 

 

(39)

 

 

V

^ск

 

 

*

 

1 - * п

V

*п Л - - Л

#

Р т ®Л

 

 

 

 

 

 

 

 

1

Следовательно,

в

/

общем

 

виде

частные производные

/ а 7 ; \

% и

1

ау'к \

 

определяются

не только

177- -д-г- К

177” I

 

дре

)

 

Уп

\

дРе 1 (Ре- Р 1),*

У ск

\

 

/(?,-/>*), <

 

 

 

 

 

 

сжимаемостью твердой фазы породы, но и начальной пористостью

пее изменением при

Сучетом равенств (38) и (39) относительную деформацию по­ ристой среды с неоднородным скелетом можно выразить систе­ мой уравнений

( У'ск ) / "" Рск<*'(р *

Р :)

1-"-п ( А-п" \ р е- Р с ) , ‘ 1 Рт ё Р Г ’

(40)

 

- ( ^ = ^ ^ - ^ - 14(1 ),-,,/^ ^

(41)

- ( т ^ ) , = 1 = ^

(Рек -

К К ) <* (А - А ) + Рт <*А-

(42)

Когда порода состоит из сцементированного упругого неод­ нородного скелета и упругих минеральных (глинистых) вклю­ чений, свободно размещенных в его норовом пространстве, при оценке объемных деформаций скелета, порового объема и твердой фазы породы следует учесть также объемные изменения включе­ ний. Деформация аналогичной среды, но с однородным скелетом была рассмотрена В. М. Добрыниным [36]. Пользуясь уравне­ ниями (40)—(42) и методом анализа, примененным В. М. Добры­ ниным, найдем дифференциальные уравнения среды с неодно­ родным скелетом и включениями в порах скелета. Обозначим соответственно объемы скелета породы с включениями Т^ск, вклю­ чений Увк, пор — Уп. Связь между указанными параметрами породы с включениями и «чистой» породы, т. е. без включений, будет следующей:

К к= ^сК; Ув = Уя - У ък=-.Уп (1-чУ , 1сп = К ( 1 - 1 ! ) ;

У ^ - У т+ Увк] УПК= Т|^П1

где т| — относительная глинистость (коэффициент заполнения пор), определяемая уравнением

 

11 — ^вк __

^*рк('1 ^п)

 

1 У'п

Р вк (1 к п ) + к п *

Здесь

Рвк — относительное

весовое содержание включений

в породе

в %.

 

 

На основании приведенных

выражений и уравнений (14),

(15)и (16) для породы с неоднородным скелетом и включениями

впорах можем записать

%^

.. 1 (

ЭУок \

^ (Р е - Р с ) +

II

Уск \

дРс )р,

 

 

 

9<Ут+Уи+ У*

 

 

(43)

 

+

*

(

др1

Х г - л . ^

 

 

1

(ЭУП\

 

 

1

/ 3 ( ^ - К вк) \

ув \ эр,

) Р1, (й(Ре—■?<)+•

\

дР1

/(Ре-Р,.), (йР1\

 

 

 

 

 

 

 

 

(44)

/

<*7Х\

 

1 /

9 (^ + ^вк) '

^ { Р е - Р д - 'г

^

Л

 

Ут \

ЭРе

у)

+ Т 7

Г*(П+У»д V дР1

^

)

ч

 

«? 1

(45)

Если воспользоваться определениями сжимаемости скелета и пор для «чистого» коллектора применительно к породе с- вклю­

чениями, то

можно

записать

 

 

 

 

 

 

 

 

»

 

 

д у с к \

1

=

4 -

 

 

 

(46)

* со! 1

8к"

 

 

д р е >

 

^СК V

ЭРе / р 1г

 

 

 

 

 

 

 

 

' Р[ 1*

 

 

 

( Э У п

\

 

_

1

1

( д У й ■

6'

1

(47)

 

1 — т]

 

к

^

\

д Р е

/ Р1,

*

 

‘Т Г

э р , .

 

Рп

1 - 4

*

 

 

/

) Р1, ~

гДе Рек и рл — сжимаемость скелета и пор породы с включениями. В уравнении (47) учитывалось условие о свободном размещении включений в поровом объеме скелета, т. е. включения не подвер­

гались сжатию при действии внешней нагрузки на скелет.

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

С помощью уравнения (44) найдем

 

 

1

( Э(Ук~Увк)\

йр1 =

у п

\

дР1

/(РГР1)’ 1

1 Г 1 (М'п\

1

Т) |^1

к п \ к п / ( р с- р 1 ) , *

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

где Р„к — сжимаемость включений.

 

 

 

 

 

 

Преобразуя

уравнение

(48),

получим

 

 

 

 

 

 

 

(

)(ре-Р<). ^

у„

(

ар" )(рГР1) / р1~

 

 

 

(1—п)<1—*п) ( “*Г)<р,-л . ‘ ~ ^

йРь

(50)

где

 

 

 

 

 

 

(*«-*/)»

 

 

 

 

 

 

 

 

 

 

 

Рв1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и - п ЬРт-

 

 

 

 

(51)

 

 

 

 

Рр-

1 — у\

 

 

 

 

Второй член уравнения (43) после преобразования будет иметь

вид

 

ЛЁЕ«\

= Л - (*?&.)

 

*,.=

 

 

 

 

 

 

 

V

УСк

) (Рв- Р/), <

Рек \

дрс

/(Р е-Р ,).

<

 

 

1

{ <МУ;+Уп + Уек) \

 

.

__П_

у; V

 

 

лп .

"

 

 

 

; (Р _Р.), щ

-

уск •

 

; (Ре-р;), , ар>1-

,

Уп

1

(

ЭУп \

.

4 - 1 ™

_ ± 1

 

 

 

_

“1" Рек

' К„

I

ар; )(рг р<), 1 Р‘ ^

Уск * ^вк V 5рг

У(ре-р.), < -

(1 -

*п) Р, *д

( ^

) (Ре. Р(.), , -

А; (1 -

Ч> ЦрРх 4Р1 ~

 

 

а;лР„„аР1~

^

 

 

) (РГР|), ,

Рт

 

(52)

Чтобы найти объемную деформацию всей твердой фазы (ске­ лета совместно с включениями), воспользуемся уравнениями (21), (46) и (47).

/<*ут\

__/а(У;+ГВк) \

_ П ( * К \

. Уж ( йУ™ \ =

V Ут Л

V Ут

Л “

\

Л -1" Ут

^ Увк

Л

1 — Аг„ (1 — -п) *

 

^

ск “

 

^ ^Ре ~

Р^ ~

 

~

1-*п(1П- п )

ЛР‘ ~ к"1 - ^ 1 - г))

=

- 1 - ^ 1 - П ) (Рск- В

Д

-

\ - ^

^ Г

ар‘-

<*»

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

в поровом пространстве среды, на основании выражений (46),

(47),

(50),

(52)

и

(53) определится системой

уравнений

 

 

-

 

 

 

-

т

 

к

,+1

<54>

 

 

 

 

 

 

 

 

 

 

(55)

 

-

( т г ) , - -т = т = * г <&. ■

< *

(Р. - й )

+

 

 

 

 

1

( 1— т1) У-р

Рт <*А*

 

 

(56)

 

 

 

1—Л-п (1— Л)

 

 

 

 

 

 

 

 

 

 

Если в

уравнениях (54)—(56)

заменить

значения

К на

/сп,

то их

можно

переписать в

следующем

виде:

 

 

кп

1 —кп

 

1

~ ( п

1 - к в

/ *к п

\

|( йкп

1^кп

\

| 1

А’п^Хр

 

Арв-р<)' *

 

1

кП

Рт й Р й

 

>

 

1

А*п|Хр

(57)

 

1

 

 

\

1

А*п

Рт й р ь

 

'(ре-р,). 1

 

 

/

 

 

 

(58)

-

 

 

 

 

~ { ^ т ) г Г ^ Г (Рск ~ АпРп) * {Р е ~ Л ) + - 4 = ^ - ^

<5 9 >

Уравнения (54)—(56) или (57)—(59) позволяют оценить де­ формацию пористой среды с помощью четырех независимых параметров — рск, рп, рт и кп.

Исследования Ф. Гассмана [93], И. Гиртсма [94], В. М. Доб­ рынина [35], Ш. Нагумо [108] и других показали наличие связи между Рек, рт и рп. Например, для породы с однородным скелетом связь между этими величинами через пористость выражается уравнением

Р ^ Р х + а д .

(60)

На основании равенства (60), путем экспериментального опре­ деления рт и рп можно оценить величину Рск. В случае неоднород­ ного скелета эта связь выражается более сложным уравнением

Рек *■“ АпРп — Рт

1

[ МпЛ

(61)

1 —

\ *Р1 /(Ре-РдЫ'

 

 

Многие экспериментальные исследования физических парамет­ ров пород проводятся в условиях изолированности порового пространства от атмосферы [5, 25, 26]. Это соответствует тому, что поры породы «закрыты», и под действием нагрузки на скелет в жидкости, находящейся в порах породы, возникает неконтро­ лируемое норовое давление. На основании уравнения (48) можно