Файл: Подводные и подземные взрывы сб. ст.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 111

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

СЕЙСМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПОДЗЕМНЫХ ЯДЕРНЫХ ВЗРЫВОВ.

Ч. I. РАСЧЕТ СЕЙСМИЧЕСКОГО СПЕКТРА1)

Р. А. Мюллер, Дж. Р. Мёрфи

Описана модель масштабного пересчета сейсмического спектра, основанная на аналитической аппроксимации функции источника ядерного взрыва, п по большому объему экспериментальных дан­ ных произведена оценка ее применимости. Результаты теории подо­ бия указывают, что решающую роль в формировании сейсмических спектров подземных ядерных взрывов играют мощность, глубина и вмещающая среда. Представленные результаты показывают, что для функции возбуждения сейсмического источника подземного ядер­ ного взрыва не выполняется подобие по корню кубическому из за­ ряда. Для взрывов на достаточно больших глубинах, обеспечиваю­ щих полный камуфлет, теория подобия предсказывает зависимость показателя степени при заряде от частоты в диапазоне значений от 0,90 на низких частотах до 0,45 на высоких частотах. Эти значения показателя степени при заряде находятся в хорошем соответствии с эмпирическими соотношениями, статистически устанавливаемыми по данным измерений.

ВВЕДЕНИЕ

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

’)

Mueller R. A., Murphy J. R., Seismic characteristics of under­

ground nuclear detonations. Part

I. Seismic spectrum scaling, Bulle­

tin of

the Seismological Society

of America, 61, № 6, 1675—1692

(1971).

 

СЕЙСМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПОДЗЕМНЫХ ВЗРЫВОВ. Ч. I 289

числительных программ, основанных на фундаменталь­ ных физических принципах [1, 3]. Несмотря на то что был достигнут значительный прогресс в этой области, многое еще остается сделать в части отыскания аде­ кватных уравнений состояния и условий хрупкого раз­ рушения для разнообразных горных материалов, в ко­ торых производятся взрывы.

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

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

10 Зак, 741


290

Р. А. МЮЛЛЕР, ДЖ. Р. МЕРФИ

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

ТЕОРИЯ

Решение сферически симметричной задачи для функ­ ции возбуждения, действующей в безграничной одно­ родной среде, было получено Шарпом [14] во временном представлении и Латтером и др. [6] в частотном выра­ жении. Решение в частотной форме легко получается из волнового уравнения и закона Гука. Предполагается, что существует окружающая точку взрыва сферическая поверхность, вне которой применима теория бесконечно малых деформаций. Радиус этой сферы будем называть упругим радиусом и обозначим его через ге1, тогда для радиальных расстояний г ^ ге\ уравнение движения приводится к одномерному волновому уравнению для потенциала смещения ср:

д ! <р

_

1

д 2(р

m

д г 2

 

с2

д /2 ’

которое имеет решение <р = ср(т) для выходящих сфери­ ческих воли, где x = i rei)/c — время запазды­ вания. Смещение Z(r,t) связано с потенциалом смеще­ ния

 

 

Z (r,

=

 

 

(2)

Закон Гука устанавливает связь между давлением

p{t) и смещением на геХ\

 

 

 

 

 

 

dZ (г, 0

L

-

2% iz (г. О

( )

P{t) = —

(А + 2ц)

дг

 

3

 

 

 

 

где

X и

р — константы Ламе. Подставляя

сначала

Z(r,i) из уравнения (2) в (3)

и применяя преобразова­

ние

Фурье,

заменяем спектр

ср(т) спектром

Z(r,i) и



СЕЙСМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПОДЗЕМНЫХ ВЗРЫВОВ. Ч. I 291

получаем

Р(to) re1

1П " .j"™(СО\|

 

 

с2

 

 

Z(co) =

/

A

(

(4)

 

4(Д.

Г

ГС

/

[0q + /со0и — P<B“J

 

 

где Z (со) — спектр

смещения

на

 

г;

р(со)— спектр

функ­

ции давления,

действующей

на

ге\\

ц — модуль

сдвига;

с — скорость

волны сжатия;

 

 

 

+

2,u)/4j.i; со — уг­

ловая частота и соо =

с/ге\.

 

 

поверхностной

 

точке

Спектр Фурье

на

удаленной

 

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

Z(co)Dist = r(co)-Z(co)r.

(5)

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

2 ,

(со)

 

П СО)

2 , ( 0 )

Z 2

(ш )

D Is t

Г (со)

Z-2 (ш )

 

 

 

 

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

2, (® )

Pi (“)

г е \ №

\ X

 

2 2(ш) D is t

р 2 ( ш )

г еиУ-\с2

 

х /

с?+

г20 2

/со“ 02 + 0

— 2Ра) ®02®2 + Р2“ 4

с2 +

гчо?

\

0

(7)

 

 

 

 

'01 + 0

— 2Pi) “oi“2 + Р2“4

Аналитическая аппроксимация формы профиля дав­ ления, действующего на упругом радиусе, может быть подобрана на основе анализа наблюдений внутри среды при ряде подземных ядерных взрывов, включающих «Салмон» [12], «Хардхэт» [11] и «Шоал» [15]. Эта

10*


292

Р. А. МЮЛЛЕР, ДЖ. Р. МЕРФИ

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

p(t) = (p0e-«* + p 0c)H(t),

(8)

где р0 + Рос = Рои— пиковое ударное давление, р0с — установившееся давление, а — постоянная затухания и

Рис. 1. Профиль давления на границе упругой зоны.

Н(t) — единичная ступенчатая функция. Спектр p(t) вы­ ражается в виде

р{ш) =

( P a s Р о с )

Р о с

(9)

а + гео

гео

Из уравнений (7) и (9) получаем окончательное со­ отношение подобия между амплитудными спектрами двух различных взрывов:

г , (со)

'~е1,<ЧЙ2 - /

(«2 +

М2)

, [

(g>2P0s, +

PQC,a l) Х

^2 (со) iD Is t

те \Н ^ \ V

(а | +

со2)

У

(со2Pos, +

Рос,а1)

 

 

„2 . _2..2

Р2СО4 +

(l — 2Рз) COiX +

C0g2

 

X

Cj “Г Г со

(10)

С22 +, г9®2

Pi®1 +

(l — 2Pi) ®oi®2 +

 

 

“ 01