Решение задачи о точном уплощении Земли для волн Рэлея тема диссертации и автореферата по ВАК РФ 04.00.22, кандидат физико-математических наук Киселев, Сергей Геннадьевич

• Киселев, Сергей Геннадьевич
• кандидат физико-математических наук
• 1999

Оглавление диссертации кандидат физико-математических наук Киселев, Сергей Геннадьевич

симметрии

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

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

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

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

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

ШЛ-формы

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

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

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

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

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

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

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

напряжений

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

разложения,

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

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

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

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

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

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

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

Заключение

Литература

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

Введение диссертации (часть автореферата) на тему «Решение задачи о точном уплощении Земли для волн Рэлея»

Введение

Поверхностные сейсмические волны лявовского и рэлеевского типов служат ценным источником информации о строении Земли. Общие сведения о её внешних областях получены, главным образом, благодаря изучению поверхностных волн. Так как распространение этих волн происходит вдоль поверхности, то убывание их энергии вследствие геометрического расхождения пропорционально 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) Для Р-БУ колебаний решена задача точного уплощения, т.е. найден метод, позволяющий рассчитать колебания сферически слоистой среды при помощи быстрых и эффективных алгоритмов, разработанных для плоско-слоистых полупространств. Ранее такая задача была решена только для 8Н колебаний.

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

3) Используя одинаковую форму уравнений Штурма - Лиувилля для плоско и сферически-слоистых сред, можно указать матричное преобразование подобия, которое переводит уравнение для сферической среды в уравнение для плоской. Тем самым заданной сферической среде сопоставляется плоская среда, т.е. осуществляется точное уплощение.

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

5) Как правило, такую систему можно решить только численно. Однако выделен класс сферически слоистых сред, для которых эта система является автономной и допускает аналитическое решение. Сравнение этих решений с приближённым уплощением Бисваса показывает те особенности, которые не учитываются э тим приближённым уплощением.

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

7) Точное решение задачи уплощения представлено в виде ряда по обратным степеням радиуса сферы. Нулевому члену этого ряда соответствует приближение Бисваса. Учёт следующего, первого, члена позволяет уточнить параметры уплощающей среды. В результате этого собственные частоты колебаний Земли, заданной в виде моделей Гутенберга, PREM, IASP91, а также модели Луны, при расчёте предлагаемым методом уплощения получаются точнее, чем при использовании уплощения Бисваса.

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

Литература

1. АкиК, Ричарде П. Количественная сейсмология. М.: Мир, 1983. Т. 1, 520 с.

2. Гантмахер Ф.Р, Теория матриц. М., Наука 1988. 552 с.

3. Гервер М.Л., Каждан Д. А. Нахождение скоростного разреза по дисперсионной кривой. Вопросы единственности // Некоторые прямые и обратные задачи сейсмологии. М.: Наука, 1968. С 78-94. (Вычисл. сейсмология; Выи.4).

4. Завадский В.В., Киселев С. Г., Макеев O.A., Маркушевич В.М. Рэлеевские волны в средах Пикериса. II. Дисперсионные свойства. Теоретические проблемы геодинамики и сейсмологии. (Вычислительная сейсмология. Выпуск 27.) 1994, стр. 158-170, М., Наука.

5. Завадский В.В., Киселев С.Г., Макеев O.A., Маркушевич В.М. Рэлеевские волны в средах Пикериса, Доклады Академии Наук, 1995, 343(6), стр. 813-817, М.,Наука.

6. Киселёв С.Г., Маркушевич В.М., Цемахман A.C. Матричная задача Штурма-Лиувилля для рэлеевских волн в цилиндрически симметричных телах. Геодинамика и прогноз землетрясений, (Вычислительная сейсмология Выпуск 26), 1993, стр. 168-176, М., Наука.

7. Киселёв С.Г., Маркушевич В.М. О разделении переменных в уравнениях рэлеевских колебаний слоистых тел, Доклады Академии Наук, 1993, том 332, N 3, стр.297-300, М„ Наука.

8. Киселёв С.Г., Маркушевич В.М. Рэлеевские колебания слоистых тел как матричная задача Штурма-Лиувилля, Доклады Академии Наук, 1994, том 335, N 1, стр.29-31, М., Наука.

9. Киселев С.Г., Кузнецов А.Н., Маркушевич В.М. Задача уплощения Земли: происхождение, методы точного решения и разложение в ряд. Теоретические проблемы в геофизике. (Вычислительная сейсмология. Выпуск 29.) 1997, стр.28-43, М., Наука.

10. Киселев С.Г., Кузнецов Ä.H., Маркушевич В.М., Цемахман A.C. Разложение на множители и форма Штурма - Лиувилля уравнений для Р-SV колебаний слоистых сред. Теоретические проблемы в геофизике. (Вычислительная сейсмология. Выпуск 29.) 1997, стр.44-69, М., Наука.

11. Киселев С.Г. Обобщение метода отражений на неоднородные среды Пикериса Теоретические проблемы в геофизике. (Вычислительная сейсмология. Выпуск 29.) 1997, стр.70-80, М., Наука.

12. Киселев С.Г., Маркушевич В.М. Представление уравнений рэлеевских волн в форме Штурма-Лиувилля Математические вопросы теории распространения волн. Выпуск 26. Записки научных семинаров ПОМИ РАН. 1997, 239, стр. 110-116, ПОМИ РАН, Санкт-Петербург.

13. Киселев С.Г., Кузнецов Ä.H., Маркушевич В.М. Исследование алгоритма точного уплощения Земли для P-SV колебаний и сравнение с приближенным уплощением Бисваса на примере разреза Гутенберга. Вычислительная сейсмология. Выпуск 31, М., Наука, в печати.

14. Кузнецов А. Н. Функция Лагранжа и разделение переменных для упругих колебаний в осесимметричной слоистой среде.//' Теоретические проблемы геодинамики и сейсмологии (Вычислительная сейсмология, вып. 27). М. Наука, 1994, стр. 171-190.

15. Левшин А.Л. Поверхностные и каналовые сейсмические волны. // М.Наука. 1973. 176 с.

16. Лидский В. Б., Нейгауз М. Г. К методу прогонки в случае самосопряжённой системы второго порядка. ЖВМ и МФ, т. 2, №1, М. 1962с. 161-165.

17. Лурье AM. Пространственные задачи теории упругости. М., 1955. 492 с.

18. Маркушевич В.М., Резников Е.Л. Исследование строения симметричной твердой среды по стоячим SH-волнам на поверхности./ЛГеоретическая и. вычислительная геофизика. №2. Результаты исследований по международным геофизическим проектам. М. Междуведомственный геофизический комитет при президиуме АН СССР. 1974, стр. 5-34.

19. Маркушевич В.М. Вынужденные гармонические колебания рэлеевского типа и матричная рассеяния. Математические методы в сейсмологии и геодинамике. (Вычислительная сейсмология., вып. 19). М.: Наука, 1986. С. 119-135.

20. Маркушевич В.М. Подстановка Пикериса и некоторые спектральные свойства задачи Рэлея. Теория и алгоритмы интерпретации геофизических данных. (Вычислительная сейсмология, вып. 20). М.: Наука, 1989. С. 117-127.

21. Молотков Л. А. Матричный метод в теории распространения волн в слоистых упругих и жидких средах.Л., 1984, 201 с.

22. Снеддон И.Н., Берри Д.С. Классическая теория упругости. М. Гос. изд. физ.-мат. лит. 1961.220 стр.

23. Шкадинская Г. В. Теория и метод расчёта поверхностных волн Рэлея в неоднородных средах. Диссертация на соискание учёной степени к. ф.-м. н. ИФЗ им. О. Ю. Шмидта АН СССР, М. 1970,189 с.

24. Шкадинская Г. В. Метод расчёта поверхностных волн Рэлея в шаре. // Алгоритмы интерпретации сейсмических данных (Вычислительная сейсмология, вы п. 5. М. Наука, 1969, стр. 178-188.

25. Aho-Zena, A. Dispersion function computations for unlimited frequency values. Geophys. J. R. astr. Soc.,58, 1979, pp.91-105.

26. Alterman Z, Jarosh H, Pekeris C.L. Propagation of Rayleigh waves in the earth, Geophys. J. 4(1961), 219-241.

27. Bhattacharya S.N. Extention of the Thomson-Haskell method to non-homogeneous spherical layers. Geophys. J. R. astr. Soc. (1976) 47, 411-444.

28. Bhattacharya S.N. Extended formulation for Rayteigh-wave computations at high frequencies in a sherical layered earth. Geophys. J. R. astr. Soc. (1986) 84, 311-329

29. Biswas N.N., Knopoff L. Exact earth-flattening calculation for Love waves 11 Bull. Seismol. Soc. Amer., 1970, V.60, 1123-1127.

30. Biswas NN. Earth-flattening procedure for the propagation of Rayleigh wave // Pure Appl. Geophys., 1972, V.96, 61-74.

31. Chapman, C.H. The Earth flattening transfonnation in body wave theory. Geophys. J. R. astr. Soc.,35,1973, pp.55-70.

32. Chapman, C.H. The turning point of elastodynamic waves. Geophys. J. R. astr. Soc., 39, 1974, pp.613-621.

33. Doornbos, D.J. Seisffloiogical algorithms, Academic Press, 1988,

34. Ewing W. M.., Jardetzky W.S., Press F. Elastic waves in layered media. McGraw-Hill Book Co., Inc. 1957. 380c.

35. Gomherg J.S., Masters T.G. Wave form modelling using locked-mode synthetic and differential seismograms: application to determination of the structure of Mexico, Geophys. J. 4(1988), 193-218.

36. Henkin G.M. Markushevich V.M. Inverse problems for the surface elastic waves in a layered half-space, in Inverse Problems, Advances in Electronics and Electron Physics, Supplement 19, 190-204, Academic Press, Ltd., 1987.

37. Hoskins, L.M. The strain of gravitating sphere of variable density and elasticity, Trans. Am. Math. Soc., January, 1-43,1920.

38. Kennett B.L.N. Seismic wave propagation in stratified media. Cambridge University Press, 1983. 342 c.

39. Knopoff, L. A matrix method for elastic wave problems. Bull. Seism. Soc. Am., 54, 1964, pp. 431-438.

40. Koo B. Y.'C., M. Katzin. An exact Earth-flattening procedure in propagation around a sphere. Journal of Research of the National Bureau of Standards. D. Radio Propagation. Vol. 64D. №1. 1960. P.61-64.

41. Markushevich V.M. The detemination of elastic parameters of a half-spase using a monochromatic vibration field at the sutface, Wave Motion, 9, 37-49, 1987.

42. Markushevich V.M., Kiselev S.G., Steblov G.M., Tsemahman A.S. Reduction of the vector elastodynamic equation to the matrix Sturm- Liouville boundary value problem. Nonlinear Waves in Solids, Proceedings of IUTAM Symposium

on Nonlinear Waves in Solids, 1995, p.201-205 ASME Press (American Society of Mehcanical Engineers)

43. Mutter G. Earth-flattening approximation for body waves derived from geometric ray theory-improvements, corrections and range of applicability. Journal of Geophysics, v.42, p.429-436,1977.

44. Pekeris C. L. Accuracy of the Earth-flattening approximation in the theory of microwave propagation. Phys. Review, vol.70, № 7,8,1946, p. 518-522.

45. Pryce M. H. L. The diffraction of radio waves by the curvature of the Earth. Advances in physics, vol. II, № 5-8, 1953, p. 67-95.

46. Takeuchi, H. and Saito, M. (1972), Seismic surface waves. Methods in Computational Physics, Vol.11 (ed. BA.Bolt), pp.217-295. Academic Press, New York,

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

Система координат и V к К к 1

1, Пл оская декартова 1 1 1 1 0 0

2. Цилиндрическая I 1 1 1 р 0 0

3. Цилиндрическая II 1 1 ч 1 0 1

4. Плоская полярная 1 я 1 1 1 0

5. Сферическая 1 ч ч со 1 1

Табл.1. Компоненты метрического тензора для выделенных систем

координат

кг к кз к\ к5 ¿7

1,2. 0 1 0 1 0 г 0 1 1

3. X г 1 1 Но Г Х+(1 "г" 0 X. Го 2 г 1 1

4. Г'О 0 1 -2 1 0 е**

5. П) х/2 0 еф -2 0 а/2 х/2

Табл.2, Выражения для коз ффициентов в различных системах

координат

/I к /з

1,2. 1 1 0

3. В. г г Го

4/ 0 X (2 (/Г1)' + 27)

5.

Табл.3. Выражения для коэффициентов и в различных системах

координат

т\ т2 т% 0*4

1 з 1 0 0 1

3. г п> ¡40 / 2(1' {я+2 /<) (3 я+5 р) \ 1 я+2 /4 ¿за. г

г v и (Х+ц) 2 гр (\+1л)2 ) г л+//

4. е2х 2/4 0 Х-

5. е2х -2 до (/Г1)' 1/2 е~2х

Табл.4, Выражения для коэффициентов m¡ в различных системах

координат

а2 €а-2(/1 + 2)2(я + 3)€3

аг (7 4-1) О (2 (у 4- З)«2 +• (13 у 4- 34)» + 17 у + 42) €г

щ -(у + 2) (2 и 4- 2 7 4- 9) е

ь2 2(й-4(й42)е2)

2 (я - 1.) (7 4- 2) е

2(7 + 2)

Табл.5. Коэффициенты разложения элемента Кц

йг2 4 (л + 2) (2 и2 4- (6 7 4-16)» 4- 9 (74-2)) е4~ (2 (7 + 3)«* + 4 (6 7 4- И)п + 33 74-50)Ое2 4-2(7 4- 1)П2

-(у 4*2)е((2(74* I)«3 + 2(117+ 5)??2 4» (417-26)?! + 3 (57*~22))€2 + 2 {-у п 4- П 4- 7 4- 9) Л)

Оо -(7 4- 2) (2 (7 4-1) л2 4- 4 (у* + 5 7 4- 3) я + 472 + 13 7-6) с2

¿2 4(я2 +6л4-8)(7 + 2)б2

¿1 4 (и + 4) (у + 2)2 6

н 2 (7 + I) (7 + 2)2

Табл.6, Коэффициенты разложения элемента Ьц

«2 -2 (и2 е2 + 2 и е2 + О)

Щ -2 (2 п + 3) (у + 2) б

а0 -у2 - 5 у - 6 / /

Ъг 2 (П - 4 (п + 2) е2)

Ьх 2 (/г - 1)(у + 2)е

Ьо 2 {у+ 2) "1

Табд.7, Коэффициенты разложения элемента Кц

«4 -2 е (2 е4 п5 + 26 е4 «4 + 116 е4 я3 + (216 е4 + 11 П е2) п2 + (144 е4 + 58 О е2 -2 О2) п + (72 е2 - 7 £1) О)

Дз -2 ((и + 2) (2 (2 у + 5) пъ + (51 у +110) л2 + (149 у + 306) и + 108 (у + 2)) е4 -((5 у + 7) л2 + (13 у + 16) п - 2 (у + 6)) Не2 - (у + 1)О2)

-(у+ 2)е(2 ((у + 7) и3 + (31 у + 94) п2 + (113 у + 300) я + 15 (7 у + 18)) е2 -(6 у + 10/2 + 9у + 21)П)

-(у + 2) ((2 (у + 3)п2 + (17 у2 + 81 у + 102) п + 35 у2 + 165у + 198) е2 -(у2 + 4 у + 3) О)

й?0 -(у + 2)2 (2 у2 + 9 у + 9)€

¿4 2 (П 4 (й + 2) е2)2

-4 (и - 1) (у + 2) е (4 (й + 2) е2 - П)

62 2 (у + 2) (((у + 2) й2 - 2 (у + 6) п + у - 14) е2 + 2 П)

¿4 4 (й - 1) (у + 2)2 е

Ьо 2 (у + 2)1

Табл.8. Коэ ффициенты разложения, э лемента Ьц

<*2 12е2 (9(V4 - 5 У3 -42 v2 - 70 V-29)е4 + 3 (3 V2 + 5 V - 3)Пе2 +Пг)

ал -12 (9 (2 V5 - 5 у4 - 41 V3 - 5 V2 + 132 V + 100) б6 + 3(4 у3 - 5 V2 - 41 у - 25) П б4 + 2 V П2 б2)

Щ (3 (4 V2 + 1) е2 - 4 П) (3 (у2 - 4) е2 + О)2

Ьг -12 V е2 (27 (у6 - 20 V5 - 198 V4 - 630 V3 - 900 V2 - 590 V - 154) е6 + 9 (7 V4 + 20 V3 - 20 V2 - 80 V - 32)Пе4 + 3 (5 V2 + 10 V - 2)02 е2 + О?)

ъх 6 V €2 (27 (4 V7 - 25 V6 - 223 V5 - 265 V4 + 692 V3 + 1510 V2 + 544 V - 200) б6 + 18 (6 V5 - 20 V4 — 173 у3 - 190 V2 + 146 V + 175) Пб4 + 9 (4 V3 - 5 v2 - 41 V - 25) П2 б2 + 4 V О3)

Ь0 —v (3 (4 V2 + 1) б2 - 4 О) (3 (v2 - 4)б2 + П)3

й\ 108 (V2 + 5 у + 5)б4

й-) 36 (v2 + 5 V + 5) б3 (3 (V2 - 4) б2 + О)2

Табл.9. Коэффициенты системы уравнений (99),

1 V2

2 (у2 + 5 V + 5)2

3 12 V2 б 2 + 3 б 2 - 4 П

4 (3 V2 б2 - 12 е2 + П)4

5 (9 V2 б2 + 30 V б2 + 18 б2 + П)2

6 675V4б6 + 43362е6 -6003Пб4+ 294П2 б2 - 4 О3 + 45 V3(183б6-4б4Л) + 3 V2 (12591 б6 - 741 Пе4 + 4 П2 е2) + 45 V(1536б6 - 151Ле4 + 4 П2 е2)

Табл. 10, Сомножители уравнения (100),

г*

У

Рис Л. Зависимость коэффициентов ус,,/ ср, ср, ср исц от безразмерного частотного параметра Ш = со2 сг /Зд2. На графике коэффициентам соответствуют цвета:

цвет синий красный зеленый жёлтый

коэффициент \1 Су / Ср V /ср СР сц

глубина,км

Рис.2. Уплощение однородной сферы. На графике обозначены:

цветом параметры среды толщиной линий начальное значение матрицы Н .... сплошностью линий метод уплощения

синим УР сплошными ТОЧНЫЙ

красным V* тонкими единичное штриховыми в ряд по 1 / а

зелёным р жирными специальное пунктирными Бисваса

о

^ 12

л

& Ю

о а о

5 8 К

<а

о 6

М (0

100

250

150 200

период, сек

Рис.3. Дисперсионные кривые для модели Гутенберга, красные - уплошение Бисваса, зелёные - уплощение с дополнительным членом - 1 /а.

период, сек

Рис.4. Область рис.3, где наиболее существенно повышение точности.

0 ^

1 Ю

О

о о,

о

о 8

04 «

И О

5 6

100 150 200 250

период, сек

Рис.5. Дисперсионные кривые для модели PREM, красные - уплощение Бисваса, зелёные - уплощение с дополнительным членом ~ 1 /а.

75

100

200

125 150 175 период,сек

Рис.6. Дисперсионные кривые для модели Луны, красные - уплощение Бисваса, зелёные - уплощение с дополнительным членом ~ I /а.

225 250