Бесплатный автореферат и диссертация по геологии на тему

Решение задачи о точном уплощении Земли для волн Рэлея

ВАК РФ 04.00.22, Геофизика

Текст научной работыДиссертация по геологии, кандидата физико-математических наук, Киселев, Сергей Геннадьевич, Москва

/

РОССИЙСКАЯ АКАДЕМИЯ НАУК

Международный институт теории прогноза землетрясений

РЕШЕНИЕ ЗАДАЧИ О ТОЧНОМ УПЛОЩЕНИИ ЗЕМЛИ ДЛЯ ВОЛН РЭЛЕЯ

Специальность 04,00,22, - "Геофизика"

Диссертация на соискание ученой степени кандидата физико-математических наук

и математической геофизики

На правах рукописи

КИСЕЛЁВ Сергей Геннадьевич

УДК 550.310: 515.984.54

Научный руководитель доктор физ.-мат. наук профессор В.М.Маркушевич

МОСКВА 1999

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Введение....,.....................,....................................................,....................5

Глава 1. Уравнения рэлеевских волн в слоистых средах....................................12

1.1. Уравнения упругих колебаний.,........................................................12

1.2. Осесзшметрвчные деформации.................................................12

1.3. Определение слоистой среды.......................................................................................13

1.4. Разделение колебаний на два типа: $Н и Р-8У....................................................14

1.5. Рэлеевские уравнения в криволинейной системе

координат........................................................................................................14

1.5.1. Компоненты тензора деформаций...................................................14

1.5.2. Закон Гука..........................................................................................................15

1.5.3. Уравнения движения....................................................................15

1.5.4. Условия Ламэ............................................................................16

1.6. Преобразование уравнений движения......................................... 16

1:7. Условия разделения неременных........................................................17

1.8. Типы слоистых сред, допускающих разделение переменнных.................................................................................................................18

1.9. Единая запись уравнений для сред различных

симметрии...............................................................................................................................................................21

Глава 2. Единая форма Штурма-Лиувилля уравнений для

слоистых тел различных симметрий......................................................25

2.1. Сведение рэлеевских уравнений к ШЛ-форме........................................25

2.2. Частичная факторизация рэлеевского оператора................................................27

Глава 3. Уплощение как преобразование уравнений

ШЛ-формы.......................................................................................................29

3 Л. Инвариантное преобразование уравнений......................................29

3.2. Уплощающая система уравнений........................................................30

3.3. Сведение уплощающей системы к четырём уравнениям..........................................................................................30

3.3.1. Вывод выражения для £21...*..........................................................31

3.3.2, Вывод выражения для Ь^г • -...............................................................32

3.4. Нахождение параметров плоской среды...........................................32

3.5. Преобразование трансформант смещений и

напряжений..........................................................................................33

Глава 4. Аналитический пример уплощения............................................35

4.1. Стационарные точки уплощающей системы.....................................................36

4.2. Параметры плоской среда......................................................................40

4.3. Область физически допустимых значений............................................................41

4.4. Преобразование смещений и напряжений.......................................42

4.5. Уплощение как преобразование подобия.........................................43

4.6. Оценка выделенной частоты.............................................................44

4.7. Обобщение на произвольную частоту.............................................................44

4.8. Уплощение Бисваса для аналитического примера......................................48

Глава 5. Методы численного уплощения...............................................52

5.1. Гладкое уплощение............................................................................52

5 Л Л. Вывод уравнений дня начальных условий....................................................53

5.1.2. Решение уравнений для начальных условий.................................55

5.1.3. Анализ дискриминанта........................,.............-....................................56

5.2. Уплощение однородной сферы.....................................................58

Глава 6. Разложение точного решения в ряд.................................................60

6 Л. Произвольная сферически-слоистая среда..................................60

6.2. Однородный непуаесонов слой...................................................61

6.2.1. Вывод уравнения при первой степени

разложения,.................................................................................................61

6.2.2. Решение уравнения......................................................................................65

6.2.3. Нахождение параметров среды..........................................................................67

Глава 7. Расчёт дисперсии методом уплощения........................................................68

7.1. Метод Бисваса.................................................................................................68

7.2. Метод разложения в ряд точного решения.....................................68

7.3. Пакет используемых программ.........................................................69

7.4. Сравнение результатов расчётов.................................................69

Заключение............................................,....................................................70

Литература..........................................................................................................72

Таблицы и графики.......................................................................................77

Введение

Поверхностные сейсмические волны лявовского и рэлеевского типов служат ценным источником информации о строении Земли. Общие сведения о её внешних областях получены, главным образом, благодаря изучению поверхностных волн. Так как распространение этих волн происходит вдоль поверхности, то убывание их энергии вследствие геометрического расхождения пропорционально 1 / г, где г ~ расстояние от источника до точки наблюдения, в отличие от объёмных, для которых эта зависимость пропордаоначьна 1 / г1. Таким образом, на достаточно далёких расстояниях от источника все наблюдаемые вблизи поверхности сейсмические явления целиком связаны с поверхностными волнами. В частности, наибольшие разрушения при не слишком глубокофокусны х землетрясениях вызываются именно поверхностными волнами.

Вопреки своему названию, поверхностные волны могут при низких частотах проникать в толщи Земли достаточно глубоко. Например, если период колебаний в них более 400 секунд, то они проникают более, чем на половину радиуса, то есть достигают земного ядра, переходя на ещё более низких частотах в собственные колебания Земли.

Поскольку эти волны проникают достаточно глубоко, то они содержат информацию не только о внешних, но и о достаточно глубоких областях. Однако прочесть эту информацию значительно более сложно, чем из объемных волн. Теория обратных задач да® поверхностных волн фактически начала разрабатываться Пикерисом и Слихтером; первый исследовал обратную сейсмическую задачу, а второй - обратную задачу для электромагнитных полей (у Слихтера рассматривалась асимптотика). Более строгий и в прикладном отношении более полезный для практики подход к обратным задача^ для поверхностных волн стая возможен после разработки теории обратных спектральных задач Гельфандом, Левитаном и Марченко в 50-60х годах. Обобщение этой теории на матричные задачи было сделано

группой исследователей под руководством Г.М.Хенкина в работах [36 Применением этой теорий к сейсмическим колебаниям типа Лява и Рэлея в последнее время занимается груша исследователей под руководством В .М.Маркушевича.

Для отработки решения обратных задач необходимо уметь моделировать рэлеевские колебания. Расчёт рэлеевских волн в некоторой модельной упругой среде является очень трудоёмкой вычислительной задачей. Для облегчения расчётов используются различные упрощения. Например, рассматривается среда, неоднородная только в одном, обычно вертикальном направлении, называемая слоистой. В дальнейшем будут рассматриваться только такие среды. Точное определение слоистой среды даётся ниже в разделе 1.3.

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

Если сферичность продолжать игнорировать и рассматривать колебания плоской среды, то будет внесена некоторая ошибка. Эту ошибку можно компенсировать следующим искусственным приёмом. Сферическую среду следует заменять плоской с "подправленными" упругими свойствами и плотностью. Иногда такое "исправление" плоской среды оказывается столь успешным, что отличие колебаний в плоской среде от сферической исчезает полностью, правда плоские колебания также необходимо подправлять, но эта поправка также точно находится. Нахождение "поправок" двух сортов, а

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

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

где щ ~ щ{г) - трансформанта смещения в упругой волне как функция г -радиус-вектора из центра сферы, д* = - упругий модуль сдвига, р3 ~ р${г) - плотность среды, а> - частота, § - волновое число. Штрих' обозначает дифференцирование по г, то есть Г ~д/¡дг.

Если в это уравнение подставить

где штрих' теперь обозначает дифференцирование но % - а Ща ¡г), а - радиус

Земли.

Уравнение (4) по форме совпадает с уравнением лявовских волн с трансформантой смещения и/ в плоско-слоистом полупространстве со свойствами /г/ и р/ и глубиной измеряемой от поверхности.

А* - МЛ?) = {а}г? у/ОЫа/г)), Ра ~ рАг) = аъ Г5 р/(Ы(а/ г»,

Щ = Щ(г) - (г/а)м/0а(<а/г)),

(3)

(2)

то обнаружится, что и/ удовлетворяет уравнению и/((02 + И/Ду- + Д/М/ « О,

(4)

Чтобы получить модуль сдвига ¡1 / подправленной плоской среды, необходимо модуль сдвига цх сферической среды умножить, как следует из

(2), на (а¡г)~3, а плотность на а"3 г5. При этом плоское решение Uf будет отличаться от сферического us, как следует из (3), на множитель а ¡г.

Таким образом, по заданным свойствам fis(r) и ps(r) сферически слоистой среды находятся из (2) такие свойства ji/(z) и p/(z) некоторой плоской, что исходное сферическое уравнение (1) принимает вид плоского (4) для некоторой функции и/, известным образом, по формуле (3), связанной с искомой us, то есть "уплощается". Появляется возможность вместо уравнения (1) исследовать более простое уравнение (4) и все результаты перенести на исходную задачу, используя преобразование (3).

Точное уплощение для лявовских волн было найдено Гервером и Кажданом

[3] и Бисвасом и Кноповым [29]. Были попытки найти подобное преобразование и для рэлеевских волн, но до сих пор оно не было найдено. Это объясняется более сложным характером рэлеевских колебаний по сравнению с лявовскими; в частности, они описываются системой двух уравнений второго порядка.

Вместо точного уплощения для рэлеевских волн предлагались различные варианты приближённых. Наиболее известен среди них вариант приближённого уплощения, предложенный Бисвасом [30]. При этом приближённом уплощении формулы, аналогичные (2) и (3), имеют вид:

Mr) = Mf(Ha/r)l A.s(r) = АДЫя/г)),

j (5)

Ps(r) = (a jrf р/(\п(а/r)),

it sir) = (г/ a) и/(ln(a / /•)), (6)

где Xs - ещё один упругий модуль, входящий наряду с ¡is в уравнения рэлеевских волн. Необходимо обратить внимание, что теперь, в рэлеевском случае, амплитуды волн us и й/ являются двухкомпонентными векторами, и обозначены жирным шрифтом.

При подстановке формул (5) и (6) в уравнения рэлеевских волн в сферической среде, получаются аналогичные уравнения для плоской среды плюс слагаемые вида /(г)/л, где /(г) - числитель слагаемого, зависящий от радиуса. Для небольших глубин г%а выполняется /(/")«а, следовательно формулы (5) и (6) определяют приближённое уплощение.

Несмотря на то, что данное уплощение справедливо для небольших глубин, с точки зрения практики эти глубины часто оказываются достаточными. Так в статье Бисваса [30] показано, что этим уплощением вполне можно пользоваться при небольшой погрешности для расчёта рэлеевских волн, которые проникают не глубже половины радиуса Земли.

Погрешности при этом оказываются не столь велики для не очень низких частот и поэтому многие расчёты, например, Гомберг и Мастере [35], предпочитают делать, используя именно уплощённую модель, а не сферическую.

Известно и много других вариантов приближённого уплощения, например Альтерман, Яром, Пикерис [26]. То обстоятельство, что в течение почти 30 лет предлагались различные варианты уплощений (приближённых) для рэлеевских волн и делались безуспешные попытки найти точное» говорит о том, что точное уплощение является актуальной задачей, решение которой даст существенные преимущества при расчётах сейсмограмм.

Эти преимущества состоят в следующем. При расчётах сейсмограмм наиболее эффективными оказываются матричные методы, такие как метод матричного нропагатора, метод отражений. Эти методы детально разработаны именно для плоско-слоистых сред, но имеют ограниченное применение для сферических. Применение матричных методов для сферических сред известно (Бхаттачария, [27,28]), но их эффективность сильно уступает аналогичным методам для плоских сред, так как требует разбиения сферической среды на слишком большое количество тонких слоёв, особенно глубинных. Применение уплощения позволило бы использовать для

сферических сред быстрые матричные методы расчёта сейсмограмм, первоначально разработанные для плоских сред.

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

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

Последовательное изложение этой теории и её применение к решению классической задачи уплощения Земли для рэлеевских волн и составляет содержание настоящей работы. В первой главе формулируется «остановка задачи, даётся определению тому классу слоистых упругих сред, которые входят в рассмотрение излагаемой теории, приводятся уравнения рэлеевских колебаний в произвольной криволинейной ортогональной системе координат и находится класс сред, для которых возможно разделение переменных в данной системе уравнений. Во второй главе разделённые уравнения приводятся к единой форме Штурма-Лиувщщя. Третья глава посвящена решению задачи уплощения с помощью единого представления уравнений в сферической и плоской средах.

Точное решение задачи уплощения для рэлеевских волн уже не имеет такого простого вида, как формулы (2) и (3) или (5) и (6). В общем случае

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

Для проверки эффективности метода точного уплощения рассчитаны дисперсионные кривые для моделей Земли: Гутенберга и PREM, а также для модели Луны.

Глава 1. Уравнения рэлеевских волн в слоистых средах

В дальнейшем интерес будут представлять рэлеевские колебания как пяоско-слоистых сред, которые естественно рассматривать в декартовой системе координат, так и сферически слоистых сред, наиболее просто описываемых в сферической системе координат, поэтому необходимо вывести уравнения этих �