_ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА_
2020 Геология Том 19, № 1
ГЕОТЕКТОНИКА И ГЕОДИНАМИКА
УДК 530.311
Геодинамическая модель глубинного строения палеосубдукционной зоны на восточной окраине Русской палеоплиты и распределение месторождений нефти и газа
С.В. Гаврилова, А.Л. Харитоновь
аИнститут физики Земли им. О.Ю. Шмидта РАН
123242, Москва, ул. Большая Грузинская, д.10. E-mail: [email protected]
^Институт земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова РАН
108840, Москва, Троицк, Калужское шоссе, 4. E-mail: [email protected] (Статья поступила в редакцию 27 августа 2019 г.)
Использована известная гипотеза M. Barazangi о том, что квазилинейные геолого-тектонические зоны представляют собой полосчатые структуры, которые параллельны палеовулканическому Уральскому горному поясу и могут находиться в пределах некоторого интервала расстояний от пояса. По величине этого интервала (~103 км) и периодическому расположению квазилинейных геолого-тектонических зон (с поперечным размером ~ 300 км) оценивается скорость палеосубдукции (~ 5 - 6 см в год) на примере некоторых районов Сибири. Альтернативой столь высокой скорости палеосубдукции может быть угол палеосубдукции значительно меньше 10°. Результаты расчета могут быть применимы при изучении древних и современных литосферных блоков.
Ключевые слова: геодинамическая модель, палеосубдукционная зона, месторождения нефти и газа.
DOI: 10.17072/psu.geol.19.1.1
Введение
Известно, что построение надежной физико-математической модели глубинного строения любого региона имеет существенное значение как в теоретическом плане -для исследования тектонического строения этого региона, так и в прикладном - для прогнозирования поисков полезных ископаемых на исследуемой территории. Постановка задачи данной статьи связана с разработкой методики оценки некоторых глубинных гидродинамических параметров мантии Земли на основе определения расстояния от вулканической островной гряды, использующей известную гипотезу о том, что квазилинейные геолого-тектонические зоны представляют собой полосчатые структуры, которые параллельны горному вулканическому поясу, расположенному у ребра мантийного суб-дукционного клина, и могут находиться в
пределах некоторого интервала расстояний от горного вулканического пояса. Первые исследования в этом направлении были сделаны зарубежными коллегами на примере суб-дукционной зоны в районе литосферной плиты Наска (Вага2ап§1 й а1., 1976), а также по несколько другой, но подобной методике оценки (Кап§, 1971). На территории России имеется несколько регионов, где по параметрам расположения островных вулканических дуг экспериментально могут быть проверены основанные на рассматриваемой гипотезе теоретические расчеты углового направления субдукционного погружения литосфер-ных или палеолитосферных плит. Одним из самых важных регионов России, где можно показать преимущество расчетов геодинамических параметров мантии в зоне палеосуб-дукционного клина для прогноза зон расположения месторождений нефти и газа, является территория Западной Сибири (рис. 1).
© Гаврилов С В., Харитонов А.Л., 2020
Рис.1. Схема полосчатой структуры расположения субмеридиональных геолого-тектонических поясов Западно-Сибирского региона (заштрихованы). Геолого-тектонические (в том числе нефтегазоносные) пояса связаны с тепломассопереносом вещества по серии глубинных литосферных разломов от центров пар конвективных противоположно вращающихся вихрей Карига, возникших в астеносфере, прилегающей к палеосубдукционному мантийному клину, расположенных под «подошвой» ли-тосферного слоя Уральского палеоокеана (Западно-Сибирской палеоплиты): 1 - тектонические границы Западно-Сибирского региона; 2 - границы Западно-Сибирской углеводородной провинции; 3 -номера геолого-тектонических зон, составляющих пояса определенного вида углеводородных месторождений (нефть или газ или конденсат): 1-11-1Х - Приуральско-Фроловско-Южно-Ямальский газовый пояс, Ш-1У-УП-Х- Каймысовско-Среднеобско-Надым-Пурский нефтяной пояс, У-У1-УШ -Васюган-Пайдугинско-Пур-Тазовский конденсатный пояс, 4 - границы геолого-тектонических зон; 5 - номера геолого-тектонических зон; 6 - контуры месторождений углеводородов; 7 - обозначение разных типов месторождений углеводородов (Порфирьев, 1982)
Существенной особенностью локальных и региональных геолого-тектонических зон Зауралья и Западной Сибири является периодичность их пространственного расположения вдоль широтного направления (Ступак,
Берзин, 2001). Так, например, вытянутые в меридиональном направлении эти периодические региональные геолого-тектонические зоны представляют собой своеобразные полосы шириной ~ 300 км, параллельные
Уральскому хребту и расположенные с периодичностью примерно 300 км (Порфирь-ев, 1982).
Распределение гравитационных аномалий в редукции Буге (Петров, 2001) поперек поля напряжений Уральской геотектонической системы и ряд геологических аргументов в работе М.Г. Когана и Р.Ю. Кукулиевой (1988) свидетельствуют, что после формирования Уральского палеовулканического пояса в пермский период эволюционного развития Земли (295-240 млн лет назад) происходило взаимное поддвигание литосферы Русской палеоплиты (РП) под литосферу Восточно-Сибирской палеоплиты (ВСП) в районе окраинного палеоокеанического бассейна - Уральского палеоокеана, примерно соответствующего современному расположению Западно-Сибирской низменности (рис. 2).
Рис. 2. Палеотектоническая реконструкция формирования Уральского палеовулканического пояса в процессе формирования палеосубдукционной зоны при палеотектонической коллизии Русской и Восточно-Сибирской палеоплит (ВСП). Зоны палеосубдукции: 1 - палеозойские; 2 - мезозойско-кайнозойские; 3 - ныне действующие; 4 - зона сжимающих напряжений при взаимном поддвиге Русской и Восточно-Сибирской литосферных палеоплит в палеозойский период эволюции Земли; 5 - зона расположения Уральского палеовулканического пояса; 6, 7 - расположение некоторых крупных месторождений полезных ископаемых (Арешев, Гаврилов и др., 1997)
Хорошо известна аналогичная современная субдукционная зона, которая в данный
период только развивается в пределах Охотского моря - окраинной зоны Тихоокеанской литосферной плиты, поддвигающейся под литосферу Азиатского континента (рис. 3) (Гайнанов, 1991).
Ад, мГял
100 -
50 -
.Модель А
-50
Курильская дупка-
. жтба
Охотское море ничёская "гряда хр.Витт
100 200 300 Ш 500 600
100 -
200 -
300 -
Ш -
500 -
Е00 -
Н. км
И7 ИЭ
Рис. 3. Плотностная модель астеносферы и зоны субдукции Тихоокеанской литосферы в переходной зоне Охотское море - Курильская гряда -Тихий океан (Гайнанов, 1991). Кривые «модель А» и «модель В» - остаточные мантийные аномалии силы тяжести - dg после исключения из наблюденной аномалии силы тяжести влияния изменения мощности слоев земной коры, полученной по сейсмическим данным; 1 - гравитационный эффект эклогитового слоя литосферы (Лс =+0.02 г/см3), 2 - гравитационный эффект погружающейся плиты литосферы (Лс =+0.05 г/см3), 3 - суммарный гравитационный эффект плотностных неоднородностей литосферы и астеносфероы, 4 - «подошва» литосферы
На рис. 3 хорошо видны после основного максимума, расположенного над осью желоба, два экспериментально зафиксированных минимума значений поля силы тяжести, свя-
занные с вертикальными зонами изменения плотности пород в астеносфере, под Охотским окраинным бассейном. По мнению некоторых исследователей, эти вертикальные зоны разуплотнения в астеносфере также связаны с конвективными ячейками Карига, которые образуют вертикальные зоны теп-ломассопереноса вещества из мантии к поверхности Земли.
Предполагая, что геолого-тектонические зоны на территории Западно-Сибирской низменности (Уральский палеоокеан) образовались над восходящими потоками нагретого вещества (процессы тепломассоперено-са), возникшими над конвективными вихрями Карига в палеосубдукционном мантийном клине, можно оценить некоторые геодинамические параметры мантии в зоне пале-осубдукции (скорость палеосубдукции, угол наклона палеомантийного клина, скорость
смещения конвективных вихрей Карига во времени) в районе палеотектонической коллизии Русской и Восточно-Сибирской лито-сферных плит в палеозойский период развития Земли, которые позволили создать большое количество местрождений нефти, газа, газоконденсата на этой территории.
В соответствии с имеющимися научными представлениями о развитии этого региона (Порфирьев, 1982; Арешев, Гаврилов и др., 1997; Павленкова и др., 2014) мы хотели показать, что конвективная зона шириной 1000 км в палеомантийном клине (рис. 4) в палеозойский период эволюции Земли была расположена на восток от Уральского палеовул-канического хребта, что соответствует образовавшимся там (на территории Западной Сибири) линейным геолого-тектоническим зонам, показанным выше на рис. 1.
Рис. 4. Схематический глубинный разрез астеносферы и зоны палеосубдукции, образовавшейся в процессе палеозойского погружения литосферы Русской палеоплиты (РП) в переходной зоне Уральского палеоокенана (территория Западно-Сибирской плиты (ЗСП): 1 - Уральская палеовулканиче-ская гряда; 2 - направление движения конвективного потока нагретого вещества (в том числе и углеводородов) из мантии к поверхности Земли по зонам напластования и по глубинным литосфер-ным разломам; 3 - глубинные литосферные разломы, уходящие корнями в астеносферу, связанные с зоной субдукции; 4 - «подошва» литосферы Уральского палеоокеана (Западно-Сибирской палеоплиты) в палеозойский р) период эволюции Земли; 5 - схематическое изображение расположения трех пар конвективных (противоположно вращающихся в паре) вихрей Карига, создающих в центре пары конвективных вихрей восходящие к поверхности Земли потоки разогретого мантийного вещества (тепломассоперенос) (в астеносфере, прилегающей к зоне субдукционного клина, образованного между «кровлей» литосферы субдуцирующей Русской палеоплиты и «подошвой» литосферы Восточно-Сибирской палеоплиты (ВСП). Горизонтальная ось (х) на данном схематическом глубинном разрезе отображает важный параметр - горизонтальное расстояние (300 км) между центрами трех пар конвективных вихрей Карига в палеосубдукционном мантийном клине. На вертикальной оси ф на данном схематическом глубинном разрезе показаны глубинные отметки расположения различных палеослоев земной коры и верхней мантии в палеозойский период эволюции Земли без мощного слоя современных осадочных пород на Русской плите, соответствующие 0, 30, 60, 90 км
По нашему мнению, конвективные вихри возникали над погруженной Русской па-леоплитой в астеносфере над мантийным клином. Конвективные вихри Карига (или просто вихри Карига), поскольку их может быть несколько, при достаточно пологой зоне субдукции должны располагаться в мантии между субдуцирующим литосфер-ным блоком восточного обрамления Русской палеоплиты и Восточно-Сибирской па-леоплиты в астеносфере под литосферой Уральского палеоокеана. Иными словами, вихри Карига обычно должны быть расположены в мантийном клине (рис. 4) между «кровлей» субдуцирующей плиты и «подошвой» плиты, под которую происходит суб-дукция. Конвективные вихри Карига в пале-осубдукционном мантийном клине обычно бывают расположены парами и закручены навстречу друг другу. Один конвективный вихрь из пары закручен по часовой стрелке, а другой - против. Между этими двумя вихрями расплавленной мантии образуется восходящий мантийный поток тепла и вещества (тепломассоперенос). Он и выносит по глубинным литосферным разломам простейшие углеводороды (Н2, СО2, СН4) в приповерхностные зоны их накопления и последующего преобразования (различного типа геолого-тектонические «ловушки») в «осадочном» слое Земли (Сефуль-Мулюков, 2012). В нашем случае, при погружении Русской па-леоплиты под Восточно-Сибирскую па-леоплиту в палеозойский период эволюции Земли, существовали три пары конвективных вихрей Карига. Причем каждая конвективная ячейка, состоявшая из пары противоположно вращающихся конвективных вихрей Карига, имела средний горизонтальный размер около 300 км. Все три пары конвективных вихрей Карига должны были быть расположены в палеомантийном клине общей протяженностью в горизонтальной проекции на поверхность Земли около 1000 км к востоку от Уральского палеовулканического пояса и образовывать соответствующую им трехсоткилометровую геолого-
тектоническую зональность на территории Западно-Сибирской палеоплиты, сохранившуюся до наших дней (рис. 1).
Рассмотрим модель палеомантийного клина между «подошвой» Западно-
Сибирской палеоплиты и верхней поверхностью («кровлей») субдуктирующей Русской палеоплиты (рис. 4), в результате относительного движения Русской палеоплиты на восток пододвигающейся под углом ß со скоростью V. Верхнюю поверхность («кровлю») субдуктирующей Русской палеоплиты в первом приближении будем считать наклонной полуплоскостью (z = x x tg ß ), а «подошву» Западно-Сибирской палеоплиты - горизонтальной полуплоскостью (z = 0) соответственно, где x - горизонтальная ось, соответствующая почти горизонтальному направлению «подошвы» Западно-Сибирской палеоплиты, а ось z направлена вертикально вниз. Материал мантийного клина считаем однородной несжимаемой жидкостью с постоянным коэффициентом вязкости
)c , равным его среднему значению. Движение вещества в мантийном клине определяется формулами (Schubert, Turcotte, Olson, 2001; Turcotte, 2002), а уравнение стационарного распределения абсолютной температуры T в мантийном клине при бесконечном числе Прандтля в приближении Бус-синеска имеет вид (Vx, Vz - компоненты скорости вдоль осей x и z) (Vx>dxxT)+(Vzxdz >T)=(x xTx)+[Tk2/(2x)c )], (1) где Tk - тензор вязких напряжений, символ д обозначает частную производную по координатам, х - коэффициент температуропроводности, Л - оператор Лапласа. Согласно (Turcotte, 2002), имеем
U = Vx = ду /dz = - [ D x arctg (x /z)] + {x x [(Cx + Dz) ] I (x2 + z2)]} - B,
W = Vz = ду /dx = - [ C x arctg (x /z)] + {z x [(Cx + Dz) ] I (x2 + z2) ] } - A, (2)
Tlk2 / (2 x )) = [(4 x )C) x (Cz - Dx)2 ] I I (x2 + z2)2, (3)
где при условии непроскальзывания при z = 0 и z = x x tg ß на границах мантийного клина
A = Vx { [(ß x (л I 2) x sin ß J Z [ ß 2 - sin2ß J I; B = V x I [(л I 2) - ß] x sin ß - [ ß x (л I 2) x cos ß J I Z [ ß 2 - sin2ß J; C = - V x { [ ß x sin ß J Z [ ß2 - sin2ß J I; D = V x I [ (ßx cos ß ) - sin ß J Z [ ß 2 -sin2ß J (4)
Отметим, что при ß < (л I 4) мощность (З) диссипативного тепловыделения максималь-
на у поверхности субдуктирующей литосферы, т. е. у нижней границы мантийного клина.
Граничные условия в уравнении (1) для температуры есть Т = Тт при г = 0 и г = х х tg 3 , где Тт = 1500°К - температура субсо-лидуса. Граничные условия Т = Тт вдали от ребра мантийного клина (при х ю) также справедливы, поскольку тепловая энергия, соответствующая температуре, отличной от Тт , на достаточном удалении от хребта передается к литосферным плитам, т. к. дисси-пативное тепловыделение при х — ю спадает до нуля.
Формула (1) показывает, что максимальная температура Ттах в мантийном клине достигается у поверхности погружающейся плиты и с точностью ~ 10% может быть представлена аналитической формулой Ттах = Тт + {(]о х^/к) X (1/[F+(ОхVхx/%)]}, (5) где ^ и О - функции угла субдукции 3, которые при (к / 9) < 3 < (к / 4) могут быть представлены приближенными интерполяционными формулами
^ = - 2.70 + 15.75 х р; О = - 2.35 х 10-3 +
1.27 х 10 х р,
(6)
где к = 4х 10 эрг/см с,К - коэффициент теплопроводности пород астеносферы, % = [к / (р-ср)] = 10"2 см2/с, р = 3.30 г/см3 - плотность пород астеносферы, ср =1.2 х 107 эрг / (г ' К) -удельная теплоемкость пород. Согласно (1), (Ттах - Тт) пропорционально среднему коэффициенту вязкости Т]с. Например, при Т]с =1021 Пас = 1022 П, р = 30°, V = 10 см/год, х = 100 км, согласно (5)-(6), величина (Ттах -Тт) ~ 400° К.
В мантийном клине возможны конвективная и термическая неустойчивости, вызываемые перепадом (Ттах - Тт) температуры и возможной (не учитываемой выше) зависимостью вязкости от температуры. Исследуя эти неустойчивости, будем считать угол р субдукции литосферы Русской палеоплиты малым, а вынужденное течение в мантийном клине близким к плоскопараллельному. Если возмущения ух , уу , н компонент скорости, р - давления и в - температуры зависят от пространственных координат х, г и времени t как
(Ух , Уу , у , р, в ) = [«(О, у (С), »(С), р(0, в (С)] х ехр [yхt + г (ххкх + ухку - tхю )] , (7)
где С= (г / х) и в нулевом приближении у и с не зависят от х, при этом ось у направлена вдоль простирания Уральского палеовулка-нического хребта и коэффициент вязкости (щ ) зависит от температуры (Т) как Т = щ х ехр С; С= (Е* + р х V*) / ( Я х Т) (8) (где Е* и V* - энергия и объем активации, Я
- универсальная газовая постоянная), то с учетом
дх = дх - (С/ х)хд( ;
дг = (1 / х)хЭс;
д2хх = д2хх - 2х (£/х) х д2х£ + 2х (С/х2 ) х дс + (С
2/ х2 )хд2сс ;
Э2гг = (1 / х2 ) х д2сс ;
пренебрегая при С<< 1 величиной С, для возмущений (11) имеем
г кхи + г ку V + м>'/х = 0, (9)
0 = - гкххр-2хщхкх2хи-щхку2 х«-кххкх + {щх[-Сх(в/ххТ)хи'+(и'/х)+1 кх х»]}'/х, (10) 0 = - г кухр - 2хщхку2 XV - т]хкххкухи + {щх[^' / х) + г кух»]}' / х, (11)
0 = - (р'/х) - (aхgхpхв) + т]х{(г кх х и'/х) -(к2х») - [1 кххи'х(^хв/ ххТ)] +
{щх[(г kyXv'/х) - к/х»]} + 2 х (щх»1)' / х2, (12) рхсрх{(ухв) - (гхсхв ) + (м>хТ Ух + (г кх х и х в ) + [(Жх в') / х]} = к(в)'' / х2 - к2(в) + [(4хщх^2х^хв ) / (х2хТ)] + 4хщх^х{[(^'х^хв ) / (ххТ)] - (и'/х) - (г кх х»)} / х, (13)
где (9)-(13) есть соответственно уравнение неразрывности, три компоненты уравнения движения и уравнение переноса тепла, а штрих обозначает производную по С. Для решения системы (9)-(13) вводим новые неизвестные функции £:
Яе,1ш ЯеДш ЯеДш .,
у1,2 = и , у3,4= V ' , у5,6 = » , уу = ] х
{иКе / [х - кх х »1т - и х £ х (вКе / ххТ)]}, у8 = т] х {и1т / [х + кх х »Ке - и х £ х (в 1т / ххТ)]},
у9 = тт х [vRe / (х - ку х »>)], ую= Т] х [vIm / (х + ку х wRe)], у\\=pRe - (2х ] х / х ), у12 = р1т -(2хщ х »'1т / х ), уШ4 = вадт, у 15,16 = % х
в адт. (14)
Уравнения (9)-(13) относительно функций (14) превращаются в систему 1-го порядка:
у': = (ххкхху6) + [(хху7) / щ] + (£хи'/Т)ху13; у'2 =
- (ххкхху5) + [(х^у8) / Т] + (£хи'/Т)ху14;
у'3 = (ххкуху6) + [(хху9) / т]; у'4 = (ххкуху5) +
[(х^_ую) / ]]; у'5 = (х хкх ху2) + (ххкуху4);
у'6 = - (ххкхху2) - (ххкуху3); у7 =
ххщх [(4 хк2х+к2у )ху1] + (3 ххх щхкххкуху3) - (хх
кххуи); у'8 = ххщх [(4 хк2х+к2у )ху2] +
(3 ххх т хкх хку ху4) + (ххкхху„);
у'9 = хх ц х [(4 хк2у+к2х)ху3] + (3 хххцхкххкуху^ -
(х*куХуп);
у'ю = хх цх [(4хк2у+к2х)ху4] + (3 ххх цхкххкуху2) + (ххкухуи);
у'п = - (ххкхху8) - (ххкуху10) - (ахрх^ххху^); у'12 = (ххкххуп) + (ххкуху9) - (ахрх^хх^_у!4);
у' 13,14 = у'15,16 /
у' 15 = [х2х(х + кхХ)ху1з] + [х2 х (ю - (кх х и) у14]
+ [(ххГ)ху5] - [(4 ху^хцх^х 4) / (рх^хТ)] +
[(4 ху7хххО) / (рх^)] + [( у 15 х ххУх% 2) / (2х()];
у\б = [х2х(х + кх2хх)хуи] - [х2 х (ю-(кх х и) у1з] + [(ххТ')хуб] - [(4хумхцх£2х£) / (рхсхТ)] + [(4 ху^хххО) / (рхср)] + [( у1бх ххУх^ 2) / ^хр)] , (15)
где подставляем Ж = ( С 2х V / ( ) = 0 при С<< 1 (V - скорость субдукции). Свободные граничные условия для уравнений (15) в терминах неизвестных (14) запишутся в виде >>5,6,7,8,9,10,13,14 = 0 при £ = 0, г = tg ( . (16) Для оценки величин у и ю примем коэффициенты в формуле (15) постоянными, равными их средним значениям, в частности ц = цу Тогда задача, представленная в формулах (15) - (16), имеет решение: у, = Л, х{[ со/"1 ] Л хС} с Л = [п/ (1в()], I = 1, 2... 16.
Найдем условия существования нетривиального решения при ку = 0. В этом случае конвекция возникает в форме вихрей Карига, а уравнения для у3,4,9д0 образуют отдельную подсистему, из которой следует, что у3,4,9д0 = 0, т. е. возмущения скорости и сдвиговые напряжения поперек направления субдукции Русской палеоплиты отсутствуют. Оставшиеся уравнения, в пренебрежении малым Ж = (С х V / ( ), приводят к следующему условию для у и ю :
{(Лх) + [х2х(у + к2х%)] - д - {[(х2х кх2) / [(х2х кх2) + Л2]2}х[(aхрхgхx3хТ ') /ц] +
[(4ху7х^х^1хЛ2 / х)]}2 + {[х2х(ю- (кх х
и)] + [(х2 хкх) х [(х2 хкх2)-Л2] х(ВхТ' + ах рхgхД)] /
[(х2х кх2) + Л2]2}2 = 0, (17)
где В = 4 х и'/Т, Д = 4х х х [ Д - (/) х С х tg
( ] / (рхер), Д = 4х цх4х [Д - (/) хСх tg р ]2 /
(р хсрхТ), и, согласно формуле (5),
Т ' = [(ТИах - Тт) / tg (] = [Ц} х V 2) / (к х tg р )] х (1 / {^ + [(х х V х О) / 2 ]}). (18)
Из условия (17), эквивалентного двум независимым уравнениям, получаем г = tgP : Уь = У (ку= 0) = 4 х цх ^х [(Д - Схг / 2)2 / (х2хрхсрхТ)] + {(рхaхgхx3хТ 'хк2х) /ц}х [(х2х к2х)+Л2]2} + 16хц}х{и'хк2ххЛ2х4х[(Д - Сх г / 2) /
[(х2х кх2)+Л2]2хрхерхТ)]} - {%х [(x2x к2х)+Л2] / x2}, (19)
о' = о (ку=0) = kxXU- {[ кхх(х2хкх2- Я2)] / [х2хкх2 + Я2)2]}х{[(Т'х^'х£ ) /Т] + 4х [р х « х g х х х (D - Сх t / 2) / р х ср ]}, (20)
где a = 3 х 10-5 1/К - коэффициент теплового расширения, g - ускорение силы тяжести, Т ' = (Ттах - Tm) / t, t = tg 0 , С = (E* + рхУ*) / (ЯхТ), E* и V* - энергия и объем активации, R - универсальная газовая постоянная, U' = {[D х (sin рх cos 0 + 0)] - [С х (sin2 0 )]} / tg 0, Л2 = Я2 + 4 х л2, Л=[л / (tg 0)], С = - V х {[ 0 х sin 0 ] / [ 0 2 - sin0]}, D = V х {[(0 х cos 0 ) -Sin 0 ] / [ 0 / - sin20 ]}, Т = (Ттах + Тт) / 2. В формуле (19) зависимость от средней вязкости возникает из-за существенной вблизи желоба термической неустойчивости, описываемой первым слагаемым в (19) и связанной с экспоненциальной зависимостью вязкости от С , т.е. от температуры.
Из формул (19)-(20) видно, что в полученном первом приближении, которым мы ограничимся, инкремент у и частота о являются функциями горизонтальной координаты х и, таким образом, должны быть плавно меняющимися функциями, что естественно при малых углах (0 ) наклона зоны субдукции. Благодаря присутствию воды в мантийном клине вязкость его материала может быть очень мала (Жарков, 2012), r¡c = 3 х 1018 Пас и менее. Примем величину кх - волнового числа - такой, что пространственный период [(2х^)/кх], разделяющий соседние восходящие конвективные потоки в мантийном клине, равен среднему наблюдаемому расстоянию между соседними геолого-тектоническими зонами, параллельными Уральскому палевулканическому хребту, изображенными на рис. 1 косой штриховкой (Порфирьев, 1982), т.е. [(2хя) / кх] = 300 км. Считая £ = 20, что соответствует энергии активации E* = 240 кДж/моль и активацион-ному объему V* = 8 см3/моль (Жарков, 2012), угол поддвигания 0 ~ 10° субдуктирующей литосферы Русской палеоплиты под литосферу Уральского палеоокеана (ЗападноСибирскую плиту) (литосферные плиты в восточном направлении субдуктируют полого, например, плита Наска под углом 20 -30° (Barazangi et al., 1976)) и подставляя приведенные выше значения других параметров,
входящих в (19), находим распределения у как функции х при различных скоростях субдукции, показанные на рис. 5. Видно, что инкремент у положителен (т.е. конвекция может проявляться) при V > 5 см/год, а при V ~ 6 см/год конвектирущая зона занимает диапазон расстояний (х2 - х1) ~103 км от горного пояса. Эта величина хорошо согласуется со средней широтной протяженностью геолого-тектонических зон Западной Сибири (рис. 1) и, следовательно, дает оценку скорости палеосубдукции в области восточнее Уральского палеовулканического пояса (современный Уральский хребет). Эта оценка практически не зависит от принятой средней
вязкости ]с , поскольку движущий конвекцию перепад температуры (Ттах - Тт), согласно формуле (5), пропорционален ]с , но,
с другой стороны, чем больше ]с , тем труднее возбудить конвекцию. Это видно из формулы (19), где три последних слагаемых, описывающих конвекцию, очень слабо или
вообще не зависят от ]с , а первое слагаемое, описывающее термическую неустойчивость,
мало при ]с < 3 х1018 Пас.
1200 1600 1800 2000 2200 *2 2600
Рис. 5. Инкременты у конвективной неустойчивости в субдукционном мантийном клине как функции горизонтального расстояния х (в км) от Уральского палеовулканического пояса при различных скоростях V (в см/год) палеосубдукции. В интервале х\ < х < х2 положительных у можно ожидать скорость конвективного потока нагретого вещества в мантийном клине Русской палеоплиты V ~ 6 см/год
Обсуждение результатов исследований
Изложенная методика может представлять некоторый интерес для геофизиков, занимающихся изучением глубинного строения мантии в зонах субдукции, т.к. развивает определенный подход к геодинамике в зонах субдукции. В этой статье обсуждаются две основные темы: поток в мантийном клине и конвективная неустойчивость в нем. Выдвинутая нами идея о том, что потоком за счет вязкости генерируется достаточно тепла, может быть обоснована следующим образом. В частности, нужно отметить, что вопрос о прогреве холодного субдуцирующего блока в нашей работе не затрагивается. Он в принципе уже решен в статье McKenzie (1969). Именно субдуцирующий блок прогревается кондуктивной теплопроводностью из вмещающей мантии и к глубине 700 км практически полностью размягчается. Мы также учитывали дополнительный диссипативный нагрев в мантийном клине, где максимум температуры достигается вблизи верхней грани субдуцирующей плиты, над верхней гранью. Без учета конвекции в мантийном клине этот дополнительный нагрев, по сравнению с работой (McKenzie, 1969), несколько сильнее бы прогревал субдуцирующий блок, но благодаря конвекции диссипативное тепло уносится вверх (там, где инкремент конвективной неустойчивости положителен, т.е. где конвекция действительно развивается). При угле поддвигания ß субдуцирующе-го блока Русской палеоплиты порядка 10° горизонтальная протяженность области, где еще возможна конвекция, составляет немного более 1000 км при скорости около 5 см/год.
Кроме того, выделенная нами система субпараллельных структур на территории Западно-Сибирской плиты объясняется как действие конвективных вихревых структур мантии вокруг зоны субдукции. Безусловно, могут быть и некоторые другие предположения на этот счет. Приведем краткую сводку публикаций на эту тему: McKenzie, 1969; Schubert, Turcotte, Olson, 2001; Turcotte, 2002; Karig, 1971. Однако подтверждением изложенного нами механизма взаимодействия конвективных мантийных вихрей, про-
явившихся в приповерхностных геологических структурах, является то, что подобные структуры были обнаружены и в других районах субдукции Мирового океана (Kar1g,
1971).
Заключение
Таким образом, можно считать, что скорость палеосубдукции Русской палеоплиты с учетом расстояния до Уральского палеовул-канического пояса может быть заметно ниже полученного значения V ~ 5 - 6 см/год только при угле субдукции ( << 10°. Мы попытались сравнить полученные результаты расчетов с независимыми данными глубинных геофизических разрезов (Петров и др., 2001; Павленкова и др., 2014), пересекающих Уральский хребет и Западную Сибирь, чтобы выяснить полученный ими угол наклона погружения (палеосубдукции) Русской плиты. Результаты таких исследований приведены на рис. 6 (Павленкова, 2014).
Рис. 6. Фрагмент глубинного геофизического разреза (по сейсмическим и гравитационным данным) по профилю «Кварц» в зоне сочленения Русской и Западно-Сибирской плит (Павленкова и др., 2014), показывающий: а - кривые измерен-
ного и расчетного гравитационного поля, представленного в редукции Буге: 1 - (кривая коричневого цвета) - наземные измерения гравитационного поля вдоль профиля «Кварц» ; 1 - (кривая серого цвета) - результаты решения прямой задачи гравиразведки по профилю «Кварц»; 2 -(кривая черного цвета) - результаты расчета сглаженного гравитационного поля по профилю «Кварц», по данным альтиметрических аномалий, измеренных на космических аппаратах. б -глубинный разрез по профилю «Кварц», на котором видно наклонное погружение литосферного слоя Русской плиты (между нижней границей земной коры (М) - Мохоровичича и границей (N1) литосферы) под Западно-Сибирскую плиту, трактуемое как остаточные фрагменты па-леомантийного клина, возникшего при палеосуб-дукционном погружении Русской палеоплиты под Восточно-Сибирскую в районе зоны расположения Уральского палеоокеана в палеозойский период эволюции Земли. На глубинном разрезе цифрами обозначена плотность горных пород (г/см3) различных слоев земной коры и верхней мантии: М- граница Мохоровичича, являющаяся нижней границей («подошвой») земной коры; граница N1 литосферного слоя мантии; граница N2 астеносферного слоя мантии; Ь, Н, Т - другие более глубокие границы слоев разного плот-ностного состава верхней мантии. На данном геофизическом разрезе вертикальный и горизонтальный масштабы разные и поэтому при рассмотрении статьи необходимо ввести некоторый коэффициент пропорциональности, позволяющий из угла (Ь), изображенного на рисунке, вычислить истинный угол наклона ()) Русской плиты в настоящий геологический период эволюции Земли. Истинный угол наклона (() Русской плиты под Западно-Сибирскую в настоящее время составляет около 8°С по данным наших расчетов с учетом косинусоидальной поправки за неортогональность направления профиля «Кварц» к линии простирания гребня Уральского хребта
Из рис. 6 хорошо видно, что экспериментально определенный по геофизическим (сейсмическим, гравитационным) данным угол наклона современной литосферы Русской плиты (слой между границей Мохоро-вичича - Ми границей - N1) составляет в среднем не более 10°, что может свидетельствовать в пользу правильности проведенных в статье расчетов угла наклона пале-осубдукционной зоны Русской палеоплиты. Этот разрез палеосубдукционной зоны (рис.
6) по структуре гравитационного поля очень похож на разрез современной Тихоокеанской субдукционной зоны в районе Охотского моря, приведенный на рис. 3 (кривая 2). Конечно, за многие миллионы лет, прошедшие с палеозойского периода активизации пале-осубдукционной зоны в зоне коллизии Русской и Восточно-Сибирской палеоплит, произошли определенные геолого-тектонические изменения в переходной зоне сочленения этих палеоплит, на территории Западно-Сибирской плиты (древняя территория, называвшаяся Уральским палеоокеаном), но тем не менее многие специалисты в области геологии считают, что на этом геофизическом разрезе (Павленкова и др., 2014) можно видеть остаточные фрагменты палеозойского палеосубдукционного погружения Русской палеоплиты. Кроме того, частота с, согласно формуле (20), привязана к средней горизонтальной скорости Uc = (D х ß ) -(Cx [(1 - ß ) / tg ß ]} = 0 в мантийном клине, и при малых углах наклона ß зоны субдукции выполняется условие (у >> с ). Это означает, что конвективные вихри Карига в целом мало смещаются в горизонтальном направлении: (с /kx ) < 2*10-2 см/год. Изложенная методика может помочь при проведении региональных поисковых и разведочных работ на углеводороды в аналогичных районах в окрестности расположения субдукционных и палеосубдукционных зон Земли.
Выводы
1. Конвективная зона, состоящая из трех пар конвективных вихрей Карига, в пале-осубдукционном мантийном клине занимает диапазон расстояний около тысячи километров от Уральского палеовулканического пояса. Эта величина хорошо согласуется со средней широтной протяженностью трех геолого-тектонических зон Западной Сибири (около 300 км каждая), следовательно, дает возможность сделать оценку скорости палеосубдукции в области к востоку от Уральской палеовулканической дуги (современный Уральский хребет).
2. Конвективные вихри Карига, поскольку их может быть несколько при достаточно пологой зоне субдукции, располагают-
ся в мантии над субдуцирующим блоком и под надстилающей литосферой. Иными словами, вихри Карига обычно расположены в мантийном клине, в нашем случае, при погружении (палеосубдукции) Русской па-леоплиты, они были расположены к востоку от Уральского палеовулканического пояса на расстоянии 300 км друг от друга.
3. Скорость палеосубдукции Русской па-леоплиты в палеозойский период эволюции Земли с учетом расстояния до Уральского па-леовулканического пояса может быть заметно ниже полученного значения V ~ 5 - 6 см/год только при угле субдукции ß << 10°.
Библиографический список
Арешев Е.Г., Гаврилов В.П., Донг Ч.Л., Зао Н., Попов О.К., Поспелов В.В., Шан Н.Т., Шнип О.А. Геология и нефтегазоносность фундамента Зондского шельфа. М.:Нефть и газ, 1997. 288 с.
Гаврилов С. В. Исследование механизма формирования островных дуг и задугового раздвигания литосферы // Геофизические исследования. 2014. Т.15, №4. С. 35-43.
Гайнанов А.Г., Пантелеев В.Л. Морская гра-виразведка. М.: Недра, 1991. 214 с.
Жарков В.Н. Физика земных недр. М.: Наука и образование, 2012. 384 с.
Коган М.Г., Кукулиева Р.Ю. Изгибная жесткость литосферы Евразии // Доклады АН СССР. 1988. Т. 301, № 1. С.69-75.
Павленкова Н.И., Павленкова Г.А. Строение земной коры и верхней мантии Северной Евразии по данным сейсмического профилирования с ядерными взрывами.М.: Геокарт-Пресс, 2014. 193 с.
Петров А.В. Обработка данных гравиметрии по геотраверсу Уралсейс с использованием компьютерной технологии «КОСКАД 3D» // Глубинное строение и геодинамика Южного Урала. Тверь: Изд-во «ГЕРС», 2001. С. 79-84.
Порфирьев В. В. Особенности глубинного строения земной коры и теоретические обоснования неорганического генезиса нефти. Киев: Наукова думка, 1982. 328 с.
Сейфуль-Мулюков Р. Нефть и газ: глубинная природа и ее прикладное значение. М.:Торус Пресс, 2012. 216 с.
Ступак В.М., Берзин Р. Г. Структурные, кинематические и динамические характеристики Южного Урала по данным дифференциальной сейсморазведки // Глубинное строение и геодинамика Южного Урала. Тверь: Изд-во «ГЕРС», 2001. С. 43-51.
Barazangi M., Isacks B.L. Spatial Distribution of the Earthquakes and Subduction of the Nazca Plate beneath South America // Geology. 1976. Vol.4. P. 686-692.
Karig D.E. Origin and development of marginal basins in the Western Pacific // Journal Geophysical Researches. 1971. Vol.76, N.11. P. 2542-2561.
McKenzie D.P. Speculation on the consequences and causes of plate motions // Geophys. J. Roy. Astron. Soc. 1969. Vol, 18. P. 1-32.
Turcotte D.L., Schubert G. Geodynamics. Cambridge: Cambridge University Press, 2002. 448 p.
Schubert G., Turcotte D.L., Olson P. Mantle Convection in the Earth and Planets. New York: Cambridge University Press, 2001. 940 p.
Geodynamic Model of the Deep Structure of a Paleo Subduction Zone on the East Edge of the Russian Paleo Plate and Distribution of the Oil and Gas Deposits
S.V. Gavrilova, A.L. Kharitonovb
aSchmidt Institute of Physics of Earth of the Russian Academy of Sciences 10 Bolshaya Gruzinskaya Str., Moscow 123242, Russia. E-mail: [email protected] bPushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation of the Russian Academy of Sciences. 4 Kaluzshskoe Highway, Troitsk, Moscow 108840, Russia. IZMIRAN. E-mail: [email protected]
Known hypothesis of M. Barazangi that quasilinear geological and tectonic zones represent the banded structures, which are parallel to the Urals paleo volcanic mountain belt was used. These zones can be in some interval of distances from a mountain paleo volcanic belt. On the size of this interval (~ 103 km) and the periodical arrangement of quasilinear geological and tectonic zones (of about ~ 300 km width), the paleo subduction speed (~ 5 - 6 cm a year) was estimated on the example of some Siberian regions. Keywords: geodynamic model; paleo subduction zone; oil and gas deposits.
References
AreshevE.G., Gavrilov V.P., Dong Ch.L., Zao N., Popov O.K., Pospelov V.V., Shan N.T., Shnip O.A. 1997. Geologiya i neftegazonosnost fundamenta Zondskogo shelfa [Geology and oil and gas potential of the Zond shelf basement]. Neft i Gas, Moskwa, p. 288. (in Russian).
Gavrilov S.V. 2014. Issledovanie mekhanizma obrazovaniya ostrovnykh dug i zadugovogo razdvi-ganiya litosfery [Investigation of the mechanism of island arc formation and the back-arc spreading of the lithosphere]. Geofizicheskie Issledovaniya, 15(4): 35-43. (in Russian)
Gaynanov A.G., Panteleev V.L. 1991. Morskaya gravirazvedka [Sea gravity survey]. Nedra, Moskwa, p. 214. (in Russian)
Zharkov V.N. 2012. Fizika zemnykh nedr [Physics of the Earth's Interiors]. Nauka i Obrazovanie, Moskwa, p. 384. (in Russian)
Kogan M.G., Kukulieva R.Yu. 1988. Izgibnaya zhestkost litosferi Evrazii [Flexural rigidity of lithosphere of Eurasia]. Dokladi AN SSSR, 301(1): 6975. (in Russian)
Pavlenkova N.I., Pavlenkova G.A. 2014. Stroenie kory i verkhney mantii Severnoy Evrazii po dannym seysmicheskogo profilirovaniya s yadernymi vzryv-ami [The Crust and the Upper Mantle Structure of the Northern Eurasia as Constrained by the Data on Seismic Profiling with Nuclear Explosions], Ros-Nedra, IFZ RAN, Geokart-Press, Moskwa, p. 193. (in Russian)
Petrov A.B. 2001. Obrabotka dannykh gravi-metrii po geotraversu Uralseys s ispolzovaniem kompiyutornoy tekhnologii Koskod-3D [Processing of gravity measurements data on a geotraverse Ural-seys with use of computer technology Koskod-3D]. In: Glubinnoe stroenie i geodinamika yuzhnogo Urala. Izd. Gers, Tver, pp. 79-84. (in Russian)
Porfiriev V.V. 1982. Osobennosti glubinnogo stroenia zemnoy kori i teoreticheskie obosnovania neorganicheskogo genezisa nefti [Features of the deep structure of Earth's crust and theoretical justifications of inorganic genesis of the petroleum]. Kiev, Naukova Dumka, p. 328. (in Russian)
Seiful-Mulyukov R. 2012. Neft i gas: glubinnaya priroda i ee prikladnoe znachenie [Petroleum and gas: Inorganic abiotic nature and its applicability]. Torus Press, Moskwa, p. 216. (in Russian)
12
C.B. raepunoe, A.M. XapumoHoe
Stupak V.M., Berzin R.G. 2001. Strukturnye, kin-ematicheskie i dinamicheskie kharakteristiki Yu-zhnogo Urala po dannym differentsialnoy seysmora-zvedki [Structural, kinematical and dynamical characteristics of Southern Urals on the data of differential seismic]. In: Glubinnoe stroenie i geodinamika yujnogo Urala. Izd. Gers, Tver, pp. 43-51. (in Russian)
Barazangi M., Isacks B.L. 1976. Spatial Distribution of the Earthquakes and Subduction of the Nazca Plate beneath South America. Geology. 4: 686-692.
Karig D.E. 1971. Origin and development of marginal basins in the Western Pacific. Journal Geophysical Researches. 76(11): 2542-2561.
McKenzie D.P. 1969. Speculation on the consequences and causes of plate motions. Geophys. J. Roy. Astron. Soc., 18: 1-32.
Turcotte D.L., Schubert G. 2002. Geodynamics. Cambridge, Cambridge University Press, p. 448.
Schubert G., Turcotte D.L., Olson P. 2001. Mantle Convection in the Earth and Planets. New York, Cambridge University Press, p. 940.