Файл: Бабинец А.Е. Гидрогеологические и геохимические особенности глубоководных отложений Черного моря.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.06.2024

Просмотров: 129

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

внешней силы F поршни, перемещаясь навстречу друг другу со скоростью V, уплотняют образец, создавая внутри его давление Р. Скорость движения поршней постоянна ѵ = const. Образец уплот­ няется вследствие фильтрации порового раствора через торцевые фильтры.

Поместим начало координат в центр образца, а ось Z направим вдоль оси цилиндра. Для описания процесса уплотнения воспользу­ емся зависимостью Дарси — Герсеванова:

U,-wt = - k - ^ - .

(3,18)

Здесь Uz — скорость фильтрации в направлении оси Z, ѵг скорость перемещения твердой фазы (скелета грунта) в направлении оси Z, е — коэффициент пористости, k — коэффициент фильтрации,

дН

 

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

гидравлический

градиент по оси Z.

 

 

 

 

 

 

 

Запишем

зависимость

Дарси — Герсеванова для элемента пло­

щади сечения

одного

из

торцов образца. В данном случае

имеем

 

 

 

Ut>0,

 

ѵ2<0, \иг\

= \ѵг\,

- ^ - < 0 .

 

 

 

 

 

Уравнение

(3,18)

примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

h

 

 

 

 

( З '

1 9

)

 

 

 

 

 

 

 

U + eU =

 

dz

 

 

 

 

 

Принимая

во

внимание,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{У =

 

V

и

дН

 

1

дР

 

 

 

 

 

 

 

 

 

dZ

 

у

ÔZ '

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дР

 

2Р-

 

 

 

 

 

 

а средний

градиент

давления - ^ - ~

- уЧ

поскольку при

Z = О

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

Р = Р( при Z = j

Р =

 

0,

уравнение

(3,19) можно записать так:

 

 

 

 

 

\

(1+е<) =

^ -

. - ^ .

 

 

(3,20)

 

 

Отсюда для коэффициента фильтрации

получим следующее вы­

ражение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*« = v t h t l î

+ e t ) .

 

 

(3,21)

 

 

 

Индексы

у

величин,

входящих

в

выражения

(3,20)

и

(3,21),

обозначают,

что они взяты для фиксированного

момента

времени.

Подробнее об их смысле будет сказано далее.

 

 

 

 

 

 

Сразу же заметим, что использование среднего градиента

давле­

ния является лишь первым приближением в решении задачи. Соглас­

но теории фильтрационного уплотнения,

давление в толще

грунта

распределяется нелинейно и зависит как

от

координаты Z, так и

от времени. Точный градиент давления

в

уплотненном

образце

здожно определить, лишь решив дифференциальное уравнение уплот-

42


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

Рассчитанные по (3,21) коэффициенты фильтрации осреднены в размерах испытываемого образца и несколько завышены против коэффициентов фильтрации осадков в условиях их естественного

Рис. 13. Автоматизированная установка для компрессионных испытаний донных осадков (блок-схема):

ТД — тензодатчик, ЯЧ — ячейка, ИГП-10 — пресс, ДП — датчик перемещения, У/77" — усилитель постоянного тока, БП — блок питания, Р — регистратор.

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

Описанная схема компрессионных испытаний образцов донных осадков была реализована на установке (рис. 13) с прибором одно­ осного сжатия типа ИГП-10. Это механический пресс с постоянной скоростью подачи штампа, оборудованный тензодатчиком нагрузок

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

ивеличина сжатия образца регистрировалась самопишущим элект­ ронным потенциометром типа Н-361. Образец помещался при помо­ щи цилиндрического пробоотборника в ячейку, подобную описанной выше. В зависимости от литологических особенностей осадка дли­ тельность сжатия образца составляла от 10—15 мин до 1—1,5 ч. Таким образом, на одной установке за 8-часовую смену можно

43

исследовать от 6 до 20 образцов. Это во много раз выше производитель­ ности обычных компрессионных приборов.

Пресс ИГП-10 способен создать нагрузку до 2000 кГс и имеет в комплекте сменные тензодатчики на диапазоны 0—100, 0—250, 0—500, 0—1000 и 0—2000 кГс. Тензодатчик 0—100 кГс обеспечива­ ет наиболее точное измерение усилий в области начальных нагрузок,

однако при этом нельзя повышать нагрузку

более 100 кГс, что не

дает возможности отжать достаточное для

дальнейших анализов

количество порового раствора.

 

 

 

В наших исследованиях использован диапазон

нагрузок

0—

500 кГс, что соответствует давлению в ячейке 0—67

кГс/см2.

При

6,3 5,3. 4,3 ДЗ 2,3 7,3 0,3 40 30 20 10 0t,MUH

Рис. 14.

Диаграмма сжатия образца (станция 1070, 660—675 см, ѵ =

= 0,26

мм/мин).

этом выделялось достаточное количество порового раствора (1 — 15 мл). Для расчета коэффициента фильтрации определяли удельный вес порового раствора, выделенного из каждого испытываемого об­ разца; методика описана в главе IV.

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

граммы записываются метки перемещения штампа пресса.

Каждая

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

скорости

мм/мин). По мере сближения штампов пресса нагрузка

на обра­

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

44


Выбирая фиксированные значения нагрузок F, рассчитываем соответствующие значения давлений Pt:

Pi = -Ç ,

(3,22)

где S — площадь сечения ячейки.

 

В формуле (3,21) для коэффициента фильтрации ѵ =

const. При­

нимая это во внимание рассчитываем постоянные коэффициенты () для выбранных ступеней нагрузки F{:

С І ~ - Щ - А .

(3,23)

Здесь А вводится для выражения коэффициента фильтрации в подходящих единицах скорости.

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

фиксированных нагрузок. Эти числа меток

обозначим /. Принимая

во внимание скорость ѵ, рассчитываем h:

 

К = hk + vi,.

(3,24)

Коэффициент пористости е можно рассчитать, используя значе­ ния конечной пористости nit равной по определению

ПІ =

-

£ - ,

(3,25)

где Ѵк — конечный объем керна,

Ѵп

— объем пор керна.

 

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

" • - - Ж -

< 3 ' 2 6 )

Здесь ЛѴ,. — приращение объема;

 

AVt = VltS.

(3,27)

Коэффициент пористости, по известному соотношению, будет

равен

 

в < = - Т ^ Г .

(3.28)

Определив для каждой ступени ht и ес,

рассчитаем коэффициент

фильтрации по формуле

 

k = yht(l+el)Cl.

(3,29)

В принципе можно было бы исходить в расчетах из начальных значений ен и hH, однако точность их определения гораздо ни­ же, чем точность определения ек и /ік . Это обусловлено тем, что

45


зачастую образцы осадков имеют весьма малую прочность и вы­ сокую влажность.

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

Рт =

(3-30)

где Р — вес конуса, h — глубина погружения, а — постоянный коэффициент, зависящий от угла конуса. Для тридцатиградусного конуса а = 1,108.

Теплопроводность. Для определения теплопроводности донных осадков использовались запарафинированные монолиты с условно ненарушенной структурой; измерения выполнены методом мгно­ венного линейного источника с использованием термопары константан — медь. Нагреватель и термопара вводились в образец путем прокалывания парафиновой стенки, которая обеспечивала изоля­ цию образца от внешней среды в процессе измерений.

Теплопроводность Т рассчитывалась по выражению для распре­ деления температуры при распространении теплового импульса от линейного мгновенного источника [152]

 

т = - щ ш г е х р

-Щ- •

( 3 ' З І )

где Т — температура,

вызываемая

мгновенным линейным

источни­

ком, Q — количество тепла, выделяемое единицей длины источника,

г — расстояние точки

измерения температуры от источника тепла,

т — время, р — плотность среды.

 

 

 

Применяя для этого выражения условие максимума

 

 

 

 

= 0 ,

 

получим

от

т=>т м а к с

 

 

 

 

 

 

 

 

k

= -4Ï—-

 

(3,32)

После подстановки

 

 

макс

 

 

k в (3,31) имеем

 

 

 

т

_

Q

. - 1

 

или

 

рСя/-2

 

 

 

Qe-1

 

 

 

 

 

 

 

 

макс

 

 

Теплопроводность

находим по формуле

 

 

 

К = kpC.

 

(3,33)

46


Подставляем сюда значения k, рС:

 

 

Qe—l

 

Qe~

 

(3,34>

 

4т„

пяг'

4я.т

Т

 

 

 

 

Q рассчитывается

 

^

макс' макс

 

по закону Джоуля — Ленца:

 

 

 

Q - = 0 , 2 4 - ^ - ,

 

 

 

(3,35)

 

 

 

 

 

 

Т а б л и ц а 3

Теплопроводность глубоководных глинистых

илов

Черного моря

 

Стан­

 

 

 

 

 

Теплопроводность,

Интервал, см

Характеристика ила

 

• 10-'

кал

ция

 

;

 

 

 

 

 

 

см • град • сек

1070

115—130

Темно-серый

однородный

 

2,40 ±0,05

1309

96—120

Черный

однородный

2,05 ±0,03

1627

180—195

Зеленовато-серый однородный

2,25 ±0,17

 

355—370

Темно-серый

»

 

1,90±0,19

1629

194—203

Зеленовато-бурый

тонкосло­

2,20 ±0,11

 

 

истый

 

 

 

 

407—420

Светло-серый

однородный

1,85±0,18

1633

70—85

Серый

однородный

2,40 ±0,06

 

220—234

Темно-серый

»

 

2,30 ±0,06

где / — сила тока, питающего нагреватель, V — напряжение, / — время действия источника тепла, / — длина нагревателя.

Таким образом, получаем формулу для расчета теплопровод­ ности

.

0,24/W

 

(3,36)

Іёпг

Т

Г *

 

 

4 е п т м а к с '

макс'

ѵ

'

В процессе измерений определялись максимальная температура Т'макс и время Тмакс (табл. 3) на установке, схема которой изобра­ жена на рис. 15.

Выводы. Полученные нами результаты определения некоторых водно-физических свойств глубоководных донных отложений Чер­ ного моря приведены в табл. 4 и 5 и в виде графиков на рис. 16—20.

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

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

47