RU2457514C1 - Способ определения предвестника цунами - Google Patents

Abstract

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

Description

Изобретение относится к области сейсмологии, а именно к способам определения предвестника цунами, а более конкретно к способам регистрации преимущественно акустических сигналов, предвестников цунами.

Известны способы определения предвестника цунами (патент RU №22082184 [1], патент RU №2066467 [2], авторское свидетельство SU №1300393 [3]), включающие формирование упругих колебаний, их регистрацию, сравнение их с опорным сигналом, определение сейсмических параметров среды.

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

Известен также способ определения опасности цунами (патент RU №22066466 [4]), включающий размещение в прибрежной зоне на глубине более 100 м групп устройств регистрации, соединение их трактом связи с наземными станциями приема и обработки сигналов путем поэтапного определения опасности цунами. На расстоянии 2-4 тысяч километров от берега устанавливают еще одну группу устройств регистрации, а группу устройств в прибрежной зоне размещают на расстоянии, обеспечивающем необходимое время на защиту охраняемого района, определяемом на основании формульной зависимости. Факт возникновения цунами устанавливают по сигналам дальних устройств, а по сигналам ближних устройств регистрации, установленных в прибрежной зоне на глубине более 100 м, определяют степень опасности волны цунами для охраняемого района.

Поэтапное определение опасности цунами обеспечивает повышение надежности прогноза цунами по сравнению с аналогами [1-3], однако размещение регистрирующих устройств на глубинах более 100 м огранивает информативность получения первичных сигналов и как следствие снижает достоверность прогноза, так как известно, что наибольшая информативность первичных сигналов наблюдается на глубинах 6-10 м от уровня прилива, вблизи берега и вдоль континентальных шельфов.

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

Для повышения информативности в способе сейсмического микрорайонирования (авторское свидетельство SU №251694 [5]), включающем размещение исследуемых и опорных пунктов наблюдений на участках с различными инженерно-геологическими условиями, регистрацию в них сейсмических колебаний от землетрясений из потенциально опасных и других очаговых зон, определение динамических параметров сейсмических колебаний и их вариаций в каждом исследуемом пункте наблюдений относительно опорных в заданном частотном диапазоне исследований, дополнительно проводят трехкомпонентную регистрацию сейсмических колебаний по ортогональной ориентированной на потенциально опасные очаговые зоны сети профилей, с расстоянием между пунктами наблюдений, не превышающим 1/3-1/4 длины волны наиболее высокочастотных сейсмических колебаний, образующих информативные вариации амплитуд, а расстояние между профилями составляет 1/3-1/4 минимального пространственного периода информативных амплитудных вариаций высокочастотного диапазона сейсмических колебаний.

Однако ввиду того, что в общем случае величина амплитуд давления сейсмических сигналов зависит от величины вызвавшего сигнал вертикального смещения дна (определяемого произведением скорости смещения на длительность импульса); волнового сопротивления вод (определяемого произведением плотности воды на скорость звука); угла преломления акустической волны, вышедшей из дна в воду, а также удаления горизонта наблюдений от дна, достоверные сигналы могут регистрироваться на высоких частотах (50…80 Гц и выше), что ограничивает применение известного способа только при размещении пунктов наблюдений на участках с однородными инженерно-геологическими условиями, что существенно снижает информативность данного способа.

Повышение информативности достигается в способе (патент RU №1787273 [6]), заключающемся в задании региональных кусочно-непрерывных профилей, ориентировании их в крест простирания изучаемых тектонических элементов, задании поперечных профилей и проведении по ним наблюдений, в котором региональные профили задают в виде пар квазипараллельных кусочно-непрерывных профилей, а поперечные - в виде пересекающихся друг с другом кусочно-непрерывных профилей, ориентируют поперечные профили по простиранию изучаемых тектонических элементов, создают вокруг этих элементов замкнутый полигон, при этом положение каждой последующей пары профилей уточняют после получения данных в предыдущей паре профилей, а расстояние между региональными профилями определяют по размерам изучаемых тектонических элементов, что повышает информативность за счет возможности изучения сложно построенных сред.

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

В известном способе сейсмического микрорайонирования (патент RU №1787276 [7]), заключающемся в возбуждении сейсмических колебаний невзрывным импульсным источником, регистрации их сейсмоприемниками, расположенными на участках с различными инженерно-геологическими условиями, определении значения скоростей поперечных волн, частотных характеристик зарегистрированных колебаний и оценке на основе этих характеристик приращения балльности, дополнительно возбуждают повышенные сейсмические колебания, по сравнению с первоначальными колебаниями, а в качестве величины, характеризующей частотную характеристику колебаний, используют величину, обратную средневзвешенному периоду в полосе частот 0,3-30 Гц, определяют приращение балльности для дополнительного возбуждения и вводят величину разности балльностей в качестве поправки за нелинейные эффекты в полученные ранее данные наблюдений, использовавшие маломощный, сейсмический источник, что повышает надежность и точность за счет более полного учета нелинейных свойств грунта.

Существенным недостатком данного способа является необходимость создания развивающегося напряжения в грунте не менее 0,1 и 5 кг/см², что в условиях морского дна является сложной технической задачей.

Известен также способ сейсмической разведки (RU №1787275 [8]), который включает возбуждение и регистрацию интерференционной системой сейсмических сигналов по системе многократного профилирования и обработку полученных данных, в котором по сейсмограммам, полученным в результате предварительных работ на участке профиля исследований, строят скоростно-угловые спектры из соотношения временной задержки для годографа от двойного времени пробега волны по нормали к отражающей границе, удаления взрыв - прибор, эффективной скорости до границы и угла наклона границы, по построенным спектрам выделяют основные сейсмические волны и проводят последующие сейсмические работы на профиле для выбранных параметров сейсмических волн интерференционной регистрирующей системой с оптимальной характеристикой направленности, параметры которой определяют из соотношения в зависимости от текущего угла, кратности интерференционной системы, опорной частоты сигнала и запаздыванием между двумя годографами для углов наклона и текущих.углов, что повышает эффективность сейсмической разведки в сложно построенных средах. Однако технический эффект данного способа может быть получен только в сухопутных условиях при отсутствии воздействия окружающей среды.

В известном способе сейсмической разведки (патент RU №1787274 [9]), включающем разделение геологического объекта на глубинные этажи, определение наибольшей частоты сейсмических сигналов, приходящих с каждого этажа, расчет для каждого глубинного этажа шага квантования по времени меньшего чем 1/4 наибольшей частоты сейсмического сигнала для i-го глубинного этажа, а по пространству меньшему или равному отношению длины волны наибольшей частоты сейсмического сигнала для i-го этажа к углу наклона фронта приходящей волны, возбуждение, прием группами сейсмоприемников, цифровую регистрацию с рассчитанными шагами квантования по времени и пространству для каждого глубинного этажа и обработку принимаемых сигналов, в котором для каждого глубинного этажа определяют длину базы группирования сейсмоприемников по выражению, в котором длина базы равна или меньше отношения длины волны наибольшей частоты сейсмического сигнала для i-го этажа к четырем синусам угла наклона фронта приходящей волны, что позволяет повысить отношение сигнал/помеха на этапе приема и повысить точность исследования при приеме сигналов.

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

В известном способе обнаружения возможности наступления катастрофических явлений (патент RU №2030769 [10]), включающем измерение параметра геофизического поля в контролируемом районе и суждение по полученным данным о возможности наступления катастрофических явлений, в котором измерения производят непрерывно, выявляют колебания измеряемого параметра и при обнаружении синусоидальных колебаний возрастающей частоты, имеющих амплитуду, статистически достоверно отличающуюся от фоновой для контролируемого района, и период от 100 до 1000000 с, судят о наличии возможности наступления катастрофических явлений, что повышает достоверность прогноза. Однако непосредственное использование этих сигналов в качестве прямых предвестников цунами осложнено наличием помех, создаваемых шумами морской среды различного происхождения, в связи с чем возникает задача выделения подводных сейсмических сигналов на фоне шумов морской среды.

Наиболее широкий спектр сигналов, последующая обработка которых позволяет получить более достоверный прогноз наступления катастрофических явлений, можно получить при использовании способа сейсмической разведки (патент RU №2030766 [11]), включающего возбуждение упругих колебаний, их регистрацию сейсмоприсмпиками, каждый из которых содержит три датчика, расположенных под углом 45 градусов к горизонтальной плоскости, и обработку полученных записей с выделением полезного сигнала, в котором одновременно возбуждают упругие волны Р и S-типа, регистрацию проводят сейсмоприемниками, каждый из которых дополнительно содержит четвертый датчик, при этом все датчики равномерно распределены по азимуту, при обработке полученных записей рассчитывают прямоугольные декартовые координаты полного вектора волнового поля в каждой точке приема путем сравнения модулей декартовых проекций, рассчитанных в каждой точке приема, с модулем полного вектора в данной точке приема, выделяют три монотипные линейно-поляризованные волны РР-, SV-, SH-типа и нелинейнополяризованную волну, которые используют в качестве полезного сигнала.

Регистрация сигналов фазы (РР), характеризующей приход волн, связанный с распространением в земной коре волн сжатия, и фазы (S), характеризующей вторичный приход волн, связанный с распространением в земной коре поперечных волн, скорость распространения которых приблизительно в два раза меньше скорости распространения продольных волн, существенно повышает надежность прогноза предвестника землетрясений. Однако в тех случаях, когда очаг землетрясения достаточно близок ко дну океана, в гидроакустических сигналах присутствует Т-фаза, третичный приход волн, скорость которого близка к скорости звука в воде (Walker D.A. and Bernard E.N. Comparison of T-Phase Spectra and Tsunami. Amplitudes for Tsunamigenic and Other Earthquakes. J. Geophys. Res., 98, № C7, p.12557-12565, 1993 [12]). Подводные же землетрясения, очаги которых достаточно близки к поверхности дна океана, могут обуславливать его значительные вертикальные подвижки, вызывающие на поверхности воды гравитационные волны, которые, в свою очередь, распространяясь в мелкой воде прибрежной зоны океана, могут вызывать цунами (Okal E.A. Seismic parameters controlling far-field tsunami amplitudes: a review. Nat. Hazards. 1, p.69-96, 1988 [13]). В данном способе, как и в известных способах [1-10], сигналы, обусловленные третичным приходом волн, не учитываются, что существенно снижает вероятность прогноза.

В известном способе определения предвестников цунами (патент RU №2292569 [14]), включающим размещение в прибрежной зоне и на удалении от прибрежной зоны групп устройств регистрации сигналов на глубинных горизонтах наблюдений, равномерно распределенных по азимуту, соединение их трактом связи с внешними станциями приема и обработки сигналов путем поэтапного определения опасности возникновения цунами с определением динамических параметров сейсмических колебаний и их вариаций в каждом исследуемом пункте наблюдений в заданном частотном диапазоне с обработкой регистрируемых сигналов в высокочастотном и низкочастотном диапазонах сейсмических колебаний с выделением фаз сигналов, характеризующих приход волн, связанный с распространением в земной коре волн сжатия и поперечных волн, в котором дополнительно регистрируют гидроакустические сигналы с Т-фазой; устройства регистрации размещены на горизонтах наблюдений, кратных 25 м, при максимальном горизонте наблюдений, равном 100 м; определяют приход акустической волны сейсмического происхождения по величине сдвига частоты рассеянного излучения, при этом посредством средств регистрации, размещенных в ближней зоне от очага землетрясения, выполняют анализ низкочастотных составляющих рассеянного сигнала, а в качестве опорных квазигармонических высокочастотных сигналов используют шумы судоходства; посредством средств регистрации, размещенных в прибрежной зоне, момент появления и направление прихода сейсмоакустических волн определяют путем узкополосной фильтрации и спектрального анализа волн на комбинационных частотах.

Совокупность отличительных признаков данного способа по сравнению с известными способами [1-13], заключающихся в регистрации гидроакустических сигналов с Т-фазой, размещении устройств регистрации сигналов на горизонтах наблюдений, кратных 25 м, при максимальном горизонте наблюдений, равном 100 м, определение прихода акустической волны сейсмического происхождения по величине сдвига частоты рассеянного излучения, при этом посредством средств регистрации, размещенных в ближней зоне от очага землетрясения, выполняют анализ низкочастотных составляющих рассеянного сигнала, а в качестве опорных квазигармонических высокочастотных сигналов используют шумы судоходства, а посредством средств регистрации, размещенных в прибрежной зоне, момент появления и направление прихода сейсмоакустических волн определяют путем узкополосной фильтрации и спектрального анализа волн на комбинационных частотах, что существенно повышает информативность за счет регистрации более широкого спектра сигналов путем выделения подводных сейсмоакустических сигналов на фоне шумов морской среды и как следствие этого надежность и достоверность прогноза вероятности возникновения цунами.

В тоже время возбуждение волн цунами землетрясениями в сжимаемой жидкости сопровождается генерацией гидроакустических полей в более широком частотном диапазоне, а их энергия может превышать энергию волн цунами. При этом низкочастотные поля (F<1 Гц), так же как и волны цунами, возбуждаются в основном за счет вертикальных подвижек дна в эпицентре землетрясения. Возбуждение высокочастотных гидроакустических полей (фаза Т) происходит на значительно большей площади и существенно зависит от рельефа дна. Поэтому высокочастотные гидроакустические поля содержат относительно меньше информации собственно о формах движения дна в эпицентре землетрясения. Наличие интенсивных низкочастотных акустических полей в очаге цунами наблюдается в диапазоне 0,05-0,4 Гц, при этом основная энергия упругих колебаний превышает энергию волны цунами примерно в 300 раз (Левин Б.В., Носов М.А. Физика цунами, М.: Янус - К, 2005. 360 с.). На очень низких частотах (ниже 0,01 Гц) вследствие пренебрежимо малой толщины слоя океана по сравнению с длиной волны происходит непосредственное возбуждение анемобарических волн за счет перепадов атмосферного давления (Рабинович А.Б. Длинные гравитационные волны в океане: захват, резонанс, излучение. СПб.: Гидрометеоиздат, 1993, 325 с.). Часть энергии микросейсм, распространяющихся под малыми углами к вертикали, рассеиваются в толще Земли в виде объемных волн в соответствии с законом 1/R². Другая часть микросейсм вследствие рефракции или отражения от нижележащих слоев возвращается к верхней границе и претерпевает повторные отражения и преобразования продольных волн в поперечные и наоборот. При этом могут образовываться поверхностные волны разных типов, которые могут распространяться на большие расстояния с малым затуханием (энергетический коэффициент затухания пропорционален 1/R). При этом образуются волны Релея, Стоунли и Лява. Скорость волн Релея всегда больше скорости звука в воде. Поэтому на достаточно высоких частотах, когда длина волны в водном слое соизмерима с глубиной океана, часть энергии волн Релея переходит в воду. Амплитуда волн при этом уменьшается. Оценки показывают, что влияние слоя воды при глубине океана 4 км начинает сказываться на частотах около 0,01 Гц. На частоте около 0,1 Гц волна, отраженная от поверхности жидкости, проходит на дно в противофазе, т.е. происходит максимальное подавление релеевской волны. Наибольшее затухание при этом претерпевает основная мода, так как ее пучность расположена на границе вода - грунт. Высшие моды затухают меньше, поскольку имеется ряд пучностей этих мод в нижележащих слоях. За счет обмена акустической энергией между жидкостью и упругим основанием при достаточной глубине океана может возникнуть и распространяться вдоль дна поверхностная волна Стоунли. При этом вдоль вертикали по обе стороны границы располагаются неоднородные затухающие волны. При глубине океана 4 км образование волн Стоунли возможно на частотах начиная примерно с 1 Гц, а на частотах выше 10 Гц ограничивающим влиянием глубины океана можно пренебречь. Скорость волны Стоунли меньше скорости волн в воде и грунте, поэтому отсутствуют потери энергии за счет "вытекающих" волн. Отсюда следует, что возможное распространение волн Стоунли вдоль морского дна на большие расстояния на высоких частотах в отличие от волн Релея. Поверхностные волны Лява представляют собой поперечные колебания с горизонтальной поляризацией. Поэтому они не могут непосредственно возбуждаться волнами, падающими на границу из водной среды или за счет перепадов анемобарического давления. Их появление в составе микросейсм связано с преобразованием волн Релея на неоднородностях земной коры, а также с сейсмической эмиссией из коры и верхней мантии.

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

Выявленных недостатков, присущих известным способам [1-14], лишен способ определения предвестников цунами (заявка RU №2010116097 от 29.04.2009 [15]), техническим результатом которого является повышение надежности и достоверности прогноза возникновения цунами за счет регистрации микросейсмических волн.

При этом технический результат достигается за счет того, что в способе определения предвестников цунами, включающем размещение в прибрежной зоне и на удалении от прибрежной зоны групп устройств регистрации сигналов на глубинных горизонтах наблюдений, равномерно распределенных по азимуту, соединение их трактом связи с внешними станциями приема и обработки сигналов путем поэтапного определения опасности возникновения цунами с определением динамических параметров сейсмических колебаний и их вариаций в каждом исследуемом пункте наблюдений в заданном частотном диапазоне с обработкой регистрируемых сигналов в высокочастотном и низкочастотном диапазонах сейсмических колебаний с выделением фаз сигналов, характеризующих приход волн, связанный с распространением в земной коре волн сжатия и поперечных волн, с регистрацией гидроакустических сигналов с РР, S и Т-фазой, и размещением устройств регистрации на горизонтах наблюдений определяют приход акустической волны сейсмического происхождения по величине сдвига частоты рассеянного излучения, при этом посредством средств регистрации, размещенных в ближней зоне от очага землетрясения, выполняют анализ низкочастотных составляющих рассеянного сигнала, а в качестве опорных квазигармонических высокочастотных сигналов используют шумы судоходства; посредством средств регистрации, размещенных в прибрежной зоне, момент появления и направление прихода сейсмоакустических волн определяют путем узкополосной фильтрации и спектрального анализа волн, в котором выделение фаз типа РР, S и Т осуществляют путем узкополосной фильтрации посредством рекурсивных фильтров Баттерфорта, при этом входную фильтрацию осуществляют посредством рекурсивных фильтров с целочисленными коэффициентами, а сигналы с частотой дискретизации 100 Гц и ниже подвергают фильтрации с коэффициентами в виде чисел с плавающей запятой; регистрацию гидроакустических сигналов осуществляют посредством широкополосных донных сейсмографов с не менее чем тремя сейсмическими каналами, при этом сигналы анализируются тремя независимыми детекторами, а сигнал обнаружения вырабатывается при совпадении сигналов тревоги по крайней мере в двух каналах из трех; спектральный анализ выполняют как объемных волн фаз РР и S, так и поверхностных волн Лява, Релея и Стоунли.

Совокупность отличительных признаков, заключающихся в том, что выделение фаз типа РР, S и Т осуществляют путем узкополосной фильтрации посредством рекурсивных фильтров Баттерфорта, при этом входную фильтрацию осуществляют посредством рекурсивных фильтров с целочисленными коэффициентами, а сигналы с частотой дискретизации 100 Гц и ниже подвергают фильтрации с коэффициентами в виде чисел с плавающей запятой; регистрацию гидроакустических сигналов осуществляют посредством широкополосных донных сейсмографов с не менее, чем тремя сейсмическими каналами, при этом сигналы анализируются тремя независимыми детекторами, а сигнал обнаружения вырабатывается при совпадении сигналов тревоги по крайней мере в двух каналах из трех; спектральный анализ выполняют как объемных волн фаз РР и S, так и поверхностных волн Лява, Релея и Стоунли, позволяет повысить надежность и достоверность прогноза возникновения цунами за счет регистрации и последующей обработки всех зарегистрированных сейсмических волн, включая и микросейсмические волны различного происхождения.

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

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

Поставленная задача решается за счет того, что в способе определения предвестника цунами, включающем размещение групп устройств регистрации сейсмических сигналов на глубинных горизонтах наблюдений в прибрежной зоне и на удалении от нее с целью поэтапного определения опасности возникновения цунами, соединение их трактом связи с внешними станциями приема и обработки сейсмических сигналов, регистрацию сейсмических сигналов, в котором устройства регистрации размещают на глубинных горизонтах наблюдений, кратных 25 м, при максимальном горизонте наблюдений, равном 100 м, равномерно распределенных по азимуту, регистрацию сейсмических сигналов выполняют с выделением фаз типа РР, З и Т, приход акустической волны сейсмического происхождения определяют по величине сдвига частоты рассеянного излучения, при этом посредством устройств регистрации, размещенных на удалении от прибрежной зоны, выполняют анализ низкочастотных составляющих рассеянного сигнала, используя в качестве опорных квазигармонических высокочастотных сигналов шумы судоходства, а посредством устройств регистрации, размещенных в прибрежной зоне, определяют момент появления и направление прихода сейсмических волн путем узкополосной фильтрации и спектрального анализа волн, выделение фаз типа РР, S и Т осуществляют путем узкополосной фильтрации посредством рекурсивных фильтров Баттерфорта, при этом входную фильтрацию осуществляют посредством рекурсивных фильтров с целочисленными коэффициентами, а сигналы с частотой дискретизации 100 Гц и ниже подвергают фильтрации с коэффициентами в виде чисел с плавающей запятой; регистрацию сейсмических сигналов осуществляют посредством широкополосных донных сейсмографов с не менее чем тремя сейсмическими каналами, при этом сигналы анализируются тремя независимьми детекторами, а сигнал обнаружения вырабатывается при совпадении сигналов тревоги по крайней мере в двух каналах из трех; спектральный анализ выполняют как объемных волн фаз РР и S, так и поверхностных волн Лява, Релея и Стоунли, в котором, в отличие от прототипа, посредством сейсмографов определяют координаты гипоцентра морского землетрясения и его магнитуду, при магнитуде землетрясения более 6, посредством группы устройств регистрации, размещенных на границе атмосфера - водная поверхность, дополнительно регистрируют скорость и направление ветра и морского волнения, влажность воздуха, атмосферное давление, барический градиент электрических разрядов в атмосфере, частоту звуковых волн в атмосфере, определяют коэффициент корреляции для измеренных значений скорости и направления ветра и морского волнения, влажности воздуха, атмосферного давления, барического градиента электрических разрядов в атмосфере, частоты звуковых волн в атмосфере, по которому оценивают возможность появления волны цунами.

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

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

Так как появление волны цунами характеризуется такими признаками, как:

- ненормально повышенная слышимость звуков в воздухе (1 Гц/с, до 60 Гц);

- резкое падение атмосферного давления в течение 6-12 часов;

- увеличение абсолютной влажности за 4 часа на 2 мм;

- электрические разряды в атмосфере имеют большой барический градиент, порядка 20-30 мбар;

- увеличение скорости ветра (Справочник штурмана. Под ред. В.Д.Шандабылова. М. Военное издательство. 1968. с.362-365), то по трассе возможного перемещения волны цунами, посредством средств измерения, размещенных на средствах регистрации (буйковых станциях) регистрируют такие параметры, как скорость и направление ветра и морского волнения, влажность воздуха, атмосферное давление, барический градиент электрических разрядов в атмосфере, частоту звуковых волн в атмосфере, определяют коэффициент корреляции для измеренных значений скорости и направления ветра и морского волнения, влажности воздуха, атмосферного давления, барического градиента электрических разрядов в атмосфере, частоты звуковых волн в атмосфере, по которому оценивают возможность появления волны цунами.

Способ реализуется следующим образом.

Как и в прототипе, размещают средства регистрации сейсмических сигналов, представляющие собой широкополосные донные сейсмографы непосредственно на границе вода - грунт в прибрежной зоне и на удалении от прибрежной зоны, а также на разных горизонтах по глубине с использованием для установки автономных донных станций, подводных обсерваторий, заякоренных платформ. Аналогом широкополосных сейсмографов являются широкополосные сейсмографы типа ЭХП-17, ЭХП-20.

На береговых станциях устанавливают широкополосные сейсмографы типа "G.Streckeisen Messgeratebau" (Швейцария) и "Guralp" (Англия), типа STS-1 и CMG-3. Выполняют регистрацию сейсмических сигналов с выделением фаз типа РР, S, и Т. При этом регистрируют сейсмические шумы на частотах 0,008-20 Гц, гидродинамический шум моря на частотах 0,01-100 Гц, давление волн цунами на дно на частотах 0,01-0,01 Гц.

Сигнал Т фазы, принятый в береговом клине, определяется в диапазоне частот 34…75 Гц при частоте квантования 160 Гц методом псевдодифференциального параболического уравнения, обеспечивающего выделение звуковых полей в двумерном неоднородном океане с переменными рельефом дна и профилем скорости звука с заданной точностью для любого диапазона углов скольжения локальных нормальных волн с учетом взаимодействия между ними. Так как наблюдаемый сигнал S(t) является суммой сигналов от последовательно возбуждаемых слоев, то, представляя сигнал в виде вектора столбца временных отсчетов и обозначая через s_i векторы-столбцы сигналов от последовательно возбуждаемых слоев, имеем S(S₁, S₂,… S_n), (a₁, a₂,…, a_n), a_i - суть амплитуды рассеивателей. В качестве решающей статистики используют сумму квадратов амплитуд, имеющую максимальное значение для сигнала ожидаемой структуры. Оценку получают методом наименьших квадратов, так как система линейных уравнений является неопределенной. Выполняя оценку для каждого момента времени, получают ее временную зависимость. Присутствие в ней максимума означает наличие в источнике ожидаемой структуры возбуждения звукового поля. При построении графика решающей статистики абсцисса глобального максимума соответствует времени прихода совокупного рассеянного сигнала.

Определяют приход акустической волны сейсмического происхождения по величине сдвига частоты рассеянного излучения, при этом посредством средств регистрации, размещенных в ближней зоне от очага землетрясения на внешних станциях приема и обработки сигналов, выполняют анализ низкочастотных составляющих рассеянного сигнала, а в качестве опорных квазигармонических высокочастотных сигналов используют шумы судоходства, а посредством средств регистрации, размещенных в прибрежной зоне, определяют момент появления и направление прихода сейсмоакустических волн путем узкополосной фильтрации и спектрального анализа. Выделение фаз типа РР, S и Т осуществляют путем узкополосной фильтрации посредством рекурсивных фильтров Баттерфорта, при этом входную фильтрацию осуществляют посредством рекурсивных фильтров с целочисленными коэффициентами, а сигналы с частотой дискретизации 100 Гц и ниже подвергают фильтрации с коэффициентами в виде чисел с плавающей запятой. Регистрацию сейсмических сигналов осуществляют посредством широкополосных донных сейсмографов с не менее чем тремя сейсмическими каналами, при этом сигналы анализируются тремя независимыми детекторами, а сигнал обнаружения вырабатывается при совпадении сигналов тревоги по крайней мере в двух каналах из трех.

Спектральный анализ выполняют как объемных волн фаз РР и S, так и поверхностных волн Лява, Релея и Стоунли.

При регистрации сейсмических сигналов на дне одним из важных направлений использования широкополосных донных сейсмографов является исследование микросейсмических шумов, возбуждаемых морскими и океаническими волнами. Микросейсмы проявляются в широком диапазоне частот и служат естественным фоном, который определяет порог чувствительности сейсмографов. При этом также выполняется регистрация характерных микросейсм с периодом около 6 секунд, также выявляются микросейсмы с периодами 20 и 100 секунд, что позволяет выделить как объемные Р и S волны, так и поверхностные волны Лява (колебания в диапазоне частот 0,0125-0,05 Гц), Релея и Стоунли (1-10 Гц).

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

Сигналы сейсмического происхождения регистрируют с помощью широкополосных сейсмографов с частотным диапазоном от тысячных долей герца до десятков герц, которые получили в настоящее время повсеместное распространение как основной инструмент сейсмологических исследований (Рыков А.В. Моделирование сейсмометра // М.: ОИФЗ РАН, 1995. 87 c. Usher M.J., R.F.Burch, C.Gurlap. Wide-band feedback seismometers // Phys. Earth Planet. Interiors. 1979. v.18, p.38-50. Wielandt E., Stein I.M. A digital very-broad band seismograph // Annales Geophysical. 1986. v.4. p.227-232). С помощью этих приборов регистрируются сигналы как местных, так и удаленных землетрясений и определяют координаты эпицентра землетрясения, путем преобразования по частоте в измерительном канале принимаемого гидроакустического сигнала.

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

Аналогом дрейфующего буя является устройство, которое состоит из корпуса цилиндрической формы, мачты с передающим устройством передачи информации, измерителя параметров ветра, измерителя параметров атмосферного давления с баропортом, датчиков температуры воздуха и воды, маячка, радиолокационного углового отражателя, модуля управления с опционным блоком GPS, блока информационной памяти, центрального модуля с контроллером, измерителя высоты волны и ориентации буя, датчика скорости и направления течения, датчиков определения солености, электропроводности, мутности, содержания кислорода, содержания ионов рН, контроллера процессов окисления/восстановления, источника питания (патент RU №2328757).

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

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

Центральный модуль с контроллером включает встроенный 8-канальный 16-разрядный АЦП типа AD 7715 с внешним входами для подключения датчиков, автономную систему контроля напряжения питания, внутренний температурный датчик на базе кремневого диодного p-n перехода, два компаратора с программируемым опорным напряжением, мультиплексор, последовательный интерфейс стандарта RS-232, три таймера, обеспечивающих измерение частоты относительно опорного кварцевого генератора, и представляет собой процессор с раздельными 14-битной шиной команд и 8-битной шиной данных. Двухступенчатый конвертор позволяет осуществлять выполнение до 35 команд в течение одного машинного цикла. Аналогом является микропроцессор типа PJC 14000.

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

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

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

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

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

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

Источники информации

1. Патент RU №22082184

2. Патент RU №2066467.

3. Авторское свидетельство SU №1300393.

4. Патент RU №22066466.

5. Авторское свидетельство SU №251694.

6. Патент RU №1787273.

7. Патент RU №1787276.

8. Патент RU №1787275.

9. Патент RU №1787274.

10. Патент RU №2030769.

11. Патент RU №2030766.

12. Walker D.A. and Bernard E.N. Comparison of T-Phase Spectra and Tsunami. Amplitudes for Tsunamigenic and Other Earthquakes. J. Geophys. Res., 98, №C7, p.12557-12565, 1993.

13. Okal E.A. Seismic parameters controlling far-field tsunami amplitudes: a review. Nat. Hazards. 1, p.69-96, 1988.

14. Патент RU №2292569.

15. Заявка RU №2010116097 от 29.04.2009.

Claims (1)

1. Способ определения предвестника цунами, включающий размещение групп устройств регистрации сейсмических сигналов на глубинных горизонтах наблюдений в прибрежной зоне и на удалении от нее с целью поэтапного определения опасности возникновения цунами, соединение их трактом связи с внешними станциями приема и обработки сейсмических сигналов, регистрацию сейсмических сигналов, в котором устройства регистрации размещают на глубинных горизонтах наблюдений, кратных 25 м, при максимальном горизонте наблюдений, равном 100 м, равномерно распределенных по азимуту, регистрацию сейсмических сигналов выполняют с выделением фаз типа РР, S и Т, приход акустической волны сейсмического происхождения определяют по величине сдвига частоты рассеянного излучения, при этом посредством устройств регистрации, размещенных на удалении от прибрежной зоны, выполняют анализ низкочастотных составляющих рассеянного сигнала, используя в качестве опорных квазигармонических высокочастотных сигналов шумы судоходства, а посредством устройств регистрации, размещенных в прибрежной зоне, определяют момент появления и направление прихода сейсмических волн путем узкополосной фильтрации и спектрального анализа волн, выделение фаз типа РР, S и Т осуществляют путем узкополосной фильтрации посредством рекурсивных фильтров Баттерфорта, при этом входную фильтрацию осуществляют посредством рекурсивных фильтров с целочисленными коэффициентами, а сигналы с частотой дискретизации 100 Гц и ниже подвергают фильтрации с коэффициентами в виде чисел с плавающей запятой, регистрацию сейсмических сигналов осуществляют посредством широкополосных донных сейсмографов с не менее чем тремя сейсмическими каналами, при этом сигналы анализируются тремя независимыми детекторами, а сигнал обнаружения вырабатывается при совпадении сигналов тревоги по крайней мере в двух каналах из трех, спектральный анализ выполняют как объемных волн фаз РР и S, так и поверхностных волн Лява, Релея и Стоунли, отличающийся тем, что посредством сейсмографов определяют координаты гипоцентра морского землетрясения и его магнитуду, при магнитуде землетрясения более 6 посредством группы устройств регистрации, размещенных на границе атмосфера - водная поверхность, дополнительно регистрируют скорость и направление ветра и морского волнения, влажность воздуха, атмосферное давление, барический градиент электрических разрядов в атмосфере, частоту звуковых волн в атмосфере, определяют коэффициент корреляции для измеренных значений скорости и направления ветра и морского волнения, влажности воздуха, атмосферного давления, барического градиента электрических разрядов в атмосфере, частоты звуковых волн в атмосфере, по которому оценивают возможность появления волны цунами.