д.г.-м.н. Кашик С.А., академик МАНВШ, проф. Исаев В.П.
Институт земной коры СО РАН
Рассматривается модель формирования метана в донных отложениях озера Байкал в результате раннедиагенетических процессов. Показано, что величина равновесного давления в системе растворенный метан-газовая фаза прямо зависит от минерализации порового раствора: чем выше последняя, тем более высокое давление должно быть приложено для сдерживания газовыделений. То есть, чем глубже зашло взаимодействие вода-осадок, тем большее количество метана в растворенной форме накапливается в поровом пространстве донных осадков, и тем большие глубины водной толщи требуются для подавления газовых выбросов.
Эмиссия углеводородных газов из донных отложений озера Байкал известна давно. Наиболее многочисленные потоки метана в виде грифонов фиксируются на Селенгинском мелководье, вблизи устья р. Селенги и залива Провал). Выполненные подсчеты показали, что суммарный дебит газовыделений достигает в данном регионе 20 млн. м3
в год [1].
Локализация выбросов СН4
в придельтовой зоне озера обусловлена, с одной стороны, выносом больших количеств органического углерода р. Селенгой в виде растительного детрита (достаточно сказать, что ежегодно ею в составе взвеси сбрасывается от 57 до 630 тыс. тонн Сорг.
, или в среднем 196 тыс. тонн). С другой стороны, тем, что внутренняя ежегодная биогенная нагрузка седиментационного потока фитопланктоном в Байкале еще более впечатлительна и по данным [2] достигает около 4000 тыс. тонн.
Быстрому захоронению растительных остатков способствуют высокие скорости осадконакопления вблизи авандельты р. Селенги. Причем точность количественных определений величин скорости накопления отложений не вызывает никаких сомнений. В возникшем всего за два дня (с 13 по 14 января 1862 г.) заливе Провал выше затопленной почвы бывшей Цаганской степи за 100 лет накопился слой ила мощностью 1,8 - 3,6 м [3]. Следовательно, осадконакопление происходило здесь со скоростью от 1,8 до 3,7 см в год. При таких темпах осадкообразования растительный материал не успевал сгорать в хорошо аэрируемой байкальской воде и почти полностью захоронялся в донных отложениях. Еще одним дополнительным источником органического углерода в зоне авандельты могли служить погребенные торфяники бывшей Цаганской низменности, достигающие полуметровой мощности, а также широко развитые до глубин 1-1,5 м заросли водной растительности.
Формирование метана происходит в процессе диагенетических преобразований осадка, когда неравновесная система вода-осадок, постепенно трансформируясь на пути к стационарному состоянию, превращается в породу.
Как было показано ранее, интенсивное возрастание концентраций СН4
в поровых растворах начинается после процессов нитрат- и сульфатредукции [4], когда исчезают последние источники кислорода.
(CH2
O)106
(NH3
)16
(H3
PO4
)+ 84.8NO3
-
= 7.2CO2
+ 98.8HCO3
-
+ 16NH4
+
+
+42.4N2
+HPO4
-
+ 49.6H2
O (денитрификация) (1)
(CH2
O)106
(NH3
)16
(H3
PO4
)+ 53SO4
2-
= 106CO2
+ 16NH3
+ 53S2
2-
+ H3
PO4
+
+ 106H2
O(сульфатредукция) (2)
В противном случае, выделяющийся в результате указанных процессов кислород будет продолжать окислять органику и генерировать углекислый газ. Метанообразование осуществляется на последней стадии раннего диагенеза в результате реакции диспропорционирования [5], в этом случае органический материал (формула Редфилда [6]) разлагается по схеме:
(CH2
O)106
(NH3
)16
(H3
PO4
) = 53CH4
+ 53CO2
+16NH3
+ H3
PO4
(3)
На внешнем крае Селенгинского мелководья в поверхностном слое осадков в среднем содержится 3% Сорг.
от терригенной части [7], то есть примерно 60%, поступающего органического вещества подвергается деструкции сразу ниже поверхности раздела вода-осадок.
Чтобы попытаться оценить количество образующегося в процессе раннего диагенеза метана, нами была построена имитационная модель взаимодействия донных отложений с байкальской водой, при этом предполагалось, что весь органический углерод утилизируется по реакциям (1-3). Таким образом, система закрывалась по отношению к источникам растворенного кислорода, то есть первоначально кислородсодержащие байкальские воды в процессе раннего диагенеза не подпитывались газами из озерного резервуара. Это вполне вероятно, так как по многочисленным наблюдениям инверсия окислительного режима на восстановительный происходит уже в самых верхних частях осадках, на глубине от первых миллиметров до первых сантиметров ниже поверхности дна.Так по данным [8] мощность окисленного слоя в районе дельты Селенги колеблется от 0,6 до 23 мм.
Моделирование проводилось с помощью программного комплекса "Селектор" [9] при температуре придонной воды (3,8о
С) и различных величинах давления, имитируя протекание процессов взаимодействия вода-осадок на различных глубинах. В расчетах использовался уточненный состав байкальской воды [10] и усредненный химический состав верхней части байкальски
Результаты моделирования показывают, что величина равновесного давления в системе растворенный метан-газовая фаза прямо зависит от минерализации порового раствора: чем выше последняя, тем более высокое давление должно быть приложено для сдерживания газовыделений. Другими словами, чем глубже зашло взаимодействие вода-осадок, тем большее количество метана в растворенной форме накапливается в поровом пространстве донных осадков, и тем большие глубины водной толщи требуются для подавления газовых выбросов. Наибольших величин общей минерализации и содержаний растворенного метана поровые растворы достигают, если в процессе диагенетических реакций перерабатывается примерно половина исходного осадка. Это, безусловно, крайний вариант, но мы намеренно представили такой вариант для большей наглядности (рис. 2). В то же время следует заметить, что взятые в расчет 3% органического углерода, как среднее содержание в донных осадках, не совсем адекватно отражают истинное количество захоранивающегося растительного материала, поскольку часть его (возможно значительная) могла быть утилизирована в постседиментационных процессах.
Рис. 2. Растворимость (S) метана (квадраты) и равновесие в системе метан-вода (треугольники) в зависимости от общей минерализации (M) порового раствора и общего давления.
Вблизи линии равновесия СН4 раствор
- СН4 газ
в осадке формируется так называемый активный слой [11], вблизи верхней границы которого начинают формироваться пузырьки газа.
Из-за возникающих градиентов концентраций и давлений растворенный газ стремится мигрировать из зон пьезомаксимумов в зоны пьезоминимумов и таким образом начинает выделяться в газовую фазу. Кроме того, граничные условия равновесного давления в системе СН4 раствор
- СН4 газ
могли меняться из-за сезонных понижений уровня озера, которое может составлять более одного метра, что в свою очередь также способствует процессам газовыделения.
На основе построенной модели метаногенерации в донных отложениях, можно приблизительно оценить возможные количества выделяющегося газа в зонах разгрузки. При давлении 21,5 атм. концентрация метана в растворе составляет 13,74 мг/л или 0,904 мл/л. В условиях падения давления до нормального, а это происходит на поверхности раздела вода-атмосфера, объем газа увеличится в 21,5 раз и составит 19,44 см3
. Такое количество газа выделяется из 1л раствора. Учитывая пористость верхнего неконсолидированного слоя байкальских осадков, которая по данным [12] колеблется в пределах 87-90%, из 1000 см3
отложений выделяется ~ 17 см3
газа. Несложные расчеты показывают, что в таком случае с одного квадратного метра поверхности дна и десятисантиметровой мощности активного слоя может выделиться 170000 см3
метана. Это сопоставимо с замерами выделений метана из керна в скважине BDP-96 [13].
По данным [11] количественно эмиссию газа из осадка можно оценить из полуэмпирического уравнения:
J(CH4
) = 139X2
/(1 - X), (2)
где J - поток газа в см3
/м2
день и Х - мольное количество газа.
Измерения показали, что в метановых выделениях на Байкале СН4
составляет в среднем 75 объемных процентов и остальная часть приходится на азот [1]. В таком случае, мольная доля метана составит 0,64, а ежедневный поток 158 см3
. Таким образом, эмиссия метана из донных отложений Байкала может продолжаться годами, что подтверждается натурными наблюдениями.
Работа выполнен при поддержке РФФИ, грант 02-305-65395, Министерства образования и науки, грант E02-9.0-50, гранта программы "Университеты России" УР09.01.011, Государственного контракта с республикой Бурятия N17-4/5-?.
Список литературы
1. Исаев В.П., Коновалова Н.Г.Б Михеев П.В./ Геология и геофизика. 2002. Т. 47. N 7. С. 638-643.
2. Вотинцев К.К., Поповская Г.И. В кн.: Круговорот вещества и энергии в озерах и водохранилищах. Листвиничное на Байкале, 1973. С. 75-77.
3. Казенкина Г.А., Ладохин Н.П. // Труды Вост.-Сиб. Геологического института. 1961. Вып. 3. С. 35-49.
4. Кашик С.А., Мазилов В.Н.//ДАН. 1991.Т. 316. N 4. С. 966-969.
5. Froelich P.N., Klinkhammer G.P., Bender M.L. at al.// Geochm. et Cosmochim. Acta.1979. V. 43. P. 1075-10-90.
6. RedfieldA.C.// Am. Sci. 1958. V. 46. P. 206-226.
7. Выхристюк Л.А. Органическое вещество донных осадков Байкала. Новосибирск: Наука, 1980. 80 с.
8. Martin P., Granina L., Martens K., Goddeeris B.// Hydrobiologia. 1998. V. 367. P.163-174.
9. Карпов И.К. Физико-химическое моделирование на ЭВМ в геохимии. Новосибирск: Наука, 1981. 247 с.
10. Кашик С.А., Карпов И.К., МазиловВ.Н.// ДАН. 1993.Т. 328. С.731-734.
11. Makhov G.A., Bazhin N.M./ /Chemosphere. 1999. V. 38. P. 1453-1459.
12. Мизандронцев И.Б. В кн.: Проблемы Байкала. Новосибирск: Наука, 1978. С. 33-46.
13. Кузьмин М.И., Калмычков Г.В., Гелетий В.Ф.. и др.// ДАН. 1998. Т. 362. С. 541-543.
14. Коллектив участников проекта "Байкал-бурение"// Геология и геофизика. 2000. Т. 41. N 1. С. 3-32.