Определение поверхности аномально напряженного включения в иерархической слоисто-блоковой среде по данным акустического мониторинга Текст научной статьи по специальности «Математика»

О.А. Хачай, А.Ю. Хачай

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

Геологическая среда является открытой системой, на которую действуют внешние и внутренние факторы. Они приводят ее к неустойчивому состоянию, которое, как правило, проявляется локально в виде зон, называемых динамически активными элементами, которые являются индикаторами потенциальных катастрофических источников. Эти объекты отличаются от вмещающей геологической среды своими структурными формами, которые часто являются формами иерархического типа. Процесс их активизации может наблюдаться с помощью мониторинга волновых полей, для математического обеспечения которого разработаны новые алгоритмы моделирования с использованием метода интегральных и интегро-дифференциальных уравнений. Разработан новый подход к интерпретации волновых полей, для определения контуров или поверхностей локально напряженных иерархических объектов. Разработан итерационный процесс решения теоретической обратной задачи для случая определения конфигураций 2D иерархических включений k-го ранга. При интерпретации результатов мониторинга необходимо использовать данные таких систем наблюдения, которые настроены на исследование иерархической структуры среды.

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

Введение

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

ISSN 0236-1493. Горный информационно-аналитический бюллетень. 2016. № 4. С. 354-365. © 2016. О.А. Хачай, А.Ю. Хачай.

УДК 550.830, 539.3

заключение о фундаментальной роли блочно-иерархического строения горных пород и массивов для объяснения существования широкой гаммы нелинейных геомеханических эффектов и возникновения сложных самоорганизующихся геосистем. Иерархическая структура характерна для многих систем, особенно для литосферы Земли, где было выделено по геофизическим исследованиям более 30 иерархических уровней от тектонических плит протяженностью в тысячи километров до отдельных минеральных зерен миллиметрового размера [1]. Таким образом, земная кора представляет собой сплошную среду, включающую в себя дискретную систему блоков и, как любой синергетиче-ский дискретный ансамбль, обладает свойствами иерархичности и самоподобия [2].

При построении математической модели реального объекта необходимо в качестве априорной информации использовать данные активного и пассивного мониторинга, получаемые в ходе текущей эксплуатации объекта. В работах [3, 4] построены алгоритмы моделирования в электромагнитном случае для 3-D неоднородности, в сейсмическом случае для 2-D неоднородности для произвольного типа источника возбуждения ^-слойной среды с иерархическим упругим включением, расположенным в /-ом слое.

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

В настоящей работе, используя подход, изложенный в работах [6, 7], разработан алгоритм получения уравнения теоретической обратной задачи для акустического поля (поперечной

акустической волны) для модели упругая иерархическая неоднородность k-го ранга, плотность которой совпадает с плотностью вмещающей среды для всех иерархических рангов, расположенная в v-ом слое упругого я-слойного полупространства.

Статья организована следующим образом. В разделе 1 изложен итерационный алгоритм решения обратной задачи для 2-D дифракции линейно поляризованной упругой поперечной волны для N-слойной упругой среды с упругим включением иерархического типа на основе уравнения теоретической обратной задачи. В разделе 2 изложены результаты и выводы.

Алгоритм решения обратной задачи 2-D дифракции линейно поляризованной упругой поперечной волны для N-слойной среды с иерархическим упругим включением

Пусть односвязная область D из евклидова пространства R2, ограниченная непрерывно дифференцируемой замкнутой кривой 3D, расположена в v-ом слое я-слойного полупространства. Пусть эта область содержит в себе K несоосных односвязных иерархических включений, ограниченных непрерывно дифференцируемыми замкнутыми кривыми dDk и простирающиеся параллельно оси ОХ. Границы li слоев П (j = 1,..., я) параллельны оси OY плоскости XOY декартовой системы координат. Ось OZ направлена вертикально вниз. Поместим начало координат на верхнюю границу поверхности 1-го слоя и совместим его с точкой, являющейся проекцией на OY точки, относительно которой область D звездная. Пусть U(y, z) —комплекснозначные дважды непрерывно дифференцируемые функции, удовлетворяющие соответственно двумерному скалярному уравнению Гельмгольца:

AU + c(M)U = - f (M) (1),

+ —•

ду2 дг2 '

c(M> = Ь; M e П \ D (' = °-""> (2).

К»;M e D,,(k = 1....Ю

Пусть функция U1(y, z) удовлетворяет уравнению:

AU1 + p(M)U = - f (M) (3)

p(M) J*M e nj \ D = (4)

lcv; M e Dk (k = 1,...K)

где

A =

Рассмотрим сначала случай при k = 1.

Для М е R2 \ D (j = 0,...п) определим:

и + (М) = Ц(М) - и\М) (5)

Функция и + (М) удовлетворяет уравнению (1). На границах раздела I слоев П. выполняются граничные условия:

]

и, =иу+1;

и + = и ++1; М е у(]■ = 1,...,п - 1); (6)

Ь1

Ц = Ц+1;

дЦу = Ь дЦУ+1; ь ди++ = ь дц++1

дп 1+1 дп 1 дп 1+1 дп

Ь дЦ = Ь 1 дп 1+1 ди]+1 ; дп М е ¡у

(7)

Ь. комплексные коэффициенты (/ = 0,...,п) и в общем случае: Ьу- Ф Ьу+1; на контуре дDk: при k = 1

иу= и+++ и1у; (8)

Функция удовлетворяет уравнению:

АЦ,+ с (М)ич=- f (М) (9)

и; — функция и+ в слое П ч D; и1 — функция [Г1 в слое D;

В области D при k = 1

иа = и++ и1а; М е D; Аиа + саиа = 0; (10)

Граничные условия на ^ = 1):

и+= и+ к дЦ- (дЩ + дЦ!) = 0 (Ш

дп дп дп

При М^ю функции Ц(М),Ц + (М),и\М)

удовлетворяют

условию излучения [8]. Алгоритм вычисления функции и1, например, приведен в работе [5].

Введем функцию О(М, М0), удовлетворяющую уравнению

АО + р(М)О = -д(М, М0) (12)

и граничным условиям (6), (7), при функция удовлетворяют условию излучения [8], при функция имеет особенность типа 1п1/ р(М, Мо):

р(м, М0) = у! (у - Уо)2 + (г - 20)2 (13)

Алгоритм вычисления функции для случая, когда область Б находится в v-ом слое приведен в работе [5]. Введем функцию О1, которая совпадает с фундаментальным решением уравнения (2) при k = 1. Применим формулу Грина [8] для пары функций П+ ^;(М е R2 \ О,М0 е г) в каждом из слоев П](j = 0,...,п). Проделаем процедуру аналогично [9]: домножим полученные выражения для каждого слоя на соответственно и сложим их почленно с учетом (1)—(4), (6)—(7). В результате получим:

2пи + (Ма) = -(% ) 1 (иV (М)

дG(M, Мо)

' ' дD дп (14)

ди+

-G(М,Мо)ц);М е Пу; Мо е П; дп

В области Б применим формулу Грина для пары функций Па(М), Ga(М, М0) . В результате получим:

0 = Г (Па(М) (M, М0) - Ga(М, М0)ПШ- (15)

дп дп

дО

Сложим выражения (14) и (15), учитывая (10)—(11), а также соотношение [5]:

0 = (- V) \ ШЦм)

двШ, Мо)

' • еп дп

ди1 —

-С(М,М0)-dl; М е п;М0 е Щ

дп

Тогда получим:

(дОа (М, М0)

(16)

2пП + (М0) = Г (П+ (М) + ПУ(М))(-

дп

-&/) ^ММ1) - ь.(П + дп)(( К С(М. м„) - <17>

/ Ь Рт Рт Рт / Ь

дп дп дп

-(уЬ! )С(М, М,))Ш,

Уравнение (17) является явным уравнением теоретической обратной задачи для двумерного скалярного уравнения Гель-мгольца в слоистой среде с однородным включением при заданных значениях граничных условий [5—7]. В результате реше-

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

Согласно [3, 10] задача о дифракции линейно поляризованной упругой поперечной волны на двумерной упругой неоднородности иерархического типа, расположенной в слое я-слойной среды в рамках описанной модели сводится к решению аналогичной задачи с учетом нижеследующих изменений. Уравнение теоретической обратной задачи (17) для скалярного уравнения Гельмгольца, к которому сводится наша задача, остается в силе:

2пи+ (м0) = | ((и; (м) + иЦмж

дD

дС(М,М0\ , /ди*. ди

дСа (М, М0)

дп

) - к (^.и )((%)оа (М, Мо) -(18)

Рт Рт / к

дп дп дп

-(Д )С(М, М,))Ш,

При этом

; К =^ Ьа =^а ; (19)

, ^I, ^а , р , р1, р а — значения упругого параметра Ламе и плотность в -ом слое, в слое, где находится точка М0 и внутри неоднородности при k = 1. Важным отличием настоящей задачи от рассмотренной ранее является то, что ра = ру при всех k. Физически это означает, что аномалия в акустическом поле создается аномалией напряженного состояния среды и может быть связана с очагом либо горного удара, либо землетрясений.

и+= к и:= п1 ;и = п1; (20),

где их — составляющая вектора смещения, отличная от нуля для выбранной модели.

G(M, М0) = GSS(М,М0); Ga(М, М) = GaaS(М, М0); дБ, dl, Я =*> ; %,=а>2

2 „2 Р^ ,2_ 2 Р^ (21)

Алгоритм вычисления функции Грина (М, М0) выписан в работе [11].

Таким образом, уравнение теоретической обратной задачи при k = 1 записывается в виде:

2nu+(M.) = J ((<(M) + uL(M))(dGSS(MM -

j dn

öDi un

-Ы)Xss(M'M°Y)■ (dU++dUnr)((X »zw,m.)-(22)

/St dn dn dn / Sai

-(K)Gss (M, M.)))dl;

Пусть k = 2, т.е. искомый объект представляет собой иерархическое включение: с упругим параметром Ламе ^ и плотностью pv, dD1 — контур внешнего включения и с упругим параметром Ламе плотностью pv, dD2 — контур внутреннего включения. Включения несоосные. Требуется восстановить оба контура. Для решения нашей задачи в выражении (19) Sa = Ss , в выражении (21):

dD = 6DX; dl = dlv G¡s = Gal; ^G) = -Щ-(Gfs);

dn dn

ka = kla\ = Ю

S al

Решив уравнение (22) относительно функции г1(ф), описывающей контур dD1, вычисляем функции: ux; u+; Ux по алгоритму решения прямой задачи [4] внутри и вне неоднородности,

s

помещенной в слоистую среду, ux — компонента упругого поля в слоистой среде в отсутствии неоднородности.

',2 U2

ux (M0) = ^ uX (M0) + ^^ Ц ^ (M)GSS(M, M0)dS1

*>„ 2n Si

Г /,^ dG

+ Ux (M)^dl; M0 eSi

C(M0) = uX(M0) + k)) JJ ux(M)Gss(M, M0)dSi

(23)

2n^(M0) si

(¡^v ) f fu

(24)

J ux(M)G-dl;Mo ¿Si

•> rm

2n^(M.) D x dn

На этом первый итерационный цикл заканчивается, и мы переходим ко второму итерационному циклу k = 2. Вычисленную функцию ux(M0) (24) обозначаем как:

<к-1); (25),

а'2. д / г^а \ д , а'2\. 1,2 1,2 "2 Ру

в выражении (19) Ъ1а = , в выражении (21):

дО = 3О2; dl = dl2,Gass = Gfs;-^ &55) =&£); = к^ = ®2 ^

дп дп с,„,

Уравнение (22) переписывается в виде:

2 пи+ (М0) = Г («(М) + иЦк-1)(М))(-

С (М, Мй)

дО2 дп

)ХШМ*1)(^^ж^ С(М.М)-

/ Ы дп дп дп / Ъа2

-(Яг )^(М, М,))Ш2, (26)

Решаем уравнение (26) относительно функции г2(ф), описывающей контур ЗБ2.

Если К = 2, то задача считается решенной, если k > 2, k = = k + 1, итерационный процесс продолжается. Мы вычисляем функции: и-1 „+(к-ц

их , их , (27)

по алгоритму решения прямой задачи внутри и вне иерархической неоднородности ранга ^1, помещенной в слоистую среду (физические параметры слоистой среды остаются неизменными) [4]. ?

Ъа( к-1)

■ка(к-1) - Я икх-1(М^ (М, (28)

2П

I ик-1(М) дСSdlk М0 е S(k_l)

2П ^а(к-1) дО(к-1)1 дп

ик-1(М0)=ит-2)шй)+) я ик-1(Мю^(М,М0^

2п №0) s(к-l)

+ (*>у-£>а) Г ик-1(М) , М ^ (29)

1 I их (М ) Я и1к-1, М0 £ к-1)

Вычисленную функцию и (М0) (29) обозначаем как:

их

0

1( -1) (30),

в выражении (19) 5а = Ск , в выражении (21):

дБ = дБк; й\ = й1к;О^ = О£(О|5) = А (о£); к22а = ^ = ш2 А-.

дп дп

Уравнение (22) переписывается в виде:

,дОа55 (М, м0)

Ъак

2пп+х(Мо) = | ((и+„(И) + и^-1)(М))(-

дБ

дп

ОИ Мо)) - 5.(£ + ^Ж ус )О^(м, Мо) -

дп дп дп / ъ а

-(/с,(М, Мо))Ы; (31)

Решаем уравнение (26) относительно функции гк(ф), описывающей контур дБ • k = k + 1. Итерационный процесс (27)—(31) продолжается до k = К.

Заключение

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

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

1. Прангишвили И. В., Пащенко Ф. Ф., Бусыгин Б. П. Системные законы и закономерности в электродинамике, природе и обществе. — М.: Наука, 2001. - 525 с.

2. Кочарян Г.Г., Спивак А. А. Динамика деформирования блочных массивов. — М.: ИКЦ «Академкнига», 2003. — 424 с.

3. Хачай О. А., Хачай А.Ю. О комплексировании сейсмических и электро-магнитных активных методов для картирования и мониторинга состояния двумерных неоднородностей в N-слойной среде // Вестник ЮУРГУ. Серия: Компьютерные технологии, управление, радиоэлектроника. — 2011. — № 2(219). — С. 49—56.

4. Хачай О. А., Хачай А.Ю. Моделирование электромагнитного и сейсмического поля в иерархически неоднородных средах // Вестник ЮУРГУ. Серия: Вычислительная математика и информатика. — 2013. — Т. 2. — № 2. — С. 48—55.

5. Хачай О. А. Математическое моделирование и интерпретация переменного электро-магнитного поля в неоднородной коре и верхней мантии Земли. Докторская диссертация. — Свердловск: ИГФ УрО РАН, 1994. — 314 с.

6. Хачай О. А. Об интерпретации двумерных переменных и трехмерных стационарных аномалий электромагнитного поля // Известия АН СССР. Физика Земли. — 1989. — № 10. — С. 50—58.

7. Хачай О. А. О решении обратной задачи для трехмерных переменных электромагнитных полей // Известия АН СССР. Физика Земли. — 1990. — № 2. — С. 55—59.

8. СтрэттонДж. Теория электромагнетизма. — M.-Л.: ОГИЗ, 1948. — 409 c.

9. Дмитриев В. И. Дифракция плоского электромагнитного поля на цилиндрических телах, расположенных в слоистых средах / Вычислительные методы и программирование. — М.: Изд-во МГУ, 1965, вып. 3. — C. 307—315.

10. Купрадзе В. Д. Граничные задачи теории колебаний и интегральные уравнения. — М.-Л.: Гос. Изд-во технико-теоретической литературы, 1950. — 280 с.

11. Хачай А. Ю. Алгоритмы для математического моделирования переменных электромагнитных и сейсмических полей в источниковом приближении. Кандидатская диссертация. — Екатеринбург: УрГУ им. А.М. Горького, 2007. —176 с. ЕИЗ

КОРОТКО ОБ АВТОРАХ

Хачай Ольга Александровна — доктор физико-математических наук, ведущий научный сотрудник, Институт геофизики УрО РАН, e-mail: olgakhachay@yandex.ru,

Хачай Андрей Юрьевич — кандидат физико-математических наук, доцент, Уральский федеральный университет, e-mail: andrey.khachay@gmail.com.

Gornyy informatsionno-analiticheskiy byulleten'. 2016. No. 4, pp. 354-365. O.A. Khachay, A.Yu. Khachay

DEFINING OF THE 2-D SURFACE OF THE ANOMALY STRESSED HIERARCHICAL OBJECT LOCATED IN THE LAYERED BLOCKED GEOLOGICAL MEDIUM USING THE DATA OF ACOUSTIC MONITORING

Geological medium is an open system which is influenced by outer and inner factors that can lead it to a unstable state. That non stability is as a rule occurred locally and these zones are named as dynamically active elements, which are indicators of potential catastrophic sources. These objects differ from the embedded geological medium by their structural forms, which often are of hierarchical type. The process of their activisation can be searched, using wave fields monitoring. For that purpose we had developed earlier new algorithms of modeling wave field propagation through the local objects with hierarchical structure. Here we had developed a new approach for interpretation the distribution of wave fields for defining the contours of these local hierarchical objects. We developed an algorithm for constructing the equation of theoretical inverse problem for 2-D linear polarized longitudinal elastic wave by excitation of the N-layered elastic medium with hierarchic elastic inclusion located in the v -th layer with the density equal to the density of the embedded layer. An iteration process of solving the inverse problem for the case of certain configurations of hierarchical inclusions 2D k-th rank is elaborated. When interpreting the results of the monitoring need to use the data of such systems that are configured to study the hierarchical structure of the medium.

Key words: hierarchic medium, seismic field, algorithms of modeling, equation of theoretical inverse problem.

AUTHORS

Khachay O.A., Doctor of Physical and Mathematical Sciences, Leading Researcher,

Institute of Geophysics of Ural Branch of Russian Academy of Sciences,

620016, Ekaterinburg, Russia, e-mail: olgakhachay@yandex.ru,

Khachay A.Yu., Candidate of Physical and Mathematical Sciences, Assistant Professor,

Ural Federal University, 620002, Ekaterinburg, Russia, e-mail: andrey.khachay@gmail.com.

REFERENCES

1. Prangishvili I. V., Pashchenko F. F., Busygin B. P. Sistemnyezakony izakonomernosti v elektrodinamike, prirode i obshchestve (System laws and rules in electrodynamics, nature and society), Moscow, Nauka, 2001, 525 p.

2. Kocharyan G. G., Spivak A. A. Dinamika deformirovaniya blochnykh massivov Dynamics of deformation of block arrays), Moscow, Akademkniga, 2003, 424 p.

3. Khachay O. A., Khachay A. Yu. Vestnik Yuzhno-Ural'skogo gosudarstvennogo uni-versiteta. Seriya: Komp'yuternye tekhnologii, upravlenie, radioelektronika. 2011, no 2(219), pp. 49-56.

4. Khachay O. A., Khachay A. Yu. Vestnik Yuzhno-Ural'skogo gosudarstvennogo univer-siteta. Seriya: Vychislitel'naya matematika i informatika. 2013, vol. 2, no 2, pp. 48-55.

5. Khachay O. A. Matematicheskoe modelirovanie i interpretatsiyaperemennogo elektro-magnitnogo polya v neodnorodnoy kore i verkhney mantii Zemli (Mathematical modeling and interpretation of alternating electromagnetic field in heterogeneous crust and mantle of the Earth. Dissertation of doctor of physics and mathematics), Doctor's thesis, Sverdlovsk, IGF UrO RAN, 1994, 314 p.

6. Khachay O. A. Izvestiya ANSSSR. Fizika Zemli. 1989, no 10, pp. 50-58.

UDC 550.830, 539.3

7. Khachay O. A. Izvestiya ANSSSR. Fizika Zemli. 1990, no 2, pp. 55-59.

8. Stretton Dzh. Teoriya elektromagnetizma (Theory of electromagnetism), Moscow-Leningrad, OGIZ, 1948, 409 p.

9. Dmitriev V. I. Vychislitel'nye metody i programmirovanie (Computational methods and programming, vyp. 3, issue 3), Moscow, Izd-vo MGU, 1965, pp. 307-315.

10. Kupradze V. D. Granichnye zadachi teorii kolebaniy i integral'nye uravneniya (Boundary problems of the theory of oscillations and integral equations), Moscow-Leningrad, Gos. izd-vo tekhniko-teoreticheskoy literatury, 1950, 280 p.

11. Khachay A. Yu. Algoritmy dlya matematicheskogo modelirovaniya peremennykh ele-ktromagnitnykh i seysmicheskikh poley v istochnikovom priblizhenii (Algorithms for mathematical modeling of alternating electromagnetic and seismic fields), Candidate's thesis, Ekaterinburg, UrGU im. A.M. Gor'kogo, 2007, 176 p.

ОТДЕЛЬНЫЕ СТАТЬИ

ГОРНОГО ИНФОРМАЦИОННО-АНАЛИТИЧЕСКОГО БЮЛЛЕТЕНЯ

(СПЕЦИАЛЬНЫЙ ВЫПУСК)

АНАЛИЗ БЕЗОПАСНОСТИ ОПАСНЫХ ПРОИЗВОДСТВЕННЫХ ОБЪЕКТОВ И ВОПРОСЫ ЭФФЕКТИВНОСТИ ЭКСПЕРТИЗЫ ПРОМЫШЛЕННОЙ БЕЗОПАСНОСТИ

Гаязов Рустам Рашитович — генеральный директор, Знаемский Юрий Анатольевич — начальник лаборатории, Кадысев Владимир Степанович — руководитель НТО, Маринин Илья Алексеевич — начальник отдела экспертизы,

000 «ЮТАР».

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

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

SAFETY ANALYSIS OF HAZARDOUS PRODUCTION FACILITIES AND THE EFFECTIVENESS OF THE INDUSTRIAL SAFETY EXPERTISE Gayazov R.R.1, General Director, Znaemskiy Yu.A.1, Head of Laboratory, Kadyshev V.S.1, Head of NTO, Marinin I.A.1, Chief of Department of Examination,

1 LLC «Utar», Russia.

Statistical information collected by different industries, confirms the provision on the impact of the status of the exploited technology to the venture background of the company, as well as the importance of joint control of industrial safety of an enterprise-owner, expert organizations and the state. In modern conditions, when the main control functions of equipment to perform enterprise together with the expert organizations are of particular importance in the examination of industrial safety.

Key words: expert examination of industrial safety (industrial safety), quality management, redundant reliability, the dangerous, the equipment operating under excess pressure, lifting facilities, gas distribution facilities and gas consumption, petrochemical industry, oil refining, explosion, fire, emission of hazardous substances, depressurization.