Физика Земли, 2023, № 3, стр. 63-77

О механизмах генерации избыточного горизонтального сжатия в континентальной коре

Ю. Л. Ребецкий *

Институт физики Земли им. О.Ю. Шмидта РАН
г. Москва, Россия

* E-mail: reb@ifz.ru

Поступила в редакцию 25.10.2022
После доработки 05.12.2022
Принята к публикации 09.12.2022

Полный текст (PDF)

Аннотация

Предлагается рассматривать процессы денудации поверхности и магматизма коры для объяснения появления в ней повышенных напряжений горизонтального сжатия – избыточных по отношению в литостатическому давлению. “Откапывание” породы приводит только к частичной ее разгрузке, за счет снятия веса вышележащих пород, если она в глубине находилась за пределом текучести. Это связано с тем, что в случае “откапывания” разгрузка идет по упругому закону. Из-за этого в породе возникают остаточные напряжения горизонтального сжатия, полученные на стадии закритического катакластического течения. Другой механизм формирования дополнительных сжимающих напряжений в коре связан с процессами вулканизма и магматизма. Подъем магмы по субвертикальным коровым разломам и системам трещин возможен только в условиях, когда во фронте распространения ее давление превышает уровень горизонтального сжатия пород. Следствием этого ниже фронта распространения магмы уровень напряжений горизонтального сжатия в породах повышается до уровня ее давления. Поскольку давление в подкоровом или внутрикоровом магматическом очаге близко к литостатическому давлению вышележащих пород, то выше фронта распространения магмы в разломе напряжения, нормальные к нему, превышают уровень вертикального сжатия. Таким образом, магматизация коры способна изменить режим напряженного состояния горизонтального растяжения на горизонтальный сдвиг.

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

ВВЕДЕНИЕ

Генезис напряжений избыточного горизонтального сжатия в континентальной коре в настоящее время наиболее часто связывают с проявлением дальнодействующего давления от границ литосферных плит [Molnar, Tapponier, 1975; Зоненшайн и др., 1978; Ма Си Юань, 1990; Nazarova, Nazarov, 2002; Ружич, 1997; Ружич и др., 2016; и др.], иногда – с планетарными напряжениями, вызванными вращением Земли или другими космическими факторами [Aвсюк, 2001; Dobretsov, Chumakov, 2001; Glukhovskii, 2005; Dobretsov, 2015; Викулин, 2008; Вихри …., 2004; и др.]. В ряде работ второй половины двадцатого века наличие таких напряжений объясняли остаточными напряжениями предыдущих эпох нагружения [Пономарев, 1969; 2008; Волох и др., 1972; Herget, 1973; Sykes, Sbar, 1973; Voigth, Pierre, 1974; Марков 1977; 1980].

Первый тип генезиса избыточного горизонтального давления прямо обусловлен движением литосферных плит, их конвергенцией в виде субдукции или коллизии. При этом источником сил, перемещающих плиты, определяют давление со стороны океанского хребта и силу тяжести от погруженной и уплотнившейся части слэба в зонах субдукции. Иногда также говорят о роли касательных напряжений, действующих на подошву океанической литосферы со стороны астеносферной мантии в направлении от зоны океанского спрединга к зоне субдукции. Но, в любом случае, первый из представленных выше источников сил принимается главным [Wilson, 1966]. С конца прошлого века эти возможные источники сил, определяющие движения тектонических плит, прямо связывают с общемантийной одноярусной или двухярусной конвекцией [Трубицын, 2012; Трубицын Рыков, 2000; Лобковский, Котелкин, 2000; Добрецов и др., 2001].

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

Появление в конце прошлого века в сейсмологии массовых определений фокальных механизмов землетрясений, определяемых из анализа волновых форм [Dziewonski et al., 1981], и развитие технологии расчета тензора сейсмического момента [Кейлис-Борок, 1950] дало возможность более достоверно оценить уровень сброшенных напряжений в очагах не только сильных, но и среднесильных и даже слабых землетрясений. Было установлено, что существует достаточно сильный разброс значений снимаемых напряжений в одном регионе. Разброс значений этих напряжений увеличивался при снижении магнитуды землетрясений, а с повышением магнитуды землетрясения (М > 6.5–7.0) наблюдалось снижение уровня сбрасываемых напряжений.

Оказалось, что наименьший уровень сбрасываемых напряжений наблюдается в зонах субдукции литосферных плит, а наибольший – во внутриконтинентальных орогенах и щитах. При этом уровень касательных напряжений, снимаемых в очаге, для землетрясений с магнитудами более 8 в зонах субдукции – менее 0.3–0.5 МПа. Существуют оценки, показывающие, что на большом статистически представительном множестве событий излученная в сейсмических волнах энергия составляет 10–20% от общей упругой энергии, снятой в очаге землетрясения [Bath, Duda, 1964; Kanamori, Anderson, 1975]. В свою очередь, эта снятая упругая энергия составляет 10–20% от девиаторной части энергии упругого деформирования массива [Ребецкий, 2007]. Используя эти обобщения и выражения для КПД землетрясения и КПД упругой разгрузки после землетрясения [Ребецкий, 2007], можно показать, что снимаемые напряжения в очаге землетрясения должны составлять 1/3–1/5 часть от максимальных касательных напряжений, действующих в массиве.

Таким образом, можно ожидать, что уровень максимальных касательных напряжений в зонах субдукции составляет 1.5–2.5 МПа, что определяет для горизонтальных сжимающих напряжений величину превышения над литостатическим давлением в 3–5 МПа. Это средний уровень дополнительных сжимающих напряжений, который на дистанции в первые тысячи километров может снизиться до 0.3–0.5 МПа (при размере силового индентора в первые сотни километров). Конечно же, его недостаточно для объяснения формирования избыточных напряжений горизонтального сжатия 30–50 МПа, наблюдаемого в горных выработках внутриконтинентальных орогенов.

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

Следует помнить, что сделанные выше оценки соотношения действующих и снимаемых напряжений не следует относить к конкретным землетрясениям. Они работают именно для большого множества событий как средние. Примером отклонений от сделанной оценки можно рассматривать землетрясение Тохоку, где сброшенные напряжения в очаге в верхней части коры (глубины до 20 км) оказались даже выше уровня действовавших здесь девиаторных напряжений [Ребецкий, Полец, 2021]. Это привело к смене геодинамического типа напряженного состояния с горизонтального сжатия на горизонтальное растяжение, которое проявляется в соответствующем типе фокальных механизмов землетрясений уже более 10 лет.

В настоящей работе будут представлены результаты прогноза избыточных напряжений горизонтального сжатия, генерация которых связывается с возможностью длительного существования в породах верхних слоев коры остаточных напряжений, обусловлена действием массовых сил тяжести (МСТ), а также с процессами вулканизма и магматизма в коре.

ОСТАТОЧНЫЕ НАПРЯЖЕНИЯ ОТ ДЕЙСТВИЯ МСТ

Напряженное состояние коры, отвечающее действию только массовых сил тяжести, имеет два глубинных уровня, характеризующих разную реакцию горных пород на нагружение. В верхней части коры (в зависимости от флюидного режима: сотни метров – 1–2 км) породы в условиях действия МСТ могут деформироваться упруго, а в более глубоких горизонтах формируются необратимые деформации. Из-за разупрочняющего действия флюида эти необратимые деформации обусловлены преодолением предела трещинной – катакластической текучести, а не истинного предела пластичности кристаллов и зерен. Заметим, что чисто упругое поведение пород не относится к зонам разломов, а отвечает внутренним межразломным блокам кристаллической коры.

Как в области чисто упругого, так и упруго-катакластического деформирования напряженное состояние горных массивов, отвечающее действию только собственного веса, определяется предположением о том, что вертикальные напряжения равны весу вышележащей колонки пород, к которому добавлены условия горизонтального стеснения:

(1)
${{\sigma }_{{zz}}} = - \rho gH,\,\,\,{{\varepsilon }_{{xx}}} = {{\varepsilon }_{{yy}}} = 0.$

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

Состояние в упругой области

В случае чисто упругого поведения квазиоднородного породного массива [Динник, 1926; Jager, 1962; Price, Cosrove, 1990] горизонтальные напряжения от МСТ определяются прямо из условия стеснения (второе в (1)) и закона Гука в виде выражений:

(2)
$\sigma _{{xx}}^{e} = \sigma _{{yy}}^{e} = - \frac{\nu }{{1 - \nu }}{{\rho }_{с}}gz.$

Здесь: $\nu $ – коэффициент Пуассона; ${{\rho }_{с}}g$ – средний удельный вес колонки вышележащих пород (произведение средней плотности на ускорение свободного падения) мощности $H$; ось $z$ направлена вертикально вниз. Касательные напряжения в принятой системе координат отсутствуют ${{\sigma }_{{ij}}} = 0$ (i, j = x, y, z). Отношение ${{{{\sigma }_{{xx}}}} \mathord{\left/ {\vphantom {{{{\sigma }_{{xx}}}} {{{\sigma }_{{zz}}}}}} \right. \kern-0em} {{{\sigma }_{{zz}}}}}$ при $\nu = 0.25$ близко к 0.333.

Из выражений (1) и (2) следует линейная зависимость от глубины максимальных касательных напряжений τ и всестороннего давления p:

(3)
$\begin{gathered} {{\tau }^{e}} = {{\left( {\sigma _{{xx}}^{e} - \sigma _{{zz}}^{{}}} \right)} \mathord{\left/ {\vphantom {{\left( {\sigma _{{xx}}^{e} - \sigma _{{zz}}^{{}}} \right)} 2}} \right. \kern-0em} 2} = \frac{{1 - 2\nu }}{{2\left( {1 - \nu } \right)}}{{p}_{{lt}}}, \\ {{p}^{e}} = {{ - \left( {2\sigma _{{xx}}^{e} + \sigma _{{zz}}^{{}}} \right)} \mathord{\left/ {\vphantom {{ - \left( {2\sigma _{{xx}}^{e} + \sigma _{{zz}}^{{}}} \right)} 3}} \right. \kern-0em} 3} = \frac{{1 + \nu }}{{3\left( {1 - \nu } \right)}}{{p}_{{lt}}}, \\ {{p}_{{lt}}} = {{\rho }_{с}}gz. \\ \end{gathered} $

Для значения $\nu = 0.25$: ${{\tau }^{e}} = 0.33{{p}_{{lt}}}$ и ${{p}^{e}} = 0.56{{p}_{{lt}}}$. Верхний индекс “e” показывает, что данные компоненты тензора напряжений являются результатом чисто упругой реакции горных пород на действие массовых сил.

Граница кровли области катакластического течения

Линейный рост с глубиной касательных напряжений (3) от действия МСТ должен привести к нарушению упругого равновесного состояния пород по двум наиболее вероятным причинам: 1) хрупкое разрушение; 2) достижение предела упругости. Существует гипотеза, что для пород коры достижение истинного пластического состояния (формирование необратимых деформаций за счет внутризернового или внутрикристаллического дислокационного течения) наблюдается на границе Мохо [Николаевский, 1979]. Поэтому нарушение упругого равновесного состояния в породах верхней и средней коры возникает из-за формирования и активизации множества разномасштабных дефектов хрупкой прочности, что приводит к появлению необратимых деформаций, вызванных катакластическим – трещинным течением (псевдопластичность – геомех.).

Для расчета положения кровли области катакластического течения от МСТ можно использовать критерий Друккер–Прагера:

(4)
$\begin{gathered} {{\tau }_{C}} = {{I}_{2}} + {{k}_{f}}\left( {{{I}_{1}} + {{p}_{{fl}}}} \right) \leqslant {{\tau }_{f}} \\ {\text{при}}\,\,\,\,{{I}_{2}} = {{2\tau } \mathord{\left/ {\vphantom {{2\tau } {\sqrt 3 }}} \right. \kern-0em} {\sqrt 3 }}\,\,\,\,{\text{и}}\,\,\,\,{{I}_{1}} = - p, \\ \end{gathered} $
где: ${{\tau }_{C}}$ – кулоновы напряжения; ${{I}_{1}}$ и ${{I}_{2}}$ – первый и второй инварианты тензора напряжений; ${{\tau }_{f}}$ – внутреннее сцепление (cohesion); ${{k}_{f}}$ – коэффициент внутреннего трения; ${{p}_{{fl}}} = \lambda {{p}_{{lt}}}$ – флюидное давление в трещинном пространстве горных пород (${{\lambda }} = 0.38$ – давление близко к гидростатическому, ${{\lambda }} \approx 1$ – давление близко к литостатике).

Используя предельное соотношение в критерии (4) в виде знака равенства, находим глубину Hp перехода породы в катакластическое течение в условиях действия только МСТ:

(5)
${{H}_{p}} = \frac{{3\left( {1 - \nu } \right){{\tau }_{f}}}}{{\left[ {\left( {1 - 2\nu } \right)\sqrt 3 - {{k}_{f}}\left( {1 + \nu } \right) + 3{{k}_{f}}\left( {1 - \nu } \right)\lambda } \right]\rho g}}.$

Для параметров прочности ${{k}_{f}} = 0.6$ (средний коэффициент внутреннего трения большинства горных пород, находящихся в условиях среднего уровня обжимающего давления), ${{\tau }_{f}} = 2.5$ МПа (прочность сцепления сплошных образцов алевролитов, аргиллитов и трещиноватых пород верхней и средней коры), удельного веса – ${{\rho }_{с}}g{\text{\;}} = 2.7$ Г/см3 и коэффициента Пуассона $\nu = 0.25$ в условиях сухой породы (${{p}_{{fl}}} = 0$), находим ${{H}_{p}} \approx 1.7$ км. При гидростатическом законе распределения флюидного давления по глубине ${{p}_{{fl}}} \approx 0.38{{p}_{{lt}}}$ – открытая до поверхности трещинная пористость, переход в катакластическое состояние происходит при ${{H}_{p}} \approx 0.3$ км. Эту глубину можно принимать за более-менее стандартную для осадочных бассейнов.

Таким образом, уже верхняя кора практически вся, за исключением первых сотен метров (возможно около километра), находится за пределом катакластической текучести.

Состояние катакластического течения

Накопление неупругих деформаций в области катакластического деформирования ($H > {{H}_{p}}$) осуществляется таким образом, чтобы уменьшить разницу между напряжениями, действующими в вертикальном и горизонтальном направлении. В этой области в горизонтальном направлении формируются необратимые деформации удлинения $\Delta \varepsilon _{{xx}}^{p} = \Delta \varepsilon _{{yy}}^{p} > 0$, а в вертикальном необратимые деформации укорочения $\Delta \varepsilon _{{zz}}^{p} < 0$. В сумме необратимые деформации изменения объема можно считать равными нулю $\Delta \varepsilon _{{xx}}^{p} + \Delta \varepsilon _{{yy}}^{p} + \Delta \varepsilon _{{zz}}^{p} = 0$, а можно, используя значения интенсивности катакластического сдвига, учитывать дилатансионное изменение объема [Николаевский, 1979].

Условия бокового стеснения (1) и необратимые деформации удлинения приводят к появлению в горизонтальном направлении дополнительных упругих деформаций укорочения ($\Delta \varepsilon _{{xx}}^{{ep}} = \Delta \varepsilon _{{yy}}^{{ep}} = - \Delta \varepsilon _{{xx}}^{p} = - \Delta \varepsilon _{{yy}}^{p}$) и, как следствие, к увеличению напряжений горизонтального сжатия за счет дополнительных напряжений $\Delta \sigma < 0$. Верхний индекс в деформациях: “p” и “еp” – определяет соответственно приращение необратимой и упругой деформации в области катакластического течения. В области закритического состояния происходит вертикальное уплотнение горных пород [Ребецкий, 2008а; 2008б; 2008в] и “взводится упругая пружина”, определяющая увеличение напряжений горизонтального сжатия для катакластической стадии деформирования ($\sigma _{{ii}}^{p}$) на величину $\Delta \sigma $ относительно напряжений чисто упругой стадии ($\sigma _{{ii}}^{e}$):

(6)
$\sigma _{{ii}}^{p} = \sigma _{{ii}}^{e} + \Delta \sigma \,\,\,\,{\text{при}}\,\,\,\,\Delta \sigma < 0,\,\,\,\,i = x,y.$

Для напряжений дополнительного горизонтального сжатия из предельного соотношения в виде (4), используя для флюидного давления гипотезу о линейном изменении с глубиной (${{p}_{{fl}}} = \lambda {{p}_{{lt}}}$), находим:

(7)
$\Delta \sigma = - \frac{{3~{{\tau }_{f}}}}{{\sqrt 3 + 2{{k}_{f}}}}\left( {\frac{z}{{{{H}_{p}}}} - 1} \right)\,\,\,\,\,{\text{при}}\,\,\,\,z > {{H}_{p}}.$

Изменения напряженного состояния при вертикальных перемещениях горных пород, сопровождающегося эрозией поверхности

В случае действия вертикальных восходящих движений и денудации поверхности или наоборот из-за денудации и компенсационного восходящего движения (далее будем этот механизм определять как “откапывание” породы), вертикальная нагрузка, действующая на элементарный объем геосреды, уменьшается. При этом разгрузка породы будет происходить по упругому закону, подобно тому, как это записано в выражениях (1) и (2).

Если порода находится на глубине выше кровли коры, испытывающей катакластическое течение (${{H}^{0}} \leqslant {{H}_{p}}$), то в ней имеет место чисто упругое напряженное состояние и разгрузка приводит только к пропорциональному уменьшению напряжений. Другая ситуация будет иметь место для пород, находившихся до начала эксгумации на глубине ${{H}^{0}} > {{H}_{p}}$. В них соотношение горизонтальных и вертикальных напряжений не отвечает пропорции ${\nu \mathord{\left/ {\vphantom {\nu {\left( {1 - \nu } \right)}}} \right. \kern-0em} {\left( {1 - \nu } \right)}}$. Поэтому “откапываемая” (поднимающаяся к поверхности) порода будет иметь уровень напряжений горизонтального сжатия больший, чем имела порода на этом же глубинном уровне до начала процесса денудации поверхности. Уровень дополнительных напряжений сжатия по отношению к напряжениям $\sigma _{{ii}}^{e}$, i = x, y (напряжения чисто упругого состояния по Диннику (2)) в этом случае определится выражением:

(8)
$d\sigma = - \frac{{3~{{\tau }_{f}}}}{{\sqrt 3 + 2{{k}_{f}}}}\left( {\frac{{{{H}^{0}}}}{{{{H}_{p}}}} - 1} \right)~.$

Здесь ${{H}^{0}}$ – наибольшая глубина, на которой находилась порода до начала денудации. Из (8) следует, что дополнительные напряжения $d\sigma $ постоянные по глубине.

Таким образом, на стадии “откапывания” породы напряжения горизонтального сжатия определяются выражением:

(9)
$\sigma _{{ii}}^{{Ue}} = \sigma _{{ii}}^{e} + d\sigma ,\,\,\,\,i = x,y.$

Здесь верхний индекс “Ue” определяет, что эти напряжения отвечают стадии упругой разгрузки. Соответственно превышение напряжений горизонтального сжатия на стадии “откапывания” над их значениями на стадии осадконакопления определится выражением:

(10)
$d\sigma - \Delta \sigma = - \frac{{3~{{\tau }_{f}}}}{{\sqrt 3 + 2{{k}_{f}}}}\frac{{{{H}^{0}} - z}}{{{{H}_{p}}}}~\,\,\,\,{\text{при}}\,\,\,z < {\text{\;}}{{H}^{0}}.$

При определенном уровне эксгумации породы $dH$ можно получить состояние, в котором напряжения горизонтального сжатия превысят сжимающие напряжения, действующие в вертикальном направлении. Используя (10) и (1)–(2), находим эту глубину:

(11)
${{H}_{{13}}} = \frac{{3{{\tau }_{f}}}}{{2\rho g{{H}_{p}}\left( {\sqrt 3 + 2{{k}_{f}}} \right) - 3{{\tau }_{f}}_{{}}}}\left( {dH - {{H}_{p}}} \right).$

Здесь в рассматриваемой постановке задачи о действии только МСТ все главные нормальные напряжения одинаковые.

Из выражения (11) следует, что в случае, когда амплитуда вертикальных перемещений породы меньше глубины ее перехода в псевдопластическое состояние ($dH < {{H}_{p}}),$ в коре не возникает областей смены геодинамического типа напряженного состояния ($\left| {{{\sigma }_{{zz}}}} \right| > \left| {{{\sigma }_{{xx}}}} \right|$). Согласно (11), если амплитуды денудации поверхности была 5 км, то глубина смены типа напряженного состояния будет около 8 км, а при амплитуде денудации 10 км смена типа состояния происходит на существенно большей глубине 17–18 км.

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

Графический способ оценки величины избыточного горизонтального сжатия

Приведенные выше выражения вполне можно заменить простыми графическими построениями, которые позволят оценить уровень дополнительного сжатия, вызываемого “откапыванием” пород для заданных градиентов напряжений горизонтального сжатия упругой и упруго-пластической стадии деформирования под действием МСТ. На рис. 1 показано изменение напряжений в элементарном объеме пород осадочного бассейна, который вначале испытывал погружение и увеличение напряжений вертикального и горизонтального сжатия (стадия осадконакопления), а затем упругую разгрузку, вызванную эксгумацией породы. При этом по мере накопления вышележащих над данным объемом осадочных пород до глубины ${{H}_{p}}$ (путь нагружения ОА) имеет место чисто упругое деформирование, а глубже (путь нагружения АB) появляются необратимые деформации, приводящие к большему градиенту напряжений горизонтального сжатия в сравнении с упругой стадией.

Рис. 1.

Диаграмма нагружения для элементарного объема в глубине осадочного бассейна в процессе осадконакопления (стрелка вниз) и эксгумации (стрелка вверх). Вертикальная координатная ось – глубина нахождения в осадочном бассейне элементарного объема. Горизонтальная координатная ось напряжения сжатия (${{\sigma }_{{xx}}}$ и ${{\sigma }_{{zz}}}$). Длинный пунктир – глубина кровли упруго-катакластического деформирования ${{H}_{e}}$ и глубина переиндексации главных напряжений ${{H}_{{13}}}$ для стадии эксгумации породы; точка С – глубина, для которой для стадии эксгумации горизонтальные напряжения сжатия совпадают с вертикальными напряжениями (${{{{\sigma }_{{xx}}}} \mathord{\left/ {\vphantom {{{{\sigma }_{{xx}}}} {{{\sigma }_{{zz}}}}}} \right. \kern-0em} {{{\sigma }_{{zz}}}}} = 1$). Путь нагружения: ОА – чисто упругое нагружение объема в пределах области $H \leqslant {{H}_{e}}$; АВ – нагружение на стадии упруго-катакластического деформирования ($H \geqslant {{H}_{e}}$); ВD – стадия разгрузки – эксгумация породы; CD – участок разгрузки, для которого горизонтальное сжатие становится больше вертикального $\left| {{{\sigma }_{{xx}}}} \right| > \left| {{{\sigma }_{{zz}}}} \right|$. Слева показана ориентация осей главных напряжений (${{\sigma }_{1}} > {{\sigma }_{3}}$ при сжатии отрицательном) на стадии осадконакопления, а справа на стадии эксгумации. Характер изменения с глубиной другого нормального напряжения ${{\sigma }_{{yy}}}$ такой же, как и ${{\sigma }_{{xx}}}$. $\sigma _{{xx}}^{{Ue}}$ – горизонтальные напряжения, которые снимаются в процессе упругой разгрузки на пути BD2; $\Delta \sigma $ – дополнительное сжатие, возникающее при погружении пород на пути OAB из-за катакластического течения; $\sigma _{{xx}}^{r}$ – остаточные напряжения сжатия в породе на полном пути OABCD (полное “откапывание” породы). Пояснения в тексте.

Считается, что порода, находившаяся на глубине в точке B испытывает “откапывание”. Если в результате эта порода окажется на поверхности, то процессу эксгумации исследуемого элементарного объема породы отвечает путь разгрузки BD. В результате в коре уровень напряжений горизонтального сжатия будут отвечать линии ВСD. Точке С отвечает глубина подошвы слоя, в котором вследствие процессов эксгумации напряжения горизонтального сжатия начинают превышать уровень сжимающих напряжений, действующих в вертикальном направлении. В точке D вертикальные напряжения в породе, доставленной с глубины точки В, становятся равными нулю, а горизонтальные – не нулевые.

Для промежуточных стадий процесса, когда порода, находившаяся на глубине точки B, еще не вышла на поверхность, изменение напряжений горизонтального сжатия отражено линиями BDiI (I, i = 1, 2, 3, 4). Видно, что только для линий с индексами i = 3, 4 на поверхности горизонтальное сжатие не нулевое. При меньшей амплитуде вертикальных перемещений (i = 1, 2) на поверхность поступали породы, находившееся выше глубины подошвы области чисто упругого деформирования (He). Поэтому в коре в этих случаях не наблюдается избыточное горизонтальное сжатие, превышающее литостатическое давление.

Таким образом, для формирования в верхней части коры напряженного состояния горизонтального сжатия необходимо, чтобы амплитуда “откапывания” породы была выше некоторого значения (10).

Выводы

Выполненный в этом разделе анализ показал, что экзогенные процессы, происходящие в областях формирующихся горных поднятий, способны внести существенный вклад в изменение напряженного состояния. Так, согласно (10) уровень дополнительных напряжений сжатия, возникающий при “откапывании” пород, прямо пропорционален прочности сцепления и линейно возрастает от разницы между наибольшей глубиной, на которой находились породы, и глубиной начала катакластической текучести пород от действия МСТ. Например, при глубине эксгумации породы осадочного бассейна в 5 км дополнительное сжатие на поверхности составит около 10 МПа (прочность сцепления около 3 МПа и ${{k}_{f}} = 0.6$). Для пород кристаллического щита, например Балтийский щит, у которого согласно оценкам, выполненным в работах [Евзеров, 2001; Сим, 2012; Сидоренко, 1958], денудация за последние 150 млн лет могла составить 10 км, дополнительные сжимающие напряжения у поверхности могут быть существенно больше и достигать уровня 70 МПа (${{\tau }_{f}} = 10$ МПа , ${{H}_{p}} = 1$ км).

Рассмотренный механизм генерации напряжений связан с проявлением остаточных напряжений первого рода по классификации Н.Н. Давиденкова [1936], которые уравновешиваются в пределах образца или структурного блока (мегаскопический масштаб осреднения).

Важно отметить, что если полагать связь эксгумации пород с вертикальными движениями, возникающими на подошве коры из-за процесса маломасштабной астеносферной конвекции [Михайлов, 1983; 1999; Мягков, Ребецкий, 2019], то с глубиной уровень дополнительного сжатия (10) должен падать, снижаясь до нуля вблизи подошвы коры. С другой стороны, существуют представления о внутрикоровых течениях на уровне средней коры, сопрягающих процессы погружения и осадконакопления во впадинах, и рост горных поднятий [Макаров, 2010; 2011].

Следует напомнить, что после работ Б. Войгхта [Voigth, 1966; Voigth, St Pierre, 1974], в которых было предложено объяснение механизма формирования избыточных напряжений горизонтального сжатия, вышла серия работ Д.Л. Туркотта и его коллег [Turcotte, 1973; 1974; Haxby, Turcotte, 1976], в которых была предпринята попытка опровергнуть эти идеи. Увлеченность идеями простого объяснения широкого спектра явлений с позиции движения литосферных плит, научный вес Д.Л. Турркота – все это привело к тому, что в Геодинамике эти работы Б. Войгхта были забыты и больше не развивались.

С другой стороны, в горном деле, где практическая необходимость обеспечения безопасности работы, отвергает возможность научного протекционизма, эти идеи были восприняты и использовались [Джагер, 1975; Гудман, 1987].

ТЕКТОНОФИЗИЧЕСКИЕ АСПЕКТЫ ФОРМИРОВАНИЯ МАГМОПРОВОДЯЩИХ РАЗЛОМОВ В ПРОЦЕССЕ МАГМАТИЗМА И ВУЛКАНИЗМА КОРЫ

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

Три базовых положения

Математические и физические модели формирования даек и силлов [Pollard et al., 1975; Pollard, Segall, 1987; Rubin, 1995; Gudmundsson, 2002, 2006; 2011; Kavanagh et al., 2006; Rivalta et al., 2015] так же, как и механизм возникновения магмопроводящих коровых разломов, базируются на требовании превышения давления магмы ($P_{m}^{{}}$) в ее фронте и внутри трещины – разлома над обжимающими напряжениями в окружающем массиве. На рис. 2 показана схема взаимоотношения обжимающих напряжений ${{\sigma }_{{xx}}}$ и давления магмы $P_{m}^{{}}$ для вертикальной трещины, развивающейся из магматического очага с внутренним давлением $P_{m}^{o}$. В самом первом приближении условие продвижения магмы по разлому может быть представлено в виде:

(12)
$P_{m}^{{}} \geqslant \sigma _{{хх}}^{{}}.$
Рис. 2.

Схема внутреннего давления в вертикальном магматическом канале – разломе, выходящем из магматического очага.

Вторым физическим положением, необходимым для анализа возможности продвижения магмы по разрывам и трещинам, является закон распределения давления магмы $P_{m}^{{}}$ в магмопроводящем разломе. Наиболее просто предположить гидростатический закон его распределения по глубине:

(13)
$P_{m}^{{}} = P_{m}^{o} - hg~{{\rho }_{m}}.$

Здесь: $h$ – вертикальная дистанция от кровли магматического очага; g – ускорение силы тяжести; ${{\rho }_{m}}$ – плотность самой магмы; $P_{m}^{o}$ – давление магмы в подкоровом очаге. Согласно (13) давление магмы в разломе на каждом глубинном уровне снижается от давления в магматическом очаге на величину веса нижележащего столба магмы.

Это соотношение отвечает стационарной стадии процесса проникновения магмы по разлому, когда движение магмы либо уже отсутствует, либо имеет квазипостоянную скорость. Для нестационарной стадии, когда идет продвижение магмы по разлому вверх, ее давление во фронте ниже гидростатического закона, представленного в (13) [Rubin, 1995]. Стационарная стадия возникнет только после выхода магмы из разлома на поверхность или после прекращения движения магмы.

Еще одно физическое положение, которое мы будем использовать в тектонофизическом анализе, связано с определением уровня давления $P_{m}^{o}$ в мантийном магматическом очаге, расположенном под подошвой коры мощностью ${{H}_{c}}$. Будем считать, что этот очаг неизометричен и представляет собой достаточно обширное по латерали тело с мощностью, много меньшей его горизонтальных размеров. В этом случае можно считать, что давление расплавленной магмы в нем близко к весу столба вышележащих пород, т.е. к литостатическому давлению пород на подошве коры ($~{{p}_{{lt}}}\left( {{{H}_{c}}} \right)$. При этом допустимы малые вариации давления.

(14)
$P_{m}^{o} = ~\,\,{{p}_{{lt}}}\left( {{{H}_{c}}} \right) = ~\,\,{{H}_{c}}~{{\rho }_{c}}g.$

Здесь ${{\rho }_{с}}$ – средняя плотность пород коры.

Геомеханические параметры магматического процесса в коре

Используя представленные выше базисные положения, далее выполним анализ эволюции напряженного состояния коры в ходе проникновения магмы вверх по магмопроводящему субвертикальному разлому. В рассмотренных в этом разделе тестовых моделях будем считать, что начальный магматический очаг представлял собой большое по площади субпластовое мантийное тело, расположенное под подошвой коры на глубине ${{H}_{c}} = 40$ км. Наши расчеты будут относиться к установившейся стадии процесса, когда скорость движения магмы постоянная.

Рассмотрены два варианта начального геодинамического режима коры, один из которых отвечает напряженному состоянию горизонтального растяжения, а другой – горизонтального сжатия. Первый тип режима можно связывать с действием только МСТ с учетом того, что породы с некоторых глубин должны достигать критического катакластического состояния (см. предыдущий раздел). Второй тип геодинамического режима можно связывать с дальнодействующим влиянием напряжений горизонтального сжатия от границ сталкивающихся плит [Molnar, Tapponnier, 1975], но он также может быть объяснен в рамках концепции, изложенной в предыдущем разделе.

Плотности пород в модели принимались одинаковыми по всей толщине коры ${{\rho }_{c}} = 2.8$ г/см3. Плотность магмы задавалась в двух вариантах: большее ${{\rho }_{m}} = 2.85$ г/см3 и меньшее ${{\rho }_{m}} = 2.75$ г/см3 плотности пород коры. Также варьировалось давление магмы в мантийном очаге $P_{m}^{o}$ от меньшего до большего, чем литостатическое давление $p_{{lt}}^{{}}\left( {{{H}_{c}}} \right)$ (±2%).

Считаем, что начальному состоянию до внедрения магмы в разлом отвечает близость значений двух главных напряжений, действующих в горизонтальном направлении $\sigma _{{xx}}^{0} \approx \sigma _{{yy}}^{0}$ ($\sigma _{{xx}}^{0}$ ортогонально, а $\sigma _{{yy}}^{0}$ параллельно простиранию разлома). Для вертикальных напряжений ${{\sigma }_{{zz}}}$ верхний индекс не ставится, т.к. считается, что эти напряжения не изменяются в ходе изучаемого процесса. На рис. 4 начальные значения горизонтального сжатия определяли средний вертикальный градиент 0.85${{{{\rho }}}_{{\text{c}}}}{\text{g}}$. Он связан с достижением породы после глубины ${{H}_{p}} = $ 520 м закритического состояния (см. выражение (5)) при параметрах прочности ${{k}_{f}} = 0.6$ (коэффициент внутреннего трения), ${{\tau }_{f}} = 6\,\,~{\text{МПа}}$. В расчетах предельных напряжений учитывалось влияние давления флюида $({{p}_{{fl}}})$ в трещинно-поровом пространстве при ${{p}_{{fl}}} = 0.6~{{p}_{{lt}}}$.

Рис. 3.

Схемы распределения конечного напряженного состояния коры при одинаковом начальном напряженном состоянии, отвечающем режиму горизонтального растяжения: (a) $P_{m}^{o} < p_{{lt}}^{{}}\left( {{{H}_{c}}} \right)$, (б) $P_{m}^{o} > p_{{lt}}^{{}}\left( {{{H}_{c}}} \right)$; I – ${{\rho }_{m}}~\,\, < ~\,\,{{\rho }_{c}}$, II – ${{\rho }_{m}} > ~\,\,{{\rho }_{c}}$. По вертикали – глубина в километрах, по горизонтали напряжения и давление в ГПа. $\sigma _{{xx}}^{0} = \sigma _{{уу}}^{0}$ и $\sigma _{{xx}}^{1} \ne \sigma _{{уу}}^{1}$ – горизонтальные напряжения начального и конечного напряженных состояний, ${{\sigma }_{{zz}}} = - {{p}_{{lt}}}$ – вертикальные напряжения начального и конечного напряженных состояний, $\sigma _{{xx}}^{e} = \sigma _{{уу}}^{e}$ – напряжения для чисто упругого состояния . Точка D определяет пересечение линии магматического давления в разломе ${{p}_{m}}$ с линией начального значения горизонтальных сжимающих напряжений $\sigma _{{xx}}^{0}$, точка S – определяет пересечение линии магматического давления в разломе ${{P}_{m}}$ с линией вертикальных напряжений ${{\sigma }_{{zz}}}$, а точки A и V на поверхности обозначают тип лавовых излияний – ареальный или с созданием вулканической постройки.

Рис. 4.

Схемы распределения начальных и конечных напряжений в коре при наличии режима горизонтального сжатия в верхней ее части и горизонтального растяжения в нижней: (a) $P_{m}^{o} = p_{{lt}}^{{}}\left( {{{H}_{c}}} \right)$, (б) $P_{m}^{o} > p_{{lt}}^{{}}\left( {{{H}_{c}}} \right)$; I – ${{\rho }_{m}}~\,\, < \,\,~{{\rho }_{c}}$, II – ${{\rho }_{m}} > ~\,\,{{\rho }_{c}}$. Горизонтальный пунктир разделяет кору на режимы горизонтального сжатия и растяжения, действовавшие до начала магматизма. См. подпись к рис. 3. Серые эллипсы показывают глубины возможного формирования горизонтальных интрузий или внутрикоровых очагов. Пояснения в тексте.

Режим горизонтального растяжения в коре

В случае, когда в начальном состоянии вертикальное сжатие в коре всюду больше горизонтального и когда плотность магмы меньше плотности пород коры (рис. 3I), давление магмы $P_{m}^{{}}$ в магмопроводящем разломе везде выше горизонтального сжатия, что и позволяет ей достигнуть поверхности. При этом давление магмы повышает уровень горизонтального сжатия вблизи разлома для конечного напряженного состояния ($\sigma _{{xx}}^{1} = - {{P}_{m}}$). Для первого варианта давления магмы в очаге $P_{m}^{o}$ (рис. 3Iа) излияние лав на поверхности будет иметь ареальный тип, а во втором (рис. 3Iб) – должны формироваться вулканические постройки. Низкий уровень давления $P_{m}^{o}$ (рис. 3Iа) приводит к тому, что только в верхней части коры (выше точки S) давление магмы в разломе превышает литостатическое давление. В этом диапазоне глубин могут возникать горизонтальные интрузии (силлы). Такое напряженное состояние для всей коры также возникает и во втором случае (рис. 3Iб).

Другому случаю, для которого плотность магмы больше плотности пород коры (рис. 3II), отвечает и больший вертикальный градиент падения давления магмы в разломе. Это приводит к тому, что доставка магмы к поверхности возможна только при повышенном давлении магмы в мантийном очаге (рис. 3IIб). При давлении магмы в очаге $P_{m}^{o}$, меньшем литостатического давления (рис. 3IIа), в верхней части коры (глубины 3–5 км – точка D) происходит закрытие разлома, что не дает возможности дальнейшего продвижения фронта магмы вверх. Еще большее повышение плотности магмы будет приводить к “запиранию” разлома в диапазоне глубин средней коры. Таким образом, распространение магмы вверх по субвертикальному магмопроводящему разлому приводит к повышению уровня горизонтального сжатия, действующих в ортогональном к разлому направлении. Можно ожидать, что по латерали зона измененного напряженного состояния будет иметь линейный размер порядка 1–2 мощностей коры повышенного уровня горизонтальных напряжений.

Повышение горизонтального сжатия в направлении, нормальном к разлому, согласно условию отсутствия общего горизонтального деформирования объема пород (1) приведет также к росту горизонтальных сжимающих напряжений, действующих параллельно ему. Но увеличение сжатия этих напряжений будет в 4 раза меньшим (для упругого деформирования при коэффициенте Пуассона 0.25). Это означает, что конечное напряженное состояние в зависимости от варьируемых в представленных расчетах параметров (${{\rho }_{m}}~\,\,{\text{и\;}}P_{m}^{o}$) может остаться в режиме горизонтального растяжения (рис. 3Ia ниже точки S, рис. 3IIа) или перейти в режим горизонтального сдвига (рис. 3Ia выше точки S, рис. 3Iб, рис. 3IIб).

Во всех этих случаях напряжения наибольшего горизонтального сжатия на участке продвижения магмы будут действовать ортогонально простиранию разлома. Такое их положение определяет возможность формирования вертикальных даек (рис. 3Ia ниже точки S, рис. 3IIа), ответвляющихся от разлома, с простиранием, ортогональным разлому. Обращаем внимание, что в случае, рассмотренном на рис. 3IIб, уровень давления магмы в разломе на всем его протяжении до точки “запирания” D ниже литостатического давления. Это означает, что для всего диапазона глубин проникновения магмы в разлом не возникают условия для формирования горизонтальных интрузий.

Сочетания режимов горизонтального сжатия и растяжения в коре

Считаем, что в коре вертикальный градиент горизонтальных напряжений такой же, как и на рис. 3, при этом существует дополнительное боковое обжатие, одинаковое по всей мощности коры. Для представленных на рис. 4 случаев это дополнительное сжатие составляет 82 МПа. В рассмотренных примерах режим горизонтального сжатия наблюдается только в верхней части коры. В нижней части действует режим горизонтального растяжения. Здесь так же, как и на рис. 3, варьировались давление магмы в мантийном очаге и плотность магмы в разломе.

В рассмотренных случаях магма имеет возможность подняться до определенных глубин коры (точка D). При этом большее давление магмы в мантийном очаге, а также меньшая плотность магмы приводит и к большему проникновению вверх по разрезу магмы в разломе. Для магмы меньшей плотности ее проникновение в разлом происходит до глубин верхней коры (6 км), а для высокоплотной магмы “запирание” разлома может происходить в средней коре.

Так же, как и в предыдущем расчете, в диапазоне глубин проникновения магмы в разлом (ниже точки D) в коре возрастает интенсивность напряжений горизонтального сжатия, нормального к простиранию разлома, до значений давления магмы. Для случая, представленного на рис. 4IIа, давление магмы везде в разломе остается меньше литостатического давления пород. Поэтому здесь не изменяется геодинамический режим начального напряженного состояния, т.е. на глубинах ниже точки D остается горизонтальное растяжение.

В случаях, представленных на рис. 4Iа, 4Iб и рис. 4IIб, давление магмы в разломе ниже точки D больше, чем литостатическое давление. Это означает, что новый уровень горизонтального сжатия здесь будет больше вертикального сжатия, т.е. на этих глубинах возникает режим горизонтального сдвига. Вследствие такого изменения напряжений от подошвы коры и до точки D от магмопродводящего разлома могут возникать горизонтальные интрузии. По сути, эти горизонтальные интрузии можно рассматривать как внутрикоровые промежуточные магматические очаги.

Таким образом, распространение магмы вверх по субвертикальному магмопроводящему разлому приводит к повышению уровня напряжений горизонтального сжатия, действующих в ортогональном к разлому направлении. Это означает, что конечное напряженное состояние ниже точки D в зависимости от варьируемых в представленных расчетах параметров (${{\rho }_{m}}{\text{\;и\;}}P_{m}^{o}$) может остаться в режиме горизонтального растяжения или перейти в режим горизонтального сдвига.

Выводы

Таким образом, после магматического излияния на поверхности и остывания магмы в разломе в его окрестности устанавливается новый режим напряженного состояния, при котором главное напряжение наибольшего сжатия будет действовать уже в горизонтальном направлении оротогонально разлому. Поскольку при этом второе из напряжений, действовавшее в горизонтальном направлении ($\sigma _{{xx}}^{1} \approx \sigma _{{xx}}^{0}$), изменяется мало, то произойдет изменение геодинамического режима горизонтального растяжения на режим горизонтального сдвига (промежуточным главным напряжением становится вертикальное напряжение).

Существует устойчивое представление о движущей силе поднятия магмы вверх по разлому-каналу, связанному с меньшей ее плотностью, чем плотность окружающих пород. При большей плотности магмы происходит более быстрое снижения ее давления в канал при уменьшении глубины. Как это видно из рисунков 3II, магма продвигается вверх и при плотности, большей чем плотность окружающей породы. При этом практически на всей дистанции давление магмы меньше литостатического. Главным здесь является то, что уровень давления в канале больше, чем уровень напряжений горизонтального сжатия. Поэтому при невысоком давлении магмы во внутрикоровом магматическом очаге магма не может дойти до поверхности и на глубине около 1 км происходит закрытие разлома. Поскольку в точке D (рис. 3IIб) давление магмы меньше вертикальных (литостатика) и горизонтальных напряжений, то это означает, что канал остается “слепым”, здесь не образуется плутонического магматического очага.

Для всех рассмотренных вариантов расчета для глубин распространения магмы вверх по разлому-каналу в окружающем его массиве резко повышается уровень горизонтального сжатия, действующего по нормали к разлому, и изменяется тип напряженного состояния на горизонтальный сдвиг. То есть, даже в том случае, когда проявления вулканизма на поверхности нет (рис. 3IIa, 3IIб, рис. 4), все равно эффект повышения напряжений горизонтального сжатия в большей части глубин коры вблизи магмопроводящего разлома имеется. В этой связи далее в работе будем говорить не только об активизации вулканизма, а более широко понимать проблему, определяя ее как активизацию магматизма в коре.

На рис. 3 и рис. 4 рассматривались случаи продвижения магмы по протяженному разлому. Однако он может также соответствовать и случаю, когда основные каналы поступления магмы изометричные (цилиндрические) в плане, но при этом активизируется группа вулканов, расположенных на одном магмопроводящем разломе. Это достаточно типичная ситуация для вулканизма зон субдукции.

Если следовать положениям теории упругости, то можно ожидать, что по латерали зона измененного напряженного состояния (смена горизонтального растяжения на горизонтальный сдвиг) будет иметь линейный размер, в 2–4 раза превышающий горизонтальную протяженность магмопроводящего разлома. В случае продвижения магмы вверх от очага по единичному изометричному каналу (цилиндрический) повысится уровень обоих главных напряжений, действовавших в горизонтальном направлении и таким образом в конечном состоянии будет режим горизонтального сжатия. При этом в области измененного напряженного состояния будет небольшой, около 3–5 радиусов вулканической постройки (10–50 км).

Выполненный анализ позволяет утверждать, что в случае наличия магматических очагов в нижней части коры или литосферы для любого случая распределения напряженных состояний по глубине коры (рис. 3, рис. 4) должно происходить продвижение магмы вверх по субвертикальным разломам, в результате которого та глубинная часть коры, которая находилась в состоянии горизонтального растяжения, переходит в состояние горизонтального сдвига или сжатия.

В данной работе рассматривался механизм проникновения магмы по субвертикальным разломам. Однако он распространяется и на случай проникновения магмы вверх по системам мелких субвертикальных трещин, что выглядит как явление перколяции – просачивания. Движущей силой здесь является давление магмы, превышающее уровень напряжений горизонтального сжатия. Необходимым условием перколяционного механизма проникновения магмы вверх по разрезу является насыщение магмой слоя пород в нижней коре (магматический очаг) до критических значений для подержания в ней уровня давления не ниже литостатического давления вышележащих пород.

ОБСУЖДЕНИЕ

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

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

В то же время, увеличение напряжений горизонтального сжатия, возникающее за счет магматизма и вулканизма, проникающих в кору по субвертикальным протяженным магмопроводящим разломам, всегда неравномерно распределено в латеральном направлении. Оно существенно больше в направлении, ортогональном разломам, по которым магма поднимается вверх. Изменения уровня горизонтального сжатия должны захватывать обязательно нижнюю и среднюю кору, и только в случае излияния магмы на поверхность (ареальные или вулканические) повышенное надлитостатическое горизонтальное сжатие устанавливается во всей коре. При этом превышение горизонтального сжатия над вертикальными напряжениями определяется давлением столба магмы в вулкане. Важно отметить, что в результате действия магматизма в коре возникает геодинамический тип напряженного состояния горизонтальный сдвиг, если в коре начальное состояние определялось только массовыми силами тяжести.

Можно вполне обоснованно предположить, что появление магматизма в коре обязательно должно будет сопровождаться формированием небольшого горного рельефа, поскольку в коре будут возникать магматические интрузии – промежуточные внутрикоровые очаги. Следствием формирования рельефа будет увеличение интенсивности экзогенных процессов и, таким образом, будут одновременно действовать оба механизма генерации напряжений повышенного горизонтального сжатия. Поскольку магматизм увеличивает горизонтальное сжатие приблизительно до величины литостатического давления (рис. 3, рис. 4), то для этого случая можно получить выражение для напряжений дополнительного горизонтального сжатия действующих ортогонально разлому:

(15)
$d\sigma = - {{\rho }_{с}}g\frac{{1 - 2\nu }}{{1 - \nu }}H.$

Оно зависит только от глубины $H$, с которой произведено перемещение к поверхности породы. Так же, как и в выражении (8), напряжение $d\sigma $ это дополнительное горизонтальное сжатие относительно упругого состояния (2).

Из выражения (15) следует, что взаимосвязное действие обоих механизмов приводит к тому, что уровень горизонтального сжатия в верхних горизонтах коры становится существенно большим в сравнении с действием только денудационного механизма. Важно отметить, что взаимодействие двух механизмов приводит к тому, что в коре уровень латерального сжатия будет зависеть от направления. Наибольшее сжатие будет действовать нормально к магмопроводящим вертикальным структурам земной коры (разломы, сутуры, линеаменты цепочек вулканов).

Рассмотрим теперь возможное влияние магматизма коры и денудационных процессов в различных их проявлениях (по отдельности или совместно) на напряженное состояние различных крупных геологических объектов.

Для внутриплатформенных орогенов альпийского типа с высокой интенсивностью эрозии склонов, в которых в период их формирования отмечается вулканизм (например, Кавказа) или плутонический магматизм (например, Тянь-Шань), объяснение наличия в них повышенного уровня напряжений горизонтального сжатия с позиций двух рассматриваемых в статье механизмов генерации выглядит вполне обоснованным. То же самое можно сказать об островных дугах зон субдукции, т.к. здесь активно идут как процессы денудации [Ермаков, 2005], так и процессы вулканизма.

В центральных частях горных поднятий типа плато, например Тибет, где скорость эрозии очень низкая, а время их жизни менее нескольких десятков миллионов лет, процесс денудации поверхности не позволяет вывести вверх породы с больших глубин. Поэтому здесь вклад в формирования повышенного горизонтального сжатия остаточных напряжений, вызванных эксгумацией пород, небольшой. Но здесь должен действовать механизм генерации напряжений, связанный с процессами магматизма и вулканизма, которые, в частности, для Тибета, имели несколько волн, концентрируясь в зонах террейнов [Chung et al., 2005]. Такая линейность зон магматизма наиболее способствует повышению уровня напряжений горизонтального сжатия именно в направлении нормальном к простиранию этих зон. Соответственно основным режимом напряженного состояния, который здесь должен возникнуть, является горизонтальный сдвиг в участках повышенного вулканизма и горизонтальное растяжение там, где вулканизм проявлен менее интенсивно, что и наблюдается [Rebetsky, Alekseev, 2014; Ребецкий, Алексеев, 2022].

Для щитов платформ, находящихся в несколько приподнятом состоянии на протяжении 2–3 сотен миллионов лет, например Балтийский щит, и где по данным работ [Евзеров, 2001; Сим, 2012; Сидоренко, 1958] имели место большие амплитуды денудации (до 10–15 км), действие механизма “откапывания” породы вполне способно объяснить высокий уровень напряжений горизонтального сжатия в верхней коре.

ЗАКЛЮЧЕНИЕ

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

Наибольший эффект в повышении уровня напряжений горизонтального сжатия должен наблюдаться при совместном действии этих механизмов, имеющее место для коры островных дуг зон субдукции и континентальных орогенов, в которых проявляется вулканическая деятельность. Даже при отсутствии проявления вулканизма наличие в коре процесса магматизма на стадиях, предшествовавших орогенезу или на стадии орогенеза, также способствует появлению повышенного уровня напряжений горизонтального сжатия относительно вертикальных напряжений в верхней части коры. При этом механизм, связанный с “откапыванием” пород и формированием остаточных напряжений, дает наибольший вклад в повышение уровня горизонтального сжатия именно в верхней коре, существенно снижаясь в средней и особенно в нижней коре. В то же время, механизм, обусловленный магматизмом, дает наибольшее увеличение горизонтального сжатия в нижней коре и несколько меньший вклад в увеличение этих напряжений в верхней коре.

Список литературы

  1. Авсюк Ю.Н. Внеземные движущие силы тектоники. Основные проблемы глобальной тектоники. М.: Научный мир. 2001. С. 437–441.

  2. Викулин А.В. Энергия и момент силы упругого ротационного поля // Геология и геофизика. 2008. Т. 49. № 6. С. 559–570.

  3. Вихри в геологических процессах / А.В. Викулин (ред.). Петропавловск-Камчатский: КГПИ. 2004. 297 с.

  4. Волох Н.П., Сашурин А.Д., Липин Я.И. Исследования остаточных напряжений в крепких горных породах. Современные проблемы механики горных пород. Л.: Наука. 1972.

  5. Давиденков Н.Н. Об остаточных напряжениях. Рентография в применении к исследованиям материалов. М.-Л.: ОНТИ. 1936.

  6. Гудман Р. Механика скальных пород. 1987. М.: Строиздат. 232 с.

  7. Джагер Ч. Механика горных пород и инженерные сооружения. М.: Мир. 1975. 255 с.

  8. Динник А.Н. О давлении горных пород и расчете крепи круглой шахты // Инженерный работник. 1926. № 3. С. 1–12.

  9. Добрецов Н.Л., Кирдяшкин А.Г., Кирдяшкин А.А. Глубинная геодинамика. 2-е изд. Новосибирск: изд-во СО РАН, филиал Гео. 2001. 409 с.

  10. Добрецов Н.Л., Чумаков И.М. Глобальные периодичности в эволюции литосферы и биосферы. Глобальные изменения природной среды и эволюция геопроцессов. Новосибирск: изд-во СО РАН. 2001. С. 11–26.

  11. Евзеров В.Я. Россыпные месторождения – уникальные образования в рыхлом покрове Балтийского щита // Литология и полезные ископаемые. 2001. № 2. С. 126–133.

  12. Ермаков В.А. Островные дуги и их роль в эволюции континентальной окраины (новый взгляд на известные факты) // Вулканология и сейсмология. 2005. № 5. С. 3−–18.

  13. Зоненшайн Л.П., Савостин Л.А., Мишарина Л.А., Солоненко Н.В. Тектоника плит Байкальской горной области и Станового хребта // Докл. АН СССР. 1978. Т. 240. № 3. С. 669–672.

  14. Кейлис-Борок В.И. Исследование источников, приближенно эквивалентных очагам землетрясений // Труды Геофиз. Ин-та АН СССР. 1950. № 9.

  15. Лобковский Л.И., Котелкин В.Д. Двухярусная термохимическая конвекция в мантии и ее геодинамические следствия // Проблемы глобальной геодинамики. М.: Геос. 2000. С. 29–53.

  16. Ма Си Юань Тектонические процессы, отраженные на карте динамики литосферы Китая. Геодинамика внутриконтинентальных горных областей / Н.А. Логачев (ред.). Новосибирск: Наука. 1990. С. 341−351.

  17. Макаров В.И., Алексеев Д.В., Баталев В.Ю., Баталева Е.А., Беляев И.В., Брагин В.Д. и др. Поддвиг Тарима под Тянь-Шань и глубинная структура зоны их сочленения: основные результаты сейсмических исследований по профилю MANAS (Кашгар Сонкель) // Геотектоника. 2010. № 2. С. 23−42.

  18. Макаров В.И., Рыбин А.К., Матюков В.Е., Пушкарев П.Ю., Щербина Ф.А. Особенности глубинной структуры депрессионных областей Центрального Тянь-Шаня // Инженерные изыскания. 2011. № 1. С. 42−51.

  19. Марков Г.А. О распространении горизонтальных тектонических напряжений в зонах поднятий земной коры // Инженерная геология. 1980. № 1. С. 20−30.

  20. Марков Г.А. Тектонические напряжения и горное давления в рудниках Хибинского массива. Л.: Наука. Ленинград. Отд. 1977. 211 с.

  21. Михайлов В.О. Математическая модель эволюции структур, образующихся в результате вертикальных движений // Изв. АН СССР. Сер. Физика Земли. 1983. № 6. С. 3−18.

  22. Михайлов В.О. Моделирование процессов растяжения и сжатия литосферы внутриплитными силами // Физика Земли. 1999. № 3. С. 71−81.

  23. Мягков Д.С., Ребецкий Ю.Л. Математическая модель формирования напряженно-деформированного состояния эпиплатформенных орогенов // Геодинамика и тектонофизика. 2019. Т. 10. № 1. С. 21−41. https://doi.org/10.5800/GT-2019-10-1-040

  24. Николаевский В.Н. Граница Мохоровичича как предельная глубина хрупко-дилатансионного состояния горных пород // Докл. АН СССР. 1979. Т 249. № 4. С. 817−820.

  25. Пономарев В.С. Потенциальная энергия гидростатического сжатия горных пород и сейсмичность. Труды Третьего Всесоюзного симп. по сейсмич. режиму. Новосибирск: Наука. 1969. Ч. 1. С. 30−55.

  26. Пономарев В.С. Энергонасыщенность геологической среды. М.: Наука. 2008. 378 с.

  27. Ребецкий Ю.Л. Механизм генерации остаточных напряжений и больших горизонтальных сжимающих напряжений в земной коре внутриплитовых орогенов. Проблемы тектонофизики. К 40-летию создания М.В. Гзовским лаборатории тектонофизики в ИФЗ РАН. М.: Изд. ИФЗ РАН. 2008а. С. 431−466.

  28. Ребецкий Ю.Л. Механизм генерации тектонических напряжений в областях больших вертикальных движений // Физическая мезомеханика. 2008б. Т. 11. № 1. С. 66−73.

  29. Ребецкий Ю.Л. О возможном механизме генерации в земной коре горизонтальных сжимающих напряжений // Докл. РАН. 2008в. Т. 423. № 4. С. 538−542.

  30. Ребецкий Ю.Л. Тектонические напряжения и прочность горных массивов. М.: Академкнига. 2007. 406 с.

  31. Ребецкий Ю.Л., Алексеев Р.С. Поле современных тектонических напряжений коры Высокой Азии и возможный механизм его формирования // Геология и геофизика. 2022. Т. 62. № 6. С. 875−897.

  32. Ребецкий Ю.Л., Полец А.Ю. Напряженное состояние в афтершоковой области землетрясения Тохоку 11.03.2011 // Вулканология и сейсмология. 2021. № 4. С. 22–44.

  33. Сидоренко А.В. Доледниковые коры выветривания Кольского полуострова. М.: изд-во АН СССР. 1958. 106 с.

  34. Cим Л.А. О связи объема осадконакопления в бассейнах обрамления с величиной эрозионного сноса с Фенноскандинавского щита в мезокайнозое. Осадочные бассейны и геологические предпосылки прогноза новых объектов, перспективных на нефть и газ. М-лы XLIV Тект.сов. М.:ГЕОС. 2012. С. 398−401.

  35. Ружич В.В. Сейсмотектоническая деструкция в земной коре Байкальской рифтовой зоны. Новосибирск: Изд-во СО РАН. 1997. 144 с.

  36. Ружич В.В., Кочарян Г.Г., Левина Е.А. Оценка геодинамического влияния зон коллизии и субдукции на сейсмотектонический режим Байкальского рифта // Геодинамика и тектонофизика. 2016. Т. 7. № 3. С. 383–406. https://doi.org/10.5800/GT-2016-7-3-0214

  37. Трубицын В.П. Реология мантии и динамическая теория тектоники литосферных плит // Физика Земли. 2012. № 6. С. 3–22.

  38. Трубицын В.П., Рыков В.В. Мантийная конвекция с плавающими континентами. Проблемы глобальной геодинамики. М.: Геос. 2000. С. 7–28.

  39. Bath M., Duda S.J. Earthquake volume, fault plane area, seismic energy, strain, deformation and related quantieties // Ann. Gefis (Rome). 1964. V. 17. P. 353−368.

  40. Dobretsov N.L. Periodicity and driving forces of volcanism // Russian Geology and Geophysics. 2015. V. 56. № 12. P. 1663–1670. https://doi.org/10.1016/j.rgg.2015.11.001

  41. Dziewonski A.M., Chou T.-A., Woodhouse J.H. Determination of earthquake source parameters from wave-form data for studies of global and regional seismicity // J. Geophys. Res. 1981.V. 86. P. 2825−2852.

  42. Glukhovskii M.Z. The rotation factor and some problems of geotectonics and comparative planetology // Geotectonics. 2005. V. 39. № 6. P. 423–436.

  43. Haxby W.F., Turcotte D.L. Stress induced by the addition or removal of overburden and associated thermal effect // Geology. 1976. № 4. P. 181−194.

  44. Herget G. Variation of rock stresses with depth at Canadion iron mine // Int J. Rock Mech. 1973. Min Sci. № 10. P. 37−51

  45. Jager J.C. Elasticity Fracture and Flow. London.: Methuen & Co. LTD. 1962. 208 p.

  46. Kanamory H., Anderson D.L. Theoretical basis of some empirical relations in seismology // Bull. Seismol. Soc. Am. 1975. V. 65. P. 1073−1095.

  47. Kavanagh J.L., Menand T., Sparks R.S.J. An experimental investigation of sill formation and propagation in layered elastic media // Earth and Planetary Science Letters. 2006. V. 245. P. 799−813.

  48. Molnar P., Tapponnier P. Cenozoic tectonics of Asia: effects of a continental collision // Science. 1975. V. 189. P. 419−426.

  49. Nazarova L.A., Nazarov L.A., Dyad’kov P.G. Mathematical modeling of kinematics of Central Asian plates // Journal of Mining Science. 2002. V. 38. № 5. P. 411–417. https://doi.org/10.1023/A:1023923329058

  50. Pollard D.D., Segall P. Theoretical displacements and stresses near fractures in rock: with applications to faults, joints, veins,dikes, and solution surfaces, surfaces / Atkinson B.K. (ed.). 1987. P. 277−350.

  51. Pollard D.D., Muller O.H., Dockstader D.R. The Form and Growth of Fingered Sheet Intrusions // Geological Society of America Bulletin. 1975. V. 86. № 3. P. 351.

  52. Price N.J., Cosrove J.W. Analysis of geological stryctyres. Cambridge Univ. Press. 1990. 495 p.

  53. Rivalta E., Taisne B., Bunger A.P., Katz R.F. A review of mechanical models of dike propagation: Schools of thought, results and future directions // Tectonophysics. 2015. V. 638. P. 1−42.

  54. Rebetsky Y.L., Alekseev R.S. The field of recent tectonic stresses in Central and South-Eastern Asia // Geodynamics & Tectonophysics. 2014. V. 5. № 1. P. 257−290. https://doi.org/10.5800/GT-2014-5-1-0127

  55. Rubin A.M. Propagation of magma-filled cracks // Annu. Rev. Earth Planet. 1995. V. 23. P. 287–336.

  56. Chung, Sun-Lin,, Chua Mei-Fei, Zhang Yuquan, Xie Yingwen, Lo Ching-Hua, Lee Tung-Yi, Lan Ching-Ying, Li Xianhua, Zhang Qi and Wang Yizhao Tibetan tectonic evolution inferred from spatial and temporal variations in post-collisional magmatism // Earth-Science Reviews. 2005. V. 68. № 3–4. P.173−196.

  57. Sykes L.R., Sbar M.L. Intraplate earthquakes, lithosphere stresses, and the driving mechanism plate tectonics // Nature. 1973. V. 245. P. 298–302.

  58. Turcotte D.L. Driving mechanisms for plate tectonics // Geofisica internac. 1973. V. 13. P. 309−315.

  59. Turcotte D.L. Membrane tectonics // Royal Astron. Soc. Geophys. J. 1974. V. 36. P. 33−42.

  60. Voigth B. Beziehugen zwischen grossen horirontalen Spannugen in Gebirgen und der Tektonik und der Abtragung // First Cong. Internat. Soc. For Rock Mechanics. Lisobon. Proc. 1966. V. 2. P. 51–56.

  61. Voigth B., St Pierre B.H.P. Stress history and rock stress // Third Cong. Internat. Soc. For Rock Mechanics. Denver. Proc. 1974. V 2. P. 580-582.

  62. Wilson J.T. Did the Atlantic close and then re-open? // Nature. 1966. V. 211. P. 676–681.

Дополнительные материалы отсутствуют.