CC BY
11УДК 519.688; 534.6.08; 534.8
Е.А. Букварёв, М.В. Мартынюк
ИСПОЛЬЗОВАНИЕ ОСОБЕННОСТЕЙ ЛИНЕЙЧАТОГО СПЕКТРА ПЕРИОДИЧЕСКОЙ ПОСЛЕДОВАТЕЛЬНОСТИ ИМПУЛЬСОВ ПРИ АВТОМАТИЗИРОВАННОЙ ОБРАБОТКЕ ЭХОСИГНАЛОВ
Нижегородский государственный технический университет им. Р.Е. Алексеева
В статье рассматривается решение проблемы надёжной качественной оценки временного положения отдельных реплик в эхосигнале полученном из акустического тракта в присутствии ревербераций.
Ключевые слова: ультразвуковое зондирование, реверберирующая среда, акустическая анизотропия.
Введение
В настоящее время достаточно широко применяются методы неразрушающего контроля с использованием акустического зондирования исследуемой среды (акустические уровнемеры, толщиномеры, измерители механических напряжений). Обычно информативным параметром, извлекаемым из принятого сигнала (эхосигнала), является время распространения. По времени распространения, при известном распределении скоростей в среде, можно определить толщину материала или уровень жидкости. Методы акустической тензометрии позволяют на основании изменения времени распространения для ненапряжённого и напряжённого состояний материала определить величину внутреннего механического напряжения.
Практически всегда исследуемый акустический тракт является реверберирующим. При наличии в исследуемой среде только одного акустического тракта (такой средой является, например, стенка трубы, выполненная из слабоанизотропной стали) интервал между импульсами (репликами) определяется временем прохождения сигналом пути: излучатель -граница раздела сред - приёмник, и последовательность импульсов в эхосигнале для подобного тракта является периодической. Наличие множества реплик (откликов, отражений) повышает информативность принятого из среды сигнала, но затрудняет его автоматическую обработку.
Проблема повышения надёжности подобных автоматизированных систем является достаточно актуальной и особенно важна в случаях, когда время работы оператора должно быть максимально сокращено, например, при контроле защитных оболочек реакторных отделений АЭС, когда оператор вынужден находиться зоне повышенной радиационной опасности.
Рассматриваемый метод позволяет повысить надёжность определения интервала повторения реплик в реверберирующем тракте и уменьшить количество необходимой априорной информации.
Постановка задачи
Пусть производится зондирование по схеме (рис. 1), содержащей цепи согласованной обработки эхосигнала. Сигнал на выходе согласованного фильтра (СФ-эхосигнал) представляет собой периодическую последовательность сжатых импульсов с высокочастотным заполнением в случае отсутствия акустической анизотропии (рис. 3, б) либо сумму периодических последовательностей, отличающихся энергией и периодами повторения. Периоды повторения отличаются не более, чем на 15%, при соответствующей акустической анизотропии. Требуется определить период следования импульсов в последовательности с макси-
© Букварёв Е.А., Мартынюк М.В., 2011.
мальной энергией и положение первого импульса в этой последовательности. Метод определения должен обеспечивать автоматизированную обработку при отношении сигнал/помеха не хуже 20 дБ, количестве реплик не менее пяти.
Рис. 1. Этапы предварительной обработки эхосигнала
Определение интервала
Определение интервала предлагается построить по схеме, показанной на рис. 2. На первой стадии производится амплитудное детектирование с помощью устройства взятия модуля, последовательно соединённого с фильтром низкой частоты ФНЧ-1. Частотная характеристика последнего выбирается с учётом подавления частоты заполнения входных импульсов не менее 40 дБ. ФНЧ-2 обеспечивает дополнительную фильтрацию огибающей с целью её последующего выравнивания. Выравнивание огибающей обеспечивает узел деления А/В (блок 7), на вход которого поступает огибающая СФ-эхосигнала после нелинейного преобразования в блоке 6, выполняющего подавление сигнала меньшего заданного порога.
Затем вычисляется модуль спектра сигнала, и вся дальнейшая обработка производится в спектральной области. Спектральная оценка относится к категории интегральных и обладает преимуществом при оценке периодических сигналов в условиях флуктуаций и шумов. Во многом это вызвано уменьшенным влиянием на результат поведения сигнала во временной области в окрестностях точки. Иными словами, локальные особенности сигнала (несколько локальных максимумов на одном пике, модуляции амплитуды реплик), затрудняющие анализ эхосигнала во временной области, слабо влияют на спектр эхосигнала.
Известно, что спектр периодической последовательности импульсов является линейчатым. Положение (частота) первой гармоники в линейчатом спектре импульсной последовательности определяется периодом (интервалом) повторения импульсов, поэтому можно оценивать интервал повторения реплик на основании положения первой гармоники в спектре ВО СФ-эхосигнала, однако прямому использованию данного метода препятствуют следующие факторы:
• наличие постоянной составляющей и высокий уровень помехи в спектре вблизи неё;
• неполное выравнивание огибающей;
• наличие 2-й периодической последовательности в случае акустической анизотропии.
Рис. 2. Этапы дополнительной обработки СФ-эхосигнала
Фактор наличия постоянной составляющей решается отбрасыванием низкочастотной монотонно убывающей (с нулевой частоты) области абсолютного спектра ВО СФ-эхосигнала.
ж)
з)
Рис. 3. Этапы обработки эхосигнала при отсутствии акустической анизотропии:
а - необработанный сигнал с выхода АЦП (сигнал на входе блока 1, рис. 1); б - эхосигнал после согласованной фильтрации и отбрасывания прямого сигнала (сигнал на выходе блока 2, рис. 1); в - огибающая СФ-эхосигнала (сигнал на выходе блока 4, рис. 2); г - кривая ослабления СФ-эхосигнала (сигнал на выходе блока 5, рис. 2); д - сигнал на выходе блока нелинейной обработки (блока 6, рис. 2); е - выровненная огибающая СФ-эхосигнала (выход блока 7, рис. 2); ж - спектр ВО СФ-эхосигнала (выход блока 9, рис. 2); з - положение обнаруженных в огибающей СФ-эхосигнала реплик
ж)
з)
Рис. 4. Этапы обработки эхосигнала при наличии акустической анизотропии:
а - необработанный сигнал с выхода АЦП (сигнал на входе блока 1, рис. 1); б - эхосигнал после согласованной фильтрации и отбрасывания прямого сигнала (сигнал на выходе блока 2, рис. 1); в - огибающая СФ-эхосигнала (сигнал на выходе блока 4, рис. 2); г - кривая ослабления СФ-эхосигнала (сигнал на выходе блока 5, рис. 2); д - сигнал на выходе блока нелинийной обработки (выход блока 6, рис. 2); е - выровненная огибающая СФ-эхосигнала (выход блока 7, рис. 2); ж - спектр ВО СФ-эхосигнала (выход блока 9, рис. 2); з - положение обнаруженных в огибающей СФ-эхосигнала реплик
Неполное выравнивание огибающей приводит к искажению формы и фазы линейчатого спектра в окрестностях каждой из гармоник, что увеличивает взаимное влияние гармоник
друг на друга и приводит к ошибке определения положения максимума первой гармоники. Поэтому желательно обеспечить максимальное расстояние между гармониками, что достигается включением в область анализа максимально возможного количества реплик (не менее пяти).
Похожим образом влияет на результат наличие акустической анизотропии. В этом случае необходимо максимально обузить гармонику в спектральной области. Для этого необходимо выполнить предварительную оценку интервала повторения реплик на основании положения 1 -й гармоники абсолютного спектра ВО СФ-эхосигнала, затем выбрать область анализа выровненной огибающей ВО СФ-эхосигнала кратной данной оценке интервала и выполнить повторную обработку блоками 8^9 (рис. 2).
Определение положения первого импульса
Предлагается также использовать информацию о фазе первой гармоники линейчатого спектра [1] ВО СФ-эхосигнала для оценки положения первой реплики эхосигнала. Поскольку на данном этапе уже известна приближённая оценка периода следования импульсов (T), повторное вычисление БПФ ВО СФ-эхосигнала выполняется на интервале, кратном периоду следования импульсов: от Tдо (n + 1)T, где n = 3, 4, 5, ... - число анализируемых импульсов. Интегрирование необходимо начинать с момента T для устранения прямого сигнала излучатель-приёмник.
Рассмотрим спектральную плотность одиночного импульса S(t) на интервале от - 0.5-до + 0.5-T. Обозначим фазу спектральной плотности одиночного импульса на частоте f 0= TЛ как ф0. Требование равенства ф0 = 0 будет соблюдаться при соблюдении следующего условия:
0 0.5T
J S(t )dt = J S (t )dt.
-0.5 • T 0
Например, если одиночный импульс симметричен относительно максимума ф0 будет равно нулю, если максимум расположен в точке t = 0. Далее, положением импульса на временной оси будем называть момент времени tS, делящий импульс на две равные по площади части.
Огибающая спектра периодической последовательности импульсов
...+ S(t - 2T) + S(t - T) + S(t) + S(t + T) + S(t + 2T) + ...
отличается от спектра одиночного импульса только масштабом (вдоль оси ординат), т.е. фаза первой гармоники на частоте f 0 линейчатого спектра периодической последовательности импульсов будет совпадать с фазой спектральной плотности одиночного импульса на частоте f 0.
Сдвиг сигнала во времени на t0 приводит к изменению фазовой характеристики спектра на величину 2л/10.
Следовательно, положение n-го импульса на временной оси будет связано с фазой ф0 первой гармоники линейчатого спектра следующим элементарным соотношением:
t = (Фр + n )т ^ 2п
где n = 0, ±1, ±2, ...
Большинство реализаций быстрого преобразования Фурье (БПФ) (в том, числе и библиотека FFTW [2], используемая в Matlab) производят интегрирование на интервале от 0 до t, а не на интервале от - 0.5-t до + 0.5^t, т.е. анализируется сигнал
S(t) + S(t + T) + S(t + 2T) + .+ S(t + (N - 1)T).
Поэтому приведённая формула останется справедливой только при соблюдении следующих условий: если сигнал представляет собой периодическую последовательность из N импульсов, то длительность сигнала, поступающего на вход БПФ, должна быть строго кратна T. Желательно также, чтобы интервал анализа не превышал NT:- в этом случае все гармоники линейчатого спектра будут иметь минимальную ширину и максимальную энергию в частотной области.
Экспериментальная проверка метода
Метод был опробован при определении времени распространения в семи стальных заготовках различной толщины (11^75 мм), имеющих различную акустическую анизотропию. Проверка проводилась для продольных и сдвиговых волн. Метод показал устойчивую работу даже при малом количестве откликов (рис. 5).
а)
от» ли*
r
ft»»" •
б)
Рис. 5. Результаты обработки эхосигнала, полученного при зондировании стальной заготовки толщиной 75 мм, при сильном поглощении в тракте
Выводы
Предложен метод автоматического определения периода следования импульсов в периодической последовательности и положение первого импульса в этой последовательности. Применяя различные методы интерполяциии, можно добиться высокой точности определения временного положения некоторой характерной точки отдельного пика (например, максимума огибающей) и, соответственно, интервалов между пиками. Однако для этого необходимо с достаточной степенью надёжности определить положение самих пиков.
Предлагаемый метод может быть использован в качестве вспомогательного для предварительной приближённой оценки положения хорошо локализованных (коротких) импульсов периодической последовательности на временной оси с последующим уточнением положения (например, максимума) с помощью различных методов интерполяции. Для слабо локализованных (длинных) импульсов (например, являющихся откликами узкополосного тракта) данный метод может использоваться в качестве основного.
Библиографический список
1. Гоноровский, И.С. Радиотехнические цепи и сигналы: учебник для вузов / И.С. Гоноров-ский. - М.: Радио и связь, 1986. - 512 с.
2. [Электронный ресурс] // www.fftw.org
Дата поступления в редакцию 28.04.2011
E.A. Bukvarev, M.V Martinuk,
USING PROPERTIES OF THE LINE SPECTRUM OF THE REPETITION PULSE TRAIN FOR AUTOMATIC ECHO SIGNALS PROCESSING
This article deal with solving of the problem robust and qualitative definition temporary situation of the separate replicas in echo signal received from acoustic tract under reverberations.
Key words: ultrasonic sounding, reverberating medium, acoustic anisotropy.
МЕХАНИКА ЖИДКОСТИ, ГАЗА И ПЛАЗМЫ
УДК 517.951
И.С. Костенко
ИССЛЕДОВАНИЕ ВОЗДЕЙСТВИЯ УДАЛЁННЫХ ЦУНАМИ НА ПОБЕРЕЖЬЕ САХАЛИНСКОЙ ОБЛАСТИ
Специальное конструкторское бюро средств автоматизации морских исследований ДВО РАН
Проведено численное моделирование распространения исторических цунами сейсмического происхождения в Тихом океане. Рассматривались события с эпицентрами удаленными от Дальневосточного побережья России. Была проведена оценка времени прихода волн от удаленных цунами, силы воздействия на побережье, рассмотрены особенности прохождения цунами через Курильские проливы. Приведен спектральный анализ записей.
Ключевые слова: цунами, инструментальные данные, численное моделирование, Сахалин.
Чилийское землетрясение - сильнейшее землетрясение в истории наблюдения, магни-туда, по разным оценкам от 9,3 до 9,5, произошло 22 мая 1960 года в 19 ч 11мин по местному времени и послужило причиной цунами. Возникшие волны достигали высоты 15-20 м у побережья Южной Америки. Цунами разрушило несколько городов и смыло десятки селений. Погибло 10 тысяч человек. Еще одно землетрясение, произошедшее 27 февраля 2010 года в 03 ч 34 мин по местному времен у берегов Чили, было магнитудой 8,8 и также послужило причиной цунами.
Эти цунами привели к катастрофическим последствиям на побережье Южной Америки. Волны цунами дошли до Дальневосточного побережья России.
Численное моделирование
Для моделирования этих событий использовался вычислительный комплекс NAMI-DANCE [2, 9, 12], рекомендованный ЮНЕСКО для расчета распространения цунами [10]. Математическая модель основана на нелинейных уравнениях мелкой воды, записанных в сферических координатах, и является модификацией известного японского кода TUNAMI. Она уже использовалась для моделирования цунами в 2004 году в Индийском океане [1] и оценки цунами риска Япономорского побережья России [3].
Пространственный шаг исходной расчетной сетки для акватории Тихого океана 2,7 мин. Такой шаг слишком велик для прохождения волны через Курильские проливы, поэтому для моделирования волны в Охотском море необходимо использование вложенных сеток. Этот метод основывается на том, что каждая сетка, находящаяся внутри большей, имеет шаг в три раза меньше. Для моделирования этих случаев цунами использовали три вложенные сетки. Шаг сетки для южной части Охотского моря и Курильских островов 90 с, для Курильских проливов 30 с (рис. 1). Батиметрии получены из одноминутной батиметрии GEBCO1 Digital Atlas.
Источники чилийских цунами 1960 и 2010 годов (рис. 1) выбраны в соответствии с решением Окада [7], использующим данные о параметрах землетрясения, которые приводятся на сайте Геофизической службы США [8]. Эпицентр наиболее сильного землетрясения, произошедшего 22 мая 1960 года, локализован в точке с координатами 76.50 з.д., 39.50 ю.ш.
© Костенко И.С., 2011.
на глубине 53 км. Длина разрыва, определяемая по афтершокам, 800 км, его ширина 200 км и смещение по разрыву 30 м. Угол разлома с меридианом (strike) 10°, угол смещения плиты вглубь от разлома (dip) 10° и вдоль разлома (slip) 91°. Эпицентр землетрясения, произошедшего 27 февраля 2010 года, севернее от предыдущего в точке с координатами 72.72° з.д. и 35.846° ю.ш. на глубине 55 км. Длина разрыва 450 км, ширина 100 км, смещение по разрыву 15 м. Угол разлома (strike) 16°, угол смещения плиты вглубь (dip) 14°, вдоль разлома (slip) 1°4°.
Рис. 1. Источники чилийских цунами 1960 и 2010 годов и расположение вложенных сеток для Охотского моря и Курильских островов
Рассчитанные начальные волны представляют собой знакопеременное смещение с максимальным поднятием в результате землетрясения 1960 года на 10.8 м и опусканием на 6 м, так что максимальная высота волны в очаге составила около 17 м. Максимальная высота волны в очаге, образовавшемся 27 февраля 2010 года, - около 6 м.
Накат волны на берег не моделировался, и в расчетах используются условия полного отражения (приближение вертикальной стенки на глубине 5 м). Цунами моделировалось на протяжении 35 ч с момента начала землетрясения.
В связи с тем, что источники землетрясения располагались недалеко друг от друга, по сравнению с расстоянием от источников землетрясения до Курильских островов и о. Сахалин, то время прохождения волны через Тихий океан составило порядка 21-22 часа от обоих источников.
Однако высоты волн, образовавшиеся от цунами 2010 года, были на порядок меньше, например, в Северо-Курильске (о. Парамушир) в 1960 году высота цунами составила 4,7 м, в п. Малокурильское (о. Шикотан) 4 м, а в 2010 году высоты волн составили 2,28 м и 1,1 м соответственно.
В 1960 году на Курильских островах при максимальном подъеме уровня воды в г. Северо-Курильск были затоплены дома на побережье, склады и причал, смыты мост и некоторые строения в Южно-Курильске. В Малокурильском были затоплены причалы и некоторые здания, мост вынесло на сушу. Суда сорваны с якорей. При обоих цунами волна проникла в Охотское море.
Проявления Чилийских цунами 1960 и 2010 годов на Курильских островах и о. Сахалин приведены в табл. 1 [11].
Таблица 1
Высоты волн цунами, проявившиеся на побережье Сахалина и Курильских островов,
вызванные землетрясениями в Чили
Наименование населенного пункта Высота цунами в 1960 году, м Высота цунами в 2010 году, м
г. Северо-Курильск (о. Парамушир) 4,7 2,28
п. Малокурильское (о. Шикотан) 4 1,1
г. Курильск (о. Итуруп) 1,0 0,1-0,15
г. Поронайск (о. Сахалин) 0,8 0,25
п. Стародубское (о. Сахалин) Измерения не проводились 0,2-0,3
г. Корсаков (о. Сахалин) 1,2 0,2
г. Холмск (о. Сахалин) 0,1 0,08
Максимальные распределения амплитуд волн цунами от обоих источников схожи по географии расположения. Максимальные волны приходятся на Курильские острова (рис. 2), которые являются естественной преградой на пути распространения волн цунами, и в Охотское море цунами заходит, уже потеряв достаточное количество энергии. В настоящее время ведется работа по вычислению этих потерь. Анализ наблюдений за высотами волн цунами в северо-западной части Тихого океана это подтверждает.
Рис. 2. Максимально распределение амплитуд Чилийского цунами 27 февраля 2010 года
Сопоставление измеренных и расчетных данных
В настоящей работе проведено сравнение численного моделирования Чилийского цунами 2010 года с натурными записями уровня моря. Для изучения особенностей проявления цунами в разных пунктах Сахалинской области были использованы записи уровня моря в г. Корсаков и п. Стародубское, предоставленные Сахалинским центром цунами и запись уровня в п. Охотское в ковше рыбозавода, которая проводилась ИМГиГ ДВО РАН. Чтобы сделать сравнительный анализ произошедшего события с данными моделирования в этих населенных пунктах при выполнении моделирования были поставлены компьютерные мареографы.
На рис. 3 представлены совмещенные записи уровня моря в г. Корсаков (о. Сахалин) по данным наблюдений и моделирования в первые часы прихода цунами. Из записи наблюдений вычтен предвычисленный прилив. Наблюдается достаточно хорошее совпадение
уровней. Цунами началось с небольшого поднятия на 5-10 см спустя 25 ч после начала землетрясения. Период колебаний длился около 4 ч. Время прихода волны при моделировании и в записи наблюдения практически совпадают.
Рис. 3. Сопоставление записей Чилийского цунами 2010 года в г. Корсаков (по данным наблюдения и моделирования)
На рис. 4 представлены совмещенные записи уровня моря в п. Стародубское (о. Сахалин) по данным наблюдений и моделирования. Из записи наблюдений вычтен предвычис-ленный прилив. Из графика видно достаточно хорошее совпадение высот и периодов первых волн записи моделирования и наблюдения, характерный период колебаний составляет около 100 мин. Волна цунами пришла спустя примерно 21 час после землетрясения, и цунами и началось с небольшого подъёма (5 см).
Уровень моря, см
30 п
-20 V-,-,-,-,-,-,-,-,-,-,—
3:00 6:00 !>:00 12Л0 15:00
Рис. 4. Сопоставление записей чилийского цунами 2010 года в п. Стародубское (по данным наблюдения и моделирования)
На рис. 5 представлены совмещенные записи уровня моря в п. Охотское (о. Сахалин) по данным наблюдений и моделирования. Из записи наблюдений вычтен предвычисленный прилив. Из графика видно достаточно хорошее по времени приходу цунами и периоду колебаний. Характерный период колебаний составляет около 100 мин. Волна цунами пришла
спустя примерно 21 ч после землетрясения, и цунами началось с небольшого подъёма (5 см). Амплитуды колебаний уровня моря модельного ряда больше, так как модельный (компьютерный) мареограф располагался в открытом море, поэтому высоты получились больше, чем у стационарного мареографа, расположенного в закрытом ковше рыбозавода.
Рис. 5. Сопоставление записей чилийского цунами 2010 года в п. Охотское (по данным наблюдения и моделирования)
Датчики измерения уровня моря располагаются в ковше порта, а компьютерные мареографы - вне его на глубине около 10 м. Этим объясняется немногим более ранний приход волны цунами по данным моделирования.
Также было проведено сравнение по времени прихода волн цунами и их высоте для населенных пунктов Курильских островов, в которых поводились наблюдения [4, 5, 6]. Анализ показывает достаточно хорошее совпадение записей цунами.
На рис. 6 представлено сравнение спектральной плотности мощности для рассчитанных записей Чилийских цунами 1960 и 2010 года в районе г. Северо-Курильска (Курильские острова). Отсюда видно, что усиление наблюдается на одних и тех же периодах: 85, 32 и 19 мин.
Рис. 6. Оценка спектральной мощности плотности для полученных модельных записей волновых процессов (г. Северо-Курильск)
Выводы
Выполнено моделирование чилийских цунами 1960 и 2010 годов с помощью вычислительного комплекса, основанного на уравнениях мелкой воды, в сферических координатах с учётом вращения Земли. Проведено сравнение рассчитанных и реально наблюдавшихся
записей уровня моря для пунктов на побережье о. Сахалин и Курильских островах. Также было получено практически полное согласие расчета времени прихода волны в населенные пункты, которые удалены на значительное расстояние от эпицентров цунами. Сравнительный анализ рассчитанных записей цунами на шельфе Курильских островов показывает, что параметры волны на удалённых пунктах задаются формой шельфа, и для случаев чилийских цунами они являются схожими.
Выполнение этой работы поддержано грантом РФФИ 09-05-00971
Библиографический список
1. Зайцев, А.И. Моделирование распространения катастрофического цунами (26 декабря 2004 г.) в Индийском океане / А.И. Зайцев [и др.]// Доклады РАН. 2005. Т. 402. № 3. С. 388-392.
2. Куркин, А.А. Модифицированный вычислительный комплекс «ЦУНАМИ» для оценки рисков, связанных с цунами / А.А. Куркин [и др.] // Известия АИН им. А.М. Прохорова. Прикладная математика и механика. 2004. Т. 9. С. 88-100.
3. Куркин, А.А. Сравнительная оценка цунами опасности япономорского побережья России на основе численного моделирования / А.А. Куркин [и др.] // Океанология. 2004. Т. 44. № 2. С.179-188.
4. Левин, Б.В. Первые результаты регистрации длинных волн в диапазоне периодов цунами в районе Курильской гряды на разнесенной сети станций / Б.В. Левин [и др.]// Доклады Академии наук. 2009. Т. 427. № 2. С. 239-244.
5. Файн, И.В. Расчет захваченных волн для района Курильской гряды // В сб.: Волновые процессы в северо-западной части Тихого океана. - Владивосток: ДВНЦ АН СССР, 1980. С. 87-92.
6. Файн, И.В. Частотные свойства Курильского шельфа // Генерация цунами и выход волны на берег. - Москва. 1984. С. 80-83.
7. Okada, Y. Surface deformation due to shear and tensile faults in a half-space // Bull. Seism. Soc. America. 1985. V. 75. P. 1135-1154.
8. Yalciner, A.C. Modeling and visualization of tsunamis: Mediterranean examples / A.C. Yalciner [at al.]// Tsunami and Nonlinear Waves (Ed: Anjan Kundu), Springer, 2007. P. 273-283.
9. сайт Геофизической службы США (http://earthquake.usgs.gov/earthquakes /world/events/1960 05 22.php)
10. http://ioc3.unesco.org/ptws/21 //(documents)//TsuModelMan-v3-ImamuraYalcinerOzyurt apr06.pdf.
11. http : //avi -nami. ce.metu.edu .tr.
12. http ://sakhmeteo .ru/company/structure/tsunamicenter
13. http : //namidance. ce .mrttu.edu .tr
Дата поступления в редакцию 29.04.2011
I.S. Kostenko
RESEARCH OF IMPACT OF REMOTE TSUNAMI AT THE SAKHALIN REGION COAST
In the present work numerical modeling distribution of historical tsunamis of a seismic origin in Pacific Ocean is lead. Events with an epicenter removed from Far East coast of Russia were considered. The estimation of arrival time of waves from the removed tsunami, forces of influence at coast has been lead, features of passage of a tsunami through Kuriles passages are consider. The spectral analysis of records in settlements of Sakhalin and Kuriles is lead.
Key words: tsunami, gauge records, numerical modeling, Sakhalin.
