CC BY
25УДК 004.415 В01:10.46960/1816-210Х 2020 3 53
Е.А. Рувинская, А.С. Епифанова, М.В. Кокоулина
ПЕРВАЯ ВЕРСИЯ ИНФОРМАЦИОННОЙ СИСТЕМЫ ОБРАБОТКИ И АНАЛИЗА
ДАННЫХ НАТУРНЫХ ИЗМЕРЕНИЙ ВНУТРЕННИХ ВОЛН В ОКЕАНЕ
Нижегородский государственный технический университет им. Р.Е. Алексеева
Создана система обработки натурных измерений внутренних волн, позволяющая проводить экспресс-анализ исходной записи посредством статистической обработки (вычисление моментов и построение гистограмм для ординат смещения изотермы/изопикны, высот, крутизны, периода и амплитуды угла волнового склона выделенных волн с возможностью выравнивания известными теоретическими законами распределения (нормальное, логнормальное, Вейбулла, Рэлея); для проверки качества аппроксимации используется критерий хи-квадрат) и построения спектра записи с возможностью его сравнения с моделью Гаррета-Манка (для удобства можно также выборочно отметить на графике частоты основных приливных гармоник (полусуточных (М2, 82, N2, К2) и суточных (Р1, К1, 01, Q1)). Пользователю доступны опции отображения спектров в линейном, логарифмическом и полулогарифмическом масштабе осей. Разработанный программный комплекс имеет графический интерфейс, созданный средствами вШ МаНаЪ, и является прототипом информационной онлайн-системы, которую планируется создать для анализа временных рядов изотермы/изопикны в контексте исследования динамики внутренних волн в океане. Тестирование функционала программного комплекса проведено на данных измерений температуры, полученных в шельфовой зоне Японского моря и предоставленных сотрудниками ТОИ ДВО РАН.
Ключевые слова: система обработки данных, внутренние волны, натурные измерения, спектр Гаррета-Манка, приливные гармоники, инерционная частота, гистограмма, выравнивание, логнормальное распределение, распределение Вейбулла, распределение Рэлея.
Введение
Исследование внутренних волн в океане - одна из наиболее популярных и актуальных задач современной геофизики, попытки решения которой предпринимаются в рамках различных моделей и концепций [1-7]. Интерес к изучению этого явления как важнейшего механизма переноса энергии в океане обусловлен, в первую очередь, практическими задачами. В последние десятилетия появилось множество новых исследований, подтверждающих тот факт, что внутренние волны не в меньшей (а часто - в большей) степени, чем поверхностные, могут влиять как на морскую инфраструктуру (воздействовать на подводные части морских гидротехнических сооружений, в том числе, посредством усиления сосредоточенных нагрузок, изгибающих и крутящих моментов и натяжения заякоривающих узлов, а также изменения рельефа дна), так и на функционирование морских экосистем (влиять на биопродуктивность за счет перемешивания водных масс и перераспределения питательных веществ, индуцировать перенос примесей и загрязнений [8]. Одним из наиболее важных прикладных аспектов в этом направлении является оценка воздействия интенсивных внутренних волн на нефтедобывающие платформы различных типов [9-16]. Поскольку внутренние волны представляют серьезную опасность для подобных объектов (известны случаи смещения нефтяных платформ на 200 м в горизонтальном и 10 м в вертикальном направлении, а также двукратного увеличения натяжения якорных цепей спар-платформ [17]), все больше работ по этой тематике встречается и в специализированных журналах, посвященных проблемам морской инженерии, а также в трудах специализированных конференций: например, [18] - о влиянии внутренних волн на буровые установки, [ 19, 20] - моделирование воздействия внутренних волн на спар-платформы, [21] - на шельфовые терминалы по производству сжиженного газа (СНГ) и процессы выгрузки, [22] - на полупогружные плавучие буровые установки
© Рувинская Е.А., Епифанова А.С., Кокоулина М.В.
(ППБУ) и др.) Таким образом, прикладные аспекты изучения динамики внутренних волн востребованы современной инженерной наукой.
В связи с прогрессом приборостроения, прикладной математики и методов дистанционного зондирования океана, с 1960-70-х гг. начинает появляться все больше и больше наблюдений интенсивных внутренних волн, В лаборатории моделирования природных и техногенных катастроф Нижегородского государственного технического университета им. Р.Е. Алексеева создан электронный атлас для работы с базой данных наблюдений внутренних волн, в которой систематизирована информация из более чем 2 800 литературных источников, описывающих их зарегистрированные проявления (https://lmnad.nntu.ru/igwatlas/) [ 23, 24]. Параллельно с накоплением океанологических данных о динамике внутренних волн встает вопрос также и о возможности экспресс-анализа записей натурных наблюдений (еще один проект, реализованный в нашей лаборатории, который может использоваться для численного моделирования динамики и экспресс-анализа параметров внутренних волн по данным океанологических атласов, описан в статьях [25, 26]). Информационно-аналитические системы, которые способны не только накапливать данные, но и предоставлять возможности для анализа этих данных, наиболее востребованы, и их появление является закономерным этапом на пути развития технологий в той или иной области знаний. Стоит отметить, что подобные системы активно разрабатываются, в том числе, в России, например, циклы работ ИВТ СО РАН и ИКИ РАН [27, 28] по обработке спутниковых изображений для решения задач мониторинга окружающей среды; работы ТОИ ДВО РАН (ресурсы http://oias.poi.dvo.ru/ и http://pacificinfo.ru [29, 30] - океанографические ИАС для Дальневосточного региона) и Института океанологии им. П.П. Ширшова РАН - ЕСИМО).
В настоящей статье обсуждается создание прототипа системы обработки и анализа натурных измерений внутренних волн, которая, в свою очередь станет первым этапом на пути создания информационной онлайн системы обработки и анализа натурных измерений внутренних волн и индуцируемых ими явлений на основе облачных вычислений с открытым доступом. Такая аналитическая система содержит программные средства для обработки данных - спектрального и статистического анализа записей натурных наблюдений внутренних волн (смещений заданного уровня плотности, либо скорости) с возможностью расширения функционала для анализа временных рядов в контексте исследования записей внутренних волн, полученных в ходе натурных наблюдений в океане.
Аналитическая система обработки натурных измерений внутренних волн
Интерфейс системы обработки и анализа натурных измерений внутренних волн разработан с помощью GUI Matlab и предоставляет пользователю возможности для статистического и спектрального анализа данных и графического отображения результатов.
Главное окно комплекса представлено на рис. 1. Оно содержит четыре кнопки и поле для отображения исходной записи. Кнопка «Загрузить данные» отвечает за выбор файла с входными данными, которые представляют собой временной ряд (колебания изопикны (линии одинаковой плотности), либо изотермы (линии одинаковой температуры) во времени).
После загрузки входных данных, при нажатии кнопки «Построить», в графическом окне отображается исходная запись и появляется возможность выбора способа ее дальнейшей обработки.
Исходная запись сглаживается путем спектральной фильтрации шумов с периодом менее 1 мин. В дальнейшем планируется расширить блок фильтрации и добавить возможность выбора фильтра, либо отказа от предобработки данных.
Для тестирования программного комплекса использовалась запись, предоставленная сотрудниками лаборатории статистической гидроакустики ТОИ ДВО РАН, полученная в ходе экспедиции 2011 г. в осенний период (с 11 октября 2011 (18:30) по 20 октября (10:45) (частота дискретизации данных - 1 с)) в шельфовой зоне Японского моря с термогирлянд [31-
34]. Из вертикального распределения температуры во времени в точке установки датчиков была выделена одна изотерма, соответствующая 4 °С, для которой проводилась апробация алгоритмов разработанной аналитической системы. На рис. 1 в графическом окне представлен временной ряд смещения колебания изотермы.
При нажатии кнопки «Спектральный анализ» вызывается новое окно, содержащее функционал для спектрального анализа (рис. 2). В графическом поле отображается спектр исходного сигнала £(/) (по умолчанию в линейном масштабе), полученный с помощью быстрого преобразования Фурье (дискретные спектральные амплитуды для дискретных значений временного ряда определяются формулой (1):
Sk = S (Л ) =
N
2га
=!>, ехрГ- ^ (] - IX* -1)
7=1 V N У
(1)
где N - количество отсчетов Sj, /к = - значения частот (положительных), ^ = 1/7^ -
частота дискретизации и - временной шаг между измерениями).
Программный комплекс позволяет изменять такие параметры, как масштаб осей (линейный, полулогарифмический и логарифмический), источник (тип) измерений (для аппроксимации спектром Гаррета-Манка) и приливные компоненты (для соотнесения пиков спектра с частотами суточных и полусуточных приливов).
Рис. 1. Главное окно
к
Еще в 1970-х гг. Гаррет и Манк, проанализировав результаты натурных наблюдений внутренних волн, сформулировали климатический спектр для «усредненных» океанских условий (модель GM) [35]. Он не зависит от волнообразующих факторов, показывая общие тенденции (не работает в районах акваторий с сильными течениями и топографическими особенностями). В экспериментальных работах при спектральном анализе волнового поля часто используется сравнение с полученным Гарретом и Манком законом, который претерпел различные модификации. Руководствуясь работой [36] (модель ОМ-76), мы используем уточненные спектральные модели вертикальных смещений для различных типов анализиру-
емых записей внутренних волн [3 7]: полученных с заякоренного буя («якорный» спектр) (2), либо с буксируемого судна («буксируемый» спектр) (3):
2 1 ( 12 — 1 2)1/2 (/, г) = - г1- --, (2)
я N(z) fJ
TSQ (k,z) =
f
rt- Jin
N ( z)
h N(z) N2(z) - fm
k~2, (3)
fm 2 N 2( z)
где все частоты f, fin и N) имеют размерность цикл/час, fin = 2 Qe sin ф - инерционная частота (Qe - частота вращения Земли, соответствующая периоду 24 ч, ф - географическая широта
точки наблюдений), N(z) =
g dp0(z) „
- частота Вяйсяля-Брента на горизонте z, g -
Ро(г) ^г
ускорение свободного падения, ро(г) - вертикальный профиль плотности морской воды, координата 2 направлена вниз, 2 = 0 соответствует поверхности моря. Плотность ро(^) была рассчитана по профилям температуры и солености на основании международного уравнения
состояния морской воды ТБОБ-Ю (http://www.teos-10.org/). Параметры г = ЕЬ2NN [м2цикл/ч]
и £ = { „ [ч/цикл/м] могут варьироваться в зависимости от океанских условий. Величины Е, -ЬИ
NN , Ь - безразмерный параметр, характеризующий энергию внутренних волн в «элементарном» столбе жидкости, масштаб частоты плавучести и вертикальный масштаб волновых движений (глубина стратифицированной жидкости) соответственно. Во многих задачах используются канонические значения этих параметров: Е = 2л10—5; Ь = 1300 [м]; NN = 3 [цикл/ч]. Однако, существуют процедуры для определения этих величин непосредственно из экспериментальных данных [38]. Параметр у - количество бароклинных мод (в этой модели обычно берется у = 6). Для выполнения аппроксимации в программном комплексе необходимо задать значение частоты Вяйсяля-Брента. Можно также изменить параметры «по умолчанию»: Ь, у, N, Е. Для спектра колебаний изотермы, полученного при измерениях с заякоренного буя, известно, что за исключением инерционной и приливных частот, он будет иметь непрерывный характер во всём диапазоне частот; спектральная плотность спадает по степенному закону с показателем степени от -1 до -4, а чаще от -2 до -3.
Одной из функций аналитической системы также является выбор приливных компонент: с помощью чекбоксов пользователь имеет возможность выбрать суточные и полусуточные приливные гармоники, которые отобразятся на графике полученного спектра, позволяя проверить наличие пиков на этих частотах (доступны М2 (лунный полусуточный, главная волна), 82 (солнечный полусуточный, главная волна), N2 (лунный полусуточный, большая эллиптическая волна), К2 (лунно-солнечный полусуточный, деклинационная волна); Р1 (солнечный суточный, главная волна), К1 (лунно-солнечный суточный, деклинационная волна), 01 (лунный суточный, главная волна), Q1 (лунный суточный, эллиптическая волна)). На рис. 2 приведен спектр температурных колебании и отмечены две приливные компоненты М2 и К1.
В главном окне (рис. 1) пользователь может выбрать также опцию «Статистический анализ», после чего будет вызвано новое окно (рис. 3), в графическом поле которого представлена гистограмма распределения ординат (по умолчанию). Гистограмма повторяемости ординат смещения изотермы позволяет получить общую картину диапазона изменчивости значений этой величины.
В разработанной системе обработки данных натурных наблюдений внутренних волн кроме ординат исходного временного ряда можно также проанализировать такие параметры, как высоты волн, период, крутизна, амплитуда волнового склона. По выбранной характеристике будет построена гистограмма распределения. В классическом смысле крутизна волны
- это отношение высоты волны к её длине, а амплитуда угла волнового склона (локальная крутизна) есть максимальное значение угла в точках профиля волны. Однако, поскольку мы анализируем запись по времени (знаем период волны, а не длину), то для определения этих характеристик необходимо ввести безразмерные переменные, определяющие оси абсцисс и ординат: t* = t/T*, z* = z/H*, где t и z - значения времени записи и ординат изотермы соответственно, а T* и H* - модальные значения для логнормальных кривых, выравнивающих гистограммы высот и периодов волн. Высота волны оценивается как средняя величина между вершиной и подошвой волны, период каждой волны определяется как временной интервал между соседними минимумами залегания изолинии. Всего из сглаженной записи (сглаживание проводилось путем спектральной фильтрации шумов с периодом менее 1 мин) смещений изотермы 4 °C, рассматриваемой как запись нерегулярного волнения, было выделено 1046 волн, волны менее 0,5 м не учитывались. Для примера на рис. 3 приведена гистограмма повторяемости высот волн в записи смещения изотермы 4°C. Отметим, что высоты всех рассматриваемых волн лежат в диапазоне 0,5 - 10 м. Средняя высота волны составляет 1.6 м, среднеквадратическое отклонение высоты 1,12 м (под графиком выводится среднее значение и среднеквадратическое отклонение группированного статистического ряда).
30
2 0
И
£ ю
51 -12 II
/ --1--* " -
0.005
0.01
Масштаб осей
®Линейный ЕПолулогарифмический ОЛогарифмический
0.015 0.02 0.025 0.03
f, цикл/ч
Аппроксимация спектром Гарретта-Манка
Измерения получены с ,.
0.035
0.04
0.045
0.05
Построить
ОЗаякоренного буя цБуксирующего судна ©Без аппроксимации
Построить
Выберите приливные компоненты
0М2 (Т = 12.42 ч) Ьк1 (Т = 23.93 ч)
□ S2 (Т = 12 ч) D01 (Т = 25.32 ч)
□ N2 (Т = 12.66 ч) DPI (т = 24.07 ч)
□ К2 (Т = 11.97 ч) DQ1 (Т = 26.87 ч)
Рис. 2. Окно спектрального анализа
Кроме этого, имеется возможность аппроксимации экспериментальных данных известными теоретическими распределениями плотности вероятности. В нашем программном комплексе для выравнивания гистограммы можно выбрать нормальное распределение, лог-нормальный закон, распределения Вейбулла и Рэлея [39]. Роль нормального распределения в теории вероятностей и математической статистике сложно переоценить: оно занимает уникальное положение в практических расчетах, при этом риск серьезных ошибок при его использовании невелик даже если натурные данные лучше описываются другими видами распределений. Это справедливо и при исследовании волновых процессов. Так, например, волновые ординаты (в том числе в контексте внутренних волн) часто неплохо описываются нормальным законом распределения [40]. Мы используем следующую форму нормального закона (4):
fi( x) =
1
(x-^)2 2o2
(4)
где в качестве оценки д выступает среднее выборочное, а в качестве оценки а - квадратный корень из выборочной дисперсии (более подробно об истории применения нормального закона и об оценивании его параметров можно прочитать, например, в [39]).
Еще одним часто используемым непрерывным и асимметричным «колоколообраз-ным» законом распределения является логнормальный закон (5):
(1п X / О2
/-(X) = е 282 , (5)
хол/ 2 л
где параметры £ и 5 (безразмерный) оцениваются из исходных статистических данных (см. более подробно [39]). Это одно из наиболее часто применяемых распределений (наравне с распределением Вейбулла и Рэлея) при обработке волновых записей и анализе распределений таких параметров, как высоты волн [41-43]. Так, в качестве примера использования лог-нормального закона при анализе параметров внутренних волн можно привести следующие результаты: в нашей статье [44] было показано, что высоты и периоды внутренних волн по данным из Японского моря удовлетворяют логнормальному закону; в статье [45], где анализировались записи из Белого и Баренцева морей, также для аппроксимации гистограммы высот волн используется логнормальная кривая.
Двухпараметрическое семейство кривых Вейбулла (которое по своей сути является степенным преобразованием экспоненциального закона) имеет вид (6):
c
c-1
/з(x) = - - e ^ , (6)
aya J
где с - параметр формы, который также называют наклоном Вейбулла; безразмерная величина (с =1 - экспоненциальный случай), а- параметр масштаба, который имеет ту же единицу измерения, что и величина x (при увеличении только этого параметра кривая растягивается вправо вдоль оси абсции, а ее высота уменьшается; при уменьшении а кривая «сжимается» влево, а ее высота увеличивается).
Частным случаем распределения Вейбулла с параметрами с = 2 и а= •Л b является распределение Рэлея:
/(x) = Xe2b2, (7)
b
где b - параметр распределения Рэлея, является модой этого распределения (имеет размерность x). Выбор перечисленных законов распределения обусловлен частотой их применения к анализу характеристик волновых процессов. Аппроксимация гистограммы производится средствами библиотеки Matlab dfittool.
Также гистограмма может быть построена без выравнивания теоретическим законом распределения. При аппроксимации в окне слева отображаются значения выборочного среднего и среднеквадратического отклонения для сгруппированного ряда, а также значения параметров кривой распределения (справа). Проверяется критерий Пирсона, и выводятся значения хи-квадрат и критическое значение, по которым можно судить насколько хорошо теоретическая кривая описывает исходные данные [39].
Для модельного примера гистограмма высот хорошо аппроксимируется теоретическим логнормальным законом распределения с плотностью вероятности (5) с параметрами д = 0.27 м, о = 0.63 м для высот волн H.
В дальнейшем планируется также добавить в комплекс блок с функционалом для прогнозирования ожидаемых высот волн на основании анализируемых записей [44, 46]. Хотя прогноз, выполненный по «коротким» записям (несколько десятков часов наблюдений) должен уточняться, в акваториях со сложными климатическими условиями он оправдан [45, 47], поскольку сложно обеспечить получение действительно длительных записей.
\c X \
Таким образом, разработанный программный комплекс позволяет проводить экспресс-анализ экспериментальных данных, представляющих собой запись изотермы (изопик-ны), в контексте исследования внутренних волн.
Рис. 3. Окно статистического анализа. Логнормальный закон распределения
Рис. 4. Окно статистического анализа. Закон распределения Вейбулла
Заключение
В настоящей работе с использованием библиотек МайаЬ разработана первая версия системы обработки и анализа натурных наблюдений внутренних волн в океане с графическим интерфейсом. Программный комплекс предоставляет возможности статистического анализа временных рядов, а именно: построение гистограммы и ее выравнивание с помощью известных теоретических законов распределения, вычисление моментов ряда, а также проверка критерия хи-квадрат. Пользователь может провести статистический анализ как орди-
нат исходной записи изотермы (изопикны), так и высот выделенных из записи волн, их периодов, крутизны и амплитуды волнового склона. Функционал вкладки спектрального анализа позволяет строить спектр записи в различных вариантах отображения, а также для сравнения наносить на график частоты основных приливных компонент. Доступна также опция сравнения спектра с моделью Гаррета-Манка заякоренного и буксируемого типа. Разработанный программный комплекс является прототипом информационной онлайн системы по обработке натурных данных внутренних волн в океане. Планируется расширить функционал системы за счет добавления различных опций фильтрации временных рядов, вкладки, предоставляющей возможности вейвлет-анализа, а также оценок воздействия на дно и сооружения в контексте динамики внутренних волн.
Представленные результаты получены при финансовой поддержке гранта Президента Российской Федерации по государственной поддержке научных исследований ведущих научных школ Российской Федерации НШ-2485.2020.5, а также стипендии Президента Российской Федерации молодым ученым и аспирантам (СП-1225.2019.5).
Библиографический список
1. Kurkina, O.E. Huge internal waves in the vicinity of the Spitsbergen island (Barents Sea) / O.E. Kurkina, T.G. Talipova // Natural Hazards and Earth System Science. - 2011. - V. 11. - № 3. -P. 981-986.
2. Zahibo, N. Strongly nonlinear steepening of long interfacial waves / N. Zahibo, A. Slunyaev, T. Talipova, E. Pelinovsky, O. Polukhina, A. Kurkin // Nonlinear Processes in Geophysics. - 2007. - V. 14. - № 3. - P. 247-256.
3. Куркина, О.Е. Динамика солитонов неинтегрируемой версии модифицированного уравнения Кортевега-де Вриза / О.Е. Куркина, А.А. Куркин, Е.А. Рувинская, Е.Н. Пелиновский, Т. Соомере // Письма в Журнал экспериментальной и теоретической физики. - 2012. - Т. 95. - № 1-2. - С. 98103.
4. Талипова, Т. Г. Моделирование динамики интенсивных внутренних волн на шельфе / Т. Г. Тали-пова, Е.Н Пелиновский, А.А. Куркин, О.Е. Куркина // Известия РАН. Физика атмосферы и океана. - 2014. - Т. 50. - № 6. - С. 714-722.
5. Rouvinskaya, E. Transformation of internal breathers in the idealised shelf sea conditions / T. Talipova, O. Kurkina, D. Tyugin, T. Soomere // Continental Shelf Research. - 2015. - V. 110. -P. 60-71.
6. Kurkina, O.E. Propagation regimes of interfacial solitary waves in a three-layer fluid / O.E. Kurkina, A.A. Kurkin, E.A. Rouvinskaya, T. Soomere // Nonlinear Processes in Geophysics. -2015. - V. 22. - № 2. - P. 117-132.
7. Kurkina, O. Propagation regimes and populations of internal waves in the Mediterranean Sea basin / O. Kurkina, E. Rouvinskaya, T. Talipova, T. Soomere // Estuarine, Coastal and Shelf Science. - 2017. - V. 185. - P. 44-54.
8. Woodson, C.B. The fate and impact of internal waves in nearshore ecosystems //Annual review of marine science. - 2018. - V. 10. - P. 421-441.
9. Chakrabarti, S. Handbook of Offshore Engineering (2-volume set) / S. Chakrabarti. - Elsevier. - 2005.
10.Cai, S. A method to estimate the forces exerted by internal solitons on cylindrical piles / S. Cai, X. Long, Z. Gan // Ocean Engineering. - 2003. - V. 30. - № 5. - P. 673-689.
11. Cai, S. Forces and torques exerted by internal solitons in shear flows on cylindrical piles / S. Cai, X. Long, S. Wang // Applied Ocean Research. - 2008. - V. 30. - № 1. - P. 72-77.
12.Osborne, A.R. Nonlinear ocean wave and the inverse scattering transform // Scattering. - Academic Press, 2002. - P. 637-666.
13. Song, Z.J. Comparisons of internal solitary wave and surface wave actions on marine structures and their responses / Z.J. Song, B. Teng, Y. Gou, L. Lu, Z.M. Shi, Y. Xiao, Y. Qu //Applied Ocean Research. -2011. - V. 33. - № 2. - P. 120-129.
14. Si, Z. A numerical simulation of shear forces and torques exerted by large-amplitude internal solitary waves on a rigid pile in South China Sea / Z. Si, Y. Zhang, Z. Fan // Applied Ocean Research. - 2012. -V. 37. - P. 127-132.
15.Xie, J. Strongly nonlinear internal soliton load on a small vertical circular cylinder in two-layer fluids / J. Xie, Y. Jian, L.Yang //Applied mathematical modelling. - 2010. - V. 34. - № 8. - P. 2089-2101.
16.Xie, J. A numerical study of the load on cylindrical piles exerted by internal solitary waves / J. Xie, J. Xu, S. Cai // Journal of Fluids and Structures. - 2011. - V. 27. - № 8. - P. 1252-1261.
17.Osborne, A.R. The influence of internal waves on deep-water drilling / A. R. Osborne, T.L. Burch, R.I. Scarlet //Journal of Petroleum Technology. - 1978. - V. 30. - №10. -P. 1497-1504.
18.Kurup, N.V. Study of nonlinear internal waves and impact on offshore drilling units / N.V. Kurup, S. Shi, Z. Shi, W. Miao, L. Jiang // International Conference on Offshore Mechanics and Arctic Engineering. - 2011. - V. 44335. - P. 831-840.
19.Sun, M. A new concept of spar in deep water and its hydrodynamic performance under internal wave / M. Sun, W. Huang // The Twenty-second International Offshore and Polar Engineering Conference. - International Society of Offshore and Polar Engineers. - 2012. - P.975-982.
20.Kurup, N.V. Numerical modeling of internal waves within a coupled analysis framework and their influence on spar platforms / N.V. Kurup, S. Shi, L. Jiang, M.N. Kim // Ocean Systems Engineering. - 2015. -V. 5. - № 4. - P. 261-277.
21.Xu, W. Internal wave soliton passage simulation during offloading / W. Xu, Y. Li, A. Voogt // ASME 2013 32nd International Conference on Ocean, Offshore and Arctic Engineering. - American Society of Mechanical Engineers Digital Collection. - 2013. - V. 55393. - P. V005T06A026.
22.Wei, J. Maximum Offset for the Emergency Disengagement and Control of the Deep water Semi-Submersible Platform in Internal Wave // Ocean Engineering Equipment & Technology. - 2017. - V. 4. -P. 29-36.
23.Епифанова, А.С. База данных наблюдений внутренних волн в Мировом океане / А.С. Епифанова, А.В. Рыбин, Т.Е. Моисеенко, О.Е. Куркина, А.А. Куркин, Д.Ю. Тюгин // Морской гидрофизический журнал. - 2019. - Т. 35. - № 4. - С. 395-403. DOI: 10.22449/0233-7584-2019-4-395-403.
24.Епифанова, А.С. О разработке цифрового атласа наблюдений внутренних волн в мировом океане / А.С. Епифанова, А.А. Куркин, О.Е. Куркина, Т.Е. Моисеенко, А.В. Рыбин // Труды НГТУ им. Р.Е. Алексеева. - 2019. - № 4. - C.17-26.
25. Тюгин, Д.Ю. Программный комплекс для численного моделирования внутренних гравитационных волн в мировом океане / Д.Ю. Тюгин, О.Е. Куркина, А.А. Куркин // Фундаментальная и прикладная гидрофизика. - 2011. - Т. 4, - № 2. - С. 32-44.
26. Тюгин, Д.Ю. Повышение производительности программного комплекса для моделирования внутренних гравитационных волн IGW Research с помощью Intel® Parallel Studio XE 2013 / Д.Ю. Тюгин, А.А. Куркин, Е.Н. Пелиновский, О.Е. Куркина // Фундаментальная и прикладная гидрофизика. - 2012. - Т. 5. - № 3. - C. 89-95.
27.Шокин, Ю.И. Информационная система приема, обработки и доступа к спутниковым данным и ее применение для решения задач мониторинга окружающей среды / Ю.И. Шокин, Н.Н. Добрецов, Е.А. Мамаш, В.А. Кихтенко, П.В. Воронина, В.В. Смирнов, Д.Л. Чубаров // Вычислительные технологии. - 2015. - Т. 20. - № 5. - С. 157-174.
28.Лупян, Е.А. Центр коллективного пользования системами архивации, обработки и анализа спутниковых данных ИКИ РАН для решения задач изучения и мониторинга окружающей среды / Е.А. Лупян, М.А. Бурцев, И.В. Балашов, С.А. Барталев, В.Ю. Ефремов, А.В. Кашницкий, А.А. Мазуров, А.М. Матвеев, О.А. Суднева, И.Г. Сычугов, В.А. Толпин, И.А. Уваров // Современные проблемы дистанционного зондирования Земли из космоса. - 2015. - Т. 12. - № 5. - С. 263-284.
29.Ростов, И.Д. Базы данных и информационные системы по океанографии ТОИ ДВО РАН для поддержки научных исследований и морской деятельности в ДВ регионе / И.Д. Ростов, А.А. Пан, В.И. Ростов и др. // Вест. ДВО РАН. - 2007. - № 4. - С. 85-95.
30.Ростов, И.Д. Информационная поддержка научных исследований и комплексного управления прибрежной зоной России в Японском море / И.Д. Ростов, Н.И. Рудых, В.И. Ростов и др. // Вест. ДВО РАН. - 2011. - № 6. - С. 89-93.
31.Ярощук, И.О. Экспериментальные исследования внутренних волн в прибрежной зоне Японского моря / И.О. Ярощук, А.П. Леонтьев, А.В. Кошелева, А.Н. Самченко, А.А. Пивоваров, Ф.Ф. Храп-ченков, А.Н. Швырев, Е.И. Ярощук // Подводные исследования и робототехника. - 2013. -№ 1(15). - С. 37-44.
32.Ярощук, И.О. Экспериментальные исследования внутренних волн в прибрежной зоне Японского моря / И.О. Ярощук, А.П. Леонтьев, А.В. Кошелева, А.Н. Самченко, А.А. Пивоваров, Ф.Ф. Храп-
ченков, А.Н. Швырев, Е.И. Ярощук // Подводные исследования и робототехника. - 2013. - № 1. -С. 37-44.
33. Yaroshchuk, I.O. On intense internal waves in the coastal zone of the Peter the Great Bay (the Sea of Japan) / I.O. Yaroshchuk, A.P. Leont'ev, A.V. Kosheleva, A.A. Pivovarov, A.N. Samchenko, D.V. Ste-panov, A.N. Shvyryov // Russian Meteorology and Hydrology. - 2016. - V. 41. - № 9. - P. 629-634. DOI: 10.3103/s1068373916090053.
34. Samchenko, A.N. Internal gravity waves in the coastal zone of the Sea of Japan according to the natural observations / A.N. Samchenko, I.O. Yaroshchuk, A.V. Kosheleva // Regional Studies in Marine Science. - 2018. - V. 18. - P. 156-160.
35.Garrett, C. Space-time scales of internal waves: A progress report / C. Garrett, W. Munk // Journal of Geophysical Research. - 1975. - V. 80. - № 3. - P. 291-297.
36.Desaubies, Y.J.F. Analytical representation of internal wave spectra / Y.J.F. Desaubies // Journal of Physical Oceanography. - 1976. - V. 6. - № 6. - P. 976-981.
37.Morozov, E.G. Oceanic Internal Tides. Observations, Analysis, and Modeling / E.G. Morozov // Springer International Publishing. 2018. - V. 10. - P. 978-3. DOI: 10.1007/978-3-319-73159-9.
38.Levine, M.D. A modification of the Garrett-Munk internal wave spectrum / M.D. Levine // Journal of physical oceanography. - 2002. - V. 32. - № 11. - P. 3166-3181.
39. Джонсон, Н.Л. Одномерные непрерывные распределения. Часть 1. / Н.Л. Джонсон, С. Коц, Н. Ба-лакришнан // Москва: Бином. Лаборатория знаний. - 2010. - C. 703.
40.Миропольский, Ю.З. Динамика внутренних гравитационных волн в океане / Ю.З. Мирополь-ский.. - Л.: Гидрометеоиздат. - 1981. - С. 302.
41.Guide to Wave Analysis and Forecasting // World Meteorological Organization (WMO): Geneva, Switzerland. - 1998. - V. 1998. - P. 101-119.
42.Tayfun, M. A. Distributions of wave steepness and surf parameter / M.A. Tayfun // Journal of waterway, port, coastal, and ocean engineering. - 2006. - V. 132. - № 1. - P. 1-9.
43.Kim, C. H. Nonlinear waves and offshore structures / C.H. Kim // World scientific. - 2008. - V. 27. - P. 12-53.
44.Kokoulina, M.V. Statistics of field data on internal waves» / M.V. Kokoulina, O.E. Kurkina, E.A. Ru-vinskaya, A.A. Kurkin, A.R. Giniyatullin // The Fourteenth MEDCOAST Congress on Coastal and Marine Sciences, Engineering, Managemеnt and Conservation. - 2019. -P. 733-743.
45.Свергун, Е.И. Оценка повторяемости интенсивных внутренних волн в Белом и Баренцевом морях по данным экспедиционных исследований / Е.И. Свергун, А.В. Зимин // Фундаментальная и прикладная гидрофизика. - 2017. - Т. 10. - № 2. - С. 13-19.
46.Кокоулина, М.В. Вероятностные характеристики интенсивных короткопериодных внутренних волн в Японском море / М.В. Кокоулина, О.Е. Куркина, Е.А. Рувинская, А.А. Куркин // Морской гидрофизический журнал. - 2020.
47.Зимин, А.В. Короткопериодные внутренние волны в шельфовых районах Белого, Баренцева и Охотского морей: оценка повторяемости экстремальных высот и динамических эффектов в придонном слое / А.В. Зимин, Е.И. Свергун // Фундаментальная и прикладная гидрофизика. -2018. -Т. 11. - № 2. - С. 66-72.
Дата поступления в редакцию :26.07.2020
E.A. Rouvinskaya, A.S. Epifanova, M.V. Kokoulina
THE FIRST VERSION OF THE INFORMATION SYSTEM FOR PROCESSING AND ANALYSIS OF DATA OF INTERNAL WAVES' MEASUREMENTS IN THE OCEAN
Nizhny Novgorod state technical university n.a. R.E. Alekseev
Purpose: The aim of the work is to create a system for processing field measurements of internal waves by means of statistical analysis (calculation of statistical moments and construction of histograms for ordinates of isotherm / isopyc-nal displacement, heights, steepness, period and amplitude of the wave slope angle of the identified waves with the possibility of approximation with the known theoretical distribution laws (normal, lognormal, Weibull, Rayleigh); the chi-square test is used to check the quality of the approximation) and the construction of the spectrum with the possibility of comparing it with the Garrett-Munk model (for convenience, one can also selectively mark the frequencies of the fundamental tidal harmonics ( semidiurnal (M2, S2, N2, K2) and diurnal (P1, K1, O1, Q1)).The user has access to options for displaying spectra in a linear, logarithmic and semi-logarithmic scale of the axes.
Design/methodology/approach: The developed software package has a graphical interface, created by GUI Matlab. The functional testing of the software package was carried out on the temperature measurement data obtained in the shelf zone of the Sea of Japan and provided by the staff of the POI FEB RAS.
Findings: It is a prototype of the online information system, which we are planning to create on its basis for the analysis of isotherm / isopycnal time series in the context of studying the dynamics of internal waves in the ocean. Originality/value: Such system can be successfully used for the express analysis of experimental data of internal waves' measurements.
Key words: data processing system, internal waves, field measurements, Garrett-Munk spectrum, tidal harmonics, inertial frequency, histogram, lognormal distribution, Weibull distribution, Rayleigh distribution.
