CC BY
14
РОССИЙСКИЙ ЖУРНАЛ НАУК О ЗЕМЛЕ, ТОМ 1, N 3, Февраль 1999
Модель формирования тектономагнитного эффекта в зоне разлома при сдвиге
И. А. Гарагаш
Институт физики Земли РАН
Аннотация. Изучено формирование тектономагнитного предвестника в процессе образования полосы сдвига. Магнитная вариация связана со взаимодействием деформационного и геомагнитного полей. Показано, что эффект существенно зависит от ориентации полосы сдвига относительно вектора магнитного поля Земли. Этот результат может быть использован для анализа тектономагнитных вариаций с целью краткосрочного прогноза землетрясений.
Введение
Появление тектономагнитных предвестников связано с изменением напряженного состояния породного массива при сейсмотектоническом процессе. Обычно считается, что формирование локальной геомагнитной вариации на земной поверхности является следствием пьезомагнитного эффекта намагниченных пород [Основы прогноза..., 1983]. Однако пьезоэффект нечуствителен к вариациям интенсивности напряжений и деформаций сдвига (вторых инвариантов соответствующих тензоров). В то же время именно сдвиговые деформации приводят к изменению формы и сопровождают подвижки по разломам.
Причиной возникновения предвестниковых тектономагнитных аномалий может служить взаимодействие деформационных и электромагнитных полей. Развитие этой новой области механики происходило в последние три десятилетия и стимулировалось проблемами, возникающими в геофизике, некоторых разделах оптики и акустики [Кпоро$” 1955; Оипкт, Еггпдеп, 1963]. В основе эффекта лежит образование электромагнитного поля в результате деформационного перемещения частиц проводящей среды в магнитном статическом поле. Полная система уравнений, описывающая процесс включает в себя систему квазистатических уравнений электродинамики медленно движущейся среды.
Разработка подхода, основанного на связанности магнитного и деформационного полей представляется перспективным. Согласно известным те-
©1998 Российский журнал наук о Земле.
Статья КЛЕ98009.
Онлайновая версия этой статьи опубликована 15 декабря 1998. 1ЖЬ: Ы^р://eos.wdcb.rssi.ru/rjes/RJE98009/RJE98009.htm
ориям подготовки очага землетрясения [Мячкин и др., 1975], сейсмическому событию предшествует процесс формирования магистрального разрыва, сопровождаемый сдвиговыми подвижками. По современным представлениям сейсмическое событие рассматривается как процесс локализации сдвиговой деформации. Такой взгляд основан на экспериментах с образцами горных пород [ Wawersik, Brace, 1971; Wawersik, Fairhurst, 1979], согласно которым деформирование горных пород сопровождается развитием исходных трещин и образованием новых. При этом эффективные упругие свойства постоянно уменьшаются и наступает разупрочнение материала, выражающееся в реализации нисходящей ветви диаграммы деформирования в условиях контролируемой деформации. После достижения пикового напряжения в характере распределения трещин наступают заметные изменения. Наблюдается появление преимущественной ориентации дефектов в результате частичного закрытия одних из них и ускоренного роста других. Происходит явление известное под названием локализации неупругой деформации [Шее, 1976], которое и завершает разрушение образца. Установленные в лабораторных условиях особенности разрушения материалов с большим числом трещин послужили основой для модели лавинно-неустойчивого трещинообразования (JIHT-модель) [Мячкин и др., 1975; Miaehkm et al., 1975], одно из основных положений которой состоит в том, что исходная неоднородная среда может быть приведена к некоторой эквивалентной, деформация которой будет будет равна собственно упругой деформации сплошного материала при заданных напряжениях и дополнительной, вызванной наличием дефектов. При этом образование магистрального разрыва (локализация деформации) связано с лавинным развитием внутренней неустойчивости среды, наступающей при критической плотности трещин.
В статье показано, что процесс локализации сдвиговой деформации создает благоприятные условия для появления вариаций магнитного поля в силу эффекта связанности деформационных и электромагнитных полей.
Вывод основных соотношений
Взаимодействие деформируемого твердого тела с магнитными полями изучалось в геофизике в связи с распространением сейсмических волн. Прежде всего исследователей интересовало как скажутся магнитные эффекты на скоростях сейсмических волн и их затухании. Исследовалось совместное решение уравнений движения магнитного тела \Nowacki, 1976]
+ = (1)
и соотношений
(2), можно найти изменения напряженности магнитного поля Hi.
Реализуем изложенный подход для условий локализации деформации в узких зонах. Согласно выражению (2) изменение поля деформации на малую величину щ должно привести к возмущению магнитного поля на величину /г8-. Поэтому представим величины Нг И и, В ВИДе СуММ
Щ = Н° + й,-, Щ = и? + щ (4)
где , II-1 - параметры основного состояния.
Аналогичным образом разделяются и напряжения,
<ГгЗ=(ТЪ + ТЧ- (5)
Подставляя (4) и (5) в уравнения (1) и (2) и собирая
величины одного порядка малости, получим
ТЧ,3 = 0> (6)
\72іїг' — <тц-
дЩ
т
= —а\і
дщ
Ж
(2)
м
связывающих между собой пространственные и временные изменения компонент вектора напряженности магнитного поля и вектора перемещений щ.
Здесь сгц - компоненты тензора напряжений, ¡1 -магнитная проницаемость, а - проводимость среды.
Уравнения (1) и (2) являются основными уравнениями магнитоупругости в том случае, когда свободные заряды в теле отсутствуют и токами смещения можно пренебречь. В отличие от обычных уравнения движения (1) содержат магнитный тензор Максвелла
^2, дЫ („О дщ о диі .
V +„И Я, Я, =0. (7)
Компоненты напряжений а&1 связаны с деформациями
6ц —
.0/2
(8)
соотношениями
(3)
который и обеспечивает влияние магнитного поля на деформации.
В результате решения системы уравнений (1), (2) был сделан вывод [СЬаёипск, 1956; Кпоро[[, 1955], что в такой постановке геофизических задач магнитные эффекты малы и учет их не представляет интереса. Однако при поиске предвестников землетрясений речь идет об обратной задаче, когда ищутся весьма малые аномальные изменения магнитного поля, вызываемые механическими движениями среды. Это значит, что мы можем отбросить слагаемое (3) и решить обычную несвязанную задачу о деформациях среды. После чего, подставляя найденные выражения для перемещений в правую часть уравнений
— G[(l6ik6jl ~\~ 2^/(1 ^6к1&{^£к1 (9)
где С - модуль сдвига, ь> - коэффициент Пуассона.
Поскольку подвижки, предшествующие основному толчку, как было показано выше носят характер локализации деформации, исследуем тектономагнит-ной эффект именно в связи с этим явлением.
Поле перемещений при образовании полосы сдвига имеет вид [Шее, 1976; Гарагаш, Николаевский, 1989]
и3 = и3(х2), и2 = 0, иі = иі(х2),
(10)
если ось х2 направить препендикулярно полосе, а оси Х\ И Жз поместить в ее плоскости (рис. 1).
При образовании зоны локализации деформации внутри нее происходит растрескивание материала, образование блочной структуры и затем катакла-стическое течение [МгасЬкгп е£ а1., 1975]. Эти процессы приводят к уменьшению эффективного модуля сдвига С в зоне локализации. Для того, чтобы это учесть предствавим модуль сдвига в виде
С — С0/[ 1 + х(2)ш(ж2)].
(П)
и
и
где G° - модуль сдвига ненарушенного материала, и>(х2) - функция неоднородности, зависящая как и перемещения только от одной координаты.
Функция \(¿) характеризует развитие во времени процесса локализации (рис. 2).
Исходя из того, что в процессе локализации поле напряжений не возмущается [Rice, 1976] и плоскость локализации совпадает с площадкой главного сдвига, из (9) с учетом (11) и (10) получим, что отличные от нуля сдвиговые компоненты возмущенного тензора деформации имеют вид
£-12 = «1,2 = (<rl2/G0)\(t)w(x2),
£23 = «3,2 = (сг^з/С0)\(І)ш(ж2).
Перемещением нормальным к зоне локализации и2 в дальнейшем пренебрежем в силу его малости по сравнению с перемещениями и2 и «з в плоскости локализации.
Полагая для определенности, что вектор напряженности магнитного поля Н° лежит в плоскости
^2) Рис. 2. Поверхность изменения безразмерного модуля сдвига G/G0 при формировании полосы сдвига.
(13)
dt <тц х dt
(16)
где о- - переменная Фурье, Л8- и ш - Фурье-образы
функций hi и и> соответственно.
, , Записав решение уравнения (16) и применив к нему
А1!«-) и составляет угол р с осью Ж1 рис. 1 , уравне- _ ^
" обратное преобразование Фурье, получим следую-
ние (7) окончательно преобразуем к виду
щее выражение для hi
где
, dhi д\ . hi,22 - — , г = 1, З
B=(H°a°2/G0)Sm/3.
(14)
(15)
2 /
А н°
L
/ , -V|
Лъ
hi =
exp -
оґі
<7/1
(17)
Анализ тектономагнитного эффекта при формировании магистрального сдвига
Воспользовавшись преобразованием Фурье, уравнение (14) перепишем в виде
ехр(*о-жі) > da
J at ycr/j.
где С1 - постоянная интегрирования.
Определим функцию \ следующим образом
\ = А[1 - ехр(-Аг)], (18)
где А - параметр, характеризующий скорость образования полосы локализации.
Функцию неоднородности примем в виде
j(x2) = ехр(—к2 х2)
(19)
где к - параметр характеризующий полосу локализации.
Производя интегрирование, получим
hi = tjfi.X
BA
к^/тт
Рис. 1. Схема полосы сдвига и ориентация магнит-77°
ного вектора И .
Г Í t 1 ^ •>!
exp 4Р) а\
сг/лХ — а2
cos ax2da —
Рис. 3. Графики изменения приращения напряженности магнитного поля hi от безразмерного времени А^ при различных значениях к (к = 0.1 м-1, а ¡л. = 10“3 с м“2 и В А = 300 нТ).
(20)
'DC
— ехр(—A t) J
cos ax2da
После упрощения уравнения (20) окончательно получим (см. Приложение)
hi = 2ВАку/к j\JМК+К [ sin 2 {л/ккх2)
X 2
кк2(2
exp
(A t + k)
%2
cos (2л/кк^ dÇ—
f Г K,k2Ç2
— COS [2у/ккХ2) І exp
(A t к
oo
X sin (2\/~K,kÇ) dcj — exp(—At) x
•f2
x smfôï/Kkxn) J exp(—k2(2)cos('2\/Kk()d(—
OO
X2
— cos (2yfJîjx2) J exp(—i;2^2) sin(2\/KkQd(
OO
tjfi.X
где К =
4 к2
- безразмерный параметр.
X
где Erf(x) = —= exp(—t2)dt - функция ошибки.
v71" J
о
Магнитные свойства горных пород могут колебаться в широких пределах. Так произведение cr/j, меняется от 10-2 до 10-23 сек/м2 [Rjevskiy, Novik, 1978]. Но учитывая, что напряженность магнитного поля Земли составляет около 50000 нтл и оценивая <т%
отношение как 10-2 — 10_3 при А = 1, получаем, что В А = 1 -f- 500 нтл. Заметим, что при А = 1 жесткость земной коры на оси полосы сдвига уменьшается в два раза по сравнению с начальными значениями. Это уменьшение зависит от концентрации повреждений и микротрещин. Согласно теоретическим [Салганик, 1973] и экспериментальным исследованиям [Paterson, 1978], параметр А может быть и больше единицы. Полагая, что полоса локализации при подготовке землетрясения может составлять от 1 до 103 м, найдем, что параметр к меняется в пределах от 1 до 10“4. Построим графики изменения hi от безразмерного времени A t при значениях к = 0.01, 0.005, 0.001 и к = 0.1 т-1, и// = 10_3 sek/m2, В А = 300 нтл (рис. 3). Как видно, тектономагнит-ный эффект довольно быстро достигает максимума в интервале 0.5 < Xt < 2 и затем медленно затухает. При выбранных реальных параметрах его длительность может составлять десятки минут, причем при увеличении ширины полосы сдвига это время возрастает. На рис. 4 построена поверхность изменения возмущения напряженности магнитного поля hi с течением времени t в направлении хп при к = 0.1 т-1, а ¡л. = 10_3 sek/m2, А = 0.04 sek-1, к = 0.001 и В А = 300 нтл.
(21)
Наибольший эффект должен достигаться на оси зоны локализации хо = 0. В этом случае формула (20) приобретает вид
Ы =
= —¡.\/тгВА\/кехр[—(Х1 + к)] |Ег/ ^*УАТ^~к^ — (22)
Рис. 4. Поверхность возмущения напряженности магнитного поля hi в окрестности полосы сдвига
при к = 0.1 м 1, (Tyti = 10 3 к = 0.001 и ВА = 300 нТ.
А = 0.04 с"
Заключение
Приложение А
Таким образом выполненный анализ позволяет сделать вывод, что процессы локализации деформации должны сопровождаться заметными вариациями вектора напряженности магнитного поля. При этом эффект существенно зависит от ориентации полосы сдвига по отношениям к вектору магнитного поля Земли и по разному проявляется на разных компонентах. Эти особенности должны учитываться при размещении магнитометров в сейсмоактивных районах. Для получения максимального сигнала надо стремиться ставить их там, где зоны разломов ориентированы под прямым углом к вектору магнитного поля Земли.
Полученное решение может быть использовано для целей прогноза землетрясений. Действительно, тот факт, что согласно уравнению (15) вектор возмущения напряженности магнитного поля /г8- лежит в плоскости локализации деформации х\ х3 и совпадает с направлением максимального сдвига, позволяет интерпретировать магнитные наблюдения в терминах теории подготовки землетрясения. В спокойные периоды ориентация вектора возмущения в пространстве должна меняться произвольно и соответствовать понятию белый шум. Однако по мере локализации процесса растрескивания и формирования магистрального разрыва следует ожидать появления предпочтительных направлений сдвига и, как следствие, упорядочения ориентации магнитного вектора.
Очевидно, что полученное решение позволяет определить направление подвижки, но не ее плоскость, поскольку фиксируемый в пространстве магнитный вектор принадлежит множеству площадок, для которых он является осью вращения. Чтобы выбрать из этих площадок единственную, надо привлечь дополнительные сведения об условиях, в которых происходит локализация. Как отмечалось выше, плоскость локализации ориентирована в направлении близком к плоскости максимальных сдвигающих напряжений. Предположим, что одна из компонент вектора максимального сдвига (в нашем случае магнитного вектора) будет совпадать с одним из главных сдвигов. Тогда возникает возможность однозначно определить ориентацию плоскости локализации в пространстве.
Благодарность. Это исследование поддержано Российским фондом фундаментальных исследований грант 98-05-65063.
Для вычисление интеграла
СО
г ехр( —7а2)
cos ax2da
(Al)
запишем дифференциальное уравнение
1"{х2) + V21 = J ехр(-уа2) соБ(ах2)(1а = о
(А2)
1^1 ( х1
Используя для решения уравнения (А1) метод вариаций постоянных интегрирования, искомый интеграл представим в виде
1{х2) =
__ ( *2
= ¿vf { sin(^2)/ехР ) cosKK-
(АЗ)
— cos(iae2) I exp ( —— ) sin(i/£)d£ ^ .
Литература
Гарагаш И. А., Николаевский В. Н., Неассоциированные законы течения и локализация пластической деформации, Успехи механики, 12, в. 1, 131-183, 1989.
Мячкин В. И., Костров Б. В., Соболев Г. А., Шами-на О. Г., Основы физики очага и предвестники землетрясений, Физика очага землетрясения, Наука, Москва, 6-29, 1975.
Основы прогноза землетрясений (современное состояние, проблемы и перспективы прогноза землетрясений), Прогноз землетрясений, (3), 220 с., 1983.
Салганик Р. Л., Механика тел с большим числом трещин, Механика твердого тела, Известия АН СССР, 4, 149158, 1973.
Chadwick, P., Elastic wave propogation in magnetic field, Ninth Int. Congr. Appl. Mech., Brussels, 143-153, 1956.
Dunkin, J. W., Eringen, A. C., On the propagation of waves in an electromagnetic elastic solid, J. Engng. Sci., 1, 461476, 1963.
Knopoff, L., The interaction between elastic wave motions and a magnetic field in electrical conductors, J. Geoph. Res., 60, 441-457, 1955.
Miachkin, V. I., Brace, W. F., Sobolev G. A., Dieterich J.H. Two models for earthquake forerunners, PAGEOPH, 113, 169-181, 1975.
Nowacki, W., Coupled fields in mechanics of solids, Theoretical and applied mechanics, Proceedings of the 14th IUTAM Congress, Delft, 396-416, 1976.
x
Paterson, M. S., Experimental rock deformation - the brittle field, Springer-Verlag, 254 p., 1978.
Rice, J. R., The localization of plastic deformation, Theoretical and applied mechanics, Proceedings of the 14th IUTAM Congress, Delft, 439-471, 1976.
Rjevskiy, V. V., Novik, G. I., Physical principles of rocks, Nedra, Moscow, 392 p., 1978.
Wawersik W. R., Brace W. P., Post failure behavior of a granite and diabase, Rock Mechanics, 3, 61-85, 1971. Wawersik W. R., Fairhurst C., A study of brittle rock fracture in a laboratory compression experiments, Int. I. Rock Mech. Mm. Set., 7, 561-575, 1979.
(Поступила в редакцию 15 ноября 1998.)
