Научная статья на тему 'Исследование зависимости ошибок в определении координат очага землетрясения от методов расчета (сфер и гиперболоидов)'

Исследование зависимости ошибок в определении координат очага землетрясения от методов расчета (сфер и гиперболоидов) Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
81
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЗЕМЛЕТРЯСЕНИЕ / EARTHQUAKE / СЕЙСМОДАТЧИК / ОШИБКА / ERROR / УРАВНЕНИЕ / EQUATION / ГИПЕРБОЛОИД / HYPERBOLOID / СФЕРА / SPHERE / ПРВ / SEISMIC SENSOR / PROBABILITY DISTRIBUTION DENSITY

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Асланов Гайдарбек Тагирбекович, Асланов Тагирбек Гайдарбекович, Курбанмагомедов Курбанмагомед Динмагомедович, Шахтарин Борис Ильич

Цель. Статья является продолжением серии статей по анализу влияния типа используемых поверхностей второго порядка на точность определения координат эпицентра и гипоцентра землетрясения. Метод. Для нахождения плотности распределения вероятностей ошибок в определении гипоцентра землетрясения используются метод сфер, комбинированный метод сфер и гиперболоидов, а также метод гиперболоидов. Методы, используемые для определения координат гипоцентра, гиперболоиды, имеют меньшие ошибки по сравнению с методом сфер. Это объясняется тем, что при определении времен пробега сейсмических волн введено допущение, что прирост ошибки одинаков для разности времен прихода сейсмической волны к двум сейсмодатчикам для методов, использующих гиперболоид и для разности времен прихода двух сейсмических волн к сейсмодатчику, для методов использующих сферу. Результат. В работе получено выражение для определения координат очага землетрясения методом гиперболоидов, а также плотности распределения вероятностей ошибок в определении гипоцентра землетрясения при расчетах по методу сфер, по комбинированному методу сфер и гиперболоидов, а также по методу гиперболоидов. Составлен график распределения ошибок в определении гипоцентра землетрясения при разных размещениях сейсмодатчиков и для различных значений ошибок в определении разностей времен пробега сейсмических волн. Вывод. Полученные зависимости имеют вид распределения ошибок, близких к распределению Коши. Всплеск в районе нуля у всех распределений получен в результате расчетов координат очага землетрясения при отсутствии ошибок в определении разности времен. Комбинированный метод гиперболоида и сферы имеет форму распределения близкую к методу гиперболы и кривую распределения в районе, близкому к нулю, схожую с методом сфер.

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

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Асланов Гайдарбек Тагирбекович, Асланов Тагирбек Гайдарбекович, Курбанмагомедов Курбанмагомед Динмагомедович, Шахтарин Борис Ильич

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

INVESTIGATION OF THE DEPENDENCE OF EARTHQUAKE FOCUS COORDINATE DETERMINATION ERRORS ON CALCULATION METHODS (SPHERES AND HYPERBOLOIDS)

Objectives This work is a continuation of a series of articles devoted to an analysis of the effect of the utilised second-order surface types on the accuracy of determining earthquake epicentre and hypocentre coordinates. Methods In order to find the density of error probability distribution during the determination of earthquake hypocentres, approaches using spheres, hyperboloids, as well as combined spheres and hyperboloids, are used. The hyperboloid-based methods used for determining hypocentre coordinates have fewer errors as compared to the sphere-based method. This is explained by the fact that when determining the travel times of seismic waves, it is assumed that the error increase is the same for the difference in the seismic wave arrival times to two seismic sensors (for methods using hyperboloid), and for the difference in the arrival times of the two seismic waves to the one seismic sensor (for methods using spheres). Results The expression for determining the earthquake focus coordinates by the hyperboloid method, as well as the density of error probability distribution during the earthquake hypocentre determination, is obtained for calculations using the sphere method, the hyperboloid method, as well as the combined sphere and hyperboloid method. A graph for the error distribution when determining the earthquake hypocentre is obtained for different locations of seismic sensors and for various error values concerning differences in the travel times of seismic waves. Conclusion The obtained dependencies have the form of an error distribution close to the Cauchy distribution. A wavelet in the zero regions for all distributions was obtained as a result of calculations of the earthquake focus coordinates in the absence of errors in the determination of the time difference. The combined method of the hyperboloid and the sphere has a distribution form close to the hyperboloid method, while the distribution curve in the region close to zero is similar to the sphere method.

Текст научной работы на тему «Исследование зависимости ошибок в определении координат очага землетрясения от методов расчета (сфер и гиперболоидов)»

Для цитирования: Асланов Г.К., Асланов Т.Г., Курбанмагомедов К.Д., Шахтарин Б.И. Исследование зависимости ошибок в определении координат очага землетрясения от методов расчета (сфер и гиперболоидов). Вестник Дагестанского государственного технического университета. Технические науки. 2017;44 (4):87-98. DOI: 10.21822/2073-6185-2017-44-4-8 7-98

For citation: Aslanov G.K., Aslanov T.G., Kurbanmagomedov K.D., Shakhtarin B.I. Investigation of the dependence of earthquake focus coordinate determination errors on calculation methods (spheres and hyperboloids). Herald of Daghestan State Technical University. Technical Sciences. 2016;44 (4):87-98. (In Russ.) D0I:10.21822/2073-6185-2017-44-4-87-98

ТЕХНИЧЕСКИЕ НАУКИ ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И УПРАВЛЕНИЕ

УДК 550.34.01

DOI: 10.21822/2073-6185-2017-44-4-87-98

ИССЛЕДОВАНИЕ ЗАВИСИМОСТИ ОШИБОК В ОПРЕДЕЛЕНИИ КООРДИНАТ ОЧАГА ЗЕМЛЕТРЯСЕНИЯ ОТ МЕТОДОВ РАСЧЕТА (СФЕР И ГИПЕРБОЛОИДОВ)

Асланов Г.К.1, Асланов Т.Г.2, Курбанмагомедов К.Д.3, Шахтарин Б.И.4

1'2Дагестанский государственный технический университет,

1-2

- 367026, г. Махачкала, пр. Имама Шамиля, 70, Россия,

3

Институт системных технологий,

ß

367015, г. Махачкала, пр. Али-Гаджи Акушинского, 21, Россия,

4Московский государственный технический университет им. Н.Э. Баумана,

4105005, г. Москва, ул. 2-я Бауманская, д. 5 стр. 1, Россия,

1-2e-mail: [email protected], 3e-mail:[email protected], 4e-mail: [email protected]

Резюме: Цель. Статья является продолжением серии статей по анализу влияния типа используемых поверхностей второго порядка на точность определения координат эпицентра и гипоцентра землетрясения. Метод. Для нахождения плотности распределения вероятностей ошибок в определении гипоцентра землетрясения используются метод сфер, комбинированный метод сфер и гиперболоидов, а также метод гиперболоидов. Методы, используемые для определения координат гипоцентра, гиперболоиды, имеют меньшие ошибки по сравнению с методом сфер. Это объясняется тем, что при определении времен пробега сейсмических волн введено допущение, что прирост ошибки одинаков для разности времен прихода сейсмической волны к двум сейсмодатчикам для методов, использующих гиперболоид и для разности времен прихода двух сейсмических волн к сейсмодатчику, для методов использующих сферу. Результат. В работе получено выражение для определения координат очага землетрясения методом гиперболоидов, а также плотности распределения вероятностей ошибок в определении гипоцентра землетрясения при расчетах по методу сфер, по комбинированному методу сфер и гиперболоидов, а также по методу гиперболоидов. Составлнг график распределения ошибок в определении гипоцентра землетрясения при разных размещениях сейсмодатчиков и для различных значений ошибок в определении разностей времен пробега сейсмических волн. Вывод. Полученные зависимости имеют вид распределения ошибок, близких к распределению Коши. Всплеск в районе нуля у всех распределений получен в результате расчетов координат очага землетрясения при отсутствии ошибок в определении разности времен. Комбинированный метод гиперболоида и сферы имеет форму распределения близкую к методу гиперболы и кривую распределения в районе, близкому к нулю, схожую с методом сфер.

Ключевые слова: землетрясение, сейсмодатчик, ошибка, уравнение, гиперболоид, сфера,

ПРВ

TECHNICAL SCIENCE COMPUTER SCIENCE, COMPUTER ENGINEERING AND MANAGEMENT

INVESTIGATION OF THE DEPENDENCE OF EARTHQUAKE FOCUS COORDINATE DETERMINATION ERRORS ON CALCULATION METHODS (SPHERES AND HYPERBOLOIDS)

1 2

Gaidarbek K.Aslanov , Tagirbek G. Aslanov ,

Kurbanmagomed D. Kurbanmagomedov 3, Boris I. Shakhtarin4

1'2Daghestan State Technical University,

1-2701. Shamilya Ave., Makhachkala 367026, Russia,

3Institute of System Technologies,

321 Ali-Haji Akushinsky Ave., Makhachkala 367015, Russia,

4Bauman Moscow State Technical University,

45/1 2-ya Baumanskaya Str., Moscow 105005, Russia,

1-2e-mail: [email protected], 3e-mail:[email protected], 4e-mail: [email protected]

Abstract. Objectives This work is a continuation of a series of articles devoted to an analysis of the effect of the utilised second-order surface types on the accuracy of determining earthquake epicentre and hypocentre coordinates. Methods In order to find the density of error probability distribution during the determination of earthquake hypocentres, approaches using spheres, hyperboloids, as well as combined spheres and hyperboloids, are used. The hyperboloid-based methods usedfor determining hypocentre coordinates have fewer errors as compared to the sphere-based method. This is explained by the fact that when determining the travel times of seismic waves, it is assumed that the error increase is the same for the difference in the seismic wave arrival times to two seismic sensors (for methods using hyperboloid), and for the difference in the arrival times of the two seismic waves to the one seismic sensor (for methods using spheres). Results The expression for determining the earthquake focus coordinates by the hyperboloid method, as well as the density of error probability distribution during the earthquake hypocentre determination, is obtained for calculations using the sphere method, the hyperboloid method, as well as the combined sphere and hyperboloid method. A graph for the error distribution when determining the earthquake hypocentre is obtained for different locations of seismic sensors and for various error values concerning differences in the travel times of seismic waves. Conclusion The obtained dependencies have the form of an error distribution close to the Cauchy distribution. A wavelet in the zero regions for all distributions was obtained as a result of calculations of the earthquake focus coordinates in the absence of errors in the determination of the time difference. The combined method of the hyperboloid and the sphere has a distribution form close to the hyperboloid method, while the distribution curve in the region close to zero is similar to the sphere method.

Keywords: earthquake, seismic sensor, error, equation, hyperboloid, sphere, probability distribution density

Введение. В настоящее время разработаны программы [6-9], позволяющие по магнитуде землетрясения, глубине его очага, плотности населения в районе землетрясения, типам застроек [10-12], времени суток и т.д. оценить людские потери [13], что позволяет оперативно решить вопрос о количестве привлекаемых к спасательным работам людских, материальных и технических ресурсов. Однако оценки потерь иногда не соответствуют действительности, в основном, из-за неверного определения глубины очага и эпицентра землетрясения [14-20].

Определение координат очага землетрясения, в основном, производится по разности времен пробега поперечной и продольной сейсмических волн на разнесенные в пространстве сейсмодатчики [21-22]. При этом основную ошибку в определении координат очага землетрясения вносит неверный выбор сейсмодатчиков [1-2, 19, 22-23].

В работе получены плотности распределения вероятностей ошибок в определении гипо-

центра землетрясения при расчетах по методу сфер [3], по комбинированному методу сфер и гиперболоидов [2 - 4], а такжепо методу гиперболоидов.

Постановка задачи. Для нахождения плотности распределения вероятностей ошибок в определении гипоцентра землетрясения используются метод сфер, комбинированный метод сфер и гиперболоидов, а также метод гиперболоидов. Во всех методах для определения гипоцентра землетрясения в качестве исходных данных используются значения скоростей поперечной и продольной сейсмических волн [21]. Кроме скоростей сейсмических волн, в методе сфер используются значения разностей времен прихода продольной и поперечной сейсмических волн на сейсмодатчики, в комбинированном методе сфер и гиперболоидов используются разности времен прихода продольной и поперечной волн на опорный сейсмодатчик и разность времен прихода продольной волны к двум разнесенным сейсмодатчикам, а в методе гиперболоидов используется разность времен прихода продольной волны к каждой паре из четырех сей-смодатчиков.

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

Методы исследования. В работе [3] получена формула, позволяющая вычислить координаты очага землетрясения методом сфер:

Г Р^у V у2у V I у2у ¡_ р2у г>2V Р^у

_ ~ Г3Л2 ~ Г2 Г3 ~ Г2Л3 т Г3 Г2 + К1Г2 ~ К3 Г2 ~ К1 Г3

= 2(У2Х3-Х2У3)

С>2 у у у 2 у2у у у 2 ¡_ у2у ¡_ г>2 у п2 у г>2 у < у _ К2Л3 ~ Л3Г2 ~ Л2Л3 ~ Л2Г3 "Г Л3Л2 к1л2 — К3Л2 — к1л3

= 2 (Х2У3 - У2 Х3 )

7 = ^Е 2-Х2-У2

где X, У, Z - координаты гипоцентра землетрясения;

X, - координаты соответственно второго и третьего сейсмодатчика;

Я - радиусы сфер, в центрах которых расположены сейсмодатчики с радиусами равными расстоянию до гипоцентра землетрясения, которые определяются формулой ( );

VI, У2 - скорости соответственно поперечной и продольной волны; II - разность времен прихода сейсмических волн к сейсмодатчику.

В свидетельстве на государственную регистрацию программы ЭВМ [4] приведена формула, позволяющая вычислить координаты очага землетрясения комбинированным методом сфер и гиперболоида:

( _ Х2(Уг - У2) - tr2У12(У1 - У2) - 2^±У2У2

2Х2(Уг - У2) , (Х3-Х)2 + У32 + Е2-Х2-Е32

2Г3

7 = ^Е2 -Х2 - У 2

где 1Г - разность времен прихода продольной сейсмической волны к двум сейсмодатчи-

кам.

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

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

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

В этом случае ставится задача, получить выражение для определения координат очага землетрясения методом гипербол.

На рис. 1 приведено пространственное расположение координат сейсмодатчиков и очага землетрясения. На рис. 1а приведено расположение четырех сейсмодатчиков Вг и эпицентра землетрясения О' с видом сверху.

На рис. 1б приведена косоугольная проекция для двух сейсмодатчиков и очага землетрясения, где буквой О обозначен очаг землетрясения.

О'

а)

О'

Е

Do

О

H

По теореме Пифагора имеем:

H2 +Е2 = R2; H2 +El = Rl)

(1) (2)

где Н - глубина очага землетрясения;

Ег - расстояния от эпицентра землетрясения О' до сейсмодатчиков Д-;

Ri - расстояния от очага землетрясения O до сейсмодатчиков Di

Вестник Дагестанского государственного технического университета. Технические науки. Том 44, №4, 2017 Heraldof Daghestan State Technical University.TechnicalSciences. Vol.44, No.4, 2017 _http://vestnik.dgtu.ru/ISSN (Print) 2073-6185ISSN (On-line) 2542-095Х_

По теореме косинусов получим:

Si + R2 - 2S±R cos a = Rf;

Si + E2 — 2SXE cos aa = El

где Si - расстояния между опорным сейсмодатчиком D0 и остальными сейсмодатчика-ми;

а - угол, образованный между очагом землетрясения и первым сейсмодатчиком с вершиной, в точке расположения опорного сейсмодатчика;

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

S?+R2-R?

cos а = —--

25!«

cos аа =

S^ -\-Е^ 2 Si_E

(3)

(4)

Подставляя уравнения (1) и (2) в уравнение(4), получим

cos аа =

Sl+R2-Rf

(5)

2 5! V«2- Н 2

Определим расстояние Я через расстояние Я и разности расстояний ДЯ1, из уравнений

(3) и (5)

s2+ Л Ri(I R - Л R i)

cos а = —---—

25!« S2+ Л R i( I R - Л R i)

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

cos аа = —-, —

I SiVR2- H 2

Разделим уравнение (6) на уравнение (7), получим

(6) (7)

cos a

cos аа

= V1-HI/RI

(8)

Произведя замену в уравнении (8) на sin е, после некоторых преобразований полу-

чим

cos а = cos аа cos е

(9)

где е - угол, образованный между эпицентром и гипоцентромземлетрясения с вершиной угла, в точке расположения опорного сейсмодатчика

Вестник Дагестанского государственного технического университета. Технические науки. Том 44, №4, 2017 Herald of Daghestan State Technical University.Technical Sciences. Vol.44, No.4, 2017 _http://vestnik.detu.ru/ISSN (Print) 2073-6185 ISSN (On-line) 2542-095Х_

Аналогично, получим уравнения и для углов ß и у

c o sß = c o sßßc o se (10)

c o s у = c o s уу c o s e (11)

Уравнение гиперболы, с фокусами в точках расположения опорного и первого сейсмо-датчиков в полярной системе координат имеет вид [2] :

S 2- (Л t 1V1 ) 2

R = -1 i 1 1}--(12)

2At1Vr1+251 cosa

где V1 - скорость продольной сейсмической волны;

At1 - разность времен прихода продольной сейсмической волны к опорному и первому сейсмодатчику.

Задача определения местоположения объекта на плоскости решена в [2], для перехода в трехмернное пространство произведем подстановку уравнений (9) - (11) в уравнения гипербол, образованных опорным и остальными тремя сейсмодатчиками, при этом, произведя замену Ati Vi на разность расстояний Si, получим

sl-äsl

R = т-1-1-: (13)

( )

R = s2- л s! (14) 2 (Л S2- S2 с о sß ßc о s s) (

n sf-Asf

R = T-3-3-Г (15)

( )

Подставим в уравнения (13) - (15) выражения ßß — aa = aß и уу — aa = ay, и разделим полученные уравнения на R, получим

г, 5?—At?+25i cosaaE

R = —-1-1--(16)

2A t1 v y

_ S|- Л t|+2 S2 c оs(aß+а а)£ H —-

2At,

_ S3 - Л t|+2 S3 c о s^T^а а)£

К —-

2At3

(17)

(18)

Перенесем в левую часть уравнения (16) расстояние до эпицентра Е, и подставим полученное выражение в уравнения (17) и (18).

2Д^Й - Si + Ml Е =-^-

25Х cos aa

_ Sf Si - A tf 5i+Sf(c О sa ß+s i n«/?tana a)( 2 A t ^ - Sf+A t f )

n —--(19)

2At2S1 v 7

_ Sf Si - A t| 5i+S3(c о say+s i nay tan a a )( 2 A t jfl - Sf +A t f)

n —-

2At3S1

Из уравнения (19) найдем угол tan аа,

(20)

_ 2RAtfSi-SfS1+Atf S1-Sf cosaß(2At1 R -Sf +Atf) . .

til П — "7 (21)

S2s i n a ß( 2 A tiR - Sf+A tf ) v 7

Подставив уравнение (21) в (20), после некоторых преобразований получим:

й = (5Х5253(53 б 1 п а// — 52 б 1 п ау) + 5Х(Д ^ 53 б 1 п ау — Д t252 б 1 п а//) +5352 б 1 п( ау — а// )

+Д ^5352 б 1 п( а/? — ау))/( 2 Д ^5253 б 1 п( ау — а/? ) +2 5Х (Д t352 б 1 п а/ — Д ^ 53 б 1 п ау)) (22)

Уравнение (22) позволяет определить расстояние от опорного сейсмодатчика до очага землетрясения. Поскольку нам известны разности времен прихода сейсмической волны к опорному и к остальным сейсмодатчикам, то можно найти расстояния от очага землетрясения и до остальных сейсмодатчиков:

йх = ДЗД + й; й2 = + Я; й3 = + Я;

Для проведения сравнительного анализа с результатами работ [3-5] необходимо определить координаты очага методом гипербол в декартовых координатах.

Для этого воспользуемся формулой, приведенной в работе [3]:

п2у у у2 у2у _1_ V . у2у ¡_ п2у п2у г>2у у _ к2*3 ~ Г3Л2 2 3 """ Г2Л3 "г г3 Г2 "г К1 г2 ~ К3Г2 ~ к1*3

= 2 ( К2 Х3 — ;

/?2у у у2 у2у I у у2 ¡_ у2 у ¡_ п2 у п2у /?2у _ л2л3 — л3г2 — л2л3 ■+■ л2г3 ■+■ л3л2 ■+■ «!а2 — к3л2 — «з^аэ

= 2 (ВД—ВД ) ;

7 = ^й 2 — X 2 — Г 2.

Для получения плотности распределения ошибки в определении очага землетрясения зададим область размерами 150 на 150 км, в которой будем менять взаимное расположение сейсмодатчиков. Примем за координаты очага землетрясения X = 75 км., У = 75 км., Z = 10 км.

На рис.2 и 3 приведены нормированные плотности распределения ошибки в определении гипоцентра землетрясения. Для получения распределения производилось определение координат очага землетрясения с изменением координат сейсмодатчиков с шагом 18,75 км, также были введены преднамеренные ошибки в определении разности времен, по разности координат очагов землетрясений с ошибкой и без нее, и была получена дисперсия.

На рис.2, кривой 1 соответствует метод сфер, кривой 2 - комбинированный метод сфер и гиперболоида, а кривой 3 - метод гипербол.

При расчете производилось изменение ошибки в определении разности времен прихода сейсмических волн в интервале от - 0.5 до 0.5 с шагом 0.25 с.

Рис.2. Плотность распределения ошибок в определении гипоцентра землетрясения тремя

методами

Fig. 2. Density of error distribution in determining the hypocentre of the earthquake by three methods

а) б)

в)

Рис. 3. Плотность распределения ошибок в определении гипоцентра землетрясения тремя методами с различными интервалами ошибок в определении разности времен прихода

сейсмических волн

Fig. 3. Error distribution density in the determination of the earthquake hypocenter by three methods with different error intervals in determining the difference in arrival times of seismic waves

На рис. 3 кривой 1 соответствует изменение ошибки в определении разности времен прихода сейсмических волн в интервале от - 0.5 до 0.5 с шагом 0.25 с, кривой 2 - изменение

ошибки в определении разности времен прихода сейсмических волн в интервале от - 1 до 1 с шагом 0.5 с, кривой 3 - изменение ошибки в определении разности времен прихода сейсмических волн в интервале от - 1.5 до 1.5 с шагом 0.75 с и кривой 4 - изменение ошибки в определении разности времен прихода сейсмических волн в интервале от - 2 до 2 с шагом 1 с.

На рис. 3 а приведено распределение ошибок в определении гипоцентров землетрясений методом сфер, на рис. 3б - для комбинированного метода сфер и гиперболоида и на рис. 3в -для метода гиперболоидов.

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

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

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

Вывод.

1. Все полученные зависимости имеют вид распределения ошибок, близкий к распределению Коши.

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

3. Комбинированный метод гиперболоида и сферы имеет форму распределения, близкую к методу гиперболы и кривую распределения в районе, близкому к нулю, схожую с методом сфер.

Библиографический список:

1. Асланов Т.Г., Магомедов Х.Д., Мусаева У.А., Тагиров Х.Ю. Влияние пространственного расположения сейсмодатчиков на точность определения гипоцентра землетрясения. Вестник Дагестанского государственного технического университета. Технические науки. 2016;43 (4):73-84. D0I:10.21822/2073-6185-2016-43-4-73-84

2. Асланов Т.Г. Определение координат очага землетрясения с использованием комбинированного метода. Вестник Дагестанского государственного технического университета. Технические науки. 2017;44 (2):118-125. DOI: 10.21822/2073-6185-2017-44-2-118-125

3. Асланов Т.Г., Даниялов М.Г., Магомедов Х.Д., Асланов Г.К. Об одном методе определения очага землетрясения с одновременным определением скоростей сейсмических волн // Труды института геологии Дагестанского научного центра РАН, Материалы. Издательство ДНЦ РАН. Махачкала 2010. - С. 54-59.

4. Определение координат гипоцентра землетрясения при использовании комбинированного метода сфер и гиперболоида: свид. о гос. рег. прогр. ЭВМ № 2017616673 Рос. Федерация / Асланов Т.Г., Давудова Д.Д., Мирзаханова Л.С.; заявитель и правообладатель ВГУЮ (РПА Минюста России). - № 2017611151; заявл. 13.02.17; опубл. 09.07.17.

5. Шахриманьян М.А., Нигметов Г.М., Сосунов И.В. Математическое моделирование как способ поддержки принятия решений в случае возникновения чрезвычайных ситуаций // Каталог «Пожарная безопасность» -2003. - С. 240-241.

6. Гитис В.Г., Ермаков Б.В. Основы пространственно-временного прогнозирования в геоинформатике / В.Г. Гитис, Б.В. Ермаков — М.: ФИЗМАТЛИТ, 2004.— 256 с.

7. Jing Z., Laurie G. B., MagalyK.. Mapping earthquake induced liquefaction surface effects from the 2011 Tohoku earthquake using satellite imagery // 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 2016, pp. 2328 - 2331

8. Rui J., Shuanggen J. Ionosphericacousitc and rayleigh waves detected by GPS following the 2005 Mw=7.2 northern California earthquake // 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 2016, pp. 3956 - 3959

9. Aixia D., Xiaoqing W., Xiaoxiang Y., Shumin W. The loss assessment method of building earthquake damage using The Remote Sensing and building grid data // 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 2016, pp. 4255 - 4258

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

10. Wei Z., Huan-Feng Sh., Chun-Lin H., Wan-Sheng P. Building damage information investigation after earthquake using single post-event PolSAR image // 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 2016, pp. 7338 - 7341

11. Hao D., Xin X., Rong G., Chao S., Haigang S. Metric learning based collapsed building extraction from post-earthquake PolSAR imagery // 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 2016, pp. 4742 - 4745

12. Liu L. B., Liu M., Wang J. Q. Electromagnetic environment comprehension for radar detection of vital signs at China National Training Center for earthquake search &amp; rescue // 2016 16th International Conference on Ground Penetrating Radar (GPR), 2016, pp. 1-4

13. Радоуцкий, В.Ю. Опасные природные процессы: учеб.пособие / В.Ю. Радоуцкий, Ю.В. Ветрова, Д.И. Ва-сюткина; под ред. В.Ю. Радоуцкого. - Белгород: Изд-во БГТУ, 2013. - 206 с.

14. Архангельский В.Т., Веденская И.А., Гайский В.Н. Руководство по производству и обработке наблюдений на сейсмических станциях СССР / Акад. наук СССР. Совет по сейсмологии. - Москва : Изд-во Акад. наук СССР, 1954

15. Добровольский И.П. Математическая теория подготовки и прогноза тектонического землетрясения. / И.П. Добровольский - М.: ФИЗМАТЛИТ, 2009. —240 с.

16. Гуревич П.С. Психология чрезвычайных ситуаций учебное пособие/ П.С, Гуревич - М.: ЮНИТИ-ДАНА, 2012. - 495 c.

17. Мкртычев О.В. Безопасность зданий и сооружений при сейсмических и аварийных воздействиях: монография/ О.В. Мкртычев - М.: МГСУ, 2010. - 152 c.

18. Асланов Т.Г., Тагиров Х.Ю., Асланов Г.К., Алимерденов В.Ш. Математическая модель расчета энергетического класса, интенсивности и магнитуды землетрясения в реальном масштабе времени // Вестник Дагестанского государственного технического университета. Технические науки. №2 (37), 2015 - Махачкала: ДГТУ, 2015. - С. 66-71

19. Асланов Т.Г. Разработка алгоритма определения координат очага землетрясения, с одновременным определением скоростей сейсмических волн // Научные труды молодых исследователей программы «Шаг в будущее». Том 8. «Профессионал». Москва. 2005. - С. 32-34.

20. RahinulH., ShoaibH., AkterS., AsadullahilG., TahiaF. K.. Earthquake monitoring and warning system // 2015 International Conference on Advances in Electrical Engineering (ICAEE), 2015, pp. 109 - 112

21. Alphonsa A., Ravi G. Earthquake early warning system by IOT using Wireless sensor networks // 2016 International Conference on Wireless Communications, Signal Processing and Networking (WiSPNET), 2016, pp. 1201 - 1205.

22. Асланов Т.Г., Алимерденов В.Ш. Определение структуры земли по статистическим данным времен прихода сейсмических волн произошедших землетрясений // Старт в будущее - 2013. Труды III научно-технической конференции молодых ученных и специалистов. Санкт-Петербург.

23. Асланов Г.К., Гаджиев М.М., Исмаилов Т.А., Магомедов Х.Д. О землетрясениях (прошлое и современность). - Махачкала, Информационно полиграфический центр ДГТУ. 2001.

References:

1. Aslanov T.G., Magomedov Kh.D., Musaeva U.A., Tagirov Kh.Yu. Vliyanieprostranstvennogo raspolozheniya seismodatchikov na tochnost' opredeleniya gipotsentra zemletryaseniya. Vestnik Dagestanskogo gosudarstven-nogo tekhnicheskogo universiteta. Tekhnicheskie nauki. 2016;43(4):73-84. DOI:10.21822/2073-6185-2016-43-4-73-84 [Aslanov T.G., Magomedov Kh.D., Musaeva U.A., TagirovKh.Yu. Influence of the spatial arrangement of seismic sensors on the accuracy of determining the hypocentre of an earthquake. Herald of Daghestan State Technical University. Technical Sciences. 2016;43(4):73-84. D0I:10.21822/2073-6185-2016-43-4-73-84 (In Russ.)]

2. Aslanov T.G. Opredelenie koordinat ochaga zemletryaseniya s ispol'zovaniem kombinirovannogo metoda. Vestnik Dagestanskogo gosudarstvennogo tekhnicheskogo universiteta. Tekhnicheskie nauki. 2017;44(2):118-125. DOI: 10.21822/2073-6185-2017-44-2-118-125 [Aslanov T.G. Definition of earthquake focus coordinates using a combined method. Herald of Daghestan State Technical University. Technical Sciences. 2017;44(2):118-125. DOI: 10.21822/2073-6185-2017-44-2-118-125 (In Russ.)]

3. Aslanov T.G., Daniyalov M.G., Magomedov Kh.D., Aslanov G.K. Ob odnom metode opredeleniya ochaga zemletryaseniya s odnovremennym opredeleniem skorostey seysmicheskikh voln. Trudy institute geologii Dage-

stanskogo nauchnogo tsentra RAN, Materialy. Makhachkala: Izdatelstvo DNTs RAN; 2010. 54-59. [Aslanov T.G., Daniyalov M.G., MagomedovKh.D., Aslanov G.K. About the method of seismic center determination with simultaneous definition of speeds of seismic waves. Proceedings of the Institute of Geology of the Dagestan Scientific Center of the Russian Academy of Sciences. Makhachkala: DSCRAS Publ.; 2010. 54-59. (In Russ.)]

4. Aslanov T.G., Davudova D.D., Mirzakhanova L.S. Opredelenie koordinat gipotsentra zemletryaseniya prii spol'zovanii kombinirovannogo metoda sfer i giperboloida: svid. o gos. reg. progr. EVM № 2017616673 Ros. Federatsiya / zayavitel' i pravoobladatel' VGUYu (RPA MinyustaRossii). № 2017611151; zayavl. 13.02.17; opubl. 09.07.17. [Aslanov T.G., Davudova D.D., Mirzakhanova L.S. Determination of the coordinates of the earthquake hypocenter using the combined method of spheres and hyperboloid: programme state registration EVM № 2017616673 Russian Federation / applicant and owner VGUYu № 2017611151; applied 13.02.17; publ. 09.07.17. (In Russ.)]

5. Shakhriman'yan M.A., Nigmetov G.M., Sosunov I.V. Matematicheskoe modelirovanie kak sposob podderzhki pri nyatiyaresheniy v sluchae vozniknoveniya chrezvychaynykh situatsiy. Katalog «Pozharnaya bezopasnost'»; 2003. S. 240-241. [Shakhriman'yan M.A., Nigmetov G.M., Sosunov I.V. Mathematical modeling as a method of desi-cion-making support in case of emergency situation arising. Catalogue «Firesafety»; 2003. P. 240-241. (inRuss.)]

6. Gitis V.G., Ermakov B.V. Osnovy prostranstvenno-vremennogo prognozirovaniya v geoinformatike. M.: FIZ-MATLIT; 2004. 256 s. [Gitis V.G., Ermakov B.V. Fundamentals of spatial-time prognosis in geoinformatics. Moscow: FIZMATLIT; 2004. 256 p. (In Russ.)]

7. Jing Z., Laurie G. B., Magaly K. Mapping earthquake induced liquefaction surface effects from the 2011 Tohoku earthquake using satellite imagery. 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). 2016. рр. 2328 - 2331.

8. Rui J., Shuanggen J. Ionosphericacousitc and rayleigh waves detected by GPS following the 2005 Mw=7.2 northern California earthquake. 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). 2016. рр. 3956 - 3959.

9. Aixia D., Xiaoqing W., Xiaoxiang Y., Shumin W. The loss assessment method of building earthquake damage using The Remote Sensing and building grid data. 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). 2016. рр. 4255 - 4258.

10. Wei Z., Huan-Feng Sh., Chun-Lin H., Wan-Sheng P. Building damage information investigation after earthquake using single post-event PolSAR image. 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). 2016. рр. 7338 - 7341.

11. Hao D., Xin X., Rong G., Chao S., Haigang S. Metric learning based collapsed building extraction from post-earthquake PolSAR imagery. 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). 2016. рр. 4742 - 4745.

12. Liu L. B., Liu M., Wang J.Q. Electromagnetic environment comprehension for radar detection of vital signs at China National Training Center for earthquake search &amp; rescue. 16th International Conference on Ground Penetrating Radar (GPR). 2016. рр. 1-4.

13. RadoutskiiV.Yu., VetrovaYu.V., Vasyutkina D.I. Opasnyeprirodnyeprotsessy: uchebnoeposobie (Pod red. V.Yu. Radoutskogo). Belgorod: Izd-vo BGTU; 2013. 206 s. [RadoutskiiV.Yu., VetrovaYu.V., Vasyutkina D.I. Dangerous natural processes: a manual (Edited by V.Yu. Radoutskii). Belgorod: Izd-vo BGTU; 2013. 206 p. (In Russ.)]

14. Arkhangel'skiy V.T., Vedenskaya I.A., Gayskiy V.N. Rukovodstvo po proizvodstvu i obrabotke nablyudeniyna seysmicheskikh stantsiyakh SSSR. Akad. nauk SSSR. Sovetposeysmologii. Moskva: Izd-voAkad. nauk SSSR; 1954. [Arkhangel'skiy V.T., Vedenskaya I.A., Gayskiy V.N. Observation performance and processing guide on the seismic stations of USSR. USSR Academy of Sciences. Seismology Board. Moscow: Izd-voAkad. nauk SSSR; 1954. (In Russ.)]

15. Dobrovol'skii I.P. Matematicheskaya teoriya podgotovki i prognoza tektonicheskogo zemletryaseniya. M.: FIZMATLIT; 2009. 240 s. [Dobrovol'skii I.P. Mathematical theory of preparation and forecast of a tectonic earthquake. M.: FIZMATLIT; 2009. 240 p. (In Russ.)]

16. Gurevich P.S. Psikhologiya chrezvychainykh situatsii. Uchebnoeposobie. M.: YuNITI-DANA; 2012. 495 s. [Gurevich P.S.Psychology of emergency situations. Tutorial. M.: YuNITI-DANA; 2012. 495 p. (in Russ.)]

17. Mkrtychev O.V. Bezopasnost' zdaniy i sooruzheniypriseysmicheskikh i avariynykhvozdeystviyakh. M.: MGSU; 2010. 152 s. [Mkrtychev O.V. Safety of buildings and constructions during seismic and emergency impacts. Moscow: MGSU; 2010. 152 p. (In Russ.)]

18. Aslanov T.G., TagirovKh.Yu., Aslanov G.K., AlimerdenovV.Sh. Matematicheskaya model' rascheta energeticheskogo klassa, intensivnosti i magnitudy zemletryaseniya v real'nom masshtabe vremeni. VestnikDage-stanskogo gosudarstvennogo tekhnicheskogo universiteta. Tekhnicheskie nauki. 2015;2(37):66-71. [Aslanov T.G., TagirovKh.Yu., Aslanov G.K., AlimerdenovV.Sh. Mathematical model for calculating the energy class, intensity and magnitude of an earthquake in real time. Herald of Daghestan State Technical University. Technical Sciences. 2015;2(37):66-71. (In Russ.)]

19. Aslanov T.G. Razrabotka algoritma opredeleniya koordinat ochaga zemletryaseniya, s odnovremennym opredele-niem skorostei seismicheskikh voln. Nauchnye trudy molodykh issledovatelei programmy «Shag v budushchee». Tom 8 «Professional». Moskva; 2005. P. 32-34. [Aslanov T.G. Development of an algorithm for determining the coordinates of the earthquake focus with simultaneous determination of seismic wave velocities. Proceedings of

young scientists of the «Step into the Future2» programme. Volume 8 «Professional». Moscow; 2005. рр. 32-34. (In Russ.)]

20. Rahinul H., Shoaib H., Akter S., Asadullahil G., Tahia F. K. Earthquake monitoring and warning system. Interna-tionalConferenceon Advances in Electrical Engineering (ICAEE). 2015. рр. 109 - 112.

21. Alphonsa A., Ravi G. Earthquake early warning system by IOT using Wireless sensor networks. International Conference on Wireless Communications, Signal Processing and Networking (WiSPNET). 2016. рр. 1201 - 1205.

22. Aslanov T.G., AlimerdenovV.Sh. Opredelenie struktury zemli po statisticheskim dannym vremen prikhoda seis-micheskikh voln proizoshedshikh zemletryasenii. Trudy III nauchno-tekhnicheskoikonferentsii mo-lodykhuchenykh i spetsialistov «Start v budushchee - 2013». Sankt-Peterburg; 2013. [Aslanov T.G., Alimerdenov V.Sh. Determination of the structure of the Earth according to the statistical data of the arrival times of seismic waves of occurredearthquakes. Proceedings of the III scientific and technical conference of young scientists and specialists «Start in the future - 2013». St. Petersburg; 2013. (In Russ.)]

23. Aslanov G.K., Gadzhiev M.M., Ismailov T.A., MagomedovKh.D. O zemletryaseniyakh (proshloe i sovremen-nost'). Makhachkala: Informatsionno-poligraficheskiitsentr DGTU; 2001. [Aslanov G.K., Gadzhiev M.M., Ismailov T.A., Magomedov Kh.D. About earthquakes (past and present). Makhachkala: Information and printing centre of DGTU; 2001. (In Russ.)]

Сведения об авторах:

Асланов Гайдарбек Тагирбекович - доктор технических наук, профессор, кафедра «Управление и информатика в технических системах и вычислительной техники».

Асланов Тагирбек Гайдарбекович - кандидат технических наук, докторант.

Курбанмагомедов Курбанмагомед Динмагомедович - кандидат технических наук, доцент, директор Института системных технологий.

Шахтарин Борис Ильич - доктор технических наук, профессор, кафедра «Автономные информационные и управляющие системы».

Information about the author.

Gaidarbek K.Aslanov - Dr.Sci. (Technical), Prof., Department «Management and Informatics in Technical Systems and Computer Engineering».

Tagirbek G. Aslanov - Cand. Sci. (Technical), Doctoral cand., Department «Management and Informatics in Technical Systems and Computer Engineering».

Kurbanmagomed D. Kurbanmagomedov - Cand. Sci. (Technical), Assoc. Prof., Director, Institute of System Technologies.

Boris I. Shakhtarin - Dr.Sci. (Technical), Prof., Department «Autonomous Information and Control Systems».

Конфликт интересов. Conflict of interest.

Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflict of interest.

Поступила в редакцию 28.09.2017. Received 28.09.2017.

Принята в печать 29.10.2017. Accepted for publication 29.10.2017.

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