Файл: Борисов, О. Г. Экструзии и связанные с ними газо-гидротермальные процессы.pdf

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

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

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

Добавлен: 18.10.2024

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

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

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

средних Y от X (регрессия),

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

Для определения формы связи между вариацией при­ знаков X и Y нами найден ряд

корреляционных уравнений (см. табл. 41)_. За независимую переменную Y x принято значе­

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

ляющим наклон линии регрес­ сии. Проверка по определению доверительной области, в кото­ рой с вероятностью (1—р )2 ле­

жит линия регрессии для выборки по базальтам, показа­ ла, что отклонение составляет не более +12,41. Тот факт, что исследуемая зависимость пред­ полагается линейной,для оцен­ ки силы (тесноты) связи моле­ но использовать выборочный коэффициент корреляции.

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

За результативный при­ знак принимаем общую газонасыщениость пород (в объе­ мах на объем породы); в каче-

 

 

3- С .С “

 

 

о

 

о <

 

 

 

J l O v

 

 

 

о

 

®

 

о

 

 

 

ч

 

 

 

 

 

 

си

 

 

 

 

 

я

 

 

 

 

 

X

 

 

 

 

 

А

 

 

 

 

га

 

 

са

 

 

S

 

 

 

 

 

к

м

 

о .

 

 

о

 

а О

1

 

са

 

Ж

 

 

та

 

 

 

 

X

^

 

 

 

 

о

о

g -H

o'

а

*

Ё $ £ ,

;>>Е

 

 

 

 

 

?

>

ь +

 

 

 

 

 

 

X

S

5

а

 

а

 

5 II

 

si'

*9-

 

Ч w

 

 

 

 

й х

о

 

 

 

*

 

 

 

°-s

 

 

 

 

« S

 

 

 

 

« п

с

>-

5 ©■

 

 

 

••

яс

КOJ

§ 5

си

 

 

 

 

|

|

X

 

А

 

* g

 

-f

8S С-

О

о ©

X

 

К

t

 

g

Л

X

 

•S-

 

са ’Ь*

 

н

 

 

 

£

>

о

 

§

>

X

>-

2

>•

н

 

и

 

I *

 

 

X

g

 

 

х

2

N-.

 

са

-

 

я g*

5

v*

Си c

> i j5-

 

 

S

J9

 

a

ca

с .

 

H

x

*

 

x

H

 

О

CO 2

 

s *T

 

 

X «0

 

L.

s 5

о

-в-о

 

(Dо.

*9*

 

 

 

.Т)

 

О

 

 

E-*

0)

К

3

о

 

С.

w

eO

га

 

X

с

 

еГ

 

 

л

 

 

 

X

 

 

 

о.

Л

 

со

 

 

а

 

et

 

 

О.

 

 

 

а

 

 

 

Е

 

2 2

Г--

to

CD

03

см

CD

 

 

 

 

 

t"-

CO

CD

oo

_

Г )

см

 

CD

CD

CO

СО

оо

CD

 

D -

 

см

СО

to

о CJ>

О

о

о о

о

c d CM —4t"-

оо

00

со

to

00

LO

 

со

СО

со

CD

00

 

 

о

со

С.М

c ”

О

О

о о о

о

 

 

 

 

1

1

1

Ю

in

CD

CDCD

to

05

oo

CD

CM

CO

h -

>

CD

00

Ю

о

 

О

о

см

о О

о о

о о о

 

1

 

 

1

 

 

 

CO

 

о

о

о

 

CO

00

00

r^.

1ГЛсо ю

m

см

о

а*

CO__ CM

 

Г"- сосм

П- to

 

1

 

 

со

i

 

 

 

 

 

 

 

 

 

 

 

tO

CD

D- , ,

00

LO

CD

CD

to

Tf*CD

т-

CD

CD

 

CD

О»

t"-

оз

 

o

о О

о

о

о

I 1 1

1

1 1 1

CO

oo

 

h-

о

—•*

 

C-)oo COсм

 

o'

l>

о О

СМ

 

 

 

о o' с>

 

1

1

1 1

1

CO

oo

см

см

С5

CO

CM

со

о

CM

to

r -i

 

сч

о

О

о О

О

о о

 

1

 

 

1

 

1

Ю

CO

CO

oo

СО

ю

05

CD

CD

CD

CM

» п

03

CD

COCD

CD

оо

CD

О

О

о О

о о

о

1 1

l

I

1

1

1

 

to

to

oo

ю

CD

 

CO

CD

CO

 

ю

 

 

CO

 

»—1

to

ю

 

 

 

о

 

о О

оо1о о

 

Ю

CO

 

h-

to

ю

со

LO

c 5CM

CM

со

 

00

r- о

о

ю

ю

о о о

о о о о

1

 

 

 

 

 

 

H H о

н

н

А

та.

со

X

a.

о.

ч

ч

VO

о

VO

<

та

та

 

 

О

о

СП

u

 

 

X

я

та

 

 

 

O i

<

ю Cf.

9*

13Г:


стве предопределяющих факторов рассмотрим содержание Н2 и СОо (в об. % от общей суммы газов). В табл. 43 сведены данные по множественной корреляции для некоторых разновидностей по­ род. Исследуемая зависимость между указанными факторами ли­ нейная и выражается уравнением

Y Xv = а04-.ахХ -f a2V .

(14)

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

Ба з а л ь т ы. Корреляционная связь выражена четко. Средняя ошибка коэффициентов корреляции не превышает максимальной ошибки (.Мтах) при принятом уровне значимости. Поэтому корреля­ цию считаем вполне достоверной. Коэффициент множественной кор­ реляции (^лт-г=0,68) показывает, что связь результативного при­ знака Y с двумя факторами-признаками (X и К) по абсолютной величине более тесная, чем с каждым фактором в отдельности (rXY

и />у). Исчисленные парные коэффициенты корреляции характери­ зуют степень тесноты связи между соответствующими двумя при­ знаками без исключения влияния третьего. По абсолютной величине они почти равны между собой (гд-г= 0,565; rvT= 0—0,659), что

заставляет предполагать равноценность влияния на суммарную газоиасыщенность (У) как фактора Л'(Н2), так и фактора К(С02). Для проверки этого предположение рассмотрим тесноту связи меж­ ду результативным признаком (У) и одним из факториальных признаков {X или V) при исключении влияния другого факториаль­

ного признака (см. табл. 43). Сопоставляя частные коэффициенты корреляции (rxY,v и rVytx)c вычисленными ранее парными коэффи­ циентами корреляции (гху и rrY) отмечаем, что основное влияние

на суммарную газонасыщенность пород оказывает С02. Влияние Н2 очень мало.

Сопоставление частных и парных коэффициентов корреляции, вычисленных для других пород, показывает следующую зависи­ мость общего газосодержапия пород от Н2 и С02: 1) прямую и рав­ ноценную (по абсолютной величине коэффициента корреляции) — для гранита и риолита; 2) обратную и равноценную — для диорита, андезита и диабаза; 3) зависимость для габбро аналогична базальту.

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

с газовой фазой.

132


Г а з о н а с ы щ е н н о с т ь п о р о д э к с т р у з и й . Небольшое количество проб по однотипным потокам и куполам с конкретных вулканов не позволяет дать статистической оценки различия в их газосодержании. Отдельные пробы показывают, что в куполах общее газосодержание в 1,5—2 раза меньше, чем в потоках, в стекле поро­ ды по сравнению с минералами газа меньше в 3—4 раза. Так, дацитовый поток вулкана Зимина содержит 1,45 V/V газа, а купол

того же состава — 0,97; андезитовый поток вулкана Шивелуч — 2,4, а купол — 1,15; стекло дацитового купола кальдеры Узон — 0,37, а все порфировые минералы — 3,5. Разброс значений по отдельным вулканам большой, но наблюдается общая тенденция к уменьшению содержания газа в экструзивной лаве.

Г Л А В А 14

ТЕОРЕТИЧЕСКИЙ СОСТАВ МАГМАТИЧЕСКИХ ГАЗОВ

В свое время А. Готье (ряд работ с 1900 по 1910 гг. цит. по А. 1. Ellis, 1957), занимающийся изучением газов, выделенных из горных пород, высказал предположение, что магматические газы близки к состоянию химического равновесия. Несколько позднее С. X. Вильсон (Wilson, 1938) пришел к выводу о термодинамической возможности существования некоторых соединений, таких как дву­ окись серы, метан, сульфид карбонила, бисульфид углерода в фумаролах и горячих источниках, имеющих различную температуру. Он показал, что отсутствие двуокиси серы и сульфида карбонила в газах горячих источников соответствует теоретическим предполо­ жениям. В то же время Е. Шепард (Shaphard, 1938) не разделял идею непосредственного применения термодинамики для объясне­ ния данных, полученных при анализе горных пород, считая, что сочетание катализаторов и антикатализаторов и действие адсорбции

вбоковых породах оказывает влияние на установление равновесия

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

Проанализировав имеющийся материал по вулканическим га­ зам и газам пород, А. Эллис (Ellis, 1957) отметил, что было бы небезынтересно установить теоретическое сочетание элементов в не­ которых системах с приблизительным составом естественных вулка­ нических и фумарольных газов и газов горячих источников. Теоре­ тические ассоциации молекул можно было бы затем сравнить с ассоциациями молекул в естественных газах. Это показало бы, на­ сколько близки естественные газы к состоянию равновесия. А. Эл­

133


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

Работы в этом направлении были продолжены. Изменение хи­ мического состава магматических газов в зависимости от изменения температуры и давления, но уже при взаимодействии газообразных компонентов с твердыми компонентами магмы, рассмотрел С. Мацуа (Matsuo, I960). Термодинамический аспект вулканических га­ зов был изложен в работе И. Шимаци (Shimazu, 1960), в которой он подтвердил правильность расчетов С. Мацуа. Работы по определе­ нию состава магматической газовой фазы или отдельных ее состав­ ляющих на основе термодинамических расчетов и обоснование са­ мих расчетов (Краускопф, 1961; Волков, Рузайкина, 1969; Каржавин, Вендилло, 1970 и др.) также подтвердили возможность при­ менения термодинамики для расчета теоретических и природных газовых систем.

Результаты исследования

Одним из наиболее удобных параметров при расчетах равно­ весных процессов принято считать константу равновесия (К). Кон­

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

Для определенных значений температуры и давления величина К может быть вычислена по значениям энтальпии АН, свободной энергии при постоянном давлении ДF или объеме AZ, и по завися­

щим от температуры теплоемкостям участвующих в реакциях ком­ понентов. Значение К через AZ выражается уравнением

A Z = —RT\r\K

(15)

При заданных Т и Р величина К имеет постоянное значение. При расчетах ДZR и KR различные исследователи используют не всегда сопоставимые термодинамические параметры {AZ, АН, S),

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

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

Из работы А. Эллиса (Ellis, 1957) наибольший интерес пред­ ставляет рассмотренная им система Н20 — Н2 — S2 при молекуляр­ ном отношении 100 : 2 : 1, которая стоит ближе всего к среднему от­ ношению в вулканических газах

S2+4H 20 =fM H 2+ 2 S 0 2,

(16)

3S2+4H 20 =pt 4H2S + 2 S 0 2.

(17)

134


Т а б л и ц а -14

г , ° к

Р. атм

■YH : 0

■Y H ss

'Ys o ,

Л'н 2

YS.

6 0 0

1

9 8 0 , 3

1 9 , 6

0 , 0 4

0 , 0 9 9 0 , 0 0 0

800

1

9 7 8 , 8

1 9, 1

0 , 5 2

1 , 6

0 , 0 0 9

1000

1

9 7 2 , 4

1 6 , 8

2 , 7

8 , 1

0 , 0 4 9

1200

1

9 5 7 , 9

1 1 , 7

7 , 5

2 2

, 8

0 , 1 0 5

14 00

1

9 3 8 , 8

5 , 3

1 3 , 9

4 1 . 9

0 , 0 8 0

12 00

10

9 6 6 , 8

1 5 , 0

4 ,5

1 3 , 7

0 , 0 4 8

1200

100

9 7 3 , 5

1 7 , 2

2 , 3

6

, 9

0 , 0 0 8

1200

ю о о

9 7 5 . 8

18,1

1 , 5

4

, 6

0 , 0 0 6

Результаты расчетов, выраженные в мольных фракциях (X ■103)

приведены в табл. 44.

При низких температурах и давлении в одну атмосферу система представляет собой почти бинарную смесь, при высоких температу­ рах эта смесь будет уже тройной, состоящей из Н20, S 0 2 и Н2. Уве­ личение давления способствует сохранению системы с наименьшим количеством молекул. При высоком давлении и температуре система почти полностью отвечает составу при 900° и давлении в одну ат­ мосферу. Если эта закономерность прослеживается и далее, то при температуре 1200° К по мере увеличения давления система должна приближаться к составу, отвечающему выделениям высокотемпера­ турных фумарол при нормальном давлении.

Расчеты для наиболее реальной природной системы с участием С02 (Н20 — С02 — Н2 — S2, при молекулярном отношении 100:10: :2: 1) для температур от 500 до 1400°К при давлении в одну ат­ мосферу показывают следующую закономерность изменения соста­ ва газов. По мере увеличения температуры уменьшается содержа­

ние С 02, H2S и увеличивается Н2, С02, S 02. Отношение,

равное 1

для H2S и Н2,

S 0 2, СО, отмечается при температурах около 1100,

1250 и 1350° К,

соответственно. Постоянство отношений

H2S : С02

выдерживается

до 800° К, затем лишь незначительно

начинает

уменьшаться. Отношение S 0 2: С02 по мере повышения температуры быстро возрастает. Приведенные данные неплохо согласуются с наблюдениями на вулканах.

Интересные термодинамические расчеты проведены для углево­ дородных газов Хибинского и Лавозерского массивов на Кольском полуострове, геологические условия залегания которых в извержен­

ных

горных

породах

и средний

изотопный

состав углерода

(6С =

—1,00)

этих газов, близкий к магматической С 02, свидетель­

ствуют об абиогенной

их природе

(Каржавин,

Вендилло, 1970).

Рассмотрена термодинамическая возможность образования метана в условиях земной коры из газовой фазы, содержащей СО и Н2 при

различных давлениях (1-—5000

атм)

и температуре

(500—1500° К).

Для упрощения расчетов использовалось уравнение

 

4СО+8Н2 ^

2Н20

+ С 0 2+ЗСН 4,

(18)

полученное суммированием ряда более простых уравнений, отража-

135