Научная статья на тему 'Сравнительный анализ результатов решения задачи о диффузионно-тепловой неустойчивости пламени'

Сравнительный анализ результатов решения задачи о диффузионно-тепловой неустойчивости пламени Текст научной статьи по специальности «Физика»

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

Аннотация научной статьи по физике, автор научной работы — Шрагер Э. Р., Васенин И. М., Сабденов К. О.

Проводится сравнение областей неустойчивого горения, полученные на основе применения различных модельных функции для скорости химической реакции. Показано сильное влияние величины энергии активации химической реакции на прогноз устойчивости горения при малых числах Льюиса Le. Дано теоретическое объяснение устойчивости горения при Le. 0.

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

Текст научной работы на тему «Сравнительный анализ результатов решения задачи о диффузионно-тепловой неустойчивости пламени»

превращение N^00,, П^КН4СЮ41 на самом деле происходит в два этапа или ромбически-кубическому превращению при 513 К предшествует структурный переход при 475 К в нестабильную фазу. В интервале 475...512 К наблюдается осцилляция акустических параметров, а в области постперехода осциллирует ТКЛР. Такое поведение указанных параметров, вероятно, обусловлено одновременным протеканием полиморфного превращения, химической реакции и

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

СПИСОК ЛИТЕРАТУРЫ

1. Парсонидж Н., Стейвли Л. Беспорядок в кристаллах. В 2-х ч. / Пер. с англ. под ред. Г.Н. Жижина. - Ч. 1. - М.: Мир, 1982. -434 с.

2. Мейер К. Физико-химическая кристаллография. - М.: Металлургия, 1972. - 480 с.

3. Беломестных В.Н., Теслева Е.П. Полиморфные превращения типа ориентационный порядок - беспорядок. Ч. I. Общие представления // Известия Томского политехнического университета. - 2003. - Т. 306. - № 4. - С. 17-22.

4. Беломестных В.Н., Теслева Е.П. Полиморфные превращения типа ориентационный порядок - беспорядок. Ч. II. Азотсодержащие ионно-молекулярные кристаллы натрия // Известия Томского политехнического университета. - 2004. - Т. 307. -№ 6. - С. 11-17.

5. Беломестных В.Н. Физико-химическая акустика кристаллов. -Томск: Изд-во ТРОЦа, 1998. - 183 с.

6. Беломестных В.Н., Гринева М.Н., Поздеева Э.В. Дилатометрические исследования перхлората аммония // Известия вузов. Физика. - 1981. - № 5. - С. 128. (Деп. в ВИНИТИ № 1784-81).

7. Болдырев В.В., Ворсина И.А., Гришакова Т.Е., Михайлов Ю.И. КР - спектроскопическое наблюдение первичной стадии термического разложения персульфатов // Доклады АН СССР. -1989. - Т. 306. - № 3. - С. 647-649.

8. Brown R.J.C., Weir R.D., Westrum E.F., (Jr). The thermodynamics is of heat capacities and related data of NH4CIO4 and ND4ClO4 // J. Chem. Phys. - 1989. - V. 91. - № 1. - P. 399-407.

9. Solymosi F. Structure and stability of salts of halogen oxyacids in the solid phase. - Budapest: Akad. Kiado, 1977. - 467 p.

10. Вишневский И.И., Скрипак В.Н. Поведение решеточной теплопроводности в области фазовых переходов // В сб.: Теплофизиче-ские свойства твердых веществ. - М.: Наука, 1973. - С. 44-48.

11. Гладков С.О. Теплопроводность пористых веществ // Перспективные материалы. - 2000. - № 5. - С. 9-12.

12. Иванов Е.Ю., Болдырев В.В. О механизме полиморфного перехода в кристаллах перхлората аммония // Доклады АН СССР. - 1979. - Т. 248. - № 4. - С. 862-865.

13. Stremme K.O. The crystal structures of the orientationally disordered cubic high-temperature phases of univalent metal perchlorates // Acta Chem. Scand. - 1974. - V. A28. - P. 515-527.

14. Беломестных В.Н., Похолков Ю.П., Ульянов В.Л., Хас-анов О.Л. Упругие и акустические свойства ионных, керамических диэлектриков и высокотемпературных сверхпроводников. - Томск: STT, 2001. - 226 с.

15. Hamada A., Yamamoto S., Fujiyoshi O. Crystal structure of NH4ClO4 in low temperature phase // J. Korean Phys. Soc. - 1998. - V. 32. - P. S152-S155.

16. Справочник химика. Изд. 3-е, испр. Т. 1, 2. - Л.: Химия, 1971.

УДК 536.46

СРАВНИТЕЛЬНЫЙ АНАЛИЗ РЕЗУЛЬТАТОВ РЕШЕНИЯ ЗАДАЧИ О ДИФФУЗИОННО-ТЕПЛОВОЙ НЕУСТОЙЧИВОСТИ ПЛАМЕНИ

Э.Р. Шрагер, И.М. Васенин, К.О. Сабденов*

Томский государственный университет *Томский политехнический университет E-mail: sabdenov@k21.phtd.tpu.ru

Проводится сравнение областей неустойчивого горения, полученные на основе применения различных модельных функции для скорости химической реакции. Показано сильное влияние величины энергии активации химической реакции на прогноз устойчивости горения при малых числах Льюиса Le. Дано теоретическое объяснение устойчивости горения при LexQ.

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

Различают два вида неустойчивости: диффузионно-тепловую и гидродинамическую. Первый из них возникает при нарушении баланса между тремя факторами: отводом тепла из зоны химической реакции в свежую смесь из горючего и окислителя, диффузией окислителя в зону химической реакции и скоростью ее протекания. Первые два из них характеризуются числом Льюиса Ье=Б/к, где Б и к -

коэффициенты диффузии и температуропроводности, третье же - параметром к=Е(ТЬ-Т0)/(2ЛТЬ2), где Е- энергия активации; ТЬ, Т0 - температура стационарного пламени и начальная температура смеси; Я - универсальная газовая постоянная.

При линейном анализе устойчивости при прочих равных условиях наиболее достоверная информация извлекается в предположении не одномерности по пространству возможных возмущений, накладываемых на изначально плоский стационарный фронт пламени. Поэтому исследования [1, 2] показали, что область устойчивости в плоскости (к ', Ье) сосредоточена в окрестности прямой Ье=1. Эта область имеет форму полуострова, выходящей из полуплоскости к <2 и неограниченно сужающейся в пределе к'-^ж. Но в [1] анализ устойчивости проводился на основе дисперсионного соотношения, полученного в [3], где зона химической реакции представляется границей разрыва. По одну сторону этой границы находится свежая смесь, по другую - продукты горения. В работе же [2] скорость реакции Ж представлялась 5-функцией Дирака, аргументом которой является температура Т. Эти два подхода соответствуют, в строгом смысле, бесконечно большой энергии активации. Несмотря на кажущуюся, на первый взгляд, эквивалентность в постановках задач [1, 2], дисперсионные соотношения получаются разными по форме. Это вызвано тем, что во второй работе предположение о квазистационарном протекании химической реакции отсутствует.

Модельная скорость химической реакции

с конечной энергией активации

Упомянутые выше модельные представления скорости химической реакции Ж значительно облегчают решение задач теории горения. Но результаты, полученные в рамках таких представлении, имеют ограниченную область применения. В этом отношений одной из самых привлекательных модельных форм является [4]

Ж = к0Ы ехр I -

Е

ЯТ„

•П(Т-Т.),

(1)

описывающая брутто-реакцию первого порядка и где к0 - предэкспоненциальный множитель; N -концентрация реагента; п - единичная функция Хевисайда, а Т, - температура воспламенения, явный вид которой можно определить, требуя совпадения выражения для скорости стационарного пламени асимптотически в пределе Т,^ТЬ с известным результатом из теории, где скорость химической реакции определяется законом Аррениуса.

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

В дальнейших аналитических исследованиях устойчивости пламени [5, 6] учитывалась конечность скорости химической реакции на основе выражения (1). Поэтому для выяснения характера поведения области устойчивого горения в зависимости от вида Ж представляет интерес детальное сравнение дисперсионных соотношений, полученных в работах [1, 2, 5, 6].

Дисперсионные соотношения

В простом случае брутто-реакции первого порядка для линейного анализа устойчивости достаточно рассматривать уравнения [2, 5, 6]

дЯ + * М.

д? дх

д 2Я

я' + * ^ дх2 + ду2 + Ж

дЬ+* дЬ = Ье № +д2ь 1+Ж„

д? дх |дх2 ду2 1 *

(2)

где безразмерное время ? измеряется в единицах к/(и„0)2, а координаты х, у - в единицах к/(и„°), где и0 - скорость плоского стационарного пламени, распространяющегося вдоль координаты х. Связь между Т, N Ж и введенными вместо них новыми величинами и, Ь, Ж, дается формулами

N - N и

Т-Т

Я = о

Т -Т

QN о

(Т-То), Ь = ■

N

Ж, = а(1 - Ь)п(и - и,), а =

коК

(иО)

ехр

- * = — ' " О '

Е

ЯТ„

Здесь полагается температура пламени равной адиабатической, поэтому TЬ=T0+QN0/cp, где Q - тепловой эффект реакции; N - начальная концентрация реагирующего вещества; ср - теплоемкость при постоянном давлении. Граничные условия, накладываемые на решения (2), отражают неизменность исходного состояния смеси и наличие термодинамического равновесного состояния продуктов горения: х ^ -ж: и = Ь = 0, х ^ +оо : йи / йх = йЬ / йх = О,

Исследование устойчивости сводится к исследованию поведения спонтанно возникающих в системе малых (в данном случае двумерных) возмущений. Если Я°(х), Ь°(х) - стационарные решения (2), то нестационарные распределения температуры и выгорания представляются как [7]

Я =Я0 +4(?,у)ад, Ь = Ь0 ,у)О(х), * = 1+ *',

где 4, 4 - деформации полей температуры Я и выгорания Ь; Е, О - неизвестные функции; *'- изменение скорости пламени. Если 4=4(1+е), е=еош1, то дальнейший анализ устойчивости пламени сводится, как и в [1-3], к постановке задачи в предположении постоянства скорости пламени, а для * получается формула зависимости скорости пламени от кривизны его фронта с неопределенным постоянным параметром е [7]:

* =-

(Ье -1)(1 + е) д24

е ду2 '

П

Дисперсионное соотношение, полученное с ис пользованием (1) имеет вид [5, 6]:

к - 2у х-У к+1)

к +1 ~7+Р Я2)]

- = 0,

(3)

_ -1 + + 4(Р + Я2) 1+,] 1+ 4Ье(Р+ ЬеЯ2) Г_ 2 , в_ 2Ье

_ -1 + -уД + 4Ье(д + 0 + ЬеЯ2) _ ^Д + 4gLe -1 Х_ 2Ье , к _ 2Ье ,

где О - инкремент нарастания возмущений с волновым числом Я.

Аналогичные выражения, полученные в [2, 3] и в принятых здесь обозначениях имеют вид:

к ,= (2«-1)(!*в-1) Le(в-а) + Le -1'

- (ГИ. Баренблатт, Я.Б. Зельдович, А.Г Истратов); (4)

к ,_а(Ъев-1) + Le(а-1)2 Le(в -а) + Le-1

- (А.П. Алдушин, С.Г. Каспарян);

1 + 1 + 4(О + Я2)

(5)

а _ -

2

ших длин волн). Потеря устойчивости с (1), как и в [1-3], имеет апериодический характер. Т.е. на границе устойчивости 0=0 и вблизи этого значения является вещественным числом.

Если учет конечности скорости химической реакции по модельной функции (1) при Le>1 не приводит к качественно новому поведению пламени, то этого нельзя сказать в случае Le<1. Здесь инкремент нарастания возмущений О на границе устойчивого горенияявляется чисто мнимой величиной: П=гю, где а - частота. Отсутствие действи-

тельной части О следует как из анализа корней (4), (5) [1, 2], так и из (3) [5, 6].

Входящие в (3) и (4, 5) параметры к', к хотя и различны, но мало отличаются по численной величине (численно к определяется через д и Le, причем вид параметра а задается выбором выражения для скорости пламени и0, которое лучше взять из [4]). Поэтому в дальнейшем считается к=к и сравнение результатов, следующих из соотношений (3-5) проводится с этим допущением.

Анализ результатов

Рассмотрим числа Le>1. Если Le отличается от единицы, то пламя оказывается неустойчивым по отношению к возмущениям с малыми волновыми числами. При этом формулы (4) и (5) дают практически одинаковые результаты (рис. 1) даже при небольших энергиях активации, исключение составляют большие но, скорее всего, нереальные значения Le. Из формулы (3) следует более широкая область неустойчивости, чем получается из (4, 5). Но с ростом значений к области неустойчивости, рассчитанные по приведенным трем дисперсионным соотношениям, становятся мало различающимися.

Учет нестационарности скорости сгорания вещества [2] в зоне химической реакции привел хотя и незначительному, но заметному увеличению области неустойчивости (кр. 2, 3). Но этот эффект оказался сильнее при более близкой к реальной форме скорости химической реакции (1): если при Я2=0 (к=5) поправка к формулам (4, 5) составляет 6,7 %, то для волновых чисел с Я2=0,2 она примерно равна 15 %.

Из результатов, приведенных на рис. 1, также следует, что при неизменном значении к рост числа Льюиса приводит к дестабилизации пламени, причем со стороны малых волновых чисел (или боль-

0,0 0,1 0,2 0,3 0,4 0,5 2 0,6

Я

а)

2,0 Ье 1,8

1,6

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

1,4

1,2

1,0

Рис. 1. Границы области неустойчивости пламени при Le>1: а) к=5 и б) к=10. Кр. 3 - формула (4), 2 - (5), 1 - (3). Область неустойчивости находится выше кр. 1-3

Согласно дисперсионным соотношениям (4) и (5) область неустойчивости зарождается на плоскости (Я2, Le) в окрестности точки Я2«0,2, Le=0 (рис. 2, а). Из соотношения (3) следует зарождение области неустойчивого горения с окрестности точки Я2«0,2, Le«0,1 в виде острова с быстро увеличивающимися размерами с ростом к.

По соотношениям (4, 5) при Le=0 устойчивое горение невозможно, всегда существуют возмущения с длинами волн (рис. 2, 3), которые со време-

нем будут нарастать. Дисперсионное соотношение (3) констатирует возможность существования устойчивого пламени при столь малых числах Льюиса, но при относительно небольших энергиях активации. Но с ростом к устойчивость пламени быстро снижается, особенно для малых чисел Le. Так, уже при к=6 нижний берег острова неустойчивости полностью касается оси Я2 (рис. 3, а).

Обращает на себя внимание столь сильное влияние на результаты учет нестационарности скорости химического превращения. Если при малых волновых числах и Le>1 формулы (4, 5) приводили к одинаковым результатам, то при Le<1 этого уже не наблюдается. Расчеты положения границы области устойчивости по формулам (4) и (5) отличаются «10 % для Я2<<1 (рис. 2), но для больших волновых чисел отличие возрастает многократно. С ростом к (или энергии активации) различие в результатах, получаемых из дисперсионных соотношений (3-5) хотя и медленно, но все же стирается. Это быстрее всего происходит в области малых волновых чисел (рис. 3).

Ье

0,16 0,14 0,12 0,10 0,08 0,06

Ч

го

0,0

0,2

0,4

0,6

а)

0,16 0,12 0,08 0,04

0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 ^2 б)

Рис. 2. Границы области неустойчивости пламени при Le<1: а) к=4,7; б) к=5. Кр. 1 - формула (4), 2 - (5), 3 - (3). Область находится ниже кр. 1, 2 и ограничена кр. 3

Изображение области неустойчивости в координатах ^е, а) показывает (рис. 3, б), что колебательные процессы на границе устойчивости могут происходить с двумя частотами (кр. 3). Это не следует из (4, 5): каждому числу Le соответствует только одна частота (кр. 1, 2). Общая тенденция роста частоты с увеличением волнового числа имеет место для всех кривых. Кроме того, имеет место резкий скачок вверх частоты а при очень малых числах Le, что особенно заметно на кр. 3.

Следуя работе [1], область устойчивого горения можно представить на плоскости (к, Le), рис. 4. Здесь подразумевается устойчивость для любых волновых чисел. Применение функции источника (1) показывает более низкую устойчивость пламени для чисел Le>0,15, по сравнению с 5-функциональным представлением. Но при числах Льюиса Le<0,15 горение после неуклонного снижения устойчивости приобретает дополнительную стабильность, в то время как кривая (пунктирная линия), рассчитанная по дисперсионному соотношению (5), показывает продолжающееся монотонное падение устойчивости пламени (нижние ветви на рис. 4).

Ье

0,35

Т-Г^-1-1-1-г

0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0

К

а)

Ье

0,35 -

0,30-

0,250,20- / ! %

0,15- * ♦ ♦

0,10- 1

0 1

3 4

6 7

б)

Рис. 3. Границы области неустойчивости (к=6) в координатах (Le, Я2) (а) и (Le, а) (б). Нумерация кривых и соответствующих им формул та же, что и на рис. 2. Области неустойчивости на рис. б) ограничены кр. 1-3 и осью а

2,0'

Le

1,5

1,0

0,5

0,0

Рис. 4. Границы области устойчивого горения в координатах к, Le. Пунктирные линии - расчет по (5), сплошные - по (3)

Формулу (3), раскладывая в ряд Тейлора по степеням 1Д и удерживая только первые два его члена с учетом связи у=а—1, перепишем в виде

_ (р +а - 1)(Ьев -1) + (Ье-1)0 Ье(в-а) + Ье -1 '

Даже если не обращать внимания на небольшое различие в параметрах и k то, как видно, она по форме не совпадает ни с формулой (4) Баренблатта-Зельдовича-Истратова, ни с формулой (5) Алдуши-на-Каспаряна. Это выглядит удивительным по той причине, что выражение для скорости реакции (1) в пределе очень большой энергии активации мало отличается от 5-функции. И только лишь в предельном случае k=ж все три дисперсионные соотношения приводят к равенству Ье(в-а)+Ье-1=0, один из корней которого 0.=Х указывает в отсутствии строгого выполнения равенства Ье=1 на неустойчивость пламени с бесконечно тонкой зоной химической реакции.

Проверка на адекватность полученных результатов

Химическая реакция в реальном пламени протекает сложным образом со многими путями (цепями) и образованием большого числа промежуточных компонентов. Представление же горения одностадийным процессом по схеме А^В является упрощением, где сложное пламя рассматривается как имеющее какое-то одно «топливо», с исчезновением которого связано выделение тепла [8]. Тогда коэффици-

енты диффузии и теплопроводности, предэкспонен-циальный множитель и энергия активации в скорости химической реакции в выражениях (1, 2) и др. имеют смысл т.н. формально-кинетических параметров, имеющих смысл усредненных или эффективных величин по всем возможным цепям и ступеням химической реакции [9]. Согласно приведенным выше результатам, эффективное число Льюиса в устойчивом пламени в газах должен быть приблизительно равным единице. Для смесей СН4-02, С0-02 оно равно 1,25 и 1,1 соответственно, и приведенные примеры близости эффективного числа Льюиса единице являются не единичными [8].

Равенство Ье=0 имеет место при безгазовом горении так называемых систем самораспространяющегося синтеза. Устойчивость пламени при таких числах Льюиса - довольно частое явление, хотя известны многочисленные факты наблюдения и неустойчивых режимов горения [10]. Детальное сравнение теоретических выводов с экспериментальными результатами не представляется возможным ввиду отсутствия надежных систематических данных по эффективной энергии активации при самораспространяющемся синтезе.

Выводы

Сравнительный анализ результатов устойчивости пламени, следующих из различных подходов к представлению характера и скорости протекания химической реакции показал, что:

• замена аррениусовской формы скорости реакции на модельные выражения сильнее всего оказывает влияние на условия устойчивого горения при числах Ье<1, чем при Ье>1;

• учет конечности скорости химической реакции не отражается на условиях реализации апериодической и колебательной потери устойчивости;

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

• учет конечности скорости химической реакции объясняет наблюдающуюся устойчивость плоского фронта пламени при числах Ье«0. Авторы выражают благодарность Р. С. Буркиной и А.М.Ти-

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

СПИСОК ЛИТЕРАТУРЫ

1. Гришин А.М., Берцун В.Н., Агранат В.М. Исследование диффузионно-тепловой неустойчивости ламинарных пламен // Доклады АН СССР. - 1977. - Т. 235. - № 3. - С. 550-553.

2. Алдушин А.П., Каспарян С.Г. О теплодиффузионной неустойчивости фронта горения // Доклады АН СССР. - 1979. - Т. 244. - № 1. - С. 67-70.

3. Баренблатт Г.И., Зельдович Я.Б., Истратов А.Г. К теории те-плодиффузионной неустойчивости ламинарного пламени // Прикладная механика и техническая физика. - 1962. - Т. 17. -№ 3. - С. 21-26.

4. Вилюнов В.Н., Дик И.Г., Зурер А.В., Ищенко А.Н. Зависимость скорости распространения теплодиффузионного пламени для широкого диапазона чисел Le // Физика горения и взрыва. - 1984. - № 5. - С. 35-42.

5. Сабденов К.О. О диффузионно-тепловой неустойчивости ламинарного пламени // Инженерно-физический журнал. -2002. - Т. 75. - № 4. - C. 73-79.

6. Сабденов К.О. Нахождение условий диффузионно-тепловой устойчивости пламени применением модельной функции скорости химической реакции // Известия Томского политехнического университета. - 2004. - Т. 307. - № 2. - С. 25-29.

7. Сабденов К.О. К вопросу нахождения постоянной Маркштей-на // Известия Томского политехнического университета. -2004. - Т. 307. - № 3. - С. 21-25.

8. Фристром Р.М., Вестенберг А.А. Структура пламени. Пер. с англ. - М.: Металлургия, 1969. - 364 с.

9. Зельдович Я.Б., Баренблатт Г.И., Либрович В.Б., Махвила-дзе Г.М. Математическая теория горения и взрыва. - М.: Наука, 1980. - 480 с.

10. Мержанов А.Г. Проблемы технологического горения // Процессы горения в химической технологии и металлургии. - Черноголовка: Изд-во ИХФ РАН, 1975. - С. 5-28.

УДК 539.3;539.215

ИЗУЧЕНИЕ ОСОБЕННОСТЕЙ ОТКЛИКА ГРАНИЦ РАЗДЕЛА В РАЗЛОМНО-БЛОКОВЫХ СРЕДАХ НА ИЗМЕНЕНИЕ ИХ СОСТОЯНИЯ И ДИНАМИЧЕСКИЕ ВОЗДЕЙСТВИЯ. Ч. 2. РЕЗУЛЬТАТЫ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ

С.В. Астафуров1, Е.В. Шилько1, А.В. Димаки1, В.В. Ружич2, В.В. Лопатин3, Б.Г. Шубин3, В.Л. Попов4, С.Г. Псахье1

'Институт физики прочности и материаловедения СО РАН. г. Томск 2Институт земной коры СО РАН. г. Иркутск 3Научно-исследовательский институт высоких напряжений при ТПУ. г. Томск 4Берлинский технический университет. г. Берлин, Германия E-mail: lopatin@hvd.tpu.ru

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

Введение

Результаты многолетнего мониторинга смещений по активным разломам и натурных экспериментов демонстрируют принципиальную возможность инициации в высоконапряженных фрагментах активных разломов аномальных по амплитуде смещений в «квазивязком» (криповом) режиме совместным применением водонасыщения и вибрационных воздействий [1, 2]. При этом в процессе трения краев разлома будет иметь место диссипация аккумулированной в земной коре упругой энергии и релаксация локальных напряжений. Полученные результаты позволили предложить идею принципиально нового способа управления режимом смещений во фрагментах сейсмоопасных тектонических разломов путем инициирования смещений в режиме сдвиговой ползучести за счет насыщения жидкостью фрагмента разлома в сочетаниями с циклическими воздействиями на него. При этом необходимо иметь ясное понимание механизмов влияния подобных комплексных воздействий, приводящих к переключению режима смещений берегов разлома. Отметим, что для детального изучения особенностей отклика зон активных разломов в настоящее время активно используется математический аппарат механики деформируемого твердого тела, а также методы компьютерного моделирования.

В течение многих лет при описании деформационных процессов в земной коре в большинстве случаев применялись различные способы осредне-

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

Многочисленные литературные данные говорят о том, что земная кора имеет иерархически организованное дискретное блочное строение [3-6], что обеспечивает наличие дополнительных степеней свободы и, как следствие, достаточно высокую деформационную способность среды в целом. В восьмидесятые годы прошлого века, после того как были опубликованы работы академика М.А. Садовского, началось интенсивное развитие исследований, направленных на изучение влияния блокового строения земной коры на различные геофизические процессы. Постепенно в науках о Земле сформировались два различных подхода к проблеме влияния нарушений сплошности на свойства мас-

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