ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
Геология
Вып. 4 (25)
ПЕТРОЛОГИЯ, ВУЛКАНОЛОГИЯ
УДК 552.321.6+553.46
Пластическое течение и реоморфическая дифференциация вещества в мантийных ультрамафитах
Д.Е. Савельев3, В.Б. Федосеевь
аИнститут геологии Уфимского научного центра РАН, 450077, Уфа, ул.К.Маркса, 16/2. E-mail: savl71@mail.ru
ьНижегородский государственный университет им. Н.И. Лобачевского, 603950, Нижний Новгород, пр. Гагарина, 23. E-mail: fedoseev@phys.unn.ru
(Статья поступила в редакцию 17 сентября 2014 г.)
Рассмотрены основные закономерности строения ультрамафитовых массивов офиолитовой ассоциации, которые принято считать фрагментами верхней мантии, выведенными на современный эрозионный срез в результате тектонических процессов. Показано, что главной тенденцией изменения состава и структуры мантийных пород является их стратификация, сопровождающаяся обособлением реологически наиболее «слабых» агрегатов - дунитовых тел, которые часто являются вмещающими породами для хромовых руд. На основе проведенного анализа сформулированы аргументы, не позволяющие принять существующие гипотезы петро- и рудогенеза в верхней мантии. Разработана оригинальная рео-морфическая модель дифференциации мантийного вещества, опирающаяся на термодинамический подход, позволяющая разрешить проблемы и противоречия, возникающие при использовании классических моделей (магматической и метасоматической).
Введение
В настоящее время считается общепринятым, что верхняя мантия Земли имеет ультраосновной состав [30 и др.], вещество ее находится в вязкопластическом состоянии и поэтому в ней могут возникать адвективные [10, 52] и/или конвективные течения [28, 29, 44, 46 и др.], построены многочисленные модели конвек-тирующей мантии (например, [15, 45 и др.]). Вместе с тем одно из следствий допущения о возможности возникновения вязких мантийных потоков остается обычно за рамками рассмотрения. Оно состоит в том, что течения возникают в среде, не являющейся истинным флюидом, а обладающей полифазным минеральным составом. Такие течения неизбежно должны вести к перестройке внутренней структуры и вещественным преобразованиям. Направленность структурных и вещественных изменений будет зависеть от конфигурации силовых полей, иниции© Савельев Д.Е., Федосеев В.Б., 2014
рующих потоки, и от физических свойств частиц, слагающих данную геологическую среду.
Одним из источников информации об особенностях состава и строения верхней мантии геологического прошлого являются ультрамафиты, обнажающиеся в составе офиолитовых ассоциаций складчатых поясов. В настоящей статье проведен анализ главных особенностей строения офи-олитовых ультрамафитов, сформулированы основные проблемы существующих генетических концепций, а также предложена оригинальная физическая модель реоморфической дифференциации вещества верхней мантии, позволяющая их разрешить.
Состав и строение ультрамафитов складчатых поясов
Прямую информацию о составе и строении верхней мантии Земли можно получить из трех основных источников: 1) керна глубоководного бурения и материала драгировок с дна Мирового океана, 2) ультраосновных ксенолитов из внутри-плитных базальтов континентов и океанических плит, 3) ультрамафитовых массивов, обнажающихся в составе офиолито-вых комплексов.
Многочисленные исследования представителей трех названных типов геологических объектов показали, что, несмотря на некоторые второстепенные особенности, их объединяет главное: ультраосновной (существенно перидотитовый) состав и постоянное присутствие высокотемпературных пластических деформаций главных породообразующих минералов [30, 60 и др.], на которые затем накладываются вторичные низкотемпературные преобразования пород (в основном сер-пентинизация). В настоящем исследовании мы ограничимся рассмотрением только высокотемпературной (мантийной) истории ультрамафитов и не будем учитывать наложенные низкотемпературные (коровые) деформации. Учитывая, что наиболее представительный материал о
составе и структуре верхней мантии дают офиолиты, в дальнейшем речь пойдет именно о них.
Первичные породы массивов являются ультраосновными тектонитами, претерпевшими высокотемпературное пластическое течение, и представляют собой тугоплавкий остаток (рестит), испытавший частичное плавление и экстракцию из него легкоплавкой базальтовой составляющей. Последняя часто присутствует либо во внутренних частях ультрамафитового разреза в виде даек габброидов, либо в краевых частях массивов формирует крупные дифференцированные тела пи-роксенит-габбрового состава. Наиболее интересны и дискуссионны особенности внутреннего строения реститового комплекса офиолитовых ультрамафитов. Выделим главные черты их строения:
Помимо указанных выше главных закономерностей строения ультрамафито-вых комплексов, особо следует отметить характер взаимоотношений между шпи-нелевыми перидотитами, составляющими основу разреза реститов, и дунитами, обычно имеющими подчиненное значение. В состав шпинелевых перидотитов помимо оливина (60-80% по объему) входит заметное количество пироксенов (ор-топироксен 15-25%, клинопироксен до 15%), тогда как дуниты являются практически мономинеральными оливиновыми породами. Анализ обширного фактического материала по офиолитовым комплексам мира показывает, что во всех случаях ду-ниты являются «более поздними», «эпигенетическими» образованиями по отношению к перидотитам [23, 31, 39, 54, 66 и др.]. Это выражается в пересечении дуни-товыми жилами полосчатых гарцбургитов, формировании сетчатого дунит-гарцбургитового комплекса, дунитовых штокверков, будинаже гарцбургитов, включенных в дунитовую матрицу [17]. По мнению авторов настоящей статьи, наблюдаемые явления объясняются не столько «полихронностью» происхождения пород, сколько различной их мобильностью в условиях мантийного пластического потока [35, 36].
Аналогичные взаимоотношения можно наблюдать на микроскопическом уровне между ортопироксеном (энстатитом) и оливином. Одним из первых на меньшую вязкость оливина по сравнению с энстати-том в процессе деформации мантийных ультрабазитов указал А.Николя с соавторами при описании альпийских массивов [71]. В цитируемой работе показано, что при одних и тех же условиях деформирования в зернах оливина происходит трансляционное скольжение в направлении смещения, а в зернах ортопироксена наблюдается разрыв. Подобные соотношения между породообразующими минералами ультрамафитов часто наблюдаются и на массивах Южного Урала [32, 36].
Многочисленные экспериментальные исследования также показали, что в широком диапазоне Р-Т-в условий вязкость агрегатов энстатита примерно на 25% выше, чем форстерита [57, 61], для оливина характерно явление гидролитического ослабления [64], а также более высокая скорость рекристаллизации отжига по сравнению с пироксенами [57].
Таким образом, рассмотренные выше главные особенности строения и состава офиолитовых ультрамафитов позволяют предположить, что основным механизмом их становления было твердофазное (рео-морфическое) течение в условиях верхней мантии. Прежде чем перейти к изложению физической модели реоморфической дифференциации мантийного вещества, рассмотрим существующие генетические представления и обозначим основные проблемы, связанные с предлагаемыми в них механизмами петрогенезиса.
Основные генетические концепции и их проблемы
Происхождение ультраосновных пород офиолитовых комплексов и связанных с ними месторождений хрома на протяжении долгих лет является предметом оживленной дискуссии. В первой половине XX столетия господствующей гипотезой образования хромититов являлась магматическая, в рамках которой они рассматривались как продукты дифференциации ультраосновной магмы (дунитовой или пери-дотитовой) [6, 17, 20, 22, 26, 42, 43 и др.]. Развитие метасоматической гипотезы [4, 23, 31 и др.], рассматривающей дунит-гарцбургитовые комплексы с хромовым оруденением как продукт метасоматиче-ских преобразований перидотитового (эн-статититового) субстрата, было вызвано главным образом невозможностью объяснить магматической дифференциацией факт «эпигенетичности» дунитов и хро-мититов (более тугоплавких) по отношению к окружающим гарцбургитам.
К разновидностям метасоматической гипотезы можно отнести и реакционно-магматическую модель образования ульт-рамафитов и хромовых руд в офиолитах, которая в последние годы имеет колоссальное влияние на специалистов, изучающих ультрамафиты во всем мире [5, 38, 65, 66, 78 и др.]. Суть гипотезы сводится к следующему. Гарцбургиты и ду-ниты офиолитовых комплексов являются не простыми тугоплавкими остатками от частичного плавления пиролита, а представляют собой результат реакции рестита с проникающими сквозь него расплавами, генерированными на большей глубине. Расплавы «расторгают» пироксены и «кристаллизуют» дополнительный оливин, что в предельном случае ведет к формированию дунитов. В некоторых интерпретациях дуниты считаются изолированными «каналами», по которым базальтовые расплавы транспортируются к поверхности [67]. При этом образование хромититов связывается исключительно с транспортировкой сквозь перидотиты бо-нинитовых расплавов, реакцией бонини-тов с перидотитами, расторжением пирок-сенов, осаждением оливина и высокохромистого шпинелида [78 и др.].
Против реакционной (метасоматиче-ской) интерпретации петрогенезиса в верхней мантии можно выдвинуть ряд аргументов.
Аргумент 1. В метасоматических гипотезах, как правило, отсутствует описание механизмов растворения - переноса -осаждения мантийного материала, в первую очередь весьма инертного хрома. За исключением реакционно-магматической все они апеллируют к гипотетическим агентам массопереноса, не указывая их конкретного состава.
Аргумент 2. В гарцбургитах, дунитах и хромититах отсутствуют вещественные доказательства образования их при активном участии флюидной фазы. Высокотемпературные «флюидные» минералы (амфиболы, флогопит) здесь обычно отсутствуют.
Аргумент 3. Отсутствуют реакционные взаимоотношения между породообразующими минералами ультрамафитов, например, между оливином и ортопироксе-ном. При метасоматическом образовании оливина на месте пироксенов в процессе замещения под воздействием флюидов такие соотношения должны наблюдаться повсеместно (везде, где сохранился энста-тит), иначе следует признать, что «оливи-низация» - мгновенный процесс, а это противоречит здравому смыслу.
Аргумент 4. Для метасоматических гипотез мантийные ультрабазиты представляют собой статичную субстанцию, сквозь которую проникают расплавы или флюиды и преобразуют ее химическим путем. Максимально, что «дозволяется» кристаллическому каркасу, так это испытывать незначительную деформацию, в результате которой из них удаляется частичный расплав. В то же время без обсуждения предполагается, что все процессы плавления, просачивания и реакции расплава с остаточными минералами происходят внутри поднимающегося мантийного вещества (mantle diapir, upwelling mantle). Механизм этого подъема и его влияние на структуру и вещественный состав мантийного субстрата обычно не рассматриваются.
Аргумент 5. В рамках метасоматиче-ских гипотез, в том числе реакционно-магматической, нельзя объяснить резкие контакты дунитов и гарцбургитов, а также увеличение содержания ортопироксена в
приконтактовых частях гарцбургитов, что очень часто наблюдается в офиолитовых массивах [35].
Аргумент 6. Если дунит - канал, то в нем должно быть достаточно много расплава вплоть до момента «замерзания» системы. Если в момент ее «замерзания» в канале много расплава, то это уже не ду-нит, т. к. при кристаллизации расплава должно происходить образование базаль-тоидной ассоциации (пироксены, плагиоклаз) и порода по составу будет уже соответствовать лерцолиту. Если же в «канале» при его «замерзании» расплав отсутствует, то может ли он считаться каналом, где протекал расплав? В реальных офиолитах, которые послужили полигоном для разработки реакционно-магматических моделей, наблюдаются именно дуниты.
Аргумент 7. В рамках реакционно-магматической гипотезы отсутствует механизм образования хромититов как геологических тел. Даже если принять точку зрения, что хромит кристаллизуется в ду-ните при «расторжении» пироксенов, остается непонятным, какие факторы заставляют его формировать отдельные тела с различающейся концентрацией - от ред-ковкрапленных до массивных руд.
Таким образом, традиционные гипотезы образования ультрамафитов складчатых поясов - магматическая и метасома-тическая (включая реакционно-магматическую) - сталкиваются с серьезными проблемами при объяснении генезиса хромовых руд в офиолитах. Альтернативой для них является признание ведущим механизмом дифференциации вещества реоморфическое (твердофазное) течение ультрамафитов. Со времени первых исследований пластического течения ультрамафитов, инициированных Григг-сом, достигнуты большие успехи как в понимании механизмов пластической деформации главных породообразующих минералов верхней мантии, так и процессов формирования офиолитовых комплексов [9, 37, 51, 54, 56, 57,58, 61, 71, 72, 63, 64, 68 и др.]. Вместе с тем, признавая
важную роль пластических деформаций в образовании структуры ультрамафитовых тел, исследователи в большинстве случаев обходят стороной вопрос о связи тектонического течения и дифференциации вещества верхней мантии, и лишь в отдельных работах делается вывод о ведущей роли реоморфизма в процессе мантийного пет-ро- и рудогенеза [3, 11, 34, 35, 36, 50, 54]. В следующем разделе мы представляем физическую модель реоморфической дифференциации мантийного материала, в основе которой лежит термодинамический подход.
Физическая модель реоморфической дифференциации
Рассматриваемые нами в этом разделе геологические структуры не являются изолированными системами. Они обмениваются с внешней средой энергией в виде упругих (сейсмических) волн, теплопере-носа, либо создают вокруг себя поля упругих напряжений, имеют сложный неоднородный химический и фазовый состав и т.п. Поэтому далее для описания подобных систем используется термодинамический подход.
Как известно, в геологических объектах широким распространением пользуется вертикальная стратификация, возникающая в условиях гравитационного поля и образующаяся при осаждении более плотных фаз (кристаллов и обломочных зерен) в менее плотной жидкости (расплаве или воде). Напротив, главной особенностью рассматриваемых нами ультрама-фитовых комплексов является неоднородность распределения компонентов по горизонтальному сечению. Поэтому необходимо выделить внешние факторы, способные создать горизонтальные градиенты распределения компонентов. Кроме гравитационного поля существует множество иных факторов, вызывающих изменения в составе и распределении компонентов в пространстве. К ним можно отнести неоднородность температуры и/или концентрации в слоях, граничащих с рассматриваемой системой, внешние источники акустического (вибрационного и сейсмического) воздействия, гидростатические поля (тектонические напряжения), вязкопластические потоки. С каждым из этих факторов можно связать возникновение соответствующих обобщенных сил, проявляющихся в виде разнообразных эффектов (Бенара, Киркендалла, Френкеля), сопровождающихся массопереносом и перераспределением компонентов в горизонтальном направлении. Ограничим рассмотрение только гидродинамическим полем, полагая, что вибрации и температурные градиенты отсутствуют или несущественны.
Гидродинамические поля возникают в породах восходящего мантийного диапи-ра. Эти поля возникают при вязкопласти-ческих течениях в твердофазных системах. Сами течения могут быть вызваны разными причинами: компрессией-декомпрессией в окружающих породах, разностью температуры в окружающих систему слоях (термоконвекция), конвекционные течения могут возникать и под действием гравитационного поля [2].
Скорости массопереноса и перераспределения компонентов в условиях вяз-копластических потоков могут существенно превышать скорость диффузионных процессов. Поэтому существование вяз-копластических течений любого происхождения можно рассматривать в числе факторов, участвующих в формировании структуры верхней мантии.
Экспериментальные наблюдения потоков в двухфазных системах разного типа показали, что частицы в них могут смещаться перпендикулярно оси потока в различных направлениях. Экспериментально установлено, что в вертикальном нисходящем или горизонтальном потоке частицы, имеющие плотность более высокую, чем плотность среды, смещаются к стенке, а частицы с меньшей плотностью - к оси потока. Концентрация дисперсной фазы может быть неоднородной и несимметричной, жидкие частицы могут образовывать пленки на стенках, пузырьки
всплывают преимущественно по центральной области вертикальных потоков, крупные частицы могут двигаться, сталкиваясь со стенками, ограничивающими поток. Частицы с нулевой плавучестью двигаются, концентрируясь на удалении и от стенок, и от оси потока, либо осциллируют внутри потока. На движение частиц оказывает влияние множество факторов, в частности, геометрия, скорость, плотность частиц и среды [7, 47].
Многочисленные обзорные работы, экспериментальные и численные исследования последних лет свидетельствуют, что проблема теоретического описания миграции частиц в потоках самого разного типа остаётся актуальной и в настоящее время [8, 13, 69, 73]. Классическое описание основано на аппарате гидродинамики [18, 25, 47], существуют работы, в которых привлекаются методы и понятия статистической физики и термодинамики [12, 13, 14, 41, 48].
Рассмотрим процессы массопереноса в стационарном гидродинамическом поле. Поместим тело массой ть в среду с плотностью р, которая движется равномерно и прямолинейно. При этом сечение потока, средняя скорость и эпюра скоростей остаются постоянными по всей длине потока. Простейшими примерами таких потоков являются течение Куэтта - ламинарное течение вязкой жидкости между двумя параллельными стенками/плитами, одна из которых двигается относительно другой, и течение Пуазейля - ламинарное течение вязкой жидкости через каналы в виде прямого кругового цилиндра или слоя между параллельными стенками/плитами под действием перепада давления.
Свяжем систему координат с неподвижной границей потока так, что ось х совпадет с направлением потока, а ось г будет нормальна границам.
Рассмотрим случай, при котором скорость движения центра тела, находящегося на расстоянии г от неподвижной стенки, равна скорости течения v(r) среды на том же уровне. Пусть за время наблюдения тело не успевает смещаться поперёк
потока и не вращается под действием обтекающих его потоков.
Система является стационарной, т. к. полная и внутренняя энергия, импульс, количество вещества, состав, масса, объемы жидкой и твёрдой фаз не зависят от времени. Полная энергия системы "тело + среда" имеет вид
МУ02 2
+и.
Б (т ) = - dW (т ) стремится изменить по-йт
ложение тела и определяет направление самопроизвольного перемещения компонентов системы (тела и среды), при которых полная энергия системы понижается.
В простейшем случае плоского сдвигового течения Куэтта v(z) = az действует сила
где масса системы М = ^ mi складывается из масс т\\, составляющих систему частей. У0 - скорость движения центра инерции системы в выбранной системе координат, т.е. относительно неподвижной стенки. Внутренняя энергия и включает кинетическую энергию относительного (по отношению к центру инерции системы) движения всех частей системы и потенциальную энергию их взаимодействия
Опишем состояние системы в зависимости от положения тела в потоке. Потенциальную энергию взаимодействия тела со средой примем равной нулю. Тогда внутренняя энергия системы включает только кинетическую энергию движения среды и тела относительно центра инерции системы, движущегося со скоростью У0. Пренебрегая локальными потоками среды в окрестности тела, связанными с тем, что один край двигается быстрее потока, а другой медленнее, представим полную энергию системы как
Н Т2 ^2
W(т)= Ар01V2 (г)сЪ + — | (рьУ (т)-ру? (г))сЪ, (2)
здесь А - параметр, по смыслу соответствующий площади горизонтального сечения системы. Первый интеграл соответствует энергии среды при отсутствии тела. Второе слагаемое равно энергии тела и учитывает энергию среды, которую заменило тело. Предполагается, что тело с плотностью рь имеет форму куба со стороной —, основание куба ориентировано параллельно стенкам потока.
Энергия системы (2) зависит от положения тела т Обобщенная сила
Б ( т ) =
а 2ттп
где а = — - скорость сдвига, ть = — рь, йг
т0 = — р0 - масса тела и масса среды, вытесненной телом, рь, р0 - плотность тела и среды соответственно. При нулевой плавучести тела (рь = р0) выталкивающая сила отсутствует. Для "легких" тел (рь < р0) сила направлена в сторону увеличения скорости потока. Для "тяжелых" (рь > р0) тел сила направлена в сторону уменьшения скорости потока, к его неподвижной границе.
Существенное отличие силы (3) от силы Архимеда состоит в том, что ¥(т) зависит от положения тела относительно границ системы. Для течения Куэтта это кажется странным, ведь при любом т разность скоростей в окрестности тела одинакова. Однако, вспоминая, что сила действует и на тело, и на среду, можно представить это следующим образом: слои, расположенные между неподвижной стенкой и телом, выполняют роль рычага. Как и сила Архимеда, это пример обобщенных или "термодинамических сил", действующих на систему в целом, а не на отдельные составляющие систему тела. Обобщенной координатой является координата центра тяжести тела т. То, что формально сила привязана к положению тела, создаёт иллюзию, что сила действует именно на тело.
Рассмотрим более сложный случай, когда тело движется относительно среды со скоростью = у(т) + Ду. Гидродинамика не предлагает простого и однозначного уравнения, описывающего выталкивающую силу в этом случае, поэтому воспользуемся уравнением (2), принимающим вид
r " & & )+Av) —I 2
^ H ^ /+L2 Р Ж+A)2 —pV (z) (4) Wr)=Ap\\V2 (z) dz+L j ^^-&—-dz+W¡x (4)
При больших отклонениях скорости тела от скорости среды пренебрежение вкладом локальных потоков, возникающих в окрестности тела ДЖ1ос, становится некорректным. Чем сильнее скорость тела отличается от скорости среды, тем больше вклад ДЖ1ос. Выталкивающая сила, действующая на движущееся в потоке Куэтта тело, имеет вид
F (r, Av) = aLm0)r
—cm^AvAWloc. (5) dr
В геологических системах при вязко-пластическом течении относительная скорость движения тел Дv в твёрдой среде не велика, и последним слагаемым мы пренебрежём. Сгруппируем оставшиеся слагаемые следующим образом:
F ( r, Av ) =
p0L3Ava (6)
силы) [70]. Есть другие выражения, например, FM = 1.96 ( b1Re p + b2) p0d3Ava
[24] или Fy = 0.052p0d3 Ava, описывающее силу в условиях конкретного эксперимента [77]. В (6), как во всех приведённых уравнениях, сила пропорциональна массе жидкости, вытесненной телом, относительной скорости тела, скорости сдвига, а направление силы может быть разным в зависимости от плотности тел.
Рассмотрим некоторые частные случаи. Для «легких» тел р0 >> pb при малой скорости относительного движения тел lAvl « d\\a\\ < r\\a\\ сила (6) приближенно
равна
F ( r, Av ):
f^L + Рь.л
VAv Ро J
Она направлена в сторону увеличения г и стремится сместить тела, опережающие поток Дv > 0, в направлении увеличения скорости потока. Если же тело отстаёт от потока Дv < 0 (например, пузыри в «падающем» потоке), то сила меняет знак при
и сравним эту силу с ранее описанными силами [25]. Выражение (6) не зависит от предыстории, поэтому не имеет отношения к силе Бассэ-Бусинеска, возникающей при ускоренном движении тел. В рассматриваемой модели тело не вращается, поэтому силы Магнуса и Саффмена, зависящие от угловой скорости, отсутствуют. Сила (6) никаким образом не зависит от вязкости среды и не является силой трения.
Описывая движение не вращающегося тела в сдвиговом потоке, часто применяют другую силу Саффмана [53, 74]. В используемых обозначениях она имеет вид
и направлена в сторону роста
Эта сила полуV ^ J
чена для неинерционной системы отсчета (привязана к движущемуся телу), что осложняет сопоставление с вариантом (6). Коэффициент Сц сложным образом зависит от скорости потока и свойств среды [55]. При изменении условий течения Ся может менять знак (направление действия
скорости потока при r > r0 и в сторону неподвижной стенки при r < r0.
Для тел большой плотности р0 << pb выражение (6) имеет вид
F (r, Av)«fl - — ] Р р0 13Ava .
^ Av ) Ро
Сила стремится сместить тяжёлые тела, опережающие поток Av > 0 (например, тяжелые тела в восходящем потоке), на Av
уровень r = —.
Различные варианты поведения тел наглядно иллюстрирует зависимость энергии системы от положения и скорости тел в потоке Куэтта (рис. 1). Система, стремясь понизить свою энергию, переносит тела в сторону минимумов энергии.
Для плоского течения Пуазейля v(z) = b(zH - z2), 0 < z < H, где b - параметр, определяющий вид эпюры скоростей, H - ширина потока между неподвижными стенками, пропуская
Рис. 1. Зависимость энергии системы от положения тел г в потоке Куэтта для тел, опережающих поток (а) и отстающих от потока (б). Рядом с кривыми указано отношение плотности тел к плотности среды ръ/р0. Относительная плотность пузырей принята равной 0, тел с нулевой плавучестью - 1. Стрелками указано направление силы, действующей на тела разной плотности
громоздкую аналитическую запись полной энергии системы, получим
К (г, Ь)=щЪ
Сила формально разбивается на составляющие, одна из которых зависит от плотности компонентов системы, обращаясь в 0 для тел с нулевой плавучестью (ръ = ро), другая от плотности не зависит и действует даже на тела с нулевой плавучестью. В потоке Пуазейля выталкивающая сила зависит и от положения тела относительно стенки г, и от размеров тела Ё. Вводя безразмерные координаты, можно получить выражение
К(д)= л
(V-Зд2 +д)+Х/ °
где д = ЫН, 0 < д < 1 для пространственной координаты и Я = Ё/Н, 0 < X < 1 для размеров тела, А = р0ЬНъЪ2 - коэффициент, учитывающий геометрические характеристики.
Выталкивающая сила (9) антисимметрична относительно середины потока. В центре потока направление силы меняется
на противоположное. Полином (2д -Зд + д) имеет экстремумы в точках
q = У2 ± л/3/6 (0.211, 0.789), симметричных относительно середины потока. Они соответствуют слоям, где вклад первого слагаемого в выталкивающую силу максимален.
Вклад второго слагаемого пропорционален размерам тела и линейно растёт при удалении от центра потока. Независимо от плотности тела эта составляющая выталкивающей силы направлена к стенкам. Для очень мелких частиц (X ^ 0) вторым слагаемым можно пренебречь.
В центре потока (д = 0.5) выталкивающая сила равна нулю. Это соответствует положению равновесия, которое в зависимости от размеров и плотности тел может оказаться устойчивым или неустойчивым (рис. 2).
Для тел большой плотности (ръ > р сила всюду направлена к ближайшей стенке (рис. 2), но максимальное по величине значение она имеет на некотором удалении от нее. На тела с нулевой плавучестью (ръ = ро) действует сила
К (д ) = А
которая тоже всюду
направлена от центра потока к стенкам.
Для тел любой плотности в окрестности неподвижных стенок (д = 0 и д = 1) выталкивающая сила не равна нулю, прижимая тела к неподвижной стенке,
Рис. 2. Зависимость энергии системы от положения тела г в потоке Пуазейля: а - скорости тела и потока равны; б - тело отстаёт от потока. Цифрами указано отношение плотности тела к плотности среды рь/р0. Относительная плотность каверн или пузырей принята равной 0, тел с нулевой плавучестью - 1. Стрелками показано направление силы, действующей на тела разной плотности
она становится «прижимающем» силои.
Для тел, плотность которых меньше плотности среды, в центральной части потока сила направлена к центру, а на периферии - к стенкам (рис. 2). Положения
около q - У2 ± >/3/6 (я < 0.211 и q > 0.789), где выталкивающая сила обращается в 0 и меняет направление, зависят от размера и плотности тела. Они соответствуют максимумам энергии (рис. 3) и являются неустойчивым равновесием.
Согласно экспериментальным наблюдениям, в двухфазных потоках частицы могут концентрироваться в положении Я « 0.1^0.2 и я« 0.8^0.9, причём центральная часть потока содержит сравнительно мало частиц [47]. Используя термодинамический подход, можно обнаружить, что эта ситуация может возникнуть в потоках Пуазейля, когда скорость "тяжёлых" тел меньше скорости потока Ау < 0.
Скопления легких компонентов у стенок (см. рис. 2,б) могут образовать слои, двигающиеся в обратном основному потоку направлении. В потоках возможны и другие слоистые структуры (рис. 1-2).
Рассмотренная выше модель, в которой единственная частица двигается прямолинейно и равномерно вместе с потоком среды, практически не выходит за рамки механики. При описании движения множества тел (дисперсной фазы) в потоке следует использовать в качестве термодинамического потенциала не внутреннюю энергию системы и, а свободную энергию Гиббса G =и + PV - TS или Гельмгольца G =и - TS, где Р - давление, V - объём, Т - температура, S - энтропия. Энтропийный вклад позволяет рассмотреть пространственное распределение тел в потоке. При этом появляется другая обобщенная сила, связанная с разностью концентрации дисперсной фазы.
Пусть рассматриваемая система имеет возможность изменять свою энергию, суммарное количество и состав компонентов не меняются во времени (закрытая термодинамическая система). Самопроизвольные процессы в механике происходят в направлении уменьшения энергии системы. Обмен энергией с окружающими геологическими слоями может происходить как на стенках, так и вдоль потока. Выделяющаяся при этом энергия может рассеиваться в виде гидродинамических, акустических или тепловых возмущений, излучаемых системой. В реальных условиях подобные взаимодействия с окружающей средой могут сопровождаться также изменением эпюр скоростей потока,
распределения гидростатического давления, рождением вихревых движений среды в окрестности переносимых потоком тел, перемещениями и поворотами самих тел.
В термодинамике самопроизвольные процессы происходят в направлении перехода в равновесное состояние. Критериями равновесного состояния являются экстремумы термодинамического потенциала (минимум свободной энергии или максимум энтропии системы и др.). Во внешнем поле одним из критериев термодинамического равновесия является условие инвариантности суммы химического и механического потенциала. Следуя работам [1, 49], назовём эту сумму механохимиче-ским потенциалом.
Концентрация дисперсных частиц в системе п(д) определяется их индивидуальными свойствами (масса, плотность, модуль упругости) и геометрической конфигурацией поля (однородное, центрально-симметричное, распределенное и др.). Если система содержит несколько компонентов, условием термодинамического равновесия является инвариантность ме-ханохимического потенциала для каждого из компонентов системы:
Ц (д)+ц (д)+р( д)^( д)+Ц П (д))=«та*, (10)
здесь д - координаты, р(д) - гидростатическое давление, У(д) - удельный (на моль или на одну частицу) объём компонента, для несжимаемых компонентов это константа, п/(д) - мольная доля компонента, ц° - стандартный химический потенциал компонента. Слагаемое 01п(пг) - это энтропийный вклад (энтропия смешения).
Для систем с вязкопластическим течением в однородном гравитационном поле условие инвариантности механохимиче-ского потенциала принимает вид
V (X, У, 2)
-У&* >
п, = к+pv
- доля частиц типа /, р^) - гидростатическое давление, у(х,у,£) - скорость потока, координата 2 направлена вверх. Для одно-компонентной несжимаемой жидкости П = 1 это выражение становится аналогом уравнения Бернулли.
Выражения типа (11) описывают гипсометрические или пространственные распределения дисперсных компонентов. Для объемного распределения в гидродинамическом и гравитационном поле из (11) следует
П ( д ) = Сг еХР
gz +
+р (q )V
Для слоя, расположенного на фиксированной высоте z = const с прямолинейным горизонтальным ламинарным потоком вдоль оси у, распределение упрощается до
v(х)2 / \\ + р( x )v
n (x) = Ce в . (12) Различия в концентрации частиц определены эпюрой скорости v(x) и свойствами компонентов. Для крупнодисперсных систем это выражение имеет смысл использовать только в окрестности минимумов внутренней энергии U (рис. 1-2) либо в виде отношения концентрации частиц разного типа. Например, для смеси частиц А и В в дисперсионной среде плотностью р0 на основе уравнения (12) получим
ПаМ - С ехр
- CAB еХР
-(рЛ -PBVB +V -vb)Po)■
здесь константа Саб определяется исходным составом системы, Уа,УВ - объём, а Ра, рВ - плотность частиц типа А и В соответственно. Это уравнение показывает тенденцию изменения состава дисперсной системы в разных частях потока.
При ( Ра + р0 Уа <( рв + р0 )УВ,
va ^ pb + po
или <
vb Ра + Po
+ Р (X, у, z) V + в ln (n (x, у, z)) - const здесь g - ускорение свободного падения, р,, V, - плотность и удельный (на моль или на одну частицу) объём компонентов, nl(h)
доля частиц А будет возрастать там, где скорость потока больше. Здесь возможны варианты:
- частицы одного размера
(рА + р0)/(рв + р0) > 1, что эквивалентно
рА > рВ, - менее плотные там, где скорость выше;
- частицы равной плотности
УА/Ув < 1 - мелкие частицы преобладают там, где скорость выше. В этом случае концентрация и крупных, и мелких частиц малой плотности по-прежнему уменьшается, а частиц большой плотности возрастает к периферии потока;
- частицы мелкие, но плотные (или крупные, но легкие) УЛ ( Рл + р0 )< Ув (Рв + р0 ) могут повести
себя неоднозначно, например, возможен случай, когда мелкие частицы большой плотности начнут преобладать в области больших скоростей над крупными частицами малой плотности.
Изложенную выше физическую модель вполне можно применить к рассматриваемым нами мантийным породам. Как было отмечено выше, основу офиолито-вых ультрамафитов составляют перидотиты, сложенные на 60-80% оливином, на 20-40% пироксенами (энстатитом (2035%) и диопсидом (до 10%)). Геолого-структурные и петрографические данные свидетельствуют о том, что наиболее мобильными являются оливиновые агрегаты, в предельном случае формирующие «слои» дунитов.
Из петрографических наблюдений мантийных ультрамафитов следует, что размеры структурных элементов деформации для сосуществующих оливина, ор-топироксена и хромита не одинаковы. Пироксен в пластическом течении участвует в виде зерен или их крупных фрагментов (1-5 мм, в среднем 1,5-2 мм; р=3,2 г/см3), оливин - в виде субзерен со средним размером около 0,5 мм (р=3,2 г/см3). Подставим в полученное неравенство значения плотности и размеров частиц (структурных элементов деформации) минералов при условии, что дисперсионной средой является оливин. Для пары оливин (01) -ортопироксен (Орх) получим
В перидотитах помимо главных минералов (Ol, Opx, Cpx) присутствует глиноземистый хромшпинелид в количестве 0,n %, образующий обычно ксеноморфные выделения. По данным экспериментального плавления шпинелевых перидотитов [62, 75], этот минерал одним из первых (вместе с пироксенами) переходит в расплав. Вместе с тем в дунитах уже встречается и часто образует скопления вплоть до месторождений высокохромистый идио-морфный хромшпинелид (хромит). Таким образом, дуниты и хромититы сложены наиболее тугоплавкой минеральной ассоциацией, но всегда обнаруживают более высокую мобильность по сравнению с перидотитами: внедряются в последние, формируют штокверки. Кроме того, в строении месторождений хрома отмечается еще одна важная закономерность: мощные дунитовые тела обычно включают рассеянное оруденение с мелкозернистым хромитом, а тела небольшой мощности среди перидотитов вмещают массивные хромовые руды, сложенные средне- и крупнозернистым хромитом.
Все указанные выше закономерности объяснимы лишь с позиций реоморфиче-ского генезиса мантийных ультрамафито-вых комплексов. Рассмотрим этот вопрос подробнее. Акцессорные хромиты дуни-тов и бедновкрапленных руд мелкие (0,1-0,8 мм при среднем 0,5 мм; р=4,5 г/см3), в более концентрированных рудных телах размер зерен увеличивается до 0,5-3 мм, а наиболее крупнозернистыми являются массивные хромиты (2-10 мм). Тогда частицы хромита (Chrt) согласно неравенству 1 =( Voi/Vchrt)<((pChrt+Pol)/(OOl +pOl))~ 1.2 должны преобладать в областях с низкими скоростями течения. С уменьшением относительного размера частиц хромита, участвующих в деформации, область распространения хромита может увеличиться.
Прогрессирующее пластическое течение уже расслоенных ультрамафитов может реализоваться в обстановке сжатия, что ведет к образованию на месте полосчатых дунит-перидотитовых комплексов
сложных по строению сетчатых и шлиро-во-полосчатых образований того же состава. Деформация локализуется в более слабых (мобильных) дунитовых слоях, приводя к их растеканию, а перидотиты (обычно гарцбургиты) играют роль относительно жестких ограничений как для дунитов, так и для хромититовых тел, которые мигрируют внутри пластичной ду-нитовой оболочки.
Поскольку дунит и гарцбургит имеют постепенные переходы и различие в их реологии определяется количественными соотношениями между оливином и орто-пироксеном, в ходе деформации происходит дальнейшее перераспределение зерен данных минералов: оливин стремится к формированию мономинеральных подвижных масс, ортопироксен также образует зоны обогащения вплоть до жил энста-тититов.
В итоге происходит будинаж компетентных гарцбургитов, наиболее крупные ослабленные зоны образуются по механизму трещин отрыва и заполняются реологически слабыми дунитами, они осложнены серией более тонких жил, образованных по механизму сколовых трещин. Формирующаяся таким образом структура соответствует картируемым на массивах гарцбургитового типа шлирово-полосчатым дунит-гарцбургитовым комплексам, в строении которых принимают участие крупные дунитовые тела, обрамленные густым штокверком маломощных дунитовых жил. Крупные обособления дунитов располагаются внутри зон интенсивного пластического течения с определенным шагом, который составляет около 10 км [33].
В хромовых рудах одновременно с ростом содержания рудных минералов увеличивается размер их зерен, а тела, сложенные густовкрапленными и массивными рудами, имеют тенденцию локализоваться вблизи компетентных гарцбургитовых экранов [17, 23]. Последнее объясняется тем, что согласно неравенству (13) при увеличении размеров структурных элементов деформации (зерен) хромита до 2
мм и сохранении таковых для оливина (0,5 мм), мы получаем 0.03125 = (Уо/Ус^д <(рап+Рог) /(рог +Ро/))~1.2, что говорит о сильном разделении в пластическом потоке данных минералов. Следствием этого в условиях сжатия является выжимание более мобильных оливиновых агрегатов в области с меньшим давлением, а залежи массивных хромититов ассоциируют с гарцбургитами, отделяясь от последних маломощной дунитовой оторочкой.
Таким образом, реоморфическая модель позволяет успешно объяснить все главные закономерности строения ультра-мафитовых массивов офиолитовой ассоциации. Надеемся, что дальнейшая разработка представленной модели может решить и более локальные проблемы, связанные с составом и строением рассматриваемых геологических объектов. К последним, в частности, относятся: 1) вариации состава и морфологии хромшпи-нелидов при переходе от лерцолитов к ду-нитам и хромовым рудам; 2) различный состав рудообразующих хромшпинелидов в различных массивах и их вариации в пределах крупных тел, а также ряд других.
Заключение
В работе обоснована принципиальная модель, способная объяснить закономерности перераспределения вещества верхней мантии в процессе становления офио-литовых комплексов и формирования рудных тел. Разработанная модель разрешает проблемы и противоречия, возникающие при допущении магматического или мета-соматического способа образования рассматриваемых геологических объектов.
В текстурах и структурах ультрамафи-тов офиолитовой ассоциации запечатлены признаки пластического течения, которое являлось одним из главных факторов пет-рогенеза и рудообразования. Стратификация мантийного вещества, сопровождавшаяся, в частности, формированием тел хромититов, происходила при твердофазном перераспределении минералов в породах, представляющих собой дисперсионную систему. Первопричиной пластического течения и структурирования мог быть декомпрессионный подъем мантийного вещества в пределах рифтогенных структур и сейсмическое воздействие. Дальнейшее усложнение внутренней структуры и более глубокое разделение вещества, вероятно, происходило в обстановке сжатия, что, возможно, сопоставимо с современными геодинамическими об-становками в мантии преддуговых бассейнов.
Особенности деформационного процесса могли определяться воздействием различных физических полей: акустического (вибрационного, сейсмического), гравитационного, гидродинамического и гидростатического. Вероятнее всего, в природных объектах мы имеем дело с результатом комплексного воздействия (суперпозиции) данных факторов. При всех вариантах перечисленных воздействий определяющим параметром является размер частиц, для гравитационного, акустического и гидродинамического полей существенным фактором является плотность, а для гидростатического поля существенным фактором становятся и модули упругости компонентов породы. Подстановка в условия (12, 13) числовых значений физических свойств компонентов мантийных ультрабазитов позволяет получить решения, подтверждающие постоянную ассоциацию в рудных телах хромита с оливином. На основе предложенной модели возможно рассмотрение обратной задачи о том, какие условия могли привести к формированию рудных тел различного размера, морфологии и качественного состава хромитовых руд.
Работа выполнена при поддержке гранта «РФФИ-Поволжье», проект «№14-05-97001 р_повольжье_а»
Библиографический список
Российский химический журнал. 1998. № 42. С. 36-46.
формирования хромитового оруденения в альпинотипных гипербазитах // Рудогенез. Миасс-Екатеринбург, 2008. С. 4-7.
хромитов Кемпирсайского массива // Магматизм, метаморфизм, металлогения Урала. Свердловск, 1963. С. 325-330.
расплавов в мантии под зонами спрединга и образование дунитов замещения: обзор проблемы // Геология и геофизика, 2009, Т.50, №9. С. 992-1012.
частицами. М.: Мир, 1975. 380 с.
твердыми частицами. М.: Физматлит, 2003. 192 с.
конвекция в тектоносфере // Геотектоника. 1990. № 5. С. 3-8.
B.И., Романовская М.А. К проблеме генезиса хромитового оруденения в гиперба-зитах альпинотипной формации // Вестник Московского университета. 2006. №6.
B.Г., Суслов П.В., Брюгманн Г. Структура каналов течения расплавов в мантии // Геотектоника. 2008. № 6, С. 25-45.
C.16-27.
A Plastic Flow and Rheomorfic Differentiation of the Mantle Ultramafic Rocks
D.E. Savelieva, V.B. Fedoseevb
a Institute of Geology, Ufa Scientific Center of Russian Academy of Sciences, 16/2, K. Marksa Str., Ufa 450077, Russia E-mail: savl71@mail.ru b Nizhniy Novgorod State National Research University named after Lo-bachevskiy, 23 Gagarin Pr., Nizhniy Novgorod, 603950, Russia E-mail: fedoseev@phys.unn.ru
In this paper, the general characteristics of morphological features of the ophiolitic ultramafic rock formations are discussed. The ultramafic rocks are the fragments of upper mantle, which were exposed on the surface due to tectonic events. It is shown that their main chemical and structural characteristic is a stratification accompanied by separation of the rheologically weakest dunite bodies usually containing the economic amount of chromite ore. Based on results of conducted analysis, we propose a new hypothesis of petro- and ore genesis in the upper mantle. Using the thermodynamic approach, we developed the rheomorfic model of the differentiation of the mantle matter. This model solves many problems inherent to currently used magmatic or me-tasomatic models.
References
perekhody i khimicheskoe ravnovesie v tsen-trobezhnykh polyakh. Gravitatsionnaya khimiya [Phase transitions and chemical equilibrium in the centrifugal fields. Gravitational Chemistry]. Rossiyskiy khimicheskiy zhurnal. 42: 36-46.
formy sosuda i samoproizvolnoe vozniknovenie tsirkulyatsii pri vrashchenii mnogokomponentnykh zhidkikh smesey [Effect of vessel shape and spontaneous appearance of circulation due to spinning of multi-component liquid mixtures]. Vestnik KGTU. pp. 101-104.
formation in alpinotype ultrabasic rocks]. In Rudogenez. Miass-Ekaterinburg, pp. 4-7.
I khromitov Kempirsayskogo massiva [About the origin of dunites of the Kempirsay massif]. Magmatizm, metamorfizm, metal-logeniya Urala. Sverdlovsk, pp. 325—330.
gaza s tverdymi chastitsami [Turbulent flowing of gas with solid particles]. Moskva, Fizmatlit, p. 192.