Численно-аналитическое моделирование волновых полей в неоднородных средах тема диссертации и автореферата по ВАК РФ 05.13.18, доктор физико-математических наук Фатьянов, Алексей Геннадьевич

  • Фатьянов, Алексей Геннадьевич
  • доктор физико-математических наукдоктор физико-математических наук
  • 2005, Новосибирск
  • Специальность ВАК РФ05.13.18
  • Количество страниц 252
Фатьянов, Алексей Геннадьевич. Численно-аналитическое моделирование волновых полей в неоднородных средах: дис. доктор физико-математических наук: 05.13.18 - Математическое моделирование, численные методы и комплексы программ. Новосибирск. 2005. 252 с.

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

Введение

Глава 1. Численное решение задачи Лэмба для неоднородных неупругих сред Больцмана с экспоненциальными функциями последействия

1.1. Постановка задач

1.2 Сведение к дифференциальным уравнениям

1.3. Конечные интегральные преобразования по пространственной переменной в задачах распространения неупругих волн

1.4. Некоторые вопросы сходимости и оценка точности метода

Глава 2. Полуаналитический метод решения прямых динамических задач для различных моделей слоистых сред и источников

2.1. Конечные интегральные преобразования по пространственной и временной переменным в задаче Лэмба для сред Больцмана с произвольными функциями последействия

2.2. Расчет полных волновых полей в неупругих средах

2.3. Методы решения краевых задач, полученных после отделения переменных

2.4. Полуаналитический метод расчёта волновых полей

2.5. Нестационарные. волновые поля в анизотропных средах с поглощением энергии

2.6. Постановки задач и некоторые вопросы реализации полуаналитического метода в средах с гравитацией

2.7. Волновые поля в слоистых пористых средах

2.8. Волновые поля от точечных источников дислокационного типа и источников конечных размеров

2.9. Расчёт однократных и монотипных волн без использования коэффициентов отражения и сравнение с лучевым методом

2.10. Волновые поля в разномасштабных средах

Глава 3. Метод расчета функции Грина в многомерно-неоднородных средах

3.1. Постановка задачи

3.2. Энергетический метод расчета функции Грина для многомерно-неоднородных моделей сред

3.3. Волновые поля в средах с криволинейными границами

3.4. Некоторые численные эксперименты и вопросы реализации метода для различных сред

Глава 4. Некоторые обратные динамические задачи

4.1. Оптимизационный метод одновременного определения скорости и декремента поглощения для акустических сред

4.1.1. Постановка задачи

4.1.2. Вычисление градиента функционала и построение итерационного процесса

4.1.3. Исследование чувствительности функционала и некоторые результаты расчётов

4.2. Обратная динамическая задача определения дислокационных параметров

4.2.1. Восстановление компонент направленной силы

4.2.2. Одновременное определение компонент тензора и входного импульса

4.2.3. Некоторые примеры численных расчётов

Глава 5. Численное моделирование волновых полей для некоторых моделей неоднородных сред 174 5.1. Технологические вопросы математического моделирования волновых полей

5.2. Волновые явления в средах с поглощением энергии

5.3. Моделирование волновых процессов в анизотропных средах с поглощением энергии

5.4. Волновые поля в сложно построенных средах

5.4.1. Моделирование вибросейсмических волновых полей

5.4.2. Вертикальное сейсмическое профилирование и межскважинное просвечивание

5.4.3. Волновые поля сложно построенных сред Сибири

Рекомендованный список диссертаций по специальности «Математическое моделирование, численные методы и комплексы программ», 05.13.18 шифр ВАК

Введение диссертации (часть автореферата) на тему «Численно-аналитическое моделирование волновых полей в неоднородных средах»

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

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

• оправдана во многих случаях, но требует усовершенствования при исследовании спектрального состава, затухания колебаний в различных типах волн и т.п. за счёт неидеальной упругости среды. В данной работе поглощение описывается в рамках наиболее общей модели среды Больцмана с упругим последействием [1].

В настоящее время существует большое количество различных методов решения прямых динамических задач сейсмики. Одним из

• наиболее популярных, является лучевой метод, предложенный в работе [79]. В последствие он развивался во многих работах.

Приведём только некоторые из них: [28], [80], [81], [82], [83], [88], [157]. Лучевой метод с учетом только нулевого члена лучевого ряда требует меньше вычислительных затрат чем другие методы. Кроме того, с помощью лучевого метода можно легко учесть вклад тех или иных сейсмических волн в формировании сложного волнового поля. Однако, с развитием высокоточной широкополосной сейсмологической аппаратуры, появились факты регистрации «нелучевых» сейсмических волн, которые не описываются нулевым членом лучевого ряда, и для их вычисления необходимо учесть последующие члены ряда. В работах [141], [143] такие «нелучевые» волны были обнаружены с помощью численно-аналитических методов моделирования сейсмических полей.

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

В настоящее время существует ряд методов расчета волновых полей в слоистых средах. Укажем на некоторые из них. Еще в 50-х годах Г.И. Петрашень ([21]) построил точное решение для слоисто-однородных моделей сред в замкнутом виде. Полученное в виде многомерных интегралов решение было исследовано асимптотическими методами. Сделано большое число физических выводов о протекании в них физических процессов [22]. После стандартных процедур преобразования физических переменных по времени и пространству посредством преобразования Фурье и Бесселя (в предположении, например, радиальной симметрии) уравнения динамической теории вязко упругости сводятся к системе обыкновенных дифференциальных уравнений, которые решаются для дискретных частот и волновых чисел. Эти уравнения решаются в основном разновидностями матричного метода ([23]-[27]) либо асимптотически ([28]-[30]). Широко распространён также "рефлективный" метод ([31]-[32]) расчёта части волнового поля. Все эти методы имеют свою область применимости. В матричном методе, например, при расчёте волновых полей для высоких частот (тонких слоёв) происходит потеря точности. Для преодоления этого приходится рассматривать матрицы более высоких порядков, чем это необходимо по физической постановке задачи [27]. Отметим, что в матричном методе эта проблема не решена принципиально [33]. Асимптотические методы при всей их компактности и элегантности имеют ограниченную лучевыми представлениями область применимости. При этом развитие их на более сложные модели (например, с поглощением энергии) требует значительных усилий. Рефлективный метод позволяет надёжно рассчитывать только часть волнового поля, соответствующую, например, объёмным волнам.

В начале 70-х годов в Вычислительном Центре СО РАН (в последствие институт Вычислительной Математики и Математической Геофизики) стали развиваться численно-аналитические методы решения задач сейсмики, основанные на расщеплении двумерных и трехмерных задач на серию независимых одномерных с помощью интегральных преобразований по горизонтальным координатам и с последующим решением их конечно-разностным методом ([6],[90]). Применение этих методов показало, что для получения приемлемых по точности результатов при численной дискретизации необходимо брать не меньше 20 точек на минимальную длину волны для расчёта волновых полей на расстояния порядка 100 длин волн, что приводило к большим объёмам оперативной памяти. Применение интегрального преобразования Лагерра по времени с последующей дискретизацией пространственных переменных позволяет снизить вычислительные затраты [91-92]. Однако уменьшение точности при расчёте сейсмических волновых полей на большие расстояния приводит к необходимости развития, например, конечно-разностных методов с высоким порядком аппроксимации [158]. Кроме того, применение преобразования Лагерра для сложно построенных сред, например для сред с последействием, приводит к громоздким выкладкам. При этом теряется смысл частоты колебаний, а тем самым наглядность при исследовании физики распространения волн в сложно построенных средах.

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

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

Во второй главе предложен полуаналитический метод решения прямых динамических задач для различных моделей слоистых сред и источников. Сначала выясняются возможности и условия применения конечных интегральных преобразований Фурье по t для решения задач из первой главы. После отделения пространственной и временной переменных получено двухпараметрическое семейство одномерных краевых задач, которые для дискретного набора параметров к, и лп можно решать численно с использованием современных компьютеров. На выбор численного метода при этом принципиальное значение имеют быстродействие и точность. Алгоритм должен быть устойчивым и мало зависеть от случайных ошибок, которые неизбежны при расчетах волновых полей на большие расстояния и для больших частот. Требуется также правильно передавать решения краевой задачи с кусочно-непрерывными коэффициентами, так как литологически реальные среды имеют блочную структуру и параметры среды, а вместе с ними коэффициенты краевых задач, рвутся на некоторых горизонтах. Кроме того, алгоритм должен быть эффективным в плане затрат компьютерных ресурсов. В силу этого, интегро-интерполяционным методом для сильно контрастных сред строится трехточечная разностная схема со вторым порядком аппроксимации, дополненная известными условиями сопряжения в точках разрыва коэффициентов. Далее для нахождения решения используется метод прогонки.

В случае слоистых сред предлагается алгоритм, позволяющий рассчитывать полные волновые поля для слоисто однородных сред Больцмана с произвольными функциями последействия. Количество слоев, их толщины, параметры среды могут быть произвольными величинами. Вычисления при этом строятся рекуррентно, по однотипным формулам, что допускает естественное распараллеливание на многопроцессорных системах. Заметим, что К. Аки и П. Ричарде ([33]) говорят о нецелесообразности перехода к потенциалам из за возникающих трудностей. В данном методе переход к потенциалам не только не усложняет задачу, но приводит к её значительному упрощению. Полуаналитический метод развит для расчёта любых линейных моделей геофизических сред. В работе рассматриваются нестационарные волновые поля в анизотропных средах с поглощением энергии. В основу теории расчёта волн в анизотропных неупругих средах взяты определяющие соотношения Вольтерра [41], учитывающие влияние упругого последействия. А именно, упругие постоянные си заменяются I интегральными операторами q: cttxscttx(t)-c9 \hIJ{t-z)x{z)dz, со в которых с'у и htJ определяют, соответственно, уровень и спектральный состав поглощения энергии. В практике оперируют с и коэффициентами поглощения различных типов волн, декрементами затухания и другими аналогичными величинами. Поэтому для моделирования неупругих эффектов разработаны методы определения с0 и hu по различным физическим параметрам поглощения. В связи с этим вначале рассматриваются спектральные характеристики поглощения для "просто" неупругих сред. В качестве модельных при определении коэффициентов поглощения выбраны одномерные уравнения распространения

82ип плоских продольных и поперечных волн: = p-^-f,

Я2 г? р-Ь., в этом случае, часто встречающаяся на практике линейная зависимость коэффициента поглощения от частоты аппроксимировалась с использованием подхода, предложенного в [45]. В этой работе отказываются от частотно независимого поглощения, но с затуханием, обеспечивающем его практическое постоянство в некотором частотном диапазоне. Как показывают многочисленные экспериментальные исследования, в ограниченном диапазоне частот для широкого класса горных пород характерна близкая к линейной зависимость коэффициентов поглощения продольных и поперечных волн от частоты. При этом исключается случай сильно водо-насыщенных пород, в которых эта зависимость близка к квадратичной. Здесь уже более правомерно рассматривать случай пористой среды, заполненной жидким флюидом, что так же рассмотрено в работе. Модель Максвелла не даёт линейной зависимости от частоты [46]. Показано, что это отклонение от линейности тем меньше, чем уже спектр входного сигнала по отношению к доминирующей частоте. Т.е. для входных импульсов с не очень широким спектром модель

Максвелла практически даёт линейность коэффициента от частоты и отсутствие дисперсии, часто наблюдаемую на практике. При введении физических параметров поглощения в анизотропной среде учитывается, что в ней характер распространения волн зависит от угла между направлением распространения и осью Oz. Показано что для замкнутого описания поглощения в рамках линейной теории наследственности требуется знать пять декрементов поглощения: ахр, д^, д^, д||5, аор. Отметим, что теорию поглощения в анизотропной среде можно строить, зная декременты (коэффициенты) поглощения, например, для нескольких направлений. Данный подход, на взгляд автора, представляется более оптимальным, так как позволяет использовать обширную справочную литературу по поглощению в упругих средах. Введение поглощения в анизотропных средах, в первом приближении, можно осуществить следующим образом. Предположим, что декременты поглощения квазипродольных и квазипоперечных волн совпадают с декрементами продольных и поперечных волн: д^ = д|и, = аор = Ар, au = д||5 = д5.

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

Далее во второй главе рассматриваются постановки задач и некоторые вопросы реализации полуаналитического метода в среде с гравитацией. В своей фундаментальной работе [49] Ляв рассчитал статические деформации и малые колебания однородного гравитирующего сжимаемого шара. Выделим работы [50]-[51], в которой охватываются единым методом исследования проблемы общей теории упругих колебаний Земли. В данной работе в рамках слоистой модели Земли в сферической системе координат рассматриваются уравнения упругих колебаний гравитирующего шара [49] с различными (в том числе тензором сейсмического момента) источниками. Решение находится с помощью аналитических преобразований после применения ортогонального векторного разложения [33]. Данный метод применён также для расчёта гравитационных волн цунами в рамках широко известной модели Подъяпольского. Отметим, что здесь удалось одновременно рассчитать как упругие, так и гравитационные волн, что не было известно раньше. Расчёт гравитационных волн требует учитывать то обстоятельство, что их скорость существенно (десятки раз) ниже скорости упругих волн в твёрдой Земле. Это приводило к необходимости измельчать пространственный шаг в конечно-разностных методах, которые ранее использовались для решения этой задачи. Всё это приводило к значительному увеличению объёма вычислительной работы, что делало эту задачу неподъёмной. Полуаналитический метод, в силу его большей эффективности, позволил сделать это. Важное значение в задачах разведочной геофизики имеет численное моделирование распространения упругих волн в слоистых средах с пластами пористых, насыщенных вязкой жидкостью пород. Для описания волновых процессов в таких средах чаще всего используются уравнения теории Био в том или ином виде. В данной работе рассматриваются уравнения Био для расширенного частотного диапазона, приведенные в [59]. Развит метод расчёта волновых полей и проведено моделирование процессов распространения пористых колебаний. Отметим, что слоистая упруго-жидкая модель является частным случаем модели Био [87], число неравных коэффициентов которой в изотропном случае сокращается до трёх [159].

Данный подход применим для рассмотрения источников дислокационного типа и источников конечных размеров. Для описания этих явлений используется понятие тензора сейсмического момента [33]. В работе решена задача моделирования волновых полей для такого источника в анизотропной среде. Прямая динамическая задача в этом случае сведена к четырём задачам типа P-SV и двум для SH волн [70]. Это даёт возможность устойчивого и точного расчёта волновых полей в средах с произвольным количеством слоёв и с произвольным отношением длин волн различных типов к мощностям слоёв для произвольных тензоров сейсмического момента. Анализ данных, проведённых многими исследователями, показывает, что источник землетрясений часто имеет протяжённые размеры по пространству. Для учёта этого явления в работе сконструирован источник конечных размеров. Показано, что он является пространственным аналогом классического источника типа центра давления. Отметим, что таким образом с помощью развитой методики можно рассчитывать волновые поля для произвольно заданного источника конечных размеров.

В задачах моделирования волновых полей в сложно построенных средах большое, а зачастую и определяющее значение имеют алгоритмы, позволяющие рассчитывать динамику отдельно взятых волн. В настоящее время единственным методом численного анализа волнового поля по частям является асимптотический лучевой метод. Его применение, однако, имеет известные ограничения. В работе впервые предлагается алгоритм расчёта однократных и монотипных волн для слоистонеоднородных сред с произвольным числом слоёв на основе специальных разложений точных решений не имеющий ограничений лучевого метода. В [77] проведено сравнение лучевого метода с полуаналитическим в случае однократных продольных и поперечных волн [78]. Показано, что отличия наблюдаются не только в точке выхода головной волны, как считалось ранее, но и в некоторой области, зависящей от длительности входного импульса. При этом как показывают численные расчёты, продольные волны лучевым методом рассчитываются точнее, чем поперечные.

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

Для построения модели разномасштабной неоднородности используется телеграфный случайный процесс [160]. Отличие данного подхода от других исследований заключается в том, что так построенная разномасштабная среда позволяет точно рассматривать среднее поле для мелко- и крупномасштабных неоднородностей. Отметим, также, что в разномасштабной среде наблюдается наличие параметрического резонанса, как и в уравнениях для микронеоднородных сред, полученных в [161].

В третьей главе предлагается метод расчета функции Грина в многомерно-неоднородных средах. Существующие методы расчёта в таких средах в общем можно разделить на аналитические и численные. Все их привести в данной работе представляется невозможным. Отметим только некоторые из них [79]-[92]. В данной главе развивается метод расчёта функции Грина для многомерно-неоднородных моделей сред. Он основан на спектральных разложениях решения по пространственным и временным частотам. Поскольку полного разделения переменных здесь уже не происходит, возникают обыкновенные матричные дифференциальные уравнения. Их решения строятся аналитически. В итоге получен алгоритм, позволяющий моделировать волновые поля для блоковой геометрии сред. Как показали численные эксперименты, после соответствующей аппроксимации, данный алгоритм позволяет также моделировать произвольные границы. Расчёт функции Грина позволяет, в случае необходимости, размножать волновые поля для произвольного количества источников без дополнительных вычислительных затрат, что актуально, например, для задач площадной сейсморазведки и решения обратных задач оптимизационными методами.

Данный метод без принципиальных изменений переносится на все постановки из главы 2. Отметим только, что в матричном виде данный подход в этом случае полностью сохраняет свою структуру. Меняются только в зависимости от рассматриваемой среды соответствующие матрицы. Поскольку сборка решения в физической области осуществляется путём суммирования независимых величин, это даёт возможность эффективной программной реализации метода на многопроцессорных вычислительных комплексах. Это обеспечивает возможность проводить моделирование для любых реальных сред с произвольными геометрией и параметрами [162]-[165].

Энергетический метод расчета функции Грина для многомерно-неоднородных моделей сред с учётом принципа взаимности позволяет учесть наличие слоистой пачки с произвольными величинами мощностей слоёв и скоростных параметров [101].

В работе теоретически показано, что число обусловленности соответствующих матриц практически не зависит от свойств среды. Это позволяет говорить о вычислительной независимости метода от структуры и свойств сложно построенной среды. Доказано, что число обусловленности зависит в первой степени от количества удерживаемых членов ряда. Отметим, что при сеточных методах решения аналогичных задач числа обусловленности возникающих при этом матриц величины порядка 0{\) (например, [75]). Это и позволяет решать матричные системы стандартными методами без потери точности. В случае произвольного изменения параметров среды по латерали (х) в каждом неоднородном слое соответствующие интегралы вычисляются с помощью процедур интерполяции и сглаживания [90]. В ряде случаев, когда, например, границы рассматриваемой области уходят на бесконечность или когда граничная поверхность имеет геометрически сингулярные точки, например, острые рёбра, можно получить несколько математически корректных решений уравнения, среди которых лишь одно верно описывает исследуемое явление [96]. Дополнительное физическое условие, необходимое для однозначного определения решения, в этом случае известно как условие на ребре [96]. Оно заключается в требовании конечности энергии поля, запасённого в любом конечном объёме в окрестности ребра. В работе доказано, что сингулярные точки для так построенного решения не являются источниками излучения.

В третьем пункте данной главы предлагается метод расчёта волновых полей в средах с криволинейными границами раздела. Такие среды можно рассчитывать методом, описанным выше. Однако для таких границ после применения интегральных преобразований решение удаётся получить в явном виде. Для решения этой задачи вводится мощность слоя. Эта величина вещественна в случае плоской границы и совпадает с обычной мощностью слоя. В случае криволинейной границы расстояние становится комплексным. Использование комплексных величин приводит к тому, что условия ограниченности решения на бесконечности и принцип излучения Зоммерфельда выполняются автоматически. При этом наличие в решении дифракционной составляющей автоматически приводит к выбору нужного решения [166]-[167]. Без принципиальных трудностей данный метод переносится на различные геофизические среды. Учёт, например, анизотропного поглощения производится в соответствие с [34], [42]. Кроме того, можно рассматривать различные распределённые и сосредоточенные источники, например, тензор сейсмического момента [69]. Отметим, что в работе [99] в рассмотрение вводится комплексная фаза, что позволило решить некоторые задачи, к которым лучевой метод в обычной форме неприменим.

Далее в третьей главе рассматриваются некоторые примеры расчёта волновых полей для сложно построенных сред. По выше приведённому алгоритму составлен комплекс программ для проведения численного моделирования волновых полей в средах, возникающих в практических задачах. Расчёты проводились на 64-х разрядных процессорах. В качестве иллюстрации возможностей предлагаемого метода рассматривается моделирование волновых полей для вулканической области. Среда, в том числе и магматическая камера, аппроксимируется с помощью произвольного количества прямоугольных блоков. Так как в итоге находится функция Грина, это даёт возможность моделирования для различных режимов входных сигналов (монохроматическом, свип-сигналов и др.). Приведён пример расчёта волновых полей для очаговой зоны живущего вулкана Эльбрус (Собисевич А.Л. и др., Катастрофические процессы, Москва, 2002). Модель характеризуется пачкой слоёв, соответствующей осадочным породам, гранитам и базальтам. Магматическая камера взята в виде квадратного включения. При этом очаговая зона является не полой, а содержит среду с более низкой скоростью.

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

На практике требуется проводить моделирование для различных геофизических сред, включая упругие, неупругие, анизотропные и т.д. Данный метод без принципиальных трудностей переносится на все постановки из главы 2. По своей структуре предлагаемый метод состоит в матричном пересчёте, в отличие от скалярного для "чисто" слоистых сред. Иначе говоря, модель среды "сидит" в элементах соответствующих матриц. Таким образом, использование матриц позволяет универсализировать структуру алгоритма для различных моделей сред. В качестве примера в работе приведён расчёт для реальной нефтяной структуры, модель которой взята из промышленного пакета "Tesoral" (Канада).

Для частных геометрий среды использование, например, вместо преобразований Фурье преобразования Бесселя позволяет проводить трёхмерные (3D) расчёты на современной вычислительной технике [102].

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

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

В последнее время в связи с быстрым развитием вычислительной техники одним из наиболее часто используемых методов решения обратных задач является метод минимизации целевого функционала [108-112]. В [113] исследована обратная задача для уравнения акустики в интегральной постановке и получена оценка скорости сходимости в среднем метода наискорейшего спуска.

В работе оптимизационным методом решена задача одновременного определения скорости и декремента поглощения для акустических сред по информации о режиме колебаний на свободной поверхности. Для решения данной задачи использован способ поглощения, описанный в [33]. Как следует из [33], подход Больцмана, если функция последействия получена на основе хорошо установленных механизмов внутреннего трения таким образом, что обеспечивает практически постоянное значение Q на сейсмических частотах, согласуется с логарифмической зависимостью фазовой скорости от частоты. На этой основе получено представление для комплексной скорости в случае модели Максвелла. Искомое решение ищется как точка минимума целевого функционала, характеризующего собой квадратичное отклонение зарегистрированного волнового поля от рассчитанного для текущей модели среды. Получены формулы для градиента целевого функционала в комплексной области и исследована его чувствительность. Это позволило построить итерационный процесс одновременного определения скорости и декремента поглощения в спектральной области. Показано, что за счёт выбора оптимальных диапазонов пространственных частот в зависимости от интервала временных частот, можно повысить чувствительность целевого функционала к низкочастотным компонентам скоростной функции и декремента поглощения, как и при восполнении информации за счёт соответствующего выбора пространственных частот в случае отсутствия поглощения [108]-[111]. Для построения точки минимума целевого функционала были использованы комбинации методов сопряжённых градиентов (модификация Флетчера-Ривса) и наискорейшего спуска [117]-[118]. Результаты численных расчётов показали, что наиболее эффективно использовать комбинацию приведённых методов. Приведены результаты численных экспериментов и выяснены основные физические факторы, влияющие на сходимость метода в целом. По зарегистрированным данным на свободной поверхности развит метод определения всех компонент точечного тензора сейсмического момента в предположении, что временная форма сигнала в источнике неизвестна. Решение искомой обратной задачи строится на основе решения прямой, приведённой выше. Доказана устойчивость решения обратной задачи. Теоретические аспекты метода рассматриваются на примере слоистой трансверсально-изотропной неупругой среды. Основные этапы алгоритма рассматриваются на примере реконструкции источника типа точечной направленной силы. Пусть на свободной поверхности z=0 на линии наблюдения в = 0о известен в некоторых пунктах приёма вектор смещения, вызванный действием сосредоточенной силы на некоторой глубине. Обратная задача формулируется следующим образом: по этой информации требуется найти угол действия силы. Решение данной обратной задачи строится оптимизационным методом в физической области. Показана положительная определённость гессиана целевого функционала. Это обеспечивает единственность и устойчивость точки глобального минимума.

Далее рассматривается задача определения всех компонент тензора сейсмического момента. Эта задача рассматривалась во многих работах. Отметим только некоторые из них [121-123], в которых среда, как правило, считается однородной. Разработанные в данной работе эффективные алгоритмы построения функции Грина в неоднородных средах позволили снять это ограничение.

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

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

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

По алгоритмам, описанным выше, разработан комплекс подпрограмм, позволяющий для сосредоточенных источников широкого класса (центр давления, произвольно ориентированная сила, тензор сейсмического момента) проводить расчёты геофизических полей различной природы в сложно построенных средах на единой программно-алгоритмической базе с учётом поглощения, анизотропии и других усложняющих задачу факторов. В [58] начала осуществляться программная реализация и создание транспортабельного пакета прикладных программ, который может использоваться как самостоятельно для расчёта полных волновых полей при моделировании геофизических процессов в реальных средах, так и в качестве подсистемы решения прямой задачи, встроенной в промышленные обрабатывающие системы. Изначально в ППП в качестве языка программирования был выбран Fortran. Это обеспечивало транспортабельность пакета и возможность в настоящее время без особых трудностей использовать математическое обеспечение, ориентированное на вычисление конкретных специальных функций. Учёт этого, позволяет повысить производительность программ, что имеет большое значение само по себе и особенно возрастает в связи с бурным развитием объектно-ориентированного программирования и его использования для визуализации результата на экране. Моделирование волновых полей проводилось также на современных процессорах (Itanium, Opteron).

На основе разработанных программ исследованы динамические особенности распространения волн в средах с поглощением энергии. Если динамика волн в упругих средах довольно хорошо исследована, то для неупругих сред исследования носят в основном качественный характер [1],[126]-[129]. Экспериментальные исследования приводят, в основном, к независимости декремента поглощения от частоты и практической линейности коэффициента поглощения от частоты [130]-[134]. Эти два факта и взяты в основу при динамическом моделировании сейсмических процессов в неупругих средах.

Приведено количественное сравнение динамических характеристик сейсмических волн для широко распространённых в геофизике моделей вязкого трения, Максвелла, Гуревича, Дерягина для сред с разными декрементами затухания. Экспериментальные данные о поглощении в реальной среде свидетельствуют о большом разбросе as/af. Это отношение может быть как больше, так и меньше единицы [138]. Отметим многочисленные исследования дисперсионных уравнений для неупругих сред в случае различных законов поглощения (например, [134]). Основное отличие данной работы состоит в приведении количественных характеристик отношения доминирующих частот поперечных и продольных волн в зависимости от различных соотношений декрементов поглощения, согласующихся с экспериментальными данными.

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

Рассмотрен случай, когда источник колебаний, например, источник типа центра давления, расположен на расстоянии d<xP от свободной поверхности. В этом случае в полупространстве регистрируются наряду с прямой волной Р отражённые РР и PS и "нелучевая" волна s\ впервые обнаруженная в [141],[143]. Численное моделирование показало, что в упругой среде при малом заглублении источника спектр S' близок к спектру входного сигнала. С заглублением источника спектр S' монотонно сдвигается в сторону низких частот. Отметим, что численное моделирование подтверждает наличие "нелучевой" волны и в неупругом случае. При этом в неупругой среде спектр "нелучевой" волны сдвинут в сторону низких частот по сравнению с упругой и при заглублении источника так же наблюдается тенденция его сдвига в сторону низких частот. Объяснение этого факта получено из сравнения с динамикой релеевских волн.

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

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

Многочисленные экспериментальные наблюдения указывают на наличие поглощения [37] и анизотропии различных типов [38]-[40]. Несомненно, что эти факторы (за исключением структурных) оказывают наиболее существенное влияние на формирование волнового поля. В работе исследованы особенности распространения различных типов волн, а так же проведен анализ анизотропных (классификация М.В. Невского [145]) неупругих эффектов [42]. Обнаружено "нелучевое" явление, что угол раствора "петли" в волновом поле квазипоперечных волн на 8°-10° больше угла раствора петли на лучевых индикатрисах соответствующих волн для неупругих анизотропных сред. Исследованы волновые поля для анизотропной упругой и неупругой среды для модели Г.И. Петрашеня ([146]), характеризующейся двумя петлями по осям координат, с различными декрементами поглощения.

Среди систем возбуждения особое место по функциональным возможностям занимают невзрывные источники [147]-[148]. В работе приведены некоторые численные расчёты вибросейсмических полей и проведён анализ полученных результатов для реалистичной сейсмической модели Быстровского полигона (г. Новосибирск). Для этого построен алгоритм и написана программа моделирования виброграмм для реального свип-сигнала. Отметим качественное совпадение реальных и теоретических виброграмм [151].

В экспериментах с мощными вибрационными сейсмоисточниками на больших базах наблюдения впервые обнаружен эффект акустосейсмической индукции [154]. На расстоянии 20км были зарегистрированы поверхностные волны, возбуждённые акустическим излучением от вибрационного сейсмического источника в атмосфере. В работе данный практический эффект был объяснён численным моделированием в рамках комбинированной модели среды: атмосфера - упругая Земля с распределённым источником из [150].

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

Рассмотрены результаты численного моделирования для сложно построенных сред Юрубчено-Тахомской зоны (правобережье Енисея). В результате численного моделирования удалось добиться хорошего совпадения результатов реального и теоретического ВСП. Разработанные компьютерные технологии моделирования волновых полей для сложно построенных реальных сред впервые позволили выяснить (совместно с ОАО " Енисейгеофизика") причины отсутствия материалов на первичных сейсмограммах.

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

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

Разработанные алгоритмы и программы внедрены в Институте геофизики СО РАН, о также использовались и используются в ОАО "Енисейгеофизика" г. Красноярск, НИРФИ г. Горький, Институт вулканологии г. Петропавловск-Камчатский. Работа написана по публикациям: [12,16,17,19,20,34], [36,42,5358,63,64,69-71,78,94,101,102,108-112,154,162-167]. По материалам, используемым в диссертации, делались доклады на Международных конференциях: "Численные методы интерпретации сейсмических данных" (Суздаль, 1980), "Точные асимптотические и стохастические методы в геофизике" (Санкт-Петербург, 1981), "Модельная оптимизация в исследовательской геофизике"(Западный Берлин, 1991), "Математические методы в геофизике" (Новосибирск, 2003), "Международная конференция по вычислительной математике" (Новосибирск, 2004)," Вычислительные методы и решение оптимизационных задач" (Новосибирск, 2004), " The International symposium on mathematical mocteling of dynamic processes in atmosphere, ocean, and solid Earth" (Novosibirsk, 2004), "International Workshop on Active Monitoring in the solid Earth Geophysics (IWAM04)" (Japan, 2004). Полученные результаты обсуждались на семинарах Института вычислительной математики и математической геофизики и института геофизики СО РАН, Института физики земли РАН," С.-Петербургского отделения математического института РАН.

Автор выражает благодарность член-корр. РАН Б.Г. Михайленко за научные консультации и поддержку при выполнении работы, академику РАН А.С. Алексееву за внимание к работе и обсуждение основных результатов.

Похожие диссертационные работы по специальности «Математическое моделирование, численные методы и комплексы программ», 05.13.18 шифр ВАК

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Фатьянов, Алексей Геннадьевич

Основные результаты диссертации заключаются в следующем:

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

2. На основе численно-аналитического метода решена задача о распространении сейсмических волн в анизотропной неупругой среде. Получены дисперсионные соотношения для анизотропных коэффициентов поглощения. Численно исследована физика процесса распространения анизотропных неупругих волн.

3. Построена модель среды с разномасштабной неоднородностью. Получено замкнутое уравнение для среднего значения сферического поля и разработан метод его решения.

4. Развит метод расчёта многомерных волновых полей для сред с произвольной геометрией. Исследовано условие на ребре. Введена комплексная мощность слоя и на этой основе получено аналитическое решение в случае произвольных криволинейных границ. Проведено моделирование волновых полей в средах с шероховатыми границами.

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

Заключение

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

1. Шемякин Е.И. Динамические задачи теории упругости и пластичности. - Курс лекций для студентов. - НГУ, Новосибирск, 1968. -336с.

2. Ляв А. Математическая теория упругости. М.: ОНТИ, 1935. -678с

3. Алексеев А.С. Некоторые законы распространения волн в неоднородной среде. //ДАН СССР, 1955, т. 103, № 6. с. 989992.

4. Филиппов И.Г., Егорычев О.А. Нестационарные колебания и дифракция волн в акустических и упругих средах. М: Машиностроение, 1977. - 340с.

5. Алексеев А.С., Михайленко Б.Г. Метод вычисления теоретических сейсмограмм для сложно построенных моделей сред. //ДАН СССР, 1978, т. 240, № 5.-е. 1062-1065.

6. Снеддон И. Преобразование Фурье. М.: ИЛ, 1965. -684с.

7. Самарский А.А. Теория разностных схем. М.: Наука, 1977.-656с.

8. Рихтмайер Р., Мортон К. Разностные методы решения краевых задач. М.: Мир, 1972. -430с.

9. Воеводин В.В. Вычислительные методы линейной алгебры. М.: Наука, 1977. - 304с.

10. Толстов Г.П. Ряды Фурье. М.: Наука, 1980. - 382с.

11. Фатьянов А.Г., Михайленко Б.Г. Численное решение задачи Лэмба для неоднородной среды Больцмана с упругим последействием. // Математические методы интерпретации геофизических наблюдений. Новосибирск: Изд. ВЦ СО АН СССР, 1979.-с. 115-160.

12. Михайленко Б.Г. Решение задачи Лэмба для неоднородного полупространства с внутренними источниками // Математические проблемы геофизики, вып. 5, ч.1, Новосибирск, 1974. с. 187-194.

13. Гуревич Г.И. Деформируемость сред и распространение сейсмических волн. М.: Наука, 1974. -483с.

14. Степанов В.В. Курс дифференциальных уравнений. М.: ГОНТИ НКТП, 1938. -375с.

15. Фатьянов А.Г. Численное решение задачи Лэмба для неупругого неоднородного полупространства. // Численные методы в интерпретации геофизических наблюдений. -Новосибирск, 1980.-с. 144-156.

16. Фатьянов А.Г. Численное решение задачи Лэмба для вязкоупругого полупространства. В кн.: Численные методы в сейсмических исследованиях, Новосибирск: Наука, 1983. -с. 60-88.

17. Фатьянов А.Г., Михайленко Б.Г. Нестационарные сейсмические волновые поля в неоднородных вязкоупругихмоделях сред. // Математические проблемы геофизики: Модели и численные методы, Новосибирск, 1984. с. 82-131.

18. Михайленко Б.Г, Фатьянов А.Г. Полуаналитический метод расчёта нестационарных волновых полей для слоисто-однородных моделей сред. // Математические методы решения прямых и обратных задач геофизики, Новосибирск, 1981.-с. 92-104.

19. Петрашень Г.И. Распространение упругих волн в слоисто-изотропных средах, разделённых параллельными плоскостями. //Уч. Зап. ЛГУ, № 162, в. 26, 1952.-е. 3-189.

20. Петрашень Г.И. Общая количественная теория отражённых и головных волн, возбуждающихся в слоистых средах с плоско-параллельными границами раздела. // Вопросы динамической теории распространения сейсмических волн. Л., 1957. с. 70-164.

21. Thomson W.T. Transmission of elastic waves through a stratified solid medium. // J. Appl. Phys., 21, 1950. pp. 89-93.

22. Haskell N.A. The dispersion of surface waves a multilayered media. // Bull. Seis. Soc. Am., 1953. pp. 17-34.

23. Kennett B.L.N. Seismic wave propagation in stratified media. // Cambridge University Press, 1983.

24. Chapman C.H., Orcutt J.A. The computation of body wave synthetic seismograms in laterally homogeneous media. // Reviews of Geophysics, 23, 2, 1985. pp. 105-163.

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

26. Бабич В.М., Булдырёв B.C., Молотков И.А. Пространственно-временной лучевой метод. Л.: Изд. Лен. Универ., 1985.

27. Keller J.B. Surface waves on waves on water of non-uniform depth. // J. Appl. Phys., vol.31, № 6, 1960. pp. 1039-1046.

28. Tygel M., Hubral P. Transient waves in layered media. // Elsevier, Amsterdam. 1987.

29. Fuchs K., Muller G. Computation of synthetic seismogram with the reflectivity method and comparison with observation. // Geophys. J. R. Astr. Soc., 23, 1979. pp. 417-433.

30. Gerveny V., Molotkov L., Psensic I. Ray method in seismology. Charles University, Praha, 1977.

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

32. Фатьянов А. Г., Михайленко Б.Г. Метод расчёта нестационарных волновых полей в неупругих слоисто-неоднородных средах. //ДАН, 1988, т.301, №4. с. 834-839.

33. Аккуратов Г.В., Дмитриев В.И. Метод расчёта поля установившихся упругих колебаний в слоистой среде //В сб. Численные методы в геофизике, М.: МГУ, вып.2, 1979. с.З-12.

34. Ананьева Л.А., Фатьянов А.Г. Моделирование волновых неупругих полей в задачах вертикального сейсмического профилирования. // Математическое моделирование в геофизике. Новосибирск: Изд! ВЦ СО РАН, 1989.-е. 104-119.

35. Коган С.Я. Сейсмическая энергия и методы её определения. М.: Наука, 1975. - 286с.

36. Чесноков Е.М. Сейсмическая анизотропия верхней мантии Земли. М.: Наука, 1977. - 324с.

37. Crampin S. A review of wave motion in an isotropic and cracked elastic media. // Wave motion. Vol.3, № 4, 1981. p. 343-391.

38. Hoenig A. Elastic module of non-randomly cracked body. // Int.J.Solid structures, vol. 15, 1979. pp. 137-154.

39. Volterra V. Acta Math, vol.35, 1912.

40. Фатьянов А.Г. Нестационарные сейсмические волновые поля в неоднородных анизотропных средах с поглощением энергии. Новосибирск. 1989. - 44с. - ( Препринт / АН СССР Сиб Отд-ние ВЦ; 857)

41. Мартынов В.Н. Волновые поля от сосредоточенных источников в трансверсально-изотропных средах. // Из. АН СССР, Сер. Физика Земли, № 11, 1986.-е. 19-26.

42. Коган С.Я. Краткий обзор теории поглощения сейсмических волн. // Изв. АН СССР, Сер. Физика Земли, № 11, 1966.-е. 1-16.

43. Азими Ш.А. и др. Импульсные и переходные характеристики сред с линейными и квадратичными законами поглощения. // Изв. АН СССР, Сер. Физика Земли, № 2, 1968. с. 42-54.

44. Гольдин С.В. Линейные преобразования сейсмических сигналов. М.: Наука, 1974. - 320с.

45. Дерягин В.В. //Журнал геофизики, № 1-2, 207, 1931.

46. Петрашень Г.И. Распространение волн в анизотропных упругих средах. -П.: Наука, 1980. -420с.

47. Love А.Е.Н. Some problems of geodynamics. Cambridge, 1931.

48. Молоденский M.C. Общая теория упругих колебаний Земли. М.: Недра, 1989. - 78с.

49. Pekeris C.L., Jarosch Н., Alterman Z. Oscillations of the Earth. // Second Intern Report, The Weitzman Institute, Reboot, Israel, 1959.

50. Галицин А.В., Жуковский А.Н. Интегральные преобразования и специальные функции в задачах теплопроводности. Киев: Наукова думка, 1976.-282с.

51. Конюх Г.В., Фатьянов А.Г. Полуаналитический метод расчёта волновых полей для слоистой сферической модели Земли. // Математическое моделирование в геофизике. -Новосибирск, 1993. с. 58-77.

52. Конюх Г.В., Фатьянов А.Г. Волновые поля в гравитирующей слоистой Земле. // ДАН, т. 6, №2, 1992. с. 264-267.

53. Фатьянов А.Г. Численное моделирование волновых полей в неоднородном неупругом шаре. Новосибирск, 1981. - 22с. - ( Препринт / АН СССР Сиб Отд-ние ВЦ; 337)

54. Fatyanov A.G., Soboleva O.N. Numerical modeling of the gravitating liquid layer on the elastic half space. // International tsunami symposium, Novosibirsk, 1990. p. 29-34.

55. Fatyanov A.G., Soboleva O.N. Numerical modeling of the complex wave field in the gravitating liquid layer on the elastic half space. // International tsunami meetings, Novosibirsk, 1989. -91p.

56. Держи H.M., Фатьянов А.Г. Пакет программ для расчёта геофизических полей на персональных ЭВМ. // Математическое моделирование в геофизике. Новосибирск. 1993.-с. 78-89.

57. Biot М.А. Theory of propagation of elastic waves in a fluid saturated porous solid. // J. Acoustical Soc. America, vol. 28, № 2, 1956.-p. 168-191.

58. Френкель Я.И. К теории сейсмических и сейсмоэлектрических явлений во влажной почве. // Изв. АН СССР. Сер. Геогр. и геофиз. Т.8, № 4, 1944. с. 133-150.

59. Уайт Д.Э Возбуждение и распространение сейсмическихволн. М.: Недра, 1986.62. Николаевский В.Н., Басниев К.С., Горбунов А.Т., Зотов Г.А. Механика насыщенных пористых сред. - М.: Недра, 1986.

60. Мирошников В.В., Фатьянов А.Г. Численное моделирование волновых полей в пористой среде. Модель Био // Математическое моделирование в геофизике, Новосибирск, 1989. с. 83-103.

61. Мирошников В.В., Фатьянов А.Г. Полуаналитический метод расчёта волновых полей в слоистых пористых средах // Математическое моделирование в геофизике, Новосибирск, 1993.-с. 27-58.

62. Deresiwicz Н., Skalak R. On uniqueness in dynamic poroelasticity. // Bull of the seismological. Soc. America, vol.53, №4, p. 783-788.

63. Plona T.J. Observation of a second bulk compression wave in a porous medium at ultrasonic frequencies // Appl. Phys. Let, vol.36, №4, 1980. p. 259-261.

64. White J.E. Seismic reflections from gas reservoirs final repot // National Science Foundation, Contract № AER-17526, 1977.

65. Блохин A.M., Доровский B.H. Проблемы математического моделирования в теории многоскоростного континуума. -Новосибирск, 1994. 186с.

66. Фатьянов А.Г. Прямые и обратные задачи для тензора сейсмического момента в слоистых средах. // ДАН, 1991, т. 317, №6, 1991.-с. 1357-1361.

67. Фатьянов А.Г. Полуаналитический метод решения прямых динамических задач в слоистых средах. // ДАН, т. 310, №2, 1990.-с. 323-327.

68. Матевосян А.Х., Фатьянов А.Г. Метод расчёта волновых полей для одной модели сейсмического объёмного источника. Новосибирск. 1994. - 20с. - ( Препринт / АН СССР Сиб Отд-ние ВЦ ;1029)

69. Смирнов В.И., Соболев С.Л. Об одном новом методе в плоской задаче упругих колебаний. // Тр. Сейсмологического института. 1932, №20. - 32с.

70. Смирнов В.И., Соболев С.Л. О применении нового метода к изучению упругих колебаний в пространстве при наличии осевой симметрии. // Тр. Сейсмологического института. -1933, №29. -49с.

71. Зволинский Н.В. Отражённые и головные волны, возникающие на плоской границе раздела. // Изв. АН СССР, Сер. Геофиз. 1957, №10. - с. 1201-1218.

72. Коновалов А.Н. Численное решение задач теории упругости в напряжениях: Учеб. Пособие. Новосибирск: Изд. НГУ, 1979.-92с.

73. Козлов Е.А. Распознавание и подавление многократных волн в сёйсморазведке. М: Недра, 1982. -248с.

74. Алексеева М.В. Математическое моделирование сейсмического поля в многослойной упругой среде (в рамках лучевого метода). Новосибирск. 1987. - 60с. - ( Препринт / АН СССР Сиб Отд-ние ВЦ;729)

75. Fatyanov A.G., Mikhailenko B.G. Numerically-analytical method for calculation of theoretical seismograms in layered-inhomogeneous inelastic media. // Geophysical data inversion methods and applications, free University of Berlin, 1989. p. 499-530.

76. Алексеев A.C., Бабич B.M., Гельчинский Б.Я. Лучевой метод вычисления интенсивности волновых фронтов. //

77. Вопросы динамической теории распространения сейсмических волн. Л.: Гостоптехиздат, №5, 1961. - с. 3-24.

78. Бабич В.М., Булдырёв B.C. Асимптотические методы в задачах дифракции коротких волн. М: Наука, 1982. -248с.

79. Popov М.М. A new method of computation of wave fields in the high-frequency approximation. Leningrad, 1981. - 20p. -(Preprint/AN SSSR)

80. Клем-Мусатов К.Д. Теория краевых волн и её применение в сейсмике. Новосибирск: Наука, 1980.

81. Keller J.B. Geometrical theory of diffraction // J. Opt. Soc. Am., vol.52, №4, 1981.-p. 175-188.

82. Kelly K.R. et al. Synthetic seismograms: a finite- difference approach // Geophysics, №41, 1976. p. 2-27.

83. Zahradnik J. Finite- difference solutions to certain diffraction problems // Stud, and Geod., vol. 19., 1975. p. 233-244.

84. Бреховских Л.М. Волны в слоистых средах. М.: Наука, 1973.-344с.

85. Молотков Л.А. Бакулин А.В. Эффективная модель слоистой упруго-жидкой среды как частный случай модели Био. // Математические вопросы теории распространения волн Санкт-Петербург: Наука, т.25, 1995. - с. 172-195.

86. Крауклис П.В., Цепелев Н.В. Распространение колебаний вблизи границы неоднородного упругого полупространства сфазовой скоростью, близкой к скорости продольных волн. // Изв. АН СССР. Физика Земли, №9, 1971. с. 28-33.

87. Михайленко Б.Г. Сейсмические поля в сложно построенных средах Новосибирск, 1988, 312с.

88. Mikhailenko B.G. Spectral Laguerre method for the approximate solution of time dependent problems // Appl. Math. Lett, vol.12, 1999.-p. 105-110.

89. Мартынов B.H., Михайленко Б.Г. Численное моделирование волновых полей в анизотропных средах. // Математические методы в геофизике. Новосибирск, ч. 1, 2003.-с. 137-141.

90. Фёдоров Ф.И. Теория упругих волн в кристаллах. М.: Наука, 1965. -388с.

91. Фатьянов А.Г., Мирошников В.В. Энергетический метод расчёта функции Грина в многомерно-неоднородных средах. //ДАН, т. 351, №2, 1996.- с. 264-266.

92. Беллман Р. Введение в теорию матриц. М.: Наука, 1976. - 352с.

93. Миттра Р., Ли С. Аналитические методы теории волноводов. М.: Мир, 1974. - 328с.

94. Вайнштейн Л.А. Теория дифракции и метод факторизации. -М.: 1966.-431с.

95. Демидович Б.П. Марон И.А. Основы вычислительной математики М.: Наука, 1970. - 664с.

96. Франк Ф., Мизес Р. Дифференциальные и интегральные уравнения математической физики. М.: Гостехиздат, 1937. -998с.

97. Николаев А.В. Сейсмика неоднородных и мутных сред. -М.: Наука, 1972.-324с.

98. Фатьянов А.Г. Волновые поля в многомерно-неоднородных средах. // Математические методы в геофизике. Новосибирск, ч. 1,2003.-е. 142-145.

99. Ананьева Л.А., Фатьянов А.Г. Численное моделирование волновых полей в средах с цилиндрическими границами разделов. // Математические проблемы геофизики; Численное исследование геофизических задач. Новосибирск, 1988. - с. 105-121.

100. Алексеев А.С. Некоторые методы и алгоритмы интерпретации данных. М.: Наука, 1967. - с. 9-84.

101. Tarantola A.A strategy for nonlinear elastic inversion of seismic reflection data // Geophysics, Vol. 51, 1986. p. 1893-1903.

102. Santosa F., Symes W. An analysis of least-squares velocity inversion // Geophysical Monograph Series, №4, SEG, 1989. -154p.

103. Алексеев A.C., Добринский В.И Некоторые вопросы практического использования обратных динамических задач сейсмики // Математические проблемы геофизики. -Новосибирск: ВЦ СО АН СССР, 1975. Вып. 6, ч. 2. - с. 7-53.

104. Алексеев А.С., Цибульчик Г.М. О связи обратных задач теории распространения волн с задачами визуализации волновых полей // ДАН СССР, т. 242, №5, 1978.- с. 1030-1033.

105. А.С. Алексеев, А.В. Авдеев, А.Г. Фатьянов, В.А. Чеверда Замкнутый цикл математического моделирования волновых процессов в вертикально-неоднородных средах (прямые иобратные задачи) // Математическое моделирование, Т. 3, №10, 1991.-с. 80-94.

106. A.S. Alekseev, A.V. Avdeev, A.G. Fatyanov, V.A. Tcheverda Wave processes in vertically-inhomogeneous media: a new strategy for a velocity inversion // J. Inverse Problems. 1993. - Vol. 9, №3.- p. 367-390.

107. А. С. Алексеев, A.B. Авдеев, А.Г. Фатьянов, B.A. Чеверда Волновые процессы в вертикально-неоднородных средах: прямые и обратные задачи. Новосибирск. 1991. - 44с. - ( Препринт / АН СССР Сиб Отд-ние ВЦ;924)

108. Карчевский А.Л. , Фатьянов А.Г. Численное решение обратной задачи для системы упругости с последействием для вертикально неоднородной среды. // Сиб. журн. вычисл. математики, т.4, №3,2001. -с.259-268.

109. С.И. Кабанихин, К.Т. Искаков Обобщённое решение обратной задачи для уравнения колебания. // ДАН, т. 375, №1, 2000.- с. 22-24.

110. Marquardt D.W. An algorithm for least-squares estimation of nonlinear parameters//J. Soc. Indust. Appl. Math. 1963. - Vol. 11.- p.431-441.

111. Бухгейм А.Л. Разностные методы решения некорректных задач. Новосибирск: Изд. ВЦ СО АН СССР, 1986. - 148с.

112. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1972. -496с.

113. Васильев Ф.П. Численные методы решения экстремальных задач. Москва: Наука, 1988. - 549с.

114. Фадеев Д.К., Фадеева В.Н. Вычислительные методы линейной алгебры. М: Физматгиз, 1963.

115. Алексеев А.С., Кабанихин С.И. Обратные задачи и новые технологии в геофизике // Математические методы в геофизике. Новосибирск, ч. 1, 2003. - с. 11 -20.

116. Марчук Г.И. Методы вычислительной математики. М.: Наука, 1977.-456с.

117. Ерохин Г.Н., Бортников П.Б. Обратная задача определения тензора сейсмического момента. // Геология и Геофизика, №4, 1987.-с. 115-123.

118. Aki К. Pt.1. Generation and propagation of G waves from the Niidata earthquake of June 16, 1964. Pt.2. Estimation of earthquake moment, released energy and stress train drop from G - wave spectrum. //Bull. Earth. Res. Inst. Tokyo. Univ., №4, 1966.

119. Dziewonski A.M., Woodhause J.H. An experiment in systematic study of global seism city centroid-moment tensor solutions for 201 moderate and large earthquakes of 1981. // J. Geophys. Res., №88, 1983.

120. Backus G., Mulkahy M. Moment tensors and other phenomenological descriptions of seismic sources. I. Continuous displacements. // Geophys. J. Roy. Astron. Soc., №46, 1976. -p.341-361.

121. Костров Б.В. Механика очага землетрясения. М.: Наука, 1975.-176с.

122. Николаев Б.Г. О распространении нестационарных возмущений в неидеально-упругих средах. // Сб.: Вопросы динамической теории распространения сейсмических волн, ЛГУ, В.З, 1959.

123. Брук С.З. Задача Лэмба для вязкоупругой полуплоскости. // Механика твёрдого тела, №3, 1972. с. 56-63.

124. Блитштейн Ю.М., Мешков С.И., Чебан В.Г., Чигарев А.В. Распространение волн в вязкоупругих средах. Кишинёв: Штиница, 1977. -206с.

125. Шемякин Е.И. Задача Лэмба для среды с упругим последействием. //Докл. АН СССР, т. 104, №2, 1955. с. 193-196.

126. Attewell Р.В., Ramana J.V. Wave attenuation and internal friction as functions of frequency in rocks. // Geophysics, 31, 1049, 1966.

127. Knopoff L., // Rev Geophysics, 2, №4, 1964.

128. Пасечник И.Н. Характеристики сейсмических волн при ядерных взрывах и землетрясениях. М.: Наука, 1972. - 192с.

129. Берзон И.С., Пасечник И.Н., Поликарпов A.M. Определение параметров затухания Р волн в мантии Земли. // Из. АН СССР, Физика Земли, №2, 1975.

130. Кондратьев O.K. Сейсмические волны в поглощающих средах. М.: Недра, 1986. - 176с.

131. Пузырёв Н.Н., Худобина Л.Н. Обзор экспериментальных исследований и некоторые теоретические положения по изучению поперечных и обменных волн. // Сб.: Экспериментальные исследования поперечных и объёмных волн, Новосибирск, 1962.

132. Пузырёв Н.Н. Интерпретация данных сейсморазведки методом отражённых волн. М.: Гостоптехиздат, 1959. -434с.

133. Лёвшин А.Л., Ратникова Л.И., Сакс М.В. О дисперсии и поглощении упругих волн в горных породах. // Вычислительная сейсмология, Москва, №13, 1981.-е. 134-142.

134. Николаев А.В Сейсмические свойства грунтов. М.: Наука, 1965.-120с.

135. Саваренский Е.Ф. Сейсмические волны. М.: Наука, 1972. -296с.

136. Гурвич И.И. Сейсмическая разведка. М.: Недра, 1970. -552с.

137. Hron F., Mikhailenko B.G. Numerical modeling of non-geometrical effects by the Alekseev-Mikhailenko method // Bulletin of the seismological Society of America, v.71, №4, 1981. p. 10111029.

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

139. Алексеев А.С., Михайленко Б.Г. "Нелучевые" эффекты в теории распространения сейсмических волн. //Докл. АН СССР, т.267, №5, 1982.-с. 1079-1084.

140. Ляховицкий Ф.М. Проблемы сейсмической анизотропии тонкослоистых сред // Вестник МГУ, Геология, №5, 1983. с.68-75.

141. Невский М.В. Квазианизотропия скоростей сейсмических волн. М.: Наука, 1974.

142. Николаев А.В. Вибрационное просвечивание метод исследования Земли. - В кн.: Проблемы вибрационного просвечивания Земли. М.: Наука, 1977. - с.5-14.

143. Гущин В.В., Докучаев В.П., Заславский Ю.М., Конюхова И.Д. О распределении мощности между различными типами излучаемых волн в полубезграничной упругой среде. В кн.: Проблемы вибрационного просвечивания Земли. М.: Наука, 1977. -с.113-118.

144. Чичинин И.С. Вибрационное излучение сейсмических волн -М.: Недра, 1984.

145. Алексеев А.С. Развитие теории и методики активного вибросейсмического мониторинга зон аккумуляции напряжений (Проект 1.3 ГНТП "Глобальные изменения природной среды и климата"), (отчёт). Новосибирск, 1981. - 100с.

146. Gupta I.N., Hartenberger R.A. // Bull. Seismol. Soc. Amer., v.71, №6, 1981.-p. 1731-1744.

147. Разин A.B. // Изв. Академии Наук. Физика Земли, №2, 1993. с.73-77.

148. А.С. Алексеев, Б.М. Глинский, С.И. Дряхлов, В.В. Ковалевский, Б.Г. Михайленко, Б.М. Пушной, А.Г. Фатьянов, М.С. Хайретдинов, М.Н. Шорохов. Эффект акустосейсмической индукции при вибросейсмическом зондировании. //Докл. РАН, т.346, №5, 1995. -с.664-667.

149. Alekseev A.S., Kovalevsky V.V. Powerful vibrators for deep interior investigation. // 60 Annu. Intern. Meeting Soc. of Exploration Geophysicists. San-Francisco, 1990. p.956-957.

150. Крауклис П.В., Крауклис J1.A. Интерференционная медленная волна в пороакустическом слое Био. // Математические вопросы теории распространения волн. -Санкт-Петербург, №28, 1999.-е. 137-148.

151. Киселёв А.П. Высшие приближения лучевого метода и "нелучевые явления" в неоднородных вязкоупругих средах // Изв. АН СССР, Физика Земли, №11, 1992. с.35-38.

152. Kosloff D., Reshef M., Loewenthal D. Elastic wave calculations by the Fourier method. Bull. Seis. Soc. Am., 1984, 74, P. 875 891.

153. Молотков Л.А. О коэффициентах извилистости пор в эффективной модели Био. // Математические вопросы теории распространения волн. Санкт-Петербург, №28, 1999. - с. 157164.

154. Кляцкин В.И. Стохастические уравнения и волны в случайно-неоднородных средах. М.: Наука, 1980. - 336с.

155. Сибиряков Б.П. Геометрия трещиноватых сред и параметрические резонансы // Геология и Геофизика, т.43, №9, 2002. с.882-887.

156. Glinsky В.М., Fatyanov A.G. Vibroseismic monitoring of active volcanos. // The International symposium on mathematical modeling of dynamic processes in atmosphere, ocean, and solid Earth. -Novosibirsk, 2004 p. 142-145.

157. Фатьянов А.Г. Моделирование волновых полей в средах со стохастическими шероховатыми границами. // МНПК "Связь-2004". Международный семинар "Вычислительные методы и решение оптимизационных задач". Новосибирск, 2004. - с. 173177.

158. Фатьянов А.Г. Волновые поля в средах с криволинейными границами. // Международная конференция по вычислительной математике. Новосибирск, ч. 2, 2004. - с. 734-739.

159. Fatyanov A.G. Mathematical modeling of wave fields in media with curvilinear boundaries // Doklady Earth Sciences, v.401a, №3, 2005. p.437-439.

160. Фатьянов А.Г. Математическое моделирование волновых полей в средах с криволинейными границами // Докл. РАН, т.401, № 4, 2005. с.529-532.

161. Исимару А. Распространение и рассеяние волн в случайно-неоднородных средах. М.: Мир, т.2, 1981. -317с.

162. Клаербоут Д.Ф. Сейсмическое изображение земных недр. -М.: Недра, 1989.-407с.1. Оорма 51. УТВЕРлЗШЭ

163. Jhzракгох^Кычиатаголышго иои-езя^б^^З^СССг'угвззргдли

164. Директор Института гоа/югии я хчкхЗившш ie?.j. CU-ciokw Cojoaa1. All СССР- \ i'.^tjjyieKceeii »•1. JDO-} 14 ^ Н-Д-ДоДтепов1. V'' -'' v" t^C wлгсг о'жЕ5£Шса11. Эв9г.

165. Научио-исолодопата.'хъаиая (отлио-иопструдторокая) рпс5ота <раора<Зопга) разработка подупнали тиче с кого метода расчётаволновых полей в слоистых средах, выполненная по п.XX НИР

166. Г!*,а шелроши передача комплекса программ; проведение консультаций по его эксплуатации и анализ результатов.

167. Достигнутые топипеокио розу.ях<таты проведены численные расчёт» для реалистичинх моделей сред; исследована (Ткзика процесса.

168. Х*вагсгивиш.1 исис-всгтодь С.п.сА.Г.^птьдпов1. Oi йи/аачтеа;01*. И. C4/VT/>. K.fJ . сна uu^t-tv^f* d esse? г с

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