Математическое моделирование землетрясений Кыргызстана и алгоритм определения обратной задачи | Статья в журнале «Молодой ученый»

Отправьте статью сегодня! Журнал выйдет 20 марта, печатный экземпляр отправим 24 марта.

Опубликовать статью в журнале

Библиографическое описание:

Алимканов, А. А. Математическое моделирование землетрясений Кыргызстана и алгоритм определения обратной задачи / А. А. Алимканов, А. Д. Сатыбаев, С. Н. Касымбеков. — Текст : непосредственный // Молодой ученый. — 2020. — № 35 (325). — С. 6-11. — URL: https://moluch.ru/archive/325/73350/ (дата обращения: 06.03.2021).



Целью настоящей работы является построение математической модели землетрясений, и разработать численный метод решения обратной задачи, построить алгоритм решения поставленной задачи.

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

Одним из стихийных бедствий человечества является землетрясение, и оно происходит ежедневно — ощутимое и неощутимое [1, 2], и происходит около 3000 раз в день [3].

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

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

На активных участках Юге Кыргызстана также расположены сейсмические станции — ЕССН, которые ведут записи магнитуды сейсмического процесса, параметры земной среды при ощущении амплитуды землетрясения.

Известно, что Центральная Азия находится под угрозой землетрясения и задачи возникновения землетрясений и их решения до сих пор остаются не доконца решенными и этим занимаются многие ведущие ученые [4].

Постановка задачи. Процесс природного явления — землетрясение описывается уравнением сейсмики. Рассмотрим обратную задачу сейсмики следующего вида [5]:

(1)

(2)

(3)

где - плотность среды, - коэффициент Ламэ, - давление среды, - дельта-функция Дирака, - тета-функция Хевисайда, - положительные постоянные числа.

Обратная задача заключается в определении функции или из задачи (1)-(3) при задании другого коэффициента уравнения.

Решение. Обратная задача (1)-(3) с применением методов выпрямления характеристик и выделения особенностей приведена к обратной задаче с данными на характеристиках, и последняя задача решена численным конечно-разностным методом, затем решена первоначальная обратная задача (1)-(3).

Смещения от сейсмического источника предоставляет с собой функция Грина динамической теории упругости [6].

Пусть единичный пульс приложен в точке и в момент времени в направлении . Тогда функция Грина будет компонентной смещения поверхности Земли.

А математическая модель смещения почвы при землетрясении задается уравнением [7]:

(4)

гдеρ — плотность среды.

Чтоб однозначно определить решение уравнения (1) необходимо должны задавать начальные и граничные условия.

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

(5)

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

Рассмотрим в упругом теле объемом два поле смещений при различных начальных условиях в момент времени (см. табл. 1):

Таблица 1

Два поля смещений при начальных условиях

п / п

Поля смещения

Объемные силы

Граничные условия

Напряжение на поверхностях

1

2

Пусть существует момент времени t 0 , до которого , в , следовательно , тогда свертка равно нулю, т. е. (6)

Пусть поле смещения u(х,t) удовлетворяют условию

(7)

х,μ — коэффициенты Ламэ, f — сила, - оператор.

Тогда существует потенциалы для , обладающие следующими потенциалами:

; (8)

; (9)

, , (10)

, , (11)

и называется соответственно и — компонентами поле смещения .

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

— модуль Юнга, связь между нормальным напряжением и продольной деформацией в стержне;

— коэффициент Пуассона, отношение поперечной и продольной деформацией стержня при его продольном расторжении или сжатии.

— коэффициент Ламэ — связь между деформацией скошения прямого угла;

— коэффициент Ламэ — связь между деформациями сжатия — расширения и нормальными напряжениями.

Приводим зависимости модулей для изотропной среды:

; (12)

; (13)

. (14)

где — модуль объемного сжатия.

Между дилатацией и составляющими внешних массовых сил X, Y, Z существует связь:

; (15)

(16)

; (17)

.

Отсутствие внешних сил следует X=Y=Z=0, а составляющие вектора определяются выражениями:

; ; , (18)

тогда из (10) — (12) следует два уравнения

; (19)

; (20)

Уравнение (15) — (16) являются волновыми уравнениями и описывают распространения продольных и поперечных волн соответственно, скорости определяются по формулами:

; ; (21)

Отметим, что

. (22)

В совокупности продольные и поперечные волны называют объемными волнами.

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

Материальный ущерб от землетрясения в Кыргызстане составляет, за последние 15 лет, составляет около 105 млн.$, а землетрясения в большинство случаях происходят в Кыргызстане, где проживают 51 % населения Кыргызстана.

Это означает актуальность, необходимость изучения землетрясений Кыргызстана.

С точки зрения практики авторы считают, что наиболее приемлемым методом решения этих задач является, конечно-разностный метод. В этом случае, конечно, необходимо установить устойчивость решения, т. е. уравнение движения идеально упругой изотропной среды (13) — (15) можно записать в виде [8,9]:

; (23)

где , , ,

Уравнение (21) в однородной среде имеет вид

, (24)

где .

Из уравнения движения однородной среды (22) в некоторых упрощениях можно получить уравнение сейсмических волн в двумерном случае, т. е. уравнение зависит от двух пространственных переменных,

, (25)

Для задания начальных и граничных условий рассмотрим широко распространенная в геофизике модель среды, состоящих из двух полупространств x>0 и x<0 и с границей на плоскости x=0. Здесь предположим, что коэффициенты Ламэ и плотность среды гладки в этих полупространствах и имеет конечный скачок при переходе из одного полупространства в другое полупространства.

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

Литература:

  1. Потапов А. Д., Ревелис И. Л., Чернышев С. Н. Землетрясения. Причины, последствия и обеспечение безопасности. М.: Инфра-М, 2017. — С.344.
  2. Никонов А. А. Землетрясение. М.: Знание. 1984. — С.192.
  3. Уломов В. И. Динамика Земной коры Средней Азии и прогноз землетрясений. Ташкент: ФАН. 1974. — С.218.
  4. Аки И., Ричардс П. Количественная сейсмология. — М.: Мир. 1983. — С.520.
  5. Сатыбаев А.Дж., Алимканов А. А., Култаев Т. Ч. Алгоритм определение одной обратной задачи сейсмики с мгновенным и шнуровым источниками. Бишкек: Известия КГТУ им. И.Раззакова. 2016, № 3 (39), часть I. с. 175–180.
  6. Романов В. Г. Обратные задачи математической физики. — М.: Наука, 1984г. — 264 с.
  7. Сатыбаев А.Дж. Численное определение коэффициента Ламэ в уравнении сейсмики //История, культура и экономика, наука юга Кыргызстана. — Ош: КУУ, 2000.- С.148–152.
  8. Кабанихин С. И. Обратные и некорректные задачи. Новосибирск: Сибирское научное издательство. 2009г. — 457с.
  9. Яхно В. Г. Обратные задачи для дифференциальных уравнений упругости. Новосибирск: Наука, Сибирское отделение, 1990г. — 304 с.
Основные термины (генерируются автоматически): обратная задача, уравнение, момент времени, плотность среды, динамическая теория упругости, изотропная среда, конечно-разностный метод, Кыргызстан, математическая модель землетрясений, однородная среда.


Ключевые слова

обратная задача, алгоритм решения, математическая модель землетрясений, уравнение сейсмики, источники данных, численный, конечно-разностный метод, анализ методов
Задать вопрос