


№ 4 (2023)
Статьи
Поправки к теории упругого изгиба тонких плит для 2D-моделей в приближении Рейснера
Аннотация
Расчеты упругих изгибных напряжений и деформаций в литосфере обычно проводятся на основе теории Кирхгофа–Лява для тонких плит. Критерием ее применимости является малость отношения толщины плиты к ее длине. В океанических плитах благодаря выталкивающей силе мантии основные деформации распределены не равномерно вдоль плиты, а сосредоточены вблизи зоны субдукции. Поэтому эффективная длина изгибающейся части плиты в несколько раз меньше реальной длины, и критерий тонкости плит частично нарушается. В работе анализируется возможность применения уравнений изгиба толстых плит. Имеющиеся вариационные теории 3D-изгиба толстых плит намного более сложные по сравнению с теорией Кирхгофа–Лява, так как требуют решения не одного дифференциального уравнения, а трех, и из-за трудоемкости имеют ограниченное применение. Поскольку в геофизических приложениях часто используются 2D-модели, то в работе детально анализируются возможности и точность теории изгиба толстых пластин для 2D-моделей. Оригинальные уравнения 3D-изгиба толстых плит Рейснера после перехода к 2D для плоской деформации и плоского напряжения выписываются в форме, аналогичной уравнениям Кирхгофа с аддитивными поправками, и дополняются явными выражениями для продольного смещения. Сравнение аналитических решений 2D-уравнений Рейснера с точными решениями показывает, что она дает поправку только для функции изгиба плиты. Но эта поправка уточняет теорию Кирхофа–Лява почти на порядок. При этом решение уравнений оказывается практически таким же простым, как и уравнений тонких плит.



Миграция аномальных смещений на флангах формирующегося очага – модель “качели”
Аннотация
Полученные в рамках геомеханического моделирования результаты по динамике напряженно деформированного состояния перед землетрясением М = 7.1 Риджкрест 2019 г. (Южная Калифорния) указывают на попеременное развитие максимальных смещений на оконечностях разрыва, что имитирует процесс “раскачивания” эпицентральной зоны данного землетрясения. На основе данных результатов с использованием имеющихся теоретических представлений о блоковой структуре земной коры с соответствующими межблоковыми связями, предложена геомеханическая модель “Качели”, основанная на явлении прерывистого скольжения (stick-slip) трех связанных между собой блоков. Параметры этой модели связаны с реологическими свойствами земной коры в области разрушения перемычки, определяющей размеры очага сейсмического события.



Сейсмотектонические изменения перед землетрясением Суматра 26.12. 2004 г., М = 9.1
Аннотация
Исследованы записи скорости вертикального перемещения земной поверхности возле широкополосных сейсмических станций, расположенных вокруг эпицентра землетрясения Суматра 26.12. 2004 г. с магнитудой 9.1. На протяжении 5 лет с 1996 по 2000 гг. ближайшая к эпицентру и удаленная от него на 1700 км станция COCO записывала устойчивый суточный ход сейсмического шума. В 2001 г. на ней появились ступенеобразные искажения уровня записи. Они продолжались вплоть до момента землетрясения. Обнаружены также превышающие суточный ход импульсы с постепенным нарастанием и последующим спадом амплитуды колебаний длительностью в несколько минут. Они возникали при спокойных метеорологических условиях и геомагнитной активности. На станциях, удаленных от эпицентра более чем на 2000 км, эти импульсы не проявились. Предполагается, что перед землетрясением происходили подвижки в геологических разломах литосферы Индийского океана.



Электросопротивление ε-Fe при высоких давлениях ступенчатого ударного сжатия
Аннотация
Проведены in situ измерения электросопротивление сжатых и разогретых в процессе ступенчатого ударного нагружения образцов железа с гексагональной плотноупакованной решеткой (ε-Fe). Построены уравнения состояния ε-Fe. Выполнено математическое моделирование полученных экспериментальных результатов в рамках гидрокода, опирающегося на разработанные уравнения состояния. Моделирование позволило реконструировать объемно-температурную зависимость удельного электросопротивления ε-Fe при давлениях ≈20–70 ГПа и температурах ≈750–950 К. Выполнены расчеты объемно-температурной зависимости коэффициента теплопроводности ε-Fe с использованием закона Видемана–Франца. Полученные результаты по электро- и теплопроводности ударно сжатого и разогретого ε-Fe сопоставлены с литературными экспериментальными и теоретическими данными для железа и кремнистого железа.



Спрединг в северной части Евразиатской литосферной плиты по данным GPS-измерений и его возможная геодинамическая интерпретация
Аннотация
На всей, от Атлантического до Тихого океана, территории Евразиатской литосферной плиты, рассматриваемой в статье, векторы скоростей горизонтальных в восточном направлении смещений GPS-станций образуют пологую, плавную, выпуклую к северу огромную дугу. В ее составе можно выделить три относительно нешироких дугообразных полосы. Самая южная из них и наиболее обеспеченная результатами измерений полоса (дугообразная полоса А в данной статье) подробно рассмотрена в статье [Шевченко и др., 2021]. Прямые, реальные результаты опубликованных геодезических измерений показали, что длина этой полосы увеличивается на 5–10 мм/год. Дуга удлиняется. В упомянутой статье были рассмотрены и отклонены пять вариантов интерпретации этого удлинения, исходящие из представления о внешнем воздействии на горные породы земной коры/литосферы названной дугообразной полосы. Наиболее вероятной причиной удлинения представляется внутренний процесс увеличения объема упомянутых пород полосы в результате привноса в них глубинными флюидами дополнительного минерального материала и его последующей кристаллизации. В данной статье приводятся результаты аналогичных геодезических измерений на территориях двух других дуг – В и С. Оказалось, что расположенная севернее других дугообразная полоса С удлиняется (и со сходной скоростью), как и находящаяся южнее других полоса А. Но расположенная между ними полоса В смещается в восточном направлении, не удлиняясь. Геодезические измерения в осевой части Срединно-Атлантического хребта, на о.Исландия позволяют предположить, что общее расширение Атлантического океана и соответствующее восточное смещение Евразиатской литосферной плиты связано c нагнетанием пластичных магматических (основного состава) клиньев в осевой, рифтовой зоне упомянутого хребта. На фоне этого общего восточного смещения литосферной плиты отдельные ее части (дугообразные полосы А и С маркируют эти части) смещаются в том же восточном направлении с нарастающей с З на В скоростью. Этот рост скорости и соответствующее удлинение двух перечисленных полос мы связываем с упомянутым привносом дополнительного минерального материала и его последующей кристаллизацией.



Анализ современных деформаций Северо-Западного Кавказа и Предкавказья по данным ГНСС
Аннотация
По данным о горизонтальных перемещениях ГНСС охарактеризовано поле современной деформации Северо-Западного Кавказа и Предкавказья. Показано, что ее скорости достаточны для того, чтобы активизация крупных разломов этого региона сопровождалась сейсмическими событиями с моментной магнитудой до 7.0. В то же время поднятие горно-складчатого сооружения со скоростью до 12 мм/год, развитие его контрастного, глубоко расчлененного рельефа необъяснимо только влиянием внешних горизонтальных напряжений. Эти процессы во многом обусловлены внутренними изостатическим силами.



Геоэлектрическая модель центральной части Северного Кавказа и его флюидонасыщение
Аннотация
В пределах Северо-Кавказского региона выполнен комплекс магнитотеллурических и сейсмических исследований на профилях протяженностью более двух тысяч километров. Проведенная ранее интерпретация магнитотеллурических наблюдений с помощью программ одномерной, двумерной инверсий и трехмерного математического моделирования, позволила построить серию разрезов и моделей, которые рассматриваются как тестовые и стартовые при построении трехмерной геоэлектрической модели региона. На тестовых моделях опробовались возможности программы трехмерной инверсии компонент тензора импеданса в методе магнитотеллурического зондирования оценивать параметры проводящих блоков в структурах Большого Кавказа и Скифской плиты. В результирующей геоэлектрической модели, построенной по результатам трехмерной инверсии всех компонент тензора импеданса, положение низкоомных блоков коррелирует с глубинными разломами, вулканическими камерами различного генезиса, сейсмически активными зонами, характеризующимися пониженной скоростью сейсмических волн и их повышенным поглощением. Удельное электрическое сопротивление низкоомных аномалий объясняется степенью их насыщения водной фракцией флюида. Максимальная его концентрация сосредоточена в пределах пересечений систем разломов, флексурно-сдвиговых зон и в активизированных тектоническими процессами глубинных разломах.



Картирование магматических комплексов по данным гидромагнитных съемок Баренцевоморского региона
Аннотация
Для изучения магматических комплексов в осадочном чехле Баренцевоморского региона проанализированы практически все доступные на сегодняшний день материалы гидромагнитных съемок исследуемого региона. С использованием современных методов интерпретации авторами были обработаны результаты наблюдений, проведенных в 1995–2017 гг. на профилях общей протяженностью более 93 тыс. километров. Установлено, что в центральной и северной частях исследуемого региона широко распространены рои даек северо-западного простирания; в южной части дайки не выявлены. Показано, что магматические образования в осадочном чехле фиксируются локальными магнитными аномалиями двух типов – линейными, источниками которых являются дайки, и мозаичными, связываемыми с локальными очагами магматизма. При глубине моря 100–500 м и мощности осадочного чехла более 8 км преимущественная глубина залегания верхних кромок тел, создающих аномалии, составляет около 1500 м. Анализ аномального магнитного поля позволил уточнить строение осадочного чехла, историю протекания магматических процессов и тектонического развития Баренцевоморского региона.



Термическая история литосферы Колтогорско-Уренгойского Грабена Западно-Сибирского бассейна в районе СКВ. СГ-6 – численная реконструкция в рамках системы моделирования плоских бассейнов GALO
Аннотация
Система моделирования бассейнов GALO используется для численной реконструкции теплового режима литосферы Западно-Сибирского бассейна в Колтогорско-Уренгойском грабене в районе сверхглубокой скважины Тюменская СГ-6. Реконструкция объясняет особенности в формировании теплового режима бассейна, которые не рассматривались в предшествующих реконструкциях района. К таким особенностям относятся аномально высокий рост зрелости органического вещества в породах юры и триаса, высокие градиенты температур, наблюдаемые в верхнем фундаменте и триасо-пермском осадочном комплексе, аномально низкие температуры пород верхних горизонтов современного осадочного разреза бассейна. Анализ временных изменений тектонического погружения бассейна используется для оценки интенсивности и продолжительности событий тепловой активации и растяжения его литосферы. Тепловое влияние силла, внедрившегося в нижней юре в верхние горизонты фундамента, объясняло высокую степень созревания органического вещества в нижнетриасовых породах. Учет резких колебаний климата в плиоцен-четвертичное время вместе с гидротермальной активностью в низах осадочного чехла в верхнем плиоцене – нижнем плейстоцене, позволили получить глубинные профили температур и отражательной способности витринита, хорошо согласующиеся с измеренными значениями.



Изменения теплопроводности пород литосферы Западно-Сибирского бассейна в районе СКВ. Тюменская СГ-6
Аннотация
Численные реконструкции термического режима литосферы Западно-Сибирского бассейна в Колтогорско-Уренгойском грабене в районе сверхглубокой скважины Тюменская СГ-6 используются для анализа распределения теплопроводности пород бассейна с глубиной. Выделяются пять глубинных интервалов, отличающихся характером изменения теплопроводности пород с глубиной: зона пермафроста, осадочный разрез ниже этой зоны, зона аномального разуплотнения пород, консолидированная кора и мантия. Рассмотрены алгоритмы расчета теплопроводности и определены основные факторы, влияющие на ее изменение с глубиной, для каждого из указанных пяти интервалов. Резкое уменьшение теплопроводности пород в низах осадочного чехла и в кровле фундамента в районе скв. СГ-6 связано с разуплотнением пород под влиянием процессов тектонического растрескивания и гидротермального выветривания. Проведенный анализ предполагает, что при измерениях теплопроводности методом оптического сканирования не соблюдаются условия стационарности процесса и что по этой причине указанный метод может завышать действительные значения теплопроводности пород.



Отклик песчаных и глинистых грунтов при слабых и сильных сейсмических воздействиях
Аннотация
Анализируется отклик песчаных и глинистых приповерхностных грунтов (представляющих классы несвязных и связных грунтов) на сейсмические воздействия различной интенсивности по данным наблюдений in situ – по записям вертикальных групп японской сети сильных движений KiK-net. Из общего числа станций (~800) для анализа выбрано по 5 станций с приповерхностными песчаными и глинистыми грунтами, представленными в верхних слоях в наиболее чистом виде. Для “песчаных” и “глинистых” станций по методике [Pavlenko, Irikura, 2003] построены и проанализированы модели поведения грунта при сильных движениях, показывающие вертикальные распределения напряжений и деформаций, вызванных землетрясениями, в грунтовых слоях. Получены близкие оценки усиления сейсмических волн в песках и глинах при слабых движениях и близкие зависимости напряжение–деформация, характеризующие поведение приповерхностных грунтов при сейсмических движениях умеренной силы. Выполнен анализ разжижения песчаных грунтов при сильных движениях (землетрясение Тохоку 2011 г., с Мw ~ 9.0). Исследовано влияние эффектов протяженных сейсмических очагов (направленность их излучения) на поведение песчаных и глинистых грунтов и усиление в них сейсмических волн. Различия в поведении песчаных и глинистых грунтов отмечены лишь при сильных движениях: в песчаных грунтах возможно разжижение при уровне грунтовых вод порядка нескольких метров от поверхности, а в глинистых грунтах разжижения нет.



Диссипативная сейсмика для оценки параметров углеводородного коллектора
Аннотация
Естественный сейсмический фон Земли и его глубинная эмиссионная компонента являются одной из форм диссипации энергии геодинамических процессов. Методы исследования сейсмического фона как сигнала порожденного открытой сложной нелинейной системой (корой Земли) могут быть объединены в одно направление – диссипативная сейсмика. В данной статье один из таких методов – термодинамический индикатор состояния горных пород используется для дистанционной оценки локальной продуктивности коллектора на территории месторождения углеводородов. Термодинамический индикатор разработан с использованием энтропии Климонтовича и дает количественную оценку локальной неравновесности пород, связанной с активностью геофизических процессов. Выявлена монотонная зависимость, близкая к линейной, между значениями термодинамического индикатора, рассчитанными по записям сейсмического фона, и суммарной толщиной продуктивных пластов в скважинах вблизи точек регистрации. При расчете термодинамического индикатора существует предельная величина скользящего временного окна, которая должна выбираться эмпирически. Из полученных результатов следует перспективность использования термодинамического индикатора для оконтуривания границ залежи и выбора наиболее перспективных точек для бурения скважин по экстремумам в поле значений индикатора.


