Физико-математическое моделирование сил динамического давления и вязких напряжений в зонах литосферной субдукции
- Авторы: Гаврилов С.В.1, Харитонов А.Л.2
-
Учреждения:
- Институт физики Земли им. О.Ю. Шмидта РАН»
- Институт земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова РАН»
- Выпуск: № 4 (2023)
- Страницы: 70-77
- Раздел: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ
- Статья опубликована: 01.12.2023
- URL: https://vestnikugrasu.org/byusu/article/view/635194
- DOI: https://doi.org/10.18822/byusu20230470-77
- ID: 635194
Цитировать
Полный текст
Аннотация
Предмет исследования: зоны активных тектонических процессов в мантии, влияющих на литосферные блоки и границы плотностных неоднородностей.
Цель работы: изучение роли сил динамического негидростатического давления и вязких напряжений, действующих в зонах литосферной субдукции посредством физико-математического моделирования.
Метод исследования: применено численное термодинамическое моделирование тектонического строения литосферы и астеносферы в зоне субдукции.
Основные результаты исследования: показано, что относительная роль сил динамического давления преобладает в зонах, характеризуемых горизонтально удлиненными астеносферными движениями в мантии, как это происходит под протяженными океаническими и континентальными литосферными плитами. Под микро-плитами и вблизи их границ роли сил динамического давления и вязких напряжений сравнимы между собой. В областях термических диапиров преобладает роль вязких напряжений.
Ключевые слова
Полный текст
Введение
Целью данной работы является изучение роли сил динамического негидростатического давления и вязких напряжений, действующих в зонах литосферной субдукции посредством физико-математического моделирования. Приведем некоторый обзор опубликованной литературы по данному направлению исследований. Проблеме физико-математического моделирования сил динамического давления и вязких напряжений в зонах литосферной субдукции посвящено достаточно много работ как отечественных [1-5, 8], так и зарубежных [9-14] ученых.
Динамическая топография в тектонически активных зонах Земли, например, в зонах субдукции, зонах рифтов и в областях подъема термических или химических диапиров может быть обусловлена вязкими напряжениями или силами негидростатического давления, причем относительная роль этих сил зависит от геометрии областей, охваченных конвективными движениями в мантии Земли. Вязкие напряжения и силы динамического давления действуют также и в земных недрах, например, подпирая снизу или подсасывая сверху субдуцирующие блоки литосферы, и преимущественное влияние вязких напряжений и динамического давления также определяются геометрическими характеристиками той части мантии, в которой происходят конвективные движения.
Результаты и обсуждение
Рассмотрим модель конвекции в прямоугольной ячейке 0 < x < L, 0 < z < d, с началом координат в основании верхней мантии на глубине d, вертикальной осью z, направленной вверх, и горизонтальной осью x вдоль основания верхней мантии. Ячейка заполнена однородной жидкостью плотностью ρ с коэффициентом вязкости η, горизонтальные границы z = 0 и z = d изотермичны с температурами T(z = 0) = T0 и T(z = d) = T1, а вертикальные границы x = 0 и x = L считаются адиабатичными, на которых ∂T/∂x = 0. Безразмерные линеаризованные уравнения, определяющие возмущения термомеханического состояния среды в ячейке при бесконечном числе Прандтля в приближении Буссинеска. имеют вид уравнений (7.3.11) – (7.3.14) в [14], которые в обозначениях настоящей работы могут быть записаны как
(1)
(2)
(3)
(4)
где сохранен член , описывающий нестационарную задачу, и знак при изменен, так как ось z направлена вверх. Уравнения (1)–(4) есть соответственно - и -компоненты уравнения движения, уравнение неразрывности и уравнение теплопереноса, в которых и – декартовы координаты, p – динамическое (негидростатическое) давление, – тензор вязких напряжений, – плотность, – ускорение силы тяжести, – удельная теплоемкость при постоянном давлении, – абсолютная температура, – коэффициент теплопроводности, – оператор Лапласа, а символ с индексом обозначает частную производную по координатам x, z и времени .
В (1)–(4) = (κ / ρcp) – коэффициент температуропроводности, и для приведения уравнений к безразмерной форме в качестве новых единиц измерения координат и выбраны мощность слоя , скорость – величина , время – , температура и ее возмущения – характерный перепад температуры = (T0 – T1) > 0, напряжения и давление – величина ( ). В (2) безразмерное число Рэлея есть
> 0, (5)
где α – коэффициент теплового расширения.
Рассматривая двумерную конвекцию в плоском горизонтальном слое первоначально покоящейся жидкости, в которой имеется вертикальный градиент температуры , с невозмущенным термомеханическим состоянием покоя с постоянным вертикальным градиентом температуры и чисто кондуктивным переносом тепла, можно искать решение уравнений (1)–(4) с экспоненциальной зависимостью от времени по закону exp(γt).
При условии свободных непроницаемых изотермических горизонтальных и адиабатических вертикальных границ ищем решение (1)–(4), при постоянных (безразмерных) < 0 и в виде:
, , , , , , , , , , , , , , , , , (6)
где все не зависящие от координат величины A, B, C, D в (6) зависят от времени t по экспоненциальному закону , а k = πd×L–1 есть (безразмерное) волновое число. Подставляя (6) в (1)–(4), находим для безразмерного инкремента
(7)
Условие возникновения конвекции = 0 дает . Эта величина достигает минимума при , и при =1, ≈ 658. Если конвекция происходит в горизонтальном слое неограниченной длины, то возникают ячейки с пространственным периодом d . В случае если , , переменны, можно для оценки инкремента конвективной неустойчивости воспользоваться формулой (7), подставив в нее средние значения , , .
Рассмотрим подробнее вывод формулы (7) из уравнений (1)–(4). Пусть начальное возмущение температуры задается в (6) как с C>0. Это означает, что возмущение температуры в левой части ячейки положительно, а в правой части – отрицательно, т. е. в левой части ячейки вещество всплывает, а в правой – опускается, и, следовательно, конвективное движение жидкости происходит по часовой стрелке. Из (3) следует B = –Ak/π. Так как Δθ = –(π2 + k2) C·cos kx·sin πz, ∂tθ = γ C cos kx sin πz, то из (4)
C = .
где при C > 0 и Tz < 0 должно быть A < 0. Подставляя выражения (6) в уравнения (1) и (2), вычитая одно из уравнений из другого и сокращая полученный результат на A, приходим к формуле (7). Из (1) находим D=Aη(k2 + π2)/k, где k = π/L, т. е. D < 0.
Согласно выражениям в верхней строке (6)
, , , ,
при C > 0, A < 0, B > 0, D <0 компоненты скорости vx и vz соответствуют движению жидкости по часовой стрелке, т. е. всплыванию жидкости в левой части ячейки и опусканию жидкости в правой части ячейки. На верхней границе ячейки (при z = 1) возмущение динамического давления `p = – D·coskx= – Aη(k2 + π2)/k·coskx. Сила давления действующая изнутри жидкости на верхнюю границу положительна в левой части ячейки (т.е. «подпирает» границу снизу) и отрицательна в правой части ячейки (т. е. «засасывает» границу вниз). Сравним силу негидростатического давления на верхней границе ячейки с вертикальной силой вязких напряжений, действующей со стороны жидкости на верхнюю границу ячейки.
Нормальная компонента тензора вязких напряжений
τzz = – 2Dk2/(π2 + k2)×coskxcosπz = –2ηAk×coskxcosπz,
и на верхней границе z = 1, cos πz = –1, τzz = 2ηAk×cos kx, т.е. при A < 0 оказывается, что τzz на верхней границе отрицательна в левой части ячейки и положительна в правой части ячейки. Так как сила, действующая со стороны жидкости на единицу обтекаемой поверхности границы с внешней нормалью ni, равна [5, формула (15.14), в которой изменен знак нормали ni]:
fi = – pni + τiknk,
то вертикальная сила, соответствующая вязким напряжениям, равна –τzz, так как направленная внутрь жидкости нормаль на верхней границе nz = –1. Следовательно, вязкие напряжения в левой части ячейки действуют на верхнюю границу, как сила «подпора» снизу (в положительном направлении оси z), а в правой части ячейки – как сила «подсоса» вниз. Сравним конвективные силы вязких напряжений и негидростатического давления, действующие на верхнюю границу ячейки. Отношение этих сил
, (8)
откуда видно, что при k = π (т. е. в изометрической ячейке с отношением сторон 1 : 1, т. е. при L = d) эти силы равны между собой. В случае, например, удлиненной ячейки, для которой k < π, на верхней границе ячейки преобладает сила возмущенного негидростатического давления. В сильно удлиненной ячейке, для которой k << π, действие вязких напряжений на верхней границе пренебрежимо мало по сравнению с действием сил возмущенного динамического давления. Соотношение (8) справедливо не только на поверхности ячейки, но и во всем ее объеме. Следует отметить, что силы негидростатического давления и вязкие напряжения в земных недрах действуют, так сказать, «однонаправленно», т. е. в одну сторону, и этот вывод не связан именно с конвективной природой движения, а приложим к движениям различной природы.
Поскольку основные результаты работы были получены на модельных примерах, то, по-видимому, необходимо привести ссылки на опубликованные нами работы [1-4] и привести графический пример о сравнении результатов наших модельных расчетов с опубликованными результатами практических геолого-геофизических исследований [6].
Рисунок 1 – Глубинный геолого-геофизический разрез квазистационарного распределения безразмерной функции тока в мантийном субдукционном клине
На рисунке 1 изображен глубинный геолого-геофизический разрез квазистационарного распределения безразмерной функции тока в мантийном субдукционном клине, расположенном на акватории Черноморской литосферной плиты, погружающейся под Русскую континентальную литосферную плиту, с учетом эффектов диссипативного нагрева и конвективной неустойчивости для неньютоновской реологии мантии и концентрации воды Cw = 3×10-1 весовых %% в горных породах мантии. Вертикальная стрелка (с) показывает восходящий конвективный поток. В верхней части рисунка приведен реальный геолого-геофизический разрез мантии [6] по профилю от полуострова Крым по направлению на северо-запад, к Русской литосферной плите. Условные обозначения: 1 – осадочные породы; 2 – дислоцированные породы молодого фундамента; 3 – поверхность молодого фундамента; 4 – поверхность дорифейского фундамента (V=6.2-6.5 км/с); 5 – «гранитный» слой; 6 – породы основного состава (V=7.0 км/с); 7 – породы коро-мантийного слоя; 8 – поверхность коро-мантийного слоя (V=7.5-7.6 км/с); 9 – граница Мохоровичича; 10 – сейсмические горизонты верхней мантии; 11 – слои с пониженной скоростью в верхней мантии по данным ГСЗ; 12 – поверхность астеносферного слоя по геотермическим данным; 13 – поверхность астеносферного слоя по данным МТЗ; 14 – глубинные тектонические разломы; 15 – очаги землетрясений.
Из рисунка 1 видно, что в расчетной модели субдукционной зоны (нижняя часть рисунка) с погружением Черноморской литосферной плиты под Русскую был правильно рассчитан угол наклона этой плиты (27о), который соответствует наклону верхней кромки этой литосферной плиты, фиксируемому на реальном геолого-геофизическом разрезе (верхняя часть рисунка). Кроме того, из рис. 1 видно, что точно рассчитано местоположение конвективной вихревой структуры квазистационарного распределения безразмерной функции тока, связанной с формированием восходящего конвективного потока (вертикальная стрелка – с). Над этой конвективной зоной, в верхней части рисунка, на реальном геолого-геофизическом разрезе наблюдается зона плавления горных пород мантии, которую мы ассоциируем с восходящим мантийным термическим диапиром.
В качестве примера рассмотрим тектонически активную окрестность зоны субдукции, и качественно сравним силу динамического (негидростатического) давления и вязкие напряжения, действующие на субдуцирующий литосферный блок и подошву динамической топографии в этой области. Так как крупномасштабные циркуляционные движения под движущейся и субдуцирующей Черноморской субокеанической литосферной плитой и Русской континентальной плитой, с которой сталкивается субокеаническая плита, происходят внутри тех частей верхней мантии, которые сильно удлиненны в горизонтальной направлении, то в рамках рассмотренной конвективной модели эти циркуляционные движения характеризуются условием d << L, или, в безразмерном виде, k << π в формуле (8). Учитывая изложенное выше и ранее опубликованные работы [1-4] это может означать, что в окрестности зон субдукции на субдуцирующие литосферные блоки и «подошву» настилающей литосферы действуют преимущественно силы динамического давления, а вязкие напряжения несущественны. Это условие (хотя и без всякого обоснования), используется в [7] в параграфе § 6.11 об угле субдукции. Можно, также, отметить, что расчет динамической топографии в [10] выполняется путем вычисления упругого изгиба верхней части коры, подпираемой снизу динамическим давлением в вязком течении, происходящим в слое нижней коры, причем этот достаточно тонкий слой очень сильно вытянут в горизонтальном направлении. При этом также не учитываются вязкие напряжения, действием которых авторы пренебрегают, не приводя каких-либо обоснований.
Однако в окрестности зон субдукции микро-плит, например, Черноморской [1], Амурской [2], Адриатической [3] и Палео-Уральской [4], движения в астеносфере оказываются примерно изометричны, и роли сил динамического давления и вязких напряжений оказываются сравнимы между собой. Этим, возможно, объясняется то, что субдукция микро-плит и малых плит происходит под достаточно малыми углами к горизонту, так как субдуцирующий блок поддерживается снизу и «подсасывается» сверху не только силами динамического давления, но и сравнимыми силами вязких напряжений. Напротив, динамическая топография, формирующаяся над восходящими термическими диапирами, относительно узкими в горизонтальном направлении, связана с астеносферными потоками, в которых k>>π в формуле (8). Следовательно, динамическая топография над диапирами оказывается обусловлена преимущественно вязкими напряжениями.
Сравнивая с другими [7-14] (и огромным количеством ссылок по ним) полученные в данной работе результаты, следует отметить, что многочисленные термомеханические модели мантии в зонах субдукции показали, что конвекция в виде поперечных валообразных вихрей и формирующихся термических диапирах никогда ранее не анализировалась, поскольку модели с чрезвычайно малым углом субдукции и достаточно большой скоростью субдукции ранее никем не исследовались.
Заключение и выводы
Показано, что относительная роль сил динамического давления и вязких напряжений, действующих в областях верхней мантии, характеризуемых астеносферными течениями в окрестности активных тектонических зон, зависит от соотношения горизонтального и вертикального масштабов течений в астеносфере. Если горизонтальный масштаб движений значительно превышает их вертикальный масштаб, то роль сил динамического давления существенно преобладает над ролью вязких напряжений, и последними можно пренебречь. Так, в окрестности зон субдукции протяженных плит при вычислении динамической топографии и угла субдукции можно пренебречь вязкими напряжениями и учитывать только силы динамического давления. В зонах субдукции микро-плит следует учитывать как динамическое давление, так и вязкие напряжения, роли которых сравнимы, чем, по-видимому, объясняются малые углы субдукции микро-плит. Динамическая топография над термическими диапирами, напротив, обязана своим происхождением преимущественно вязким напряжениям. Новизна настоящей работы заключается в том, что в окрестности зон субдукции литосферных плит движения в астеносфере оказываются примерно изометричны, и роли сил динамического давления и вязких напряжений оказываются сравнимы между собой. Кроме того, новизна результатов данной работы может быть обусловлена еще и тем, что динамическая топография, формирующаяся над восходящими термическими диапирами, относительно узкими в горизонтальном направлении, связана с астеносферными потоками, в которых k >> π в формуле (8). Полученные в настоящей работе результаты могут быть использованы при решении некоторых задач теории упругости литосферных плит.
Об авторах
Сергей Владиленович Гаврилов
Институт физики Земли им. О.Ю. Шмидта РАН»
Автор, ответственный за переписку.
Email: gavrilov@ifz.ru
доктор физико-математических наук, главный научный сотрудник лаборатории 102
Россия, МоскваАндрей Леонидович Харитонов
Институт земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова РАН»
Email: gavrilov@ifz.ru
кандидат физико-математических наук, ведущий научный сотрудник лаборатории главного магнитного поля Земли
Россия, МоскваСписок литературы
- Гаврилов, С. В. Геотермодинамическая модель предполагаемой палеозоны литосферной субдукции в районе Черноморской впадины и ее связь с металлогенической зонально-стью Крыма и Кавказа / С. В. Гаврилов, А. Л. Харитонов. – Текст : непосредственный // Региональная геология и металлогения. – 2021. – № 87. – С. 4-16. doi: 10.52349/0869-7892-2021-87-04-16.
- Гаврилов, С. В. О субдукции Амурской микроплиты и конвективном механизме выноса диссипативного тепла и углеводородов из мантийного клина в Охотском море к востоку от острова Сахалин / С. В. Гаврилов, А. Л. Харитонов. – Текст : непосредственный // Вестник Академии наук Республики Башкортостан. – 2022. – Т. 42. – № 1(105). – С. 5-12. doi: 10.24412/1728-5283_2022_1_5-12.
- Гаврилов, С. В. Исследование величины о формировании аномального теплового пото-ка в бассейне Паннония и зоне Вардар при субдукции Адриатической плиты под Евро-азиатскую плиту / С. В. Гаврилов, А. Л. Харитонов. – Текст : непосредственный // International Journal of Professional Science. – 2021. – № 9. – С. 27-39. doi: 10.54092/25421085_2021_9_27.
- Гаврилов, С. В. Моделирование глубинного геодинамического строения зоны субдук-ции Русской палеоплиты под литосферу Уральского палеоокеана и связанное с субдук-цией распределение месторождений углеводородов / С. В. Гаврилов, А. Л. Харитонов. – Текст : непосредственный // Уральский геологический журнал. – 2021. – № 5(143). – С. 3–19.
- Ландау, Л. Д. Гидродинамика / Л. Д. Ландау, Е. М. Лифшиц. – М. : Наука. ГИФМЛ, 1986. – 736 с. – ISBN 1ЯБМ5-9221-0121-8. – Текст : непосредственный.
- Соллогуб, В. Б. Литосфера Украины / В. Б. Соллогуб. – Киев : Наукова думка, 1986. – 184 с. ISBN отсутствует. – Текст : непосредственный.
- Теркотт, Д. Л. Геодинамика / Д. Л. Теркотт, Дж. Шуберт. – М. : Мир, 1985. – 732 с. – Текст : непосредственный.
- Трубицын, В. П. Численная модель образования совокупности литосферных плит и их прохождения через границу 660 км / В. П. Трубицын, А. П. Трубицын / Физика Земли. – 2014. – № 6. – С. 138–147.
- Billen M., Hirth G. Newtonian versus non-Newtonian Upper Mantle Viscosity: Implications for Subduction Initiation // Geophys. Res. Lett. – 2005. – V.32. – (L19304). doi: 10.1029/2005GL023458.
- Clark M.K., Bush J.W.M., Royden L.H. Dynamic topography produced by lower crustal flow against rheological strength heterogeneities bordering the Tibetan Plateau // Geophys. J. International. – 2005. – V.162. – PP. 575–590.
- Gerya T.V. Future directions in subduction modeling // Journal of Geodynamics. – 2011. – V. 52. – PP. 344–378.
- Karig D. E. Origin and development of marginal basins in the Western Pacific // Journal Geo-physical Researches. – 1971. – V. 76. – № 11. – PP. 2542-2561. https://doi.org/10.1029/JB076i011p02542
- Miyashiro A. Metamorphism and related magmatism in plate tectonics // Am. Journ. Sci. – 1972. – V. 272. – PP. 629-656.
- Schubert G., Turcotte D.L., Olson P. Mantle Convection in the Earth and Planets. – New York: Cambridge University Press, 2001. – 940 p. ISBN 9780511612879. – DOI: https://doi.org/10.1017/CBO9780511612879
Дополнительные файлы
![](/img/style/loading.gif)