Научная статья на тему 'Магнитные индукционные эффекты, сопровождающие подводные землетрясения'

Магнитные индукционные эффекты, сопровождающие подводные землетрясения Текст научной статьи по специальности «Физика»

CC BY
156
35
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЗЕМЛЕТРЯСЕНИЯ / ОКЕАН / МАГНИТНАЯ ИНДУКЦИЯ / КВАЗИСТАТИЧЕСКОЕ ПРИБЛИЖЕНИЕ / САМОИНДУКЦИЯ / МАГНИТОИНДУКЦИОННОЕ ПРИБЛИЖЕНИЕ / EARTHQUAKE / OCEAN / MAGNETIC INDUCTION / QUASI-STATIC APPROXIMATION / SELF-INDUCTION

Аннотация научной статьи по физике, автор научной работы — Смагин Виктор Павлович, Сёмкин Сергей Викторович, Савченко Валерий Нестерович

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

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

Magnetic induction effects caused by an underwater earthquake

An analytical and numerical study of magnetic inductive effects caused by hydroacoustic waves arising from underwater earthquakes is carried out. Sound waves propagating in sea water lead to the generation of electric currents due to the action of the Lorentz force (ponderomotive force) in the geomagnetic field. A magnetic field induced by such currents is considered in the case when a sound wave propagates over long distances from a source in an oceanic waveguide and an essentially important role is played by the waveguide character of its motion in the ocean. The case in which the observation of an acoustic wave field and the magnetic field effects induced by it occur in the vicinity of a nonstationary seismoacoustic point source of frequency that is turned on at a certain instant of time has been studied in detail. Such a source can serve as a model of an underwater eruption or an earthquake with a localized focus. The magnetic field generated by sound waves of both stationary and non-stationary sources can be affected by several factors, including the influence of the bottom and the free surface, as well as the vertical stratification of the marine environment, both on the sound field itself and on the magnetic field induced by it. Calculations of the induced magnetic field were carried out in the quasistatic and magnetoinductive approximations (which takes into account the phenomenon of self-induction of the magnetic field), which in some situations seem unjustified. The space-time structure of the induction effects of the magnetic field for several models of acoustic waves excited by earthquakes is calculated and analyzed in detail, among which the simplest model is a point source in an infinite homogeneous oceanic medium. In the case of a non-stationary source of acoustic waves, calculations have shown that the dynamics of the variation with time of the magnetic field can be divided into two stages-transition and stationary regimes. The amplitude of the induced magnetic field in the stationary mode falls off in inverse proportion to the square of the distance to the source (while the amplitude of the oscillations of acoustic pressure and velocity decreases more slowly inversely proportional to the distance). In a number of cases, there has been established the existence of a magnetic precursor for the arrival of an acoustic wave. Taking into account the self-induction in the magnetoinduction approximation in the majority of the cases considered leads to a decrease in the amplitude of the field oscillations in the steady-state regime in comparison with the quasistatic solution.

Текст научной работы на тему «Магнитные индукционные эффекты, сопровождающие подводные землетрясения»

Технические науки

УДК 550.373 001 dx.doi.org/10.24866/VVSU/2073-3984/2017-4/169-187

В.П. Смагин1, С.В. Сёмкин2

Владивостокский государственный университет экономики и сервиса Владивосток. Россия

В.Н. Савченко3

Дальневосточный федеральный университет Владивосток. Россия

Магнитные индукционные эффекты, сопровождающие подводные землетрясения

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

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

1 Виктор Павлович Смагин - д-р физ.-мат. наук, зав. лаб. геофизических полей; 690014, Россия, Владивосток, ул. Гоголя, 41; e-mail: Li15@rambler.ru.

2 Сергей Викторович Сёмкин - канд. физ.-мат. наук, доцент кафедры информационных технологий и систем; 690014, Россия, Владивосток, ул. Гоголя, 41; e-mail: Li15@rambler.ru.

3 Валерий Нестерович Савченко - д-р физ.-мат. наук, кафедра общей и экспериментальной физики Школы естественных наук; e-mail: vanes.sav@yandex.ru.

V.P. Smagin, S.V. Semkin

Vladivostok State University of Economics and Service Vladivostok. Russia

V.N. Savchenko

Far Eastern Federal University Vladivostok. Russia

Magnetic induction effects caused by an underwater earthquake

An analytical and numerical study of magnetic inductive effects caused by hydroacoustic waves arising from underwater earthquakes is carried out. Sound waves propagating in sea water lead to the generation of electric currents due to the action of the Lorentz force (ponderomotive force) in the geomagnetic field. A magnetic field induced by such currents is considered in the case when a sound wave propagates over long distances from a source in an oceanic waveguide and an essentially important role is played by the waveguide character of its motion in the ocean. The case in which the observation of an acoustic wave field and the magnetic field effects induced by it occur in the vicinity of a nonstationary seismoacoustic point source of frequency that is turned on at a certain instant of time has been studied in detail. Such a source can serve as a model of an underwater eruption or an earthquake with a localized focus. The magnetic field generated by sound waves of both stationary and non-stationary sources can be affected by several factors, including the influence of the bottom and the free surface, as well as the vertical stratification of the marine environment, both on the sound field itself and on the magnetic field induced by it. Calculations of the induced magnetic field were carried out in the quasistatic and magnetoinductive approximations (which takes into account the phenomenon of self-induction of the magnetic field), which in some situations seem unjustified. The spacetime structure of the induction effects of the magnetic field for several models of acoustic waves excited by earthquakes is calculated and analyzed in detail, among which the simplest model is a point source in an infinite homogeneous oceanic medium. In the case of a non-stationary source of acoustic waves, calculations have shown that the dynamics of the variation with time of the magnetic field can be divided into two stages-transition and stationary regimes. The amplitude of the induced magnetic field in the stationary mode falls off in inverse proportion to the square of the distance to the source (while the amplitude of the oscillations of acoustic pressure and velocity decreases more slowly - inversely proportional to the distance). In a number of cases, there has been established the existence of a magnetic precursor for the arrival of an acoustic wave. Taking into account the self-induction in the magnetoinduction approximation in the majority of the cases considered leads to a decrease in the amplitude of the field oscillations in the steady-state regime in comparison with the quasistatic solution.

Keywords: earthquake, ocean, magnetic induction, quasi-static approximation, self-induction. Введение

Многие процессы в Мировом океане (океаносфере), такие как подводные землетрясения, извержения вулканов, подводные взрывы и т.д., сопровождаются генерацией звуковых волн в толще океанских вод. Многочисленные наблюдения [6] свидетельствуют о том, что подводные землетрясения, способные вызывать гигантские волны, цунами, имеют магнитуду не менее семи баллов. Такие землетрясения одновременно являются и мощным источником короткопериодных (0,5-1,5 с) звуковых волн в воде [6; 11]. На распространение длинных акустических волн в океане определяющее влияние оказывает существование дна и свободной поверхности, что может приводить к волноводному распространению звука. Подводный волновод образуется слоем воды в океане, ограниченным сверху свободной поверхностью

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

Лучи, удерживаемые подводным волноводом (или каналом, как часто называют естественные волноводы), не доходят ни до дна, где они могли бы частично проникнуть в грунт, ни до волнующейся поверхности, где они могли бы испытать рассеяние; поглощение же в воде сравнительно мало, и поэтому звук в волноводе распространяется на весьма большое расстояние с малым ослаблением. В качестве примера укажем, что звук взрыва тротилового заряда весом всего 2,7 кг был обнаружен гидрофоном, расположенным в канале на расстоянии 5750 км от взрыва. На преодоление этого расстояния звук затратил более одного часа. Пришедший звук резко отличался от короткого импульса самого взрыва: он растянулся на целую минуту, что соответствует в пространстве протяженности звуковой волны около 90 км. Такое расплывание сигнала характерно для волноводного распространения импульсов: оно вызвано дисперсией скорости звука в волноводе. Мы увидим, что без искажения в волноводе могут распространяться только некоторые типы гармонических волн. Дисперсия приводит к тому, что убывание звукового давления вдоль волновода в негармонической волне происходит не просто соответственно степени расхождения волны, так как волна не только расходится в стороны, но и растягивается вдоль направления распространения. Здесь приходится встречаться с очень сложными законами, а всю задачу о распространении негармонического сигнала приходится решать методом Фурье, разлагая сигнал на его гармонические составляющие.

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

Поскольку эти волны распространяются в проводящей среде в геомагнитном поле, они индуцируют электрические токи, магнитные поля которых могут быть измерены не только в воде, но и над её поверхностью, где нет самих акустических волн [3; 4]. Расчету таких магнитных индуцированных полей и посвящена данная работа.

1. индукционные сейсмомагнитные эффекты нормальных акустических

мод

Задавая тем или иным образом поле акустических скоростей v(r), порожденное подводным землетрясением, согласно закону Ома, можно рассчитать индуцированные электрические токи в морской воде ] = <г[ V, Е], обусловленные силами Лоренца, где и - электропроводность морской воды, Е - вектор индукции геомагнитного поля (везде далее полужирным шрифтом обозначены векторы), а зная электрические токи - найти индуцированное магнитное поле по закону Био-Савара-Лапласа, как многие, и мы в том числе, делали неоднократно ранее [8; 9]:

= F],(r-г')]^. (1Л)

4П |г - г'|3

В монографии [8], а затем в работе [10] было рассчитано индуцированное электромагнитное поле гармонической акустической волны в безграничном океане. Здесь мы продолжим исследование и рассмотрим модель океана как акустического волновода постоянной глубины H, ограниченного сверху свободной поверхностью, а снизу упругим полупространством, представляющим донное основание океано-сферы. в районах постоянной глубины океан можно рассматривать как волновод, ограниченный дном и свободной поверхностью воды. Для низких частот можно пренебрегать неровностями дна и неровностью свободной поверхности, вызванной морским волнением, и считать границы волновода плоскими. Кроме того, можно пренебрегать и неоднородностью среды, вызываемой изменением температуры и гидростатического давления с глубиной. Практически, если при данной частоте возможно распространение лишь нескольких первых номеров нормальных волн, то море можно рассматривать как однородный плоскопараллельный слой, лежащий на упругом полупространстве - морском грунте. Морской грунт вообще - упругое твердое тело, неоднородное по глубине. Найти нормальные волны в волноводе, ограниченном таким упругим телом, весьма сложно. Но некоторые основные черты моря как волновода можно представить себе, упрощая задачу: аппроксимируя грунт жидким однородным полупространством с некоторыми эффективными значениями плотности и сжимаемости (так называемая модель «твердой жидкости» [2]). Пренебрегая стратификацией морской воды, будем считать ее однородной средой, характеризующейся плотностью р0 и скоростью звука c0. Будем полагать, что снизу океан ограничен однородной изотропной средой с плотностью р1 и скоростью звука c1 и что в этой океанской среде распространяются только продольные акустические волны.

вычисленное в рамках этой модели акустическое давление плоской волны частотой с, распространяющейся в направлении оси x без диссипации энергии в грунт, равно [2]:

P = -iP0ék>H sin kz (H - z) ■ eik'x. (1.2)

Оси х и y системы координат горизонтальны, ось z направлена вертикально вверх, начало координат находится на дне. Волновые числа kx и ky связаны соотношением kX + kZ = k2, k = с / c0, при этом kz является решением уравнения

k H 21 -1 + tV(l -n2)k2 -kl

kzH =-n + arctg--, (1.3)

2 mkz

где l = 1,2,3... - номер нормальной моды, m = p1 /p0 , n = c0 /c1, n < 1 .

Используя (1.2) и равенство v =—1— grad P, вычислим компоненты скорости

ipc

волны: Г0

P éKH

vz = ^^kz coskz (H - z) • eikxX, Po®

v x =-iPo^^kx sin kz (H - z) . ek*x, (1.4)

Po®

v, = 0

Подставляя (1.4) в (1.1), получим связь между компонентами Bi индуцированного поля B в точке с координатами (0,0, h) и компонентами F- геомагнитного поля F там же, которую запишем в матричном виде:

г Gydx'dz'

B = D F , где D1 =-I —2-—г, (1.5)

Bi DiF j 1 2П x'2 + (h - z')2

Gj =

(h - z) 0 - (h - z)

0 - v ** ' + v z (h - z') 0

v vz* ' 0 - vx* ' y

(1.6)

Таким образом, величина и направление индуцированного поля определяются величиной и направлением F и ненулевыми компонентами матрицы Dj . Выполнив двойное интегрирование, этим ненулевым компонентам придадим такой вид: при значениях h > H (точка наблюдения находится над поверхностью моря):

Dn = А (кхe-kx (h-H) - e-kxh (kx cos kzH - kz sin kzH)) k

D13 = Al\(kze~K(h-H) -e-k'h(kx sinkzH + kz coskzH)), (1.7) k

D31 = D1U D33 = D13' D22 = D11 + D33 ;

а при значениях h < H (точка наблюдения находится под поверхностью моря, в толще морских вод):

Ai = Ak3 ((н - - ■ e-M (kx cos kzH - kz sin kzH) - 2kz sin kz (H - h))

Аз = -АЩ-к - e~k'h (kx sin kzH + kz cos kzH) + 2 kz cos kz(H-h)), (L8)

Ai = - e~k'h (kx cos kzH - kz sin kzH) + 2 kx cos kz(H-h)),

Аз = e~M (kx sin kzH + kz cos kzH) + 2kx sin kz(H-h)),

D22 = Dn + D33, A =V°aP°e х

2РоС

о

При выборе значения h = Н формулы (1.7) и (1.8) совпадают и приводят к следующим значениям коэффициентов матрицы Оу :

Ai = А^З(кх - e-kxH(kx coskzH - kz sinkzH)), k . ik„

Л k3

D13 = A^kf (kz - e-kxH(kx sinkzH + kz coskzH)), (1.9)

k

D31 = Dn, D33 = Di3, D22 = Dn + D33.

Коэффициенты D.., вычисленные по формулам (1.9), зависят от частоты

а = c0k и номера l нормальной моды. Если рассмотреть первую нормальную моду, которая возбуждается при волновых числах, определяемых неравенством

k > kKP =--, п то численный расчет показывает, что коэффициенты D

КР 2H 'J

(и, следовательно, компоненты индуцированного магнитного поля) убывают при росте к, имея максимум при к = к . Эти максимальные значения таковы:

nmax 2АН 2ч3/2л , 1 / ПП \\

Ai =-(1 - n ) (1 + -exP(—i—г))'

п n 2л11 - n

2iAH 2Л3/2П , пп (L10)

Аз =-(1 -n ) (1 -nexP(--]—r))-

п 2V1 - n2

MqOPQ hf

Величина индукционного магнитного поля Винд имеет порядок Bt

инд " инд

Росо

где Р0 - акустическое давление в точке наблюдения. Подставляя в эту формулу численные значения физических характеристик моря, получим Винд ~ 10-13 Р0.

Таким образом, чтобы оценить возможность измерения Винд , нужно знать порядок акустического давления Р0, соответствующего наблюдаемым акустическим волнам. В работе [11] приведены данные акустических измерений на расстояниях порядка нескольких сотен километров от места землетрясения. Согласно этим данным амплитуды акустических сигналов лежали в диапазоне 10-60 Дб, что

соответствует акустическим давлениям от 2,Ы0-3 до 2,Ы02 Па. Используя эти значения, получим, что величина Винд будет иметь порядок от 10-16 до 10-11 Тл, что на верхних границах значений вполне измеримо современными экспериментальными средствами.

Рис. 1. Зависимости индуцированного магнитного поля (в долях геомагнитного поля Fo) (просуммированного по нормальным модам) от частоты звуковой волны

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

Проанализируем поведение компонент магнитного поля в области частот, где

существует только одна (первая) нормальная мода. Это диапазон частот от у1 до

3v1, где v1 =

1

H Wi

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

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

близких к v1 (рис. 1). Введем углы а и в , определяющие ориентацию вектора F геомагнитного поля следующим образом:

Fx = F0 sin а cos в, Fy = F0 sin а sin в, Fz = F0 cos а. Подставляя эти выражения в (1.5), вычислим амплитуды компонент индуцированного магнитного поля как по координатным осям (Вх, By, Bz ), так и по

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

c

0

2

направлению горизонтальной составляющей геомагнитного поля, которую обозначим как Бг:

Bx = F0tJ D2 (h) + (D2 (h) cos2 p-D223 (h)) sin2 a,

(111)

By = F0^D121 (h) + D323 (h) sinasinp, Bz = FoyjD323 (h) + (D322 (h) cos2 в - D323 (h)) sin2 a,

BT = F0^ D22(h)sin2 a + (D13(h)cosacos в + D33 sin a sin2 в)2-

Перейдем к анализу возможности определить направления распространения звуковой волны по индуцированному магнитному полю. Над поверхностью воды у матрицы D у есть только два независимых элемента. Эти два элемента (в качестве которых можно взять, например, Dn и D13), а также угол в между направлением волны и горизонтальной компонентой геомагнитного поля нужно считать неизвестными, которые должны быть определены из результатов измерений. Собственно говоря, возможность такого определения и есть возможность определить направление акустической волны. Выразим указанные выше три неизвестных элемента через амплитуды индуцированного магнитного поля по трем взаимно перпендикулярным направлениям. Этими направлениями будут: первое - вертикальное (Bz), второе - по направлению горизонтальной компоненты геомагнитного поля (BT) и третье - по перпендикулярному к этим двум направлениям - Bn . Используя выражения (1.8) и (1.11), запишем квадраты этих трех компонент магнитного поля как функции от известных нам величин Dn , D13, в и а - угла между вертикалью и направлением геомагнитного поля. имеем следующий результат:

| Bz |2 = Ц sin2 a cos2 в + D2 cos2 а, | BT |2 = Ц sin2a + Ц2 (cos a cos в + sinasin2в)2, (1.12)

\ I2 = Цз(

| Bn |2 = Ц23 (cos a sin в - sin a sin в cos в)2-

I B |2

Исключая из системы (1.12) Dn и D13 и вводя обозначения a =—и

2 I Bn.

1 2 I z I

b = | ^ , численные значения которых определяются по данным натурных

,2 _ IB I

|B„.

измерений, получим уравнение, из которого определяется угол в :

(a2 cos2 p- b2 )(tg a cos p-1) = ((cos p + tg a sin2 p) cos p +1). (1.13)

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

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

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

Дальнейший анализ поведения компонент индуцированного поля в зависимости от параметров а , в , h и v мы провели численно. Расчеты позволяют сделать следующие выводы:

- все компоненты поля уменьшаются с увеличением высоты над поверхностью воды при любых значениях остальных параметров;

- вблизи морской поверхности зависимость от частоты всех компонент поля является монотонно убывающей при любой ориентации геомагнитного поля;

- вертикальная компонента индуцированного поля Bz монотонно убывает (в 4-6 раз) с ростом частоты при любых значениях остальных параметров. Максимальное значение Bz, которое достигается при частоте v1, может быть вычислено по формулам (1.10) - (1.11);

- зависимость компоненты Bz от z = h / H имеет максимум при некотором zm. При H = 400 м, v = 1,42 (близко к v1 ) и в = 0, zm монотонно увеличивается от 0,31 при а = 0 до 0,69 при а = 90. Все профили Bz (z) при выборе в = const для различных углов а пересекаются в одной точке zc (рис. 2). Такое поведение Bz (z) объясняется формулами (1.11), если hc - корень уравнения D13 (h) cos в = ±D33 (h), то при значении zc = hc / H расчетное поле Bz (z) не зависит от угла а ;

- зависимость компоненты Bz от частоты примерно аналогична частотной зависимости Bz, хотя её убывание не такое быстрое. На больших глубинах, частотная зависимость BT теряет монотонность изменения;

- зависимость BT от глубины моря имеет достаточно сложный вид (рис. 3), причем различный для разных интервалов остальных параметров. Вблизи поверхности BT ( z ) ведет себя так, что при а, меньших некоторого значения (около 20°), Bz растет с глубиной, а при больших значениях а - уменьшается;

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

Б2, по направлению горизонтальной компоненты геомагнитного поля Бт и перпендикулярную к первым двум указанным - Бп;

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

Рис. 2. Зависимость компоненты Б2 индуцированного магнитного поля (в долях геомагнитного поля Fo) от £. Частота - 1,42 Гц, в = 40°, для кривых 1-4, а равно 0, 30,

60 и 90° соответственно

Рис. 3. Зависимость компоненты Бт индуцированного магнитного поля (в долях геомагнитного поля F0 ) от £. Частота - 1,42 Гц, в = 0°, для кривых 1-4, а равно 10, 30,

60 и 80° соответственно

Учет самоиндукции для расчета сейсмомагнитного эффекта поля (магни-тоиндукционное приближение)

Решение (1.1) получено в квазистатическом приближении - без учета каких-либо динамических эффектов. Однако есть основания полагать, что при инфразвуко-вых частотах эти эффекты могут оказаться существенными. Рассмотрим задачу определения индуцированного акустическими звуковыми волнами в волноводе магнитного поля с учетом самоиндукции, т.е. в рамках так называемого магнито-индукционного приближения. Вклад токов смещения Максвелла на этих частотах пренебрежимо мал, поэтому учет динамических эффектов сводится только к учету самоиндукции. Уравнение для определения эффектов самоиндукции имеет тип уравнения теплопроводности [9]:

дВ

^ "у2В = у> Г1, (1.14)

в записи которого за эффект самоиндукции отвечает первое слагаемое приведенного уравнения. В случае, когда скорость V является гармонической функцией времени с частотой со, решение (1.14) принято записывать также в виде гармонической функции с той же самой частотой. В этом случае вместо уравнения типа теплопроводности (1.14) получим уравнение Гельмгольца Ки — V2 и = / с коэффициентом к2 = :

¡^0асоБ - У2Б = ц0аю\ [ V, Г], (1.15)

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

Как известно из курсов математической физики [7], классическое решение уравнения Гельмгольца К и — V и = / с нулевыми на бесконечности граничными условиями, имеет вид:

-кК

г f e

u = — I--dr', где R = r - r',

4П R

так что решение (1.15) можно представить как:

_ rot[y, F]e-к

B _ "4tj r dr . (116)

Полученное интегральное выражение преобразуем далее, используя тождество rot(y^) = (prot у + [Vq>, у] и теорему Стокса J rot AdV = J [dS, A], в результате

чего имеем:

№ г[[ЧО, Р],(г - г')] - г' + \уКг-г1 , В(г) = ^^--1—>_-¿г . (1.17)

4п ^ |г — г'|

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

Однако, как мы сейчас покажем, можно сохранить формулу (1.1), но иначе представлять в ней выражение для скорости акустической волны. Запишем последнюю формулу (1.17) в виде следующего интеграла:

B

• (118)

где символом Vr обозначен дифференциальный оператор Гамильтона, действующий только на радиальную координату вектора r. Во-первых, отметим, что дивергенция поля B , рассчитанного по формуле (1.18), равна нулю, как только и может быть для магнитного поля (из-за вихревого характера магнитного поля). В этом можно убедиться непосредственным вычислением divB , используя тождество div[a, b] = (b,rot a) - (a,rot b) и учитывая, что вектор [v(r'), F] не зависит от координаты r, а rotgrad (e~kR / R) = 0. Далее вычислим ротор вектора B, представленного выражением (1.18), используя известные векторные тождества из теории векторного поля:

rot[a, b] = a div b - b div a + (b, V)a - (a, V)b ,

grad(a, b) = [a, rot b] + [b, rot a] + (a, V)b + (b, V)a, и, снова учитывая, что дифференцирование по координате r не распространяется на векторное произведение [v(r'), F]

, а также то, что всегда rotgrad р = 0, получим: rot B = [v(r'), F]V2eR-dr' + ^JV([v(r'), F], V QRR)dr'' (U9)

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

kR -KR

V2 -— = -4п8( R) + k2-—, R R

окончательно получим, помня, что div B = 0, следующий результат для ротора поля:

2 г рп -— rot B = F][dr'. (1.20)

0 4я R

Таким образом, индуцированное магнитное поле с учетом самоиндукции в магнитоиндукционном приближении, можно рассчитать по закону Био-Савара-Ла-пласа, как было в формуле (1.1), если вместо вектора скорости v(r') использовать выражение:

кК ; ve, V - ■ (121)

2. Магнитные индукционные эффекты от нестационарного сейсмоакусти -ческого точечного источника в глубоком море

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

индуцированное магнитное поле в квазистатическом приближении

Индуцированное магнитное поле B может быть найдено из решения уравнения (1.14), где, как обычно, v - акустическая скорость, F - геомагнитное поле, с - проводимость морской воды. Первое слагаемое в левой части (1.14) описывает эффект самоиндукции при исследовании индукционных эффектов генерации магнитного поля. Квазистатическое приближение заключается в отбрасывании этого слагаемого. В результате для магнитного поля B получается уравнение Пуассона, решение которого (известно как закон Био-Савара-Лапласа) имеет следующий вид:

B(r) = [ >■ - ■">] . (2 1)

считая движение воды сферически симметричным, введем декартовую и сферическую системы координат, направив ось Z по геомагнитному полю F , а оси х и y таким образом, чтобы точка наблюдения с радиус-вектором r лежала в плоскости xz . Тогда декартовые компоненты вектора v(r', t), записанные как функции от сферических координат точки интегрирования r', имеют вид:

vх = v(r', t) sin ff cos p', vy = v(r', t) sin ff sin p ', vz = v(r', t) cos & , а компоненты вектора r - r ', разницы радиус-векторов, имеют следующий вид:

(г - г ') х = r sinff-r' sinff cosp', (г - г ' ) y =-r' sinff sinp',

(г - г ') z = r cos ff - r' cos ff .

Компоненты двойного векторного произведения в числителе подынтегрального выражения (2.1), также необходимые нам для расчетов, будут такими: [[v,F],(r - r')] = -Fv(rt) sin ff cosp '(r cos ff- r'cos ff), [[v, F],(r - r')] = -F v(rt) sin ff sinp '(r cos ff- r' cos ff), [[v, F], (r - r')] = F v(r', t)(r sin в sin ff cos p '-r 'sin2 ff),

а квадрат модуля разности Ir - r 'I будет равен

|r - r'| = r2 + r '2 - 2rr '(sin в sin ff cos^'+ cos fleos ff) Подставляя все эти выражения в (2.1), получим компоненты индуцированного

поля в следующем виде:

ju0aF г v(r ', t )r '2 (r cos в- r' cos ff) sin2 ff cos p' dr'dffdp'

By = 0, 0

в = -

Bz =

4n J (r2 + r'2 -2rr'(sinesinecosp'+ cosecose'))3/2 ju0aF r v(rt)r'2(rsinesinecosp '-r'sin2 e')sin ffdr'dffdp' 4n J (r2 + r'2 -2rr'(sinesinecosp'+ cosecose'))3/2

(2.2)

Эти формулы позволяют определить компоненты магнитного поля в произвольной точке пространства. Интегрирование по 0' и ф в этих интегралах можно провести в общем виде, как это показано в курсе современной электродинамики [1], выразив результат через эллиптические интегралы Лежандра. В случае же, когда точка наблюдения лежит на оси г (д = 0 ), результат интегрирования по 0' и ф может быть выражен в элементарных функциях. В этом случае находим, что горизонтальная компонента Вх = 0, а для вертикальной компоненты Вг получим:

в = ', t )г ,ъаг']- *тЪшв'

z 2 Jn Jn(r2 + r'2 - 2rr' cos fffn

Учитывая, что

sin3 ffdff 4

(r2 + r'2 - 2rr' cos Q')3 3(max(r, r'))

получим

Bz = —

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

2¿u0aF r v(rt)r'

f V(r ' , 3 dr', R(t) = ct. (2.3)

J fmavíV i- ^

3 0 (max(r, r'))

Чтобы вычислить этот последний интеграл, необходимо определить аналитический вид для скорости, стоящей в подынтегральном выражении. Для выбранного нами нестационарного точечного источника частоты с, включающегося в нулевой момент времени, акустическое давление и колебательная скорость определяются выражениями [9]:

0, t < 0, t > г / с -ехр(г^(? - г / с)), 0 < t < г / с'

P(r, t) = \ Po, . , (2.4)

v(r, t) =

0, t < 0, t > r / с

^^(1 - —)exp(c(t - r /с)), 0 < t < r /с, (2.5)

p0c r cor

TL

3

r

где с - скорость звука в морской воде. Введем безразмерные переменные х = а/' / с, г = ахг /с , т = а^ и параметр р = а / а$, где ах = /и0ас2. Тогда, подставив (2.5) в (2.3), получим интересующее нас магнитное поле в квазистатическом приближении (обозначаемом подстрочным индексом Q) для моментов времени т > 0 в форме следующих интегралов:

Вд (г,т) =

- А

24л ег

¡рт т

3 -3

2 0

(х2 - IX /р)е ¡рхёх, т < г

2Тл егрт

3

-¡рх

, (2.6)

| (х - ¡х/ р)е lpхdх + г | (1 - г /(рх))-ёх

\ о

т > г

Ег0 Р0 ¿и0а

где А =--¡=— . Вычислив интегралы в (2.6), получим окончательно:

р0 сы п

Бв (2,Т) =

- 2 а4л

- 2 а4л

ге'

(рг)

-(1 - (1 + грт- (рт)2/3)е-рТ), т< г

(2.7)

(рг)

3(1 - (1 + грг)е~'р + грг-е-грт), т> г

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

2л/П

\Бг

= Л-

( Р? У

рг^п^рг)) . (2.8)

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

индукционные магнитные эффекты с учетом самоиндукции (магнитоин-дукционное приближение)

Для определения динамических эффектов магнитного поля с учетом самоиндукции необходимо решить уравнение (1.14). Для безграничной среды и нулевых начальных условий решение этого уравнения имеет следующий вид [7]:

В(г, г) = (м°а)Г-' Г ( Г[Иг'> О,ГМГ - г')]ехр

16-\/ п 0 (г г)

5/2 '

(г'

I ,|2 Л |г - г| №

4(г - г')

(г',

(2.9)

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

т

z

2

еу с 7

Б(г,г) = - А ^ [ № - 7 ^2)е, ^ 0 Р

(2.10)

где

у = (г - г + х, Ч = Гф-Х , Ч2 = Гф — ,

0 0 х ,

а функция Ф имеет вид:

Ф =

г{ехр(-

(г + ж)2

-) + ехр(-

(^ - ж)2

-) + +

2( у - ж)

[ехр(-

(г + ж)2

О - ехР(-

(г - ж)2

'4 4(у - ж)' ^ 4(у - ж)' жг " х 4 4(у - ж)' ^ 4(у - ж)'

Из выражения (2.10) видно (на это указывает множитель е'рт), что при больших значениях т колебания магнитного поля становятся гармоническими с амплитудой

\В\ = -

I ^ 2

I(^ -у

0 Р

(2.11)

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

критерий применимости квазистатического приближения [8]: с > с, или р > 1.

В таблице представлено отношение 5

К

для различных значений р и

г . Видно, что при больших значениях р это отношение действительно близко к 1. Отношение 5 \Вп\ для различных значений р и г

I I т—/ | ^ т—ю

\ 0.1 0.25 0.5 1 3 4 5 7 10 50

1 0.800 0.683 0.585 0.516 0.566 0.639 0.748 1.046 0.887 0.979

2 0.666 0.530 0.456 0.450 0.841 1.047 0.880 1.048 1.033 0.981

3 0.567 0.439 0.396 0.439 0.934 0.904 1.003 1.019 0.959 0.981

4 0.489 0.378 0.362 0.461 0.885 0.934 1.018 0.992 1.010 0.981

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

5 0.429 0.335 0.341 0.520 0.969 1.007 0.970 0.973 0.988 0.982

7 0.339 0.279 0.325 0.737 0.977 0.973 0.964 0.969 1.005 0.984

10 0.254 0.230 0.357 0.645 0.915 0.986 0.974 1.009 0.989 0.986

, от \Вп\

при нару-

Из этой же таблицы видно, насколько отличается В

I т

шении условия р > 1.

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

у

1

2.5

Рис. 4. Зависимости от Т модулей (в долях геомагнитного поля) B(z,r) (кривая 1) и

^^ ^,т) (кривая 2) при р = 1 и 7 = 1

Для исследования эффектов магнитного поля в переходном режиме как в квазистатическом приближении, так и с учетом самоиндукции мы провели численный расчет ЯеВ(г,т) по формуле (2.10) для различных z и р и сравнили результат с расчетом Яе Бд (х,т) по формуле (2.7). Результаты исследования сформулированы ниже.

Заключение

Проанализировав пространственно-временную структуру магнитного поля, индуцированного нестационарной сейсмоакустической волной, полученного как в квазистатическом приближении, так и с учетом самоиндукции, мы пришли к следующим выводам:

1. Динамику изменения со временем магнитного поля можно разделить на два этапа - переходной и стационарный режимы. Амплитуда индуцированного магнитного поля в стационарном режиме уменьшается обратно пропорционально квадрату расстояния до источника (в то время как амплитуда колебаний акустического давления и скорости уменьшается медленнее - обратно пропорционально первой степени расстояния).

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

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

4. Во всех решениях есть эффект магнитного предвестника - появление магнитного поля в точке наблюдения ранее прихода акустического сигнала в эту точку (рис. 5).

0.025

0.02

О 015

0.01

0.005

-0.005

-0.01

-0.015

-0.02

-0.025

4

6

10

12

14

16

Рис. 5. Зависимости от т величины (в единицах А) индуцированного магнитного поля в квазистатическом приближении (кривая 2) и с учетом самоиндукции (кривая 1) при р = 5 и z = 5. При т = 5 акустический сигнал приходит в точку наблюдения

5. На расстояниях до источника порядка с / «100г квазистатическое решение не дает выраженного переходного режима ни при каких частотах в диапазоне от 0,5 до 10^ . В индукционном решении переходной режим становится заметным начиная с частот порядка и с ростом частоты становится все более выраженным.

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

1. Батыгин В.В., Топтыгин И.Н. Современная электродинамика. Часть 1. Макроскопическая теория. М., 2003. 736 с.

2. Исакович М.А. Общая акустика. М.: Наука, 1973. 496 с.

3. Конторович В.М. О магнитогидродинамическом эффекте в океане // Докл. АН СССР. Сер. геофиз. 1961. № 3. С. 576-579.

4. Конторович В.М. Магнитная гидродинамика океана // Вопросы магнитной гидродинамики и динамики плазмы. Рига: Изд-во АН Латв. ССР, 1962. С. 171-177.

5. Морс Ф., Фешбах Г. Методы теоретической физики. Т. 2. М.: ИИЛ, 1960. 886 с.

6. Мурти Т.С. Сейсмические морские волны цунами. Ленинград: Гидрометеоиздат, 1981. 448 с.

7. Положий Г.Н. Уравнения математической физики. М.: Высш. шк., 1964. 560 с.

8. Савченко В.Н., Смагин В.П., Фонарев Г.А. Вопросы морской электродинамики. Владивосток: Изд-во ВГУЭС, 1999. 208 с.

9. Сёмкин CB., Смагин B.H, Савченко B.H. Морские электромагнитные поля. Часть III. Инфразвуковой океанический волновод, подводные взрывы. СПб.: Изд-во Санкт-Петербургского Политехн. ун-ва, 2010. 130 с.

10. Сёмкин CB., Смагин B.H, Савченко B.H., Генерация возмущений магнитного поля при подводном взрыве II Известия РАН. Физика атмосферы и океана. 2010. Т. 46, № 1. С. 138-141.

11. У Тон Иль, Табояков А.А., Храмушин B.H. Методические разработки прогноза сильных землетрясений гидроакустическими методами (использование сейсмогидроаку-стической локации местоположения очагов землетрясений) II Морские исследования и технологии изучения природы Мирового океана. Bbm. 1. Bладивосток: ДBO РАН, 2005. 192 с.

транслитерация

1. Batygin V.V., Toptygin I.N. Sovremennaya elektrodinamika, part 1, Makroskopicheskaya teoriya, M., 2003, 736 p.

2. Isakovich M.A. Obshchaya akustika, M.: «Nauka», 1973, 496 p.

3. Kontorovich V.M. O magnitogidrodinamicheskom effekte v okeane, Dokl. ANSSSR. Ser. geofiz, 1961, No 3, pp. 576-579.

4. Kontorovich V.M. Magnitnaya gidrodinamika okeana, Voprosy magnitnoi gidrodinamiki i dinamiki plazmy, Riga: Izd-vo AN Latv. SSR, 1962, pp. 171-177.

5. Mors F., Feshbakh G. Metody teoreticheskoi fiziki, vol. 2, M.: IIL, 1960, 886 p.

6. Murti T.S. Seismicheskie morskie volny tsunami, Leningrad: «Gidrometeoizdat», 1981, 448 p.

7. Polozhii G.N. Uravneniya matematicheskoi fiziki, M., «Vysshaya shkola», 1964, 560 p.

8. Savchenko V.N., Smagin V.P., Fonarev G.A. Voprosy morskoi elektrodinamiki, Vladivostok, VGUES, 1999, 208 p.

9. Semkin S.V., Smagin V.P., Savchenko V.N.. Morskie elektromagnitnye polya, Chast' III, Infrazvukovoi okeanicheskii volnovod, podvodnye vzryvy. SPb, Izd-vo SPb. Politekhn. un-va, 2010, 130 p.

10. Semkin S.V., Smagin V.P., Savchenko V.N., Generatsiya vozmushchenii magnitnogo polya pri podvodnom vzryve, Izvestiya RAN. Fizika atmosfery i okeana, 2010, vol. 46, No 1, pp. 138-141.

11. U Ton Il', Taboyakov A.A., Khramushin V.N. Metodicheskie razrabotki prognoza sil'nykh zemletryasenii gidroakusticheskimi metodami (ispol'zovanie seismogidroakusticheskoi lokatsii mestopolozheniya ochagov zemletryasenii), Morskie issledovaniya i tekhnologii izucheniya prirody Mirovogo okeana, iss. 1, Vladivostok: DVO RAN, 2005, 192 p.

© В.П. Смагин, 2017 © С.В. Сёмкин, 2017 © В.Н. Савченко, 2017

Для цитирования: Смагин В.П., Сёмкин С.В., Савченко В.Н. Магнитные индукционные эффекты, сопровождающие подводные землетрясения // Территория новых возможностей. Вестник Владивостокского государственного университета экономики и сервиса. 2017. Т. 9. N° 4. С. 169-187. For citation: Smagin V.P., Semkin S.V., Savchenko V.N. Magnetic induction effects caused by an underwater earthquake, The Territory of New Opportunities. The Herald of Vladivostok State University of Economics and Service, 2017, Vol. 9, No 4, pp. 169-187. DOI dx.doi.org/10.24866/VVSU/2073-3984/2017-4/169-187 Дата поступления: 11.09.2017.

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