Научная статья на тему 'Решение обратной задачи гидрофизического мониторинга на примере чилийских цунами'

Решение обратной задачи гидрофизического мониторинга на примере чилийских цунами Текст научной статьи по специальности «Математика»

CC BY
239
29
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИСТОЧНИК ЦУНАМИ / ГИДРОФИЗИЧЕСКИЙ МОНИТОРИНГ / ОБРАТНАЯ ЗАДАЧА / ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ / ОБРАБОТКА ДАННЫХ НАБЛЮДЕНИЙ / TSUNAMI SOURCE / HYDROPHYSICAL MONITORING / INVERSE PROBLEM / NUMERICAL SIMULATION / OBSERVATION DATA PROCESSING

Аннотация научной статьи по математике, автор научной работы — Быков Артём Александрович, Курако Михаил Александрович, Симонов Константин Васильевич

Исследование посвящено решению задачи восстановления параметров об источнике цунами на основе данных гидрофизического мониторинга с применением асимптотических формул С.Ю. Доброхотова и численного моделирования цунами в рамках вычислительной технологии MOST. Представлены математическая постановка задачи и вычислительная технология моделирования волн цунами в океане для оценки параметров изучаемых источников цунами. Проведен ряд числительных экспериментов на примере чилийских цунамигенных землетрясений и цунами в 2010-2015 гг

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Быков Артём Александрович, Курако Михаил Александрович, Симонов Константин Васильевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Solution of the hydrophysical monitoring inverse problem in the context of Chilean tsunamis

This research is devoted to the solution of the tsunami source parameters recovery problem, based on the hydrophysical monitoring data with the using of S.Y. Dobrokhotov’s asymptotic formulae and numerical simulation within the MOST computing technology. A mathematical formulation of the problem and computational technology of the ocean tsunami waves modelling for parameters estimation of the studied tsunami sources are presented. A series of numerical experiments for the Chilean earthquakes and tsunamis (2010 2015) is carried ou

Текст научной работы на тему «Решение обратной задачи гидрофизического мониторинга на примере чилийских цунами»

Заключение

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

Литература

1. Cristani M., Cuel R. A Survey on Ontology Creation Methodologies // In Int. J. on Semantic Web & Information Systems. 2005. № 1 (2). С. 48-68.

2. Jones D., Bench-Capon T. and Visser P. Methodologies for Ontology Development. URL: http://www.iet.com/Projects/RKF/SME/methodologies-for-ontology-development.pdf

3. Артемьева И.Л. Сложно структурированные предметные области. Построение многоуровневых онтологий // Информационные технологии. 2009. № 1. С.16-21.

4. Artemieva I.L. Multilevel modular chemistry ontology: structure and management // First Russia Pacific Conf. on Computer Technology and App. 6 - 9 September, 2010. ISBN: 978-09803267-3-4. С. 12-17.

5. Щукин Е.Д., Перцов А.В., Амелина Е.А. Коллоидная химия: Учеб. для университетов и химико-технолог. Вузов / - 3-е изд., перераб. и доп. - М.: Высш. шк., 2004. 445 с.

6. Артемьева И.Л., Клещев А. С. Необогащенные системы логических соотношений. Часть 1. - Научно-техническая информация. 2000. № 7. С. 18-28.

7. Артемьева И.Л., Рештаненко Н.В. Модульная модель онтологии органической химии // Информатика и системы управления. 2004. № 2 (08). С. 64-68.

8. Артемьева И.Л., Цветников В.А. Фрагмент онтологии физической химии и его модель // Исследовано в России. 2002. № 3. С. 454-474.

Ontology model of disperse systems and their classes

Irene Leonidovna Artemieva, DSc, Professor Natalya Valentinovna Ryabchenko, PhD, Researcher

There are two modules of the developed mathematical ontology model of colloid chemistry describing disperse system ontology model in the article. The developed model will be used in creation of intelligent modeling system that integrates ontologies, knowledge and data from different chemistry section.

Keywords: Mathematical modeling; ontology model; development of knowledge-based systems.

УДК 519.6

РЕШЕНИЕ ОБРАТНОЙ ЗАДАЧИ ГИДРОФИЗИЧЕСКОГО МОНИТОРИНГА НА ПРИМЕРЕ ЧИЛИЙСКИХ ЦУНАМИ

Артём Александрович Быков, аспирант, Тел.: 8 913 190 77 05, e-mail: [email protected] Институт космических и информационных технологий, Сибирский федеральный университет http://ikit.sfu-kras.ru

Михаил Александрович Курако, ст. преподаватель, Тел.: 8 923 285 67 12, e-mail: [email protected] Институт космических и информационных технологий, Сибирский федеральный университет http://ikit.sfu-kras.ru Константин Васильевич Симонов, профессор, д.т.н, Тел.: 8 913 595 49 02, e-mail: [email protected] Институт вычислительного моделирования СО РАН http://icm.krasn.ru

Исследование посвящено решению задачи восстановления параметров об источнике цунами на основе данных гидрофизического мониторинга с применением асимптотических формул С.Ю. Доброхотова и численного моделирования цунами в рамках вычислительной технологии MOST. Представлены математическая постановка задачи и вычислительная технология моделирования волн цунами в океане для оценки параметров изучаемых источников цунами. Проведен ряд числительных экспериментов на примере чилийских цунамигенных землетрясений и цунами в 2010-2015 гг.

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

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

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

В настоящее время системы предупреждения о цунами созданы и функционируют в Российской Федерации и многих странах Тихоокеанского региона: США, Японии, Канаде, Филиппинах, Чили, Французской Полинезии, Австралии и др. Принятие решений об объявлении отмене тревоги об угрозе цунами осуществляется на основе магнитудно-географического критерия и результатов обработки данных инструментальных и визуальных наблюдений за уровнем моря.

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

Предлагается для эффективного решения задачи об источнике цунами и проведения численного моделирования разработка вычислительной методики на основе асимптотических формул и соответствующего алгоритма С.Ю. Доброхотова [1; 2; 3; 4; 5],

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

Показано, что решение задачи об источнике цунами на основе асимптотических формул С.Ю. Доброхотова [1; 2; 3; 4; 5], позволяет выйти на новый уровень и решения задачи оценки опасности цунами. Важным элементом вычислительной методики также является обработка и анализ данных наблюдений о цунами с регистрирующих БЛЯТ-станций. Оценка параметров источника цунами на основе действующих БЛЯТ-станций предполагает оценку эффективности их расстановки и, соответственно, качества имеющихся данных в ближней зоне от изучаемого источника цунами.

Таким образом, на основе исследований в работах С.Ю. Доброхотова [1; 2; 3; 4; 5] предлагается вычислительный инструмент (методика) для повышения своевременности и надежности предупреждения об угрозах цунами, путем решения обратной задачи гидрофизического мониторинга и повышения точности оценок параметров источника цунами.

Следуя указанным исследованиям [1-5], при математическом моделировании распространения волн цунами в океане возникает следующая задача:

^ + Шу(С2и) = 0, — + Уг = 0, С(х) = ,

ы ы

( х ^

х = (х,Х2) Е Я\ Со 4=0 = 0 0)

Здесь г(х, t) - возвышение свободной поверхности жидкости, х)- глубина бассейна,g - ускорение свободного падения, г0(у)-заданная функция, убывающая на

бесконечности быстрее, чем^у ,5 > 1 и локализованная в области характеристического

размера / .

Предполагается, что значение /и мало по сравнению с характеристической длиной интервала, на котором меняется глубина дна и расстоянием до точки наблюдения. В работах [2; 3; 4; 5] получены формулы для асимптотических решений линеаризованной системы уравнений мелкой воды (1). В частности, построенные в [2; 3; 4; 5] асимптотические формулы являются основой для быстрого аналитико-численного алгоритма решения задачи (1). Показано, что с течением времени решение распадается на две части - вихревую и волновую - описывает распространение волны, локализованной в окрестности замкнутой кривой (фронта волны) на плоскости Я2.

Точка х из окрестности регулярной части волнового фронта г) определяется с помощью двух координат: ф( х, t) и у( х, t), где ф( х, t) определяется из условия, что вектор у = х - X(ф, t) ортогонален вектору скорости фронта ГЦ в точке X(ф, t), т.е.

(х - X(ф, t), Xф(ф, t)) = 0. (2)

В работе [3 ] показано, что при t > 0 в некоторой окрестности волнового фронта ГЦ , не зависящей от /, и вне некоторой окрестности фокальных точек, асимптотическая формула для возвышения свободной поверхности жидкости г представляется формулой:

_ _ (£ ^, х)

г( х, г) =

С

С0 • Яе

с(( (ф, t),))

е 2 • ^

х)Аф)

V / J

ф=ф( ,

+ 01 / |. (3)

Здесь r¡(t, х) - возвышение свободной поверхности жидкости, Хф(ф, t) - точка фронта, для которой считаем возвышение, С(х, t) - скорость распространения волны, C0 - скорость распространения волны в начальный момент времени, / - малый параметр, m - индекс Морса, п(ф) = (cosф,sinф), ф е [0,2^].

Функция S ((, х), описывающая расстояние от точки наблюдения до фронта волны, определяется следующим образом [3]:

S ((, х) = .

Ж"

]В(х(ф((, х), t))

• y.

(4)

В формуле (4) у - расстояние от точки х, в которой вычисляется возвышение свободной поверхности жидкости до фронта. Функция ¥ (г,ф), описывающая характер движения волны, имеет вид [3]:

F (г,ф) = -

a • e

in¡ 4

2л/2(jb; cos2 ф + b22 sin2 ф - iz

!/2 •

(5)

Особый интерес представляют области в окрестности точек самопересечения фронта, т.к. в них необходимо считать суперпозицию двух возвышений свободной поверхности жидкости, которые относятся к разным участкам фронта. Соответствующая такой области асимптотическая формула будет иметь следующий вид [3]:

Г( х, t) = ^

к (ф, о

Co

• Re

• F

f Sj ((, х)

и

, пф )|

фj=ф)(t,

+ OI/2 |. (6)

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

Выполнен ряд численных экспериментов по моделированию поведения волн цунами в ближней зоне от изучаемых источников на основе данных регистрации волн соответствующей системой БЛЕТ-станций. В численных экспериментах изучались данные о цунамигенных землетрясениях и цунами в районе побережья Чили в 2010 г., в 2014 г. и в 2015 г. Краткая сводка данных о цунамигенных землетрясениях и ближайших БЛЯТ-станциях, зарегистрировавших цунами, представлена в таблице.

Таблица

Данные об источниках землетрясений и DART-станциях

2

Дата события Координаты эпицентра Магнитуда землетрясения DART-станции

27.02.2010 35.846° ю. ш., 72.72° з. д. 8.8 32412 (17.980° ю.ш., 86.330° з.д)

1.04.2014 4.736° ю.ш., 152.558° в.д 8.2 32401 (20.473° ю.ш., 73.429° з.д), 32402 (26.743° ю.ш., 73.983° з.д)

16.09.2015 31.570° ю.ш., 71.654° з.д. 8.3 32401 (20.473° ю.ш., 73.429° з.д), 32402 (26.743° ю.ш., 73.983° з.д)

Данные о местоположении источников цунамигенных землетрясений приведены на рисунке 1. На рисунке 2 показаны мареограммы, записанные регистрирующими DART-станциями во время изменения возвышения свободной поверхности океана для чилийских цунами 2010 (DART 32412), 2014 и 2015 гг. (DART-станции 32401 и 32402), соответственно. На рисунке 3 представлены мареограммы, полученные в результате численных экспериментов, выполненных при помощи расчетов с использованием асимптотических формул для чилийских цунами 2010, 2014 и 2015 гг., соответственно. Расчеты проводились от источника цунамигенного землетрясения к ближайшим DART-станциям (32412, 32401 и 32402). На рисунке 4 приведены расчетные карты-схемы распределения максимальной амплитуды волны в ближней зоне для источников чилийских цунами 2010, 2014 и 2015 гг., соответственно, полученные в рамках численных экспериментов, выполненных при помощи программного комплекса MOST.

USGS Community I nternel Intensity Map OFFSHORE MAULE, CHILE Feb Z7 »10 «334:14 local ».90«S 72.7J27W MS 9 Qeph: JSton laurtOlOHan

USGS Community Internet Intensity Map OFFSHORE TARAPACA, CHILE Apr I 2014 07:+$:47PM local l*.i0475 70.7»!W MS.2 Depth: 25km ID:us

|29Я respcnsr ¡n99ZlPwtii5aneii»aii< 00 km

75'W 70'W 65 *W

ШМЯГТ 1 IHII IV V VI VII VII

шш W.A Ughl ModvraMl Sbang V.ry tb«4> Sav... Vloftnll ЕЛт

И-« «п» w» vt^ujhi I тк IWUHIHil и«„| У.ИМ#

USGS Community Internet Intensity Map OFFSHORE COQUIMBO, CHILE

Sep 16 2013 06:34:33 PM local 3I.5695S 7I.6S43W M8.3 Dcplh 29 km ID:us20003k7a

Рис. 1. Расположение цунамигенных землетрясении для чилийских цунами 2010, 2014 и 2015 гг.,

соответственно

Рис. 2. Мареограммы с записями чилийских цунами с DART-станций

Рис. 3. Расчетные мареограммы для чилийских цунами (2010, 2014, 2015 гг.)

Рис. 4. Расчетные карты распределения максимальной амплитуды для чилийских цунами (2010, 2014, 2015 гг.)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

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

Авторы считают, что проведенные исследования чилийских цунами (2010, 2014, 2015 гг.) в рамках разрабатываемой вычислительной методики на основе асимптотических формул С.Ю. Доброхотова и вычислительной технологии MOST показали новые возможности для быстрого решения обратной задачи гидрофизического мониторинга. При этом сравнительный анализ расчетных и натурных мареограмм позволяет скорректировать местоположение и форму первоначального варианта источника цунами. Таким образом, предлагаемый подход предполагает предварительную вычислительную процедуру быстрой коррекции параметров источника цунами с целью повышения точности и надежности оценки опасности цунами на основе существующих вычислительных комплексов численного моделирования цунами в реальных морских акваториях.

Литература

1. Dobrokhotov S.Yu. Explicit asymptotics for tsunami waves in framework of the piston model / S. Yu. Dobrokhotov, S.Ya. Sekerzh-Zenkovich, B. Tirozzi, B. Volkov //Russ. Journ. Earth Sciences. - Moscow, 2006. № 8. P. 1-12.

2. Dobrokhotov S.Yu. Asymptotic description of tsunami waves in a frame of the piston model: the general constructions a explicitly solvable models / S. Yu. Dobrokhotov, S.Ya. Sekerzh-Zenkovich, B. Tirozzi, B. Volkov // Fundamental and Applied Geophysics. - Saint-Petersburg, 2009. № 2. P. 15-29.

3. Dobrokhotov S.Yu. Asymptotic solutions of the linear shallow-water equations with localized initial data / S.Yu. Dobrokhotov, R. Nekrasov, B. Tirozzi // Journal of Engineering Mathematics. -Moscow, 2011. № 69. P. 225-242.

4. Доброхотов С.Ю. Асимптотические решения двумерного модельного волнового уравнения с вырождающейся скоростью и локализованными начальными данными / С.Ю. Доброхотов, В.Е. Назайкинский, Б. Тироцци. - Москва : МФТИ, Алгебра и анализ, 2010. 230 с.

5. Dobrokhotov, S.Yu. Localized wave and vortical solutions to linear hyperbolic system and their application to linear shallow water equations / S.Yu. Dobrokhotov, A.I. Shafarevich, B. Tirozzi // Russian Journal of Mathematical Physics. 2008. № 2. P. 192-221.

Solution of the hydrophysical monitoring inverse problem in the context of Chilean tsunamis

Artyom Alexandrovich Bykov, Postgraduate, SFU Mikhail Alexandrovich Kurako, Senior Teacher, SFU Konstantin Vasilyevich Simonov, Professor, ICMSB RAS

This research is devoted to the solution of the tsunami source parameters recovery problem, based on the hydrophysical monitoring data with the using of S.Y. Dobrokhotov's asymptotic formulae and numerical simulation within the MOST computing technology. A mathematical formulation of the problem and computational technology of the ocean tsunami waves modelling for parameters estimation of the studied tsunami sources are presented. A series of numerical experiments for the Chilean earthquakes and tsunamis (2010 - 2015) is carried out.

Keywords - tsunami source, hydrophysical monitoring, inverse problem, numerical simulation, observation data processing.

УДК: 004.04

ЗАБОЛЕВАЕМОСТЬ ТРУДЯЩИХСЯ ПО ОТРАСЛЯМ ЭКОНОМИКИ

(НА ПРИМЕРЕ г. НОВОКУЗНЕЦКА)

Анна Егоровна Власенко, к.т.н., инженер-математик Тел.: 8 961-716-52-60, e-mail: [email protected]

Кустовой медицинский информационно-аналитический центр, г. Новокузнецк

http://nmiac.ru

Наталья Михайловна Жилина, д.т.н., доцент, профессор кафедры медицинской кибернетики и информатики Тел.: 8 913 323 31 49, e-mail: [email protected] Новокузнецкий государственный институт усовершенствования врачей Минздрава РФ

http://ngiuv.nkz.ru

Проанализирована заболеваемость с временной утратой трудоспособности (ВУТ) на примере крупного промышленного центра Сибири. Выявлены ее особенности в зависимости от отраслей производственной и непроизводственной сферы экономики, рассмотрены наиболее значимые причины ВУТ.

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

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

и решении стратегических задач рефор- _

мирования здравоохранения [1; 2]. щ[ л^ЩЁк

Большое практическое значение ' ^ ^ Жилина имеет анализ заболеваемости с временной утратой трудоспособности (ВУТ), ибо позволяет судить об уровне, динамике временной нетрудоспособности, ее причинах и определять конкретные меры по ее снижению.

Цель: провести сравнительный анализ показателей заболеваемости с ВУТ (на примере г. Новокузнецка) по отдельным отраслям экономики с выявлением наиболее значимых причин временной утраты трудоспособности.

i Надоели баннеры? Вы всегда можете отключить рекламу.