Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Подводные и подземные взрывы сб. ст

.pdf
Скачиваний:
27
Добавлен:
22.10.2023
Размер:
15.49 Mб
Скачать

Р.Мбар

Рис. 14. Измеренные и рассчитанные пиковые давления в ударной волне при взрыве «Бенхам».

О , измерения пикового давления; I - расчет при помощи водо-подпорной модели; 2—расчет при помощи водо-подпорной модели и фазового перехода в твердом компоненте туфа.

282 t. p. вУтКовиЧ

В В Е Д Е Н И Е Ф А З О В О Г О П Е Р Е Х О Д А

Фазовый переход, подобный тому, который наблю­ дал Аренс, возможем в твердом компоненте горной по­ роды, т. е. кривая разгрузки для твердого компонента горной породы от давлений 400 кбар идет резко вниз от адиабаты Гюгонио, а при давлении около 50 кбар

Р и с. 15. Измеренный и рассчитанный радиус полости при взрыве

«Б е н х а м » .

/—измеренный радиус полости; -------- расчет при помощи водо-подпорной модели и фазового перехода;---- — расчет без водо-подпорного эффекта и

фазового перехода.

эта кривая выходит на плотность, соответствующую состоянию до ударного сжатия (это видно на рис. 11). В соответствии с этим были подправлены адиабаты раз­ грузки для двухкомпонентной системы и проведены рас­ четы по программе SOC. Результаты расчетов приве­ дены на рис. 13 и 14 совместно с результатами расчета без фазового перехода. Как видно, экспериментальные данные и расчеты с учетом фазового перехода нахо­ дятся в значительно лучшем согласии. На рис. 15 пока­ зано изменение размера полости во времени, получен­ ное при помощи двух расчетов: в одном использовалась водо-подпорная модель с учетом фазового перехода, а

ВЛИЯНИЕ ВОДЫ НА ЭФФЕКТЫ ПОДЗЕМНЫХ ВЗРЫВОВ

283

в другом — не

учитывалось ни то, ни другое;

там

же

для сравнения

приведен измеренный радиус

полости.

На рис. 16 приведен профиль напряжений в зависимо­ сти от расстояния через 50 мс после взрыва, получен­ ный в этих двух расчетах.

Рис. 16. Рассчитанный импульс радиальных напряжений через 50 мс после взрыва «Бенхам».

-------- расчет при помощи водо-подпорного эффекта и фазового перехода;

------ расчет без водо-подпорного эффекта или фазового перехода.

Так как расчеты для взрывов «Хардхэт» и «Скунер» проводились первоначально без учета фазового пере­ хода, было принято, что и в этих случаях происходит фазовый переход, подобный тому, который учитывался в расчетах для взрыва «Бенхам», и были внесены соот­ ветствующие поправки в кривые разгрузки. С учетом этих изменений расчеты были повторены; оказалось, что лишь небольшое различие наблюдается между расче­ тами, выполненными по водо-подпорной модели с уче­ том и без учета фазовых переходов. На рис. 17 приведен

Р и с . 17. Рассчитанный импульс радиальных напряжений

через

12 мс после взрыва «Хардхэт».

 

-------- расчет без фазового перехода в твердом компоненте гранита;-------

расчет

с фазовым переходом.

 

Рис. 18. Рассчитанная скорость движения свободной поверхности при взрыве «Скунер».

По осп ординат: скорость свободной поверхности, м /с ;

------ расчет с фазовым

переходом в твердом компоненте туфа;--------

расчет

без фазового перехода

ВЛИЯНИЕ ВОДЫ НА ЭФФЕКТЫ ПОДЗЕМНЫХ ВЗРЫВОВ

285

график распределения радиального напряжения с рас­ стоянием через 12 мс после взрыва «Хардхэт», а на рис. 18 — рассчитанная скорость свободной поверхности для взрыва «Скунер».

ВЫВОДЫ И ЗАКЛЮЧЕНИЯ

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

Эта модель была проверена путем проведения рас­ четов по одномерной программе SOC и сравнения ре­ зультатов расчетов с экспериментальными данными, полученными в горных породах с содержанием воды от 1 до 25% по весу. Для гранита при взрыве «Хардхэт» с низким содержанием воДы различие между расчетами при помощи водо-подпорной модели и без нее оказалось относительно небольшим. Однако результаты, получен­ ные при помощи водо-подпорной модели, находятся в лучшем согласии с измерениями пиковых радиальных напряжений в области порядка - 1 кбар. Прекрасное со­ гласие между рассчитанными и измеренными скоро­ стями свободной поверхности было получено для туфа при взрыве «Скунер» с высоким содержанием воды в противоположность расчетам без введения водо-подпор­ ной модели, которые дали на 30% заниженные зна­ чения.

Для взрыва «Бенхам» с энергией 1 Мт наблюдалась значительная разница между рассчитанными и измерен­ ными временами прихода и пиковыми давлениями в об­ ласти порядка 100 кбар. Измерения указывали, что пиковые давления падают с расстоянием значительно быстрее, чем это получено из расчетов.- Дальнейшие

286 Т. Р. БУТКОВИЧ

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

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

форме импульса давления в области

порядка 1 кбар и

в кривой изменения размера полости.

непосредственно

Влияние воды в горной породе,

окружающей ядерный взрыв, зависит от ее количества. Даже небольшое количество воды может повлиять на развитие ядерного взрыва, что не объясняется лишь уменьшением прочности горной породы. Во всех прове­ денных расчетах как при помощи водо-подпорной мо­ дели, так и без нее предполагалось, что прочность гор­ ной породы одна и та же. Для гранита взрыва «Хардхэт» с содержанием воды 0,5% по весу пиковое напря­ жение на 120 м от центра взрыва (рис. 8) было на 15% выше в случае водо-подпорной модели. Для туфа при

взрыве «Бенхам», где принималось, что

содержание

воды по весу составляет

13,4%, пиковые

напряжения

на 150 м от центра взрыва

(рис. 16)

в случае водо-под­

порной модели были приблизительно

на 50% выше.

Фазовые изменения, которые происходят в силикат­ ных горных породах, сжатых ударной волной до давле­ ний, превышающих 100 кбар, существенны при опреде­ лении времени прихода и пикового давления в цент­

ВЛИЯНИЕ ВОДЫ НА ЭФФЕКТЫ ПОДЗЕМНЫХ ВЗРЫВОВ

287

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

Список литературы

1.Ahrens Т. J., Rosenberg J. Т., Shock metaniorphism: experiments on quartz and plagioclase, in «Shock Metamorphism of Natural Materials», Mono Book Co., Baltimore, Mr., 1968.

2.Ahrens T. J., Petersen C. F„ Rosenberg J. T„ Shock compression feldspars, J. Geophys. Res., 74, 2727—2746 (1969).

3.Butkovich T. R., Calculation of the shock wave from an under­

ground nuclear explosion in granite, J. Geophys. Res.,

70, 885—

892 (1965).

materials,

4 Butkovich T. R., The gas equation of state of natural

Lawrence Radiat. Lab. Rep., UCRL-14729, Livermore,

California,

1967.

 

5.Cherry J. T., Petersen F. L., Numerical simulation of stress wave propagation from underground nuclear explosions, Lawrence Ra­ diat. Lab. Rep. UCRL-72216, Livermore, California, 1970.

6.Heusinkveld M., Flolzer F., Method of continuous shock front

position

measurement, Rev. Sci. Instrum., 35, 1105

(1964).

7. Higgins

G. I-L, Butkovich T.

R.,

Effect

of water

content, yield,

medium,

and depth of burst

on

cavity

radii, Lawrence Radiat.

Lab. Rep. UCRL-50203, Livermore,

California, 1967.

 

8.ITolzer A., Mesurements and calculations of peak shock wave pa­ rameters from underground nuclear detonation, Lawrence Radiat. Lab. Rep. UCRL-7805, Livermore, California, 1964.

9.

Kennedy G. C„

Wasserburg G. J., Heard IT. C., Newton R. C.,

 

The upper

three-phase

region

in the system SiC>2 —H20,

Amer.

10.

J. Sci., 260,

501

(1962).

The tensor code, Lawrence Radiat.

Lab.

Maenchen G., Sack S.,

11.

Rep., UCRL-7316, Livermore, California, 1963.

solids

Rice M. H., McQueen R. G., Walsh J. M., Compression of

12.

by strong

shock

waves, Solid State Phys., 6, 40 (1958).

PV

Stephens D. R.,

Louis H., Lilley E. M., Loading — unloading

 

curves

for

tuff

at

the

Nevada test site, Lawrence Radiat.

Lab.

13.

Rep. UCRL-50554,

Livermore, California, 1969.

 

Terhune

R. W.,

Stubbs

T. F.,

Cherry J. T., Nuclear cratering on

 

a digital computer, Lawrence Rad. Lab. Rep., UCRL-72032, Li­

 

vermore, California, 1970.

 

 

14.Tewes H. A., Results of the Scooner excavation experiment, Law­ rence Radiat. Lab. Rep., UCRL-71832, Livermore, California, 1970.

15.Von Thiel M., Shock wave data, Lawrence Radiat. Lab. Rep. UCRL-50108, 3 vol., Livermore, California, 1966.

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

Ч. 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) и

Соседние файлы в папке книги из ГПНТБ