Математическое моделирование чрезвычайных ситуаций
Основные параметры, характеризующие землетрясения. Схема землетрясения и вызываемых им волн. Скорость распространения продольных сейсмических волн. Этапы построения математической модели опасного явления. Проверка адекватности математической модели.
Рубрика | Математика |
Вид | контрольная работа |
Язык | русский |
Дата добавления | 18.02.2020 |
Размер файла | 201,3 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Содержание
1. Перечислить и объяснить основные параметры, характеризующие землетрясения
2. Перечислить и объяснить основные этапы построения математической модели опасного явления
3. Практическое задание
Список литературы
1. Перечислить и объяснить основные параметры, характеризующие землетрясения
Землетрясения - это подземные толчки и колебания земной поверхности, возникающие в результате внезапных смещений и разрывов в земной коре.
Основной причиной наиболее разрушительных землетрясений является столкновение или иное взаимодействие тектонических плит.
Земля имеет слоистое строение: внешний твердый слой - земная кора имеет толщину от 6 км (под дном океанов) до 70 км (в горных районах континентов). Ниже коры до глубины приблизительно 2900 км простирается мантия. Центральную часть Земли от глубины приблизительно 2900 км до центра занимает ядро. Твердая оболочка земли (литосфера), включающая в себя кору и верхний слой мантии, состоит из семи-восьми гигантских плит, на которых находятся континенты и несколько десятков более мелких плит, которые перемещаются по слою сильно разогретых пластичных пород.
Толщина тектонических плит меняется от 10-80 км под океанами и 200-300 км на континентах. Скорость смещения плит - несколько сантиметров в год. Причиной этого смещения или драгера плит является вращение Земли и другие геофизические и космические факторы. Плиты раздвигаются, сдвигаются, наползают друг на друга и т.д. При взаимодействии плит возникают их упругие деформации, которые при достижении предела прочности горных пород приводят к их разрывам и разрушениям, что и проявляется в виде землетрясений.
Около 90 % землетрясений приходится на Тихоокеанский пояс, окружающий Тихий океан; 5-6 % приходится на Альпийский пояс, который тянется от Средиземного моря через Турцию, Иран, Северную Индию в Китай. Третья зона землетрясений - срединно-океанические хребты Атлантического и Индийского океанов. В России в сейсмоопасные зоны попадают около 20 % всей территории.
Источник сейсмических волн занимает определенный объем и имеет протяженность от 1 до 100 и более км. Но при теоретическом рассмотрении источник землетрясения представляется точкой, которая называется гипоцентром (фокусом).
Проекция гипоцентра на поверхность Земли вверх по вертикали называется эпицентром (рис 1). Местоположение гипо- и эпицентра определяют по времени прихода сейсмических волн.
Рис. 1. Схема землетрясения и вызываемых им волн: 1 - гипоцентр; 2 - эпицентр; 3 - Р-волны; 4 - S-волны; 5 - волны Лява; 6 - волны Релея
Землетрясение можно определить как колебания земной поверхности при прохождении волн от подземного источника энергии.
Различают несколько видов сейсмических волн. В толще горной породы от источника землетрясения распространяются первичные продольные волны (волны давления или Р-волны) и поперечные волны (волны смещения или S-волны).
Скорости распространения этих волн зависят от механических свойств горных пород и грунтов, через которые они проходят (табл. 1). При этом скорость продольных волн в среднем в 1,7 раза выше скорости поперечных волн. Объемные продольные и поперечные волны достигнув поверхности Земли вызывают поверхностные волны двух видов: волны Лява, в которых грунт смещается из стороны в сторону в горизонтальной плоскости перпендикулярно к направлению распространения и волны Релея, в которых частицы грунта перемещаются по замкнутым траекториям вверх и назад по отношению к направлению движения волны. Поверхностные волны распространяются примерно в два раза медленнее, чем поперечные.
Таблица 1 - Скорость распространения продольных сейсмических волн
Тип горной породы или грунта |
Скорость распространения Vp, км/с |
|
Граниты, базальты |
6,9 |
|
Осадочные породы |
6,1 |
|
Песчаники, известняки |
1,5 - 5,6 |
|
Полускальные грунты |
1,4 - 3,6 |
|
Крупнообломочные грунты |
1,1- 2,1 |
|
Песок |
0,7 - 1,6 |
|
Глина, суглинок |
0,5 - 1,5 |
|
Насыпные грунты |
0,2 - 0,5 |
Воздействиям продольных, поперечных и поверхностных волн соответствуют три фазы землетрясения: Р-, S- и L-фазы. Действие продольных волн (Р-фаза) похоже на удар воздушной волны. Самой разрушительной является вторая S-фаза, связанная с приходом поперечных волн, которые смещают грунт из стороны в сторону как по вертикали, так и по горизонтали.
Амплитуда колебаний земной поверхности А, мкм на расстоянии R км от эпицентра определяется следующей приближенной формулой:
Z = 10М-1,32?lgR
где М - магнитуда землетрясения, которая является условной единицей его полной энергии.
Автор понятия магнитуда Ч. Рихтер (1935) определял её как число пропорциональное десятичному логарифму амплитуды наиболее сильной волны, выпаженной в мкм, на расстоянии 100 км от эпицентра землетрясения.
Энергия, выделяющаяся при землетрясении, определяется формулой:
Е = 104,8+1,5М, Дж
Магнитуда землетрясений по шкале Рихтера может изменяться от 1 до 9. Самое сильное Ассамское землетрясение (северо-восточная Индия, 1952) имело магнитуду 8,7. Выделенная в нем громадная энергия Е ~ 1018 Дж эквивалентна энергии, выделенной при взрыве 2,2 млрд т тротила.
В среднем в мире ежегодно происходит одно землетрясение с магнитудой 8 и выше.
Интенсивность землетрясения - это мера величины сотрясения грунта, которая определяется степенью разрушения зданий и конструкций, а также характером изменения земной поверхности.
Интенсивность землетрясения в зависимости от расстояния различна, а его магнитуда всегда одна.
В РФ, как и в ряде европейских стран, используется 12-бальная международная шкала интенсивности землетрясения МSК-64. Условно землетрясения подразделяются на слабые (1-3 балла), умеренные (4 балла), сильные (5-7 баллов), разрушительные - 8 баллов, опустошительные - 9 баллов, уничтожительные - 10 баллов; катастрофические - 11 баллов; 12 - баллов сильная катастрофа.
Интенсивность землетрясения J определяется по формуле:
где расстояние от эпицентра R и глубина h измеряются в км.
Из формулы следует выражение для максимальной интенсивности землетрясения, которая достигается в эпицентре (R = 0):
J0=3+1,5V-3,51gh
Если данные о глубине гипоцентра отсутствует, то интенсивность землетрясения можно оценить по формуле:
J0 ? 1,5 (M ? 1).
Расстояние, на котором возможно возникновение колебаний заданной интенсивности, определяется соотношением
Реальная интенсивность землетрясения, учитывающая влияние грунта под застройкой и в окружающей местности, определяется соотношением:
где ?Jз, ?Jм - изменение балльности землетрясения, учитывающие типы грунта под зданиями и в окружающей местности (по сравнению с гранитом).
2. Перечислить и объяснить основные этапы построения математической модели опасного явления
Характерные особенности чрезвычайных ситуаций (ЧС), такие как внезапность возникновения, быстрота развития, неполнота и неопределенность исходной информации, разнообразие и цепной характер последствий затрудняют использование для их изучения традиционных эмпирических методов.
В связи с этим для анализа и прогнозирования чрезвычайных ситуаций все шире применяется математическое моделирование, которое является во многих случаях единственно допустимым, как например, при экспертизе особо опасных природных и техногенных явлений.
Математической моделью ЧС называется система соотношений, уравнений, неравенств, геометрических понятий и т.д., которые в математической форме отображают, воспроизводят или имитируют наиболее важные особенности и свойства реальных опасных явлений с целью анализа и прогнозирования их возникновения, развития и последствий.
Особенности математической модели во многом определяются типом моделируемой ЧС. Все ЧС можно разделить на природные, техногенные и социально-политические.
К природным ЧС относятся такие стихийные бедствия, как землетрясения, извержения вулканов, цунами, наводнения, ураганы, лавины, оползни, засухи, лесные пожары и др.
Техногенные (технологические) ЧС связаны с авариями на энергетических и промышленных объектах, а также транспортные катастрофы, которые сопровождаются взрывами, пожарами, химическим и радиоактивным заражением территорий.
К социально-политическим ЧС относятся войны, пограничные конфликты, терроризм, диверсии, саботаж.
К комбинированным природно-техногенным и природно- социальным ЧС относятся просадки грунтов, эпидемии, эпизоотии (инфекционные заболевания животных), эпифитотии (инфекционные болезни сельскохозяйственных культур) и др.
Все перечисленные выше ЧС могут быть исследованы методами математического моделирования. Замена реальной ЧС ее воображаемым виртуальным образом - математической моделью дает возможность безболезненно, сравнительно быстро и с минимальными затратами исследовать все мыслимые сценарии возникновения и развития ЧС, а также прогнозировать ее последствия.
Создание математической модели ЧС включает в себя несколько этапов. Начальным этапом является содержательное описание ЧС, которое составляется на основе всех имеющихся о ней знаний, результатов натурных обследований сходных ситуаций, консультаций с экспертами, изучения справочной и специальной литературы.
На втором этапе выполняется формализация содержательного описания модели, математическая постановка задачи с указанием всех необходимых исходных данных и искомых величин.
На третьем этапе формализованная схема ЧС должна быть преобразована в ее математическую модель. Для этого всю имеющуюся информацию необходимо выразить с помощью соотношений, неравенств, уравнений, алгоритмов. Уравнения, входящие в модель, дополняются начальными и граничными условиями, а также неравенствами, определяющими область допустимых значений вычисляемых величин.
На четвертом этапе, исследуется сама модель. Путем проведения многовариантных расчетов изучаются свойства модели и ее поведение при различных условиях.
На следующем этапе модель применяется к описанию реальных ЧС.
Путем сопоставления результатов вычислительных экспериментов с имеющимися опытными данными выполняется идентификация или уточнение параметров модели, ее тестирование, отлаадка и проверка адекватности.
После того, как адекватность модели, т.е. ее достаточное соответствие реальности, установлена, начинается использование модели для анализа и прогнозирования ЧС, происходящих в реальных условиях.
Схема построения математической модели приведена на рис. 2.
Необходимым условием получения достаточно точной и надежной математической модели ЧС является проверка ее адекватности.
Предположим, что математическая модель поражающего воздействия у некоторого фактора ЧС имеет вид:
где х - интенсивность этого фактора, и1, и2, …, ир - р параметров модели, которые предварительно задаются или определятся по опытным данным методом наименьших квадратов.
Рис. 2. Схема построения математической модели
Проверка адекватности математической модели осуществляется путем сравнения модельных (расчетных) значений yxi с эмпирическими значениями yi, найденными на m различных уровнях независимой переменной xi, где i = 1, 2, …, m. землетрясение математический модель
Простейшей предварительной мерой соответствия математической модели реальной ситуации является относительное среднее квадратическое отклонение расчетных данных от опытных:
где
-
среднее значение опытных данных,
р - число параметров модели.
Удовлетворительным можно считать значение ?~ ? 0,1? 0,2 .
Более полную оценку адекватности математической модели можно получить с помощью критерия Вилкоксона-Манна-Уитни. Для проверки гипотезы о статистической однородности выборок модельных yxi и эмпирических значений исследуемой характеристики ЧС из этих выборок составляется общий вариационный ряд и подсчитывается величина
где Ri - ранги, т.е. порядковые номера экспериментальных значений yi в общем вариационном ряде.
При объемах выборок m > 8 распределение случайной величины S близко к нормальному распределению с параметрами:
При уровне значимости б = 0,05 критическая область для гипотезы об однородности теоретической и экспериментальных выборок определяется неравенствами:
Если значение S не попадает в критическую область, то обе выборки однородны, т.е. принадлежат одной и той же генеральной совокупности и поэтому математическую модель следует признать адекватной.
Если условия позволяют на каждом из m уровней независимой переменной x выполнить по n параллельных опытов (наблюдений), результаты которых образуют матрицу наблюдений (yij), то появляется возможность выполнить дисперсионный анализ расчетных и эмпирических данных. Для этого на каждом уровне независимой переменной xi вычислим групповые средние значения yi:
по которым найдем общее среднее значение y :
3. Практическое задание
С помощью программы MS Excel путем аппроксимации данных таблицы построить математическую модель изменения площади лесного пожара. Аппроксимацию произвести линейной и квадратичной функциями. Выбрать функцию, наилучшим образом аппроксимирующую исходные данные. Спрогнозировать на основании полученной аналитической зависимости ожидаемое значение площади пожара или глубины затопления через 0,5 часа от крайнего значения указанного в таблице.
t, час. |
0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
0,6 |
0,7 |
|
S, км 2 |
0,344 |
0,656 |
0,898 |
1,045 |
1,07 |
1,162 |
1,172 |
1,264 |
Решение.
Введем таблицу значений в лист MS Excel и построим точечный график. Аппроксимируем полученный график линейной и квадратичной функциями.
Как видно из графика наилучшим образом данные аппроксимирует полином второй степени. Его R2 наибольший.
Используя его, находим прогнозное значение через 0,5 часа от крайнего значения указанного в таблице.
t, час. |
0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
0,6 |
0,7 |
1,2 |
|
S, км 2 |
0,344 |
0,656 |
0,898 |
1,045 |
1,07 |
1,162 |
1,172 |
1,264 |
0,483 |
Список литературы
1. Методология научных исследований: учебник для бакалавриата и магистратуры / Н.А. Горелов, Д.В. Круглов. - М.: Издательство Юрайт, 2015. - 290 с.
2. Методология научных исследований: учебник для магистратуры / М.С. Мокий, А.Л. Никифоров, В.С. Мокий; под ред. М.С. Мокия. - М.: Издательство Юрайт, 2015. - 255 с.
3. Основы моделирования чрезвычайных ситуаций: учебное пособие / В.Г. Шаптала, В.Ю. Радоуцкий, В.В. Шаптала; под общей ред. В.Г. Шапталы. - Белгород: Изд-во БГТУ, 2010. - 166 с.
4. Теория планирования эксперимента и анализ статистических данных: учеб. Пособие для магистров / Н.И. Сидняев. - 2-е изд., перераб. и доп. - М.: Издательство Юрайт, 2015. - 495 с.
5. Управление рисками, системный анализ и моделирование. В 2 т. Т. 1: учебник и практикум для бакалавриата и магистратуры / П.Г. Белов. - М. : Издательство Юрайт, 2015. - 460 с.
6. Управление рисками, системный анализ и моделирование. В 2 т. Т. 2: учебник и практикум для бакалавриата и магистратуры / П.Г. Белов. - М. : Издательство Юрайт, 2015. - 272 с..
7. Цой О.М. Математическое моделирование чрезвычайных ситуаций природного характера на юге Дальнего Востока / Под научн. ред. д.физ.-мат.н. И.В. Тросникова (Гидрометцентр России); МЧС России. М.: ФГБУ ВНИИ ГОЧС (ФЦ), 2012. - 192 с.
Размещено на Allbest.ru
...Подобные документы
Основные положения теории математического моделирования. Структура математической модели. Линейные и нелинейные деформационные процессы в твердых телах. Методика исследования математической модели сваи сложной конфигурации методом конечных элементов.
курсовая работа [997,2 K], добавлен 21.01.2014Создание математической модели движения шарика, подброшенного вертикально вверх, от начала падения до удара о землю. Компьютерная реализация математической модели в среде электронных таблиц. Определение влияния изменения скорости на дальность падения.
контрольная работа [1,7 M], добавлен 09.03.2016Построение математической модели технологического процесса напыления резисторов методами полного и дробного факторного эксперимента. Составление матрицы планирования. Рандомизация и проверка воспроизводимости. Оценка коэффициентов уравнения регрессии.
курсовая работа [694,5 K], добавлен 27.12.2021Принципы и этапы построения математической модели движения неуправляемого двухколесного велосипеда. Условия устойчивого движения. Вопрос гироскопической стабилизации движения. Модель движения велосипеда с гиростабилизатором в системе Matlab (simulink).
статья [924,5 K], добавлен 30.10.2015Математическое моделирование динамики биологических видов (популяций) Т. Мальтусом. Параметры и основное уравнение модели "хищник-жертва", ее практическое применение. Качественное исследование элементарной и обобщенной модификаций модели В. Вольтерра.
курсовая работа [158,1 K], добавлен 22.04.2011Операторы преобразования переменных, классы, способы построения и особенности структурных моделей систем управления. Линейные и нелинейные модели и характеристики систем управления, модели вход-выход, построение их временных и частотных характеристик.
учебное пособие [509,3 K], добавлен 23.12.2009Изучение методов решения уравнений математической физики, которые используются для расчётов распространения тепла, концентрации, волн. Решение уравнения теплопроводности интегро-интерполяционным методом (методом баланса), который применим во всех случаях.
курсовая работа [269,2 K], добавлен 15.11.2010Вводные понятия. Классификация моделей. Классификация объектов (систем) по их способности использовать информацию. Этапы создания модели. Понятие о жизненном цикле систем. Модели прогнозирования.
реферат [36,6 K], добавлен 13.12.2003История возникновения и развития математической логики как раздела математики, изучающего математические обозначения и формальные системы. Применение математической логики в технике и криптографии. Взаимосвязь программирования и математической логики.
контрольная работа [50,4 K], добавлен 10.10.2014Решение дифференциальных уравнений математической модели системы с гасителем и без гасителя. Статический расчет виброизоляции. Определение собственных частот системы, построение амплитудно-частотных характеристик и зависимости перемещений от времени.
контрольная работа [1,6 M], добавлен 22.12.2014Определения оптимизации схемы планирования эксперимента при работе со швейной машиной. Расчёт коэффициентов уравнения регрессии и выделение значимых коэффициентов прочности ткани и растяжения между лапкой и иглой. Проверка гипотезы адекватности модели.
курсовая работа [1,2 M], добавлен 30.12.2014Сущность математического моделирования. Аналитические и имитационные математические модели. Геометрический, кинематический и силовой анализы механизмов подъемно-навесных устройств. Расчет на устойчивость мобильного сельскохозяйственного агрегата.
курсовая работа [636,8 K], добавлен 18.12.2015Изучение физического процесса как объекта моделирования. Описание констант и параметров, переменных, используемых в физическом процессе. Схема алгоритма математической модели, обеспечивающая вычисление заданных зависимостей физического процесса.
курсовая работа [434,5 K], добавлен 21.05.2022Рассмотрение статических и динамических характеристик машины. Выбор математической модели систем электроприводов. Расчет параметров двигателя постоянного тока. Аппроксимация полученной переходной характеристики элементарными динамическими звеньями.
курсовая работа [833,3 K], добавлен 18.04.2014Моделирование входного заданного сигнала, построение графика, амплитудного и фазового спектра. Моделирование шума с законом распределения вероятностей Рэлея, оценка дисперсии отсчетов шума и проверка адекватности модели шума по критерию Пирсона.
курсовая работа [2,3 M], добавлен 25.11.2011Понятие и оценка необходимости в статистической обработке психологических данных. Методика и основные этапы математической обработки полученных данных, его критерии и параметры: признаки и переменные, шкалы измерения, анализ и оценка уровня значимости.
презентация [443,1 K], добавлен 28.02.2014Изучение актуальной задачи математического моделирования в биологии. Исследование модифицированной модели Лотки-Вольтерра типа конкуренция хищника за жертву. Проведение линеаризации исходной системы. Решение системы нелинейных дифференциальных уравнений.
контрольная работа [239,6 K], добавлен 20.04.2016Проектирование математической модели. Описание игры в крестики-нолики. Модель логической игры на основе булевой алгебры. Цифровые электронные устройства и разработка их математической модели. Игровой пульт, игровой контроллер, строка игрового поля.
курсовая работа [128,6 K], добавлен 28.06.2011Порядок преобразования исходных данных и построения математической модели оптимального плана доставки газет. Выбор метода решения и основные этапы его реализации. Принципы освоения и практического применения оптимизационного пакета прикладных программ.
курсовая работа [235,0 K], добавлен 25.03.2017Основные этапы обработки данных натуральных наблюдений методом математической статистики. Оценка полученных результатов, их использование при принятии управленческих решений в области охраны природы и природопользования. Проверка статистических гипотез.
практическая работа [132,1 K], добавлен 24.05.2013