CC BY
40
УДК 622.81:24 Л.В. Винокуров
СПОСОБ ОПРЕДЕЛЕНИЯ КООРДИНАТ ПОДЗЕМНОГО ЯДЕРНОГО ВЗРЫВА
Семинар № 4
звестно, что при подземном ядерном взрыве образуется, в первом приближении, сейсмическая волна в виде сферы с центральной симметрией, при этом для определения координат подземного ядерного взрыва использовалась методика, разработанная сейсмологами.
«К сожалению «обычные» точности лоцирования очагов землетрясений низки» - указывает академик РАЕН Шебалин
Н.В. Ошибка достигает 10-100 км и более.
Для быстрого реагирования и приборного определения координат и сейсмической энергии подземного ядерного взрыва требуется создание портативных легко устанавливаемых на местности технических средств. В настоящее время мной выполнены теоретические исследования и конструктивные разработки по созданию таких средств.
В основе теоретических исследований лежит теорема о вращающейся прямой вокруг центра равностороннего треугольника. Ниже приводится ее формулировка.
Если некоторая прямая иУ, первоначально совпадающая с биссектрисой АО равностороннего треугольника АВС со стороной I вращается в плоскости треугольника вокруг его центра О (см. рис.
1), то проекции 8Ь 82, 83 сторон треугольника на эту прямую изменяются по закону синуса, т.е.
51 = I 8ш(60° + р), 52 = 18Шр,
53 = 18т(60 - р) ( )
Доказательство. Перемещая прямую ИУ параллельно своему первоначальному по-
ложению до угла А, получим два прямоугольных треугольника АСС’ и АВВ’, из которых для проекций получим значения
А'С' = 81 = Ми( 60° + ф)
А'В' =83 = ^ш(60° - ф).
Перемещая прямую ИУ до угла С, получим треугольник ВСВ', из которого получим проекцию В'С' = 82 = /БШф,
что и требовалось доказать.
Из этой теоремы следует, что, если проекции сторон равностороннего треугольника определены, то положение прямой ИУ относительно заданного начала отсчета АО определяется однозначно путем использования обратных тригонометрических функций.
Для удобства использования теоремы в целях определения направления (пеленга) на источник излучения сейсмоимпульса целесообразно принять вместо прямой ИУ нормаль к фронту сейсмоволны, проходящую через центр треугольника, а вместо проекций, соответственно, временные задержки 11, 12, 13.
В общем случае Т = 8/У, (2)
где У - скорость распространения продольной сейсмоволны в объеме треугольника АВС; 8 - путь пройденный сигналом,
равный I; Т - время прохождения пути, равное I /У.
Тогда для задержек времени получим /1 = Т sin(600 + р), /2 = Т sinр,
13 = Т sin(600 -р) (3)
Если задержки времени измерены, то пеленг (угол р ) однозначно определяется из соотношений (3). Однако для большей точности определения пеленга целесообразно использовать меньшую задержку, так как при изменении аргумента в окрестности нуля синусоидальная функция имеет большую скорость роста и хорошую линейность.
На основе теоремы и приведенных рассуждений был разработан алгоритм (рис.
2) для ЭВМ, который позволяет по сигналу определять пеленг ядерного взрыва в реальном времени с помощью системы технических средств.
Технические средства представляют собой систему, включающую равносторонний треугольник из пенопласта, сторона может иметь размеры от 100 мм до 500 мм, в вершинах треугольника установлены сейсмоприемники. Треугольник устанавливается на коренных породах горизонтально.
Система (рис. 3) содержит три канала приема сейсмосигнала от взрыва. Каждый канал состоит из последовательно соеди-
Рис. 1. Определение проекций сторон треугольника
ненных сейсмоприемника 1, усилителя 2, компаратора 3, триггера 4 и счетчика 5. Система также содержит схему совпадения 6, включенную параллельно триггерам 4, генератор тактовых импульсов 7, подключенный к счетчикам
5, мультиплексор 8, подключенный к выходам счетчиков 5 и к входу ЭВМ 9, измерительный блок 10, включенный к выходу усилителя 2 и входу
ЭВМ 9.
Система работает следующим образом.
Предположим, при подземном ядерном взрыве с направления А сейсмический сигнал достигает сейсмоприемника 1 среднего канала, усиливается усилителем 2 и поступает на компаратор 3. Если сигнал по уровню больше порога срабатывания компаратора 3, то он проходит и опрокидывает триггер 4, при этом включается счетчик 5 и считает импульсы, вырабатываемые генератором тактовых импульсов 7. Одновременно с триггера 4 сигнал подается на схему 6 совпадения. Через некоторое время сигнал от очага взрыва достигает второго сейсмоприемника 1 (левый канал). При этом в левом канале, как и в предыдущем случае, срабатывают аналогичные элементы системы, и счетчик 5 левого канала тоже начинает считать импульсы, поступающие с генератора 7. С приходом сигнала от взрыва к третьему сейсмоприемнику 1 срабатывают все элементы этого канала. При срабатывании триггера 4 этого канала на входе схемы совпадения 6 одновременно присутствуют три сигнала, схема совпадения 6 срабатывает. При этом от ее выходного сигнала триггеры возвращаются в исходное со-
( начало)
™Я[—,
/ ььо л 7
Дг.у.гз/
нет
Ч>=іго->-+0^1 п
Ч^ЬЄОг
-а’гс^іп
іг
\/
£
Р
/Р€
^=24 °гг п -~
\ /
СП | (
*Ьг
4>=а,гс5'ип у
\/
240 + кггс$.і.п ^г
V
16
Ч=іго-
•
-с&сьт -=^
л/.
ьыьод
РеЬЬЛЬША1
Рис. 2. Алгоритм определения пеленга на источник сейсмоакустического импульса
стояние, счетчики 5 прекращают считать минает результат. С выхода усилителя 2
импульсы, а ЭВМ одного из каналов сигнал поступает на
9 через мультиплексор 8 считывает пока- вход измерительного блока 10, который
зания счетчиков 5, вычисляет пеленг ис- измеряет величину амплитуды и период
точника сейсмических колебаний и запо- по принципу аналого-цифрового преобра-
Рис. 3. Блок-схема системы определения координат
зования и передает результат в ЭВМ 9 для регистрации. Имеется возможность вычислить энергию взрыва.
Высокая точность определения пеленга на плоскости будет в том случае, когда точка излучения сигнала и равносторонний треугольник с сейсмоприемниками размещены в одной плоскости. С увеличением угла между направлением на очаг взрыва и плоскостью треугольника увеличивается ошибка в определении пеленга, которая при 200 достигает 6 %.
При исследовании возможностей определения координат гипоцентра в реальном времени проблема была решена дальнейшими разработками теории для полупространства.
Было выяснено, что с увеличением угла места в все проекции сторон равностороннего треугольника на вращающуюся прямую ИУ уменьшаются. При этом найдена
закономерность, что уменьшение одной малой временной задержки влечет уменьшение значений пеленга (рх, а уменьшение другой - к увеличению пеленга (р2
на такую же величину. Отсюда следует, что действительный (точный) пеленг на эпицентр может быть определен по формуле
9 = (ф\ + Ф2)/2 (4)
Таким образом, имеется возможность определить пеленг на эпицентр взрыва, не взирая на ориентацию треугольника относительно местоположения взрыва. Для этого необходимо выполнить следующие операции (см. рис. 4).
1. Определить скорость продольной волны Ур горного массива под треугольником.
где Тс - время в секундах отставания поперечной волны S от продольной волны Р. 9. Определить дальность до эпицентра
Йэц
Лэц = Лгц ■ c0sР. (10)
10. Определить глубину гипоцентра Z Z = h = dгц ■ sin в. (11)
Разработанная система может использоваться в быту физическими лицами для регистрации, например, форшоков и землетрясений. В случае регистрации форшоков с близкими по значению координатами многократно за установленный период времени, то можно предположить, что в этом квадрате возможно динамическое явление - афтершок.
Интерфейс, описанный в данной статье, может быть использован совместно с персональной ЭВМ.
2. Определить скорость поперечной волны Vs (можно из соотношения Vp * 1,7Vs).
3. Измерить временные задержки при регистрации взрыва tj и t2
t1 = T sin (p1; t2 = T sin(60 - (p2) (5)
4. Определить пеленг
p1 = arcsinT.; p2 = 600 - arcsinT (6) 5. Определить уточненный пеленг
P = (P+P2)/2.
6. Определить временные задержки для земной поверхности (для плоскости размещения треугольника)
t| = T sin p; t2 = T sin(600 - p) . (7)
7. Определить угол места в гипоцентра
СПИСОК ЛИТЕРАТУРЫ
1. Ракишев Б.Р., Винокуров Л.В. Пеленгация источников возмущения в массиве горных пород. - Алматы, НИЦ «Былым», 2002. - 236 с.
2. Винокуров Л.В. Автореферат на соис-
кание ученои степени доктора технических наук «Геоинформационные системы определения
сейсмоакустических источников в массиве горных пород». - Алматы, 2002.
Коротко об авторах
