ГИАБ. Горный информационно-аналитический бюллетень / MIAB. Mining Informational and Analytical Bulletin, 2020;(8):40-51 ОРИГИНАЛЬНАЯ СТАТЬЯ / ORIGINAL PAPER
УДК 622.235 DOI: 10.25018/0236-1493-2020-8-0-40-51
АВТОМАТИЗИРОВАННОЕ ОПРЕДЕЛЕНИЕ КООРДИНАТНОЙ СЕТКИ ВЗОРВАННОГО БЛОКА МАССИВА ПОРОД
Б.Р. Ракишев1, З.Б. Ракишева2, А.А. Орынбай1
1 Казахский национальный исследовательский технический университет имени К.И. Сатпаева (Satbayev Unuversity), Алматы, Казахстан, e-mail: b.rakishev@mail.ru 2 Казахский национальный университет имени аль-Фараби, Алматы, Казахстан
Аннотация: Даны определения понятий «координатная сетка взрываемого блока», «координатная сетка взорванного блока уступа». Описан автоматизированный метод определения узловых точек координатной сетки взорванного блока. Он базируется на учете зависимостей, связывающих начальные параметры взрываемого блока массива пород, такие, как его высота и сетка скважин с конечным расположением зафиксированных точек взорванного блока массива пород. Определяющими факторами координатной сетки взорванного блока массива пород являются характерные высота и ширина развала, которые зависят от физико-механических свойств взрываемых пород, физико-химических характеристик применяемого ВВ и параметров расположения зарядов в массиве пород. Значительное влияние на положение узловых точек координатной сетки взорванного блока оказывает коэффициент разрыхления пород по горизонтальным и наклонным слоям. Выявлена тенденция увеличения этого показателя от тыловой части развала к его периферии как в горизонтальном направлении, так и по вертикали. На основе совместного учета всех факторов, влияющих на формирование развала пород, впервые разработан аналитический метод определения узловых точек и самой координатной сетки взорванного блока. Создана компьютерная программа автоматизированного определения координатной сетки взорванного блока, учитывающая влияние физико-механических свойств взрываемых пород, физико-химических характеристик применяемого ВВ и параметров расположения зарядов в массиве пород на координаты узловых точек искомой сетки. С помощью этой программы оперативно и достаточно точно могут быть найдены положения узловых точек координатной сетки взорванного блока при различных параметрах и условиях взрывания массивов пород. Установленная таким образом координатная сетка взорванных блоков служит инструментом для определения местоположения различных элементов массива пород в развале взорванных пород.
Ключевые слова: координатная сетка взорванного блока, узловые точки, высота и ширина развала пород, горизонтальные и наклонные слои пород, коэффициент разрыхления пород, автоматизированное определение координатной сетки.
Благодарность: Работа выполнена в рамках целевой программы 2018/BR05235618 «Модернизация технологий и производств в горнодобывающей и горноперерабатывающей отраслях Республики Казахстан» при финансовой поддержке Министерства науки и образовании РК в 2018-2020 гг.
Для цитирования: Ракишев Б.Р., Ракишева З.Б., Орынбай А.А. Автоматизированное определение координатной сетки взорванного блока массива пород // Горный информационно-аналитический бюллетень. - 2020. - №> 8. - С. 40-51. DOI: 10.25018/0236-1493-20208-0-40-51.
© Б.Р. Ракишев, З.Б. Ракишева, А.А. Орынбай. 2020.
Computer-aided creation of coordinate grid for blasted rock block
B.R. Rakishev1, Z.B. Rakisheva2, A.A. Orynbay1
1 K.I. Satpayev Kazakh National Research Technical University (Satbayev Unuversity), Almaty, Kazakhstan, e-mail: b.rakishev@mail.ru 2 Al-Farabi Kazakh National University, Almaty, Kazakhstan
Abstract: The notions of 'coordinate grid of blast block' and 'coordinate grid of blasted block' are defined. The computer-aided determination of the coordinate grid points for a blasted rock block is described. The starting point is using relationships of the initial parameters of a blast rock block such as height and the final pattern of the determined points in the blasted rock block. The governing factors of a blasted block coordinate grid are the broken-rock disintegration height and width dependent on physical and mechanical properties of rocks, physicochem-ical characteristics of explosives and on blasting pattern design. Considerable influence on arrangement of the blasted block coordinate grid points is exerted by the degree of fragmentation of rocks in horizontal and inclined layers. It is found that the degree of fragmentation tends to increase from the rear to the periphery of the broken rock pile, both horizontally and vertically. Based on the integration of all determinants of a broken rock pile after blasting, the analytical method is newly developed to determine the points and create the coordinate grid of a blasted block. The new computer program for the automatic creation of a blasted block coordinate grid takes into account the influences exerted by physical and mechanical properties of rocks, phys-icochemical characteristics of explosives and blasting pattern design. The program allows fast and sufficiently accurate position fixing of a blasted block coordinate grid points at different parameters and conditions of blasting. The resultant coordinate grids of blasted blocks are the tool for locating different elements inside broken-rock disintegration.
Key words: blasted block coordinate grid, grid points, broken-rock disintegration height and width, horizontal and inclined rock layers, degree of disintegration, computer-aided creation of coordinate grid.
Acknowledgements: The study was supported by the Ministry of Education and Science of the Republic of Kazakhstan in 2018-2020 within the framework of Target Program 2018/ BR05235618: Technological and Production Modernization in the Mining and Processing Industry in the Republic of Kazakhstan.
For citation: Rakishev B. R., Rakisheva Z. B., Orynbay A. A. Computer-aided creation of coordinate grid for blasted rock block. MIAB. Mining Inf. Anal. Bull. 2020;(8):40-51. [In Russ]. DOI: 10.25018/0236-1493-2020-8-0-40-51.
Введение
Прогнозирование размещения разнородных пород уступа в развале взорванных пород имеет исключительно важное значение в горном деле. Оно актуально при отработке сложно-структурных скальных руд, при отгрузке пород различной кусковатости из различных зон развала. Установление внутренней
структуры взорванных пород позволит обоснованно выбрать рациональные технологии буровзрывных, выемочно-погрузочных работ и повысит эффективность открытой разработки многокомпонентных, сложно-структурных месторождений.
Формирование развала взорванных пород подчиняется закономерностям,
вытекающим из модели поэтапного разрушения массива горных пород взрывом зарядов взрывчатого вещества (ВВ) [1, 2]. Согласно этой модели, основное разрушение скальных горных пород происходит в первой стадии взрыва под действием волн напряжений, возбуждаемых в среде взрывом заряда ВВ, и самих газообразных продуктов взрыва.
На втором этапе (во второй стадии) взрыва вследствие влияния свободной поверхности нарушается осесимметрич-ное развитие взрывной полости и зон разрушения. Оставшиеся во взрывной полости продукты взрыва доразрушают отбиваемую часть массива и сообщают ей механическое движение в сторону свободной поверхности. Вся раздробленная, выброшенная продуктами взрыва горная масса в пространстве занимает положение, определяемое конкретным значением начальной скорости отбиваемой части уступа. Она зависит от физико-механических свойств пород взрываемого блока, физико-химических характеристик применяемого ВВ и параметров расположения зарядов ВВ во взрываемом массиве.
На третьем этапе (в третьей стадии) взрыва совершается оседание выброшенных пород раздробленной части уступа в поле силы тяжести и образуется развал взорванной горной массы. Фиксированные элементы взрываемого блока в развале занимают положения, предопределяемые исходными данными конкретного массового взрыва. Такой подход к механизму разрушения горных пород взрывом находит применение в работах ученых ближнего [3-13] и дальнего зарубежья [14 — 20].
Формирование развала горных пород подробно рассмотрены в работах [2, 3 ]. Для аналитического определения местоположения отдельных элементов взрываемого блока уступа в развале введены понятия «координатная сетка
взрываемого блока», «координатная сетка взорванного блока». Разработан аналитический метод определения узловых точек и самой координатной сетки взорванного блока.
Геометрические характеристики
элементов взорванного блока
уступа
Для создания координатных сеток взрываемый блок целесообразно разбить на определенное число частей виртуальными плоскостями, проведенными параллельно откосу уступа и горизонту. Например, при однорядном расположении скважин взрываемый блок указанными плоскостями достаточно разбить на 16 равных частей (элементов) (см. рис. 1, а). При двухрядном расположении скважин взрываемый блок разбивается на 24 части (элемента), т.е. к рассмотренному случаю добавляются два наклонных равновеликих слоя толщиной, равной половине расстояния между рядами скважин и т.д. (см. рис. 2, а).
Следы горизонтальных разбивающих плоскостей в разрезе уступа представляют собой горизонтальные окон-туривающие линии элементов уступа. В их число входят также следы кровли и подошвы уступа. Аналогично следы наклонных разбивающих плоскостей представляют собой наклонные окон-туривающие линии рассматриваемых элементов уступа. Они включают в свой состав линию откоса уступа и линию отрыва взрываемого блока от массива пород (см. рис. 1, а, 2, а).
Ранее проведенными исследованиями [2, 3] выявлены закономерности размещения фиксированных элементов взрываемого блока в развале и установлены положения оконтуривающих линий при различных способах взрывания. Из их анализа следует, что размещение рассматриваемых элементов уступа в развале полностью определя-
Рис. 1. Схема размещения элементов уступа во взрываемом (а) и взорванном блоках (б) при однорядном взрывании
Fig. 1. Arrangement of elements of bench in blast block (a) and blasted block (b) in blasting in single row of boreholes
ются положением их оконтуривающих линий во взорванной горной массе.
Совокупность взаимно пересекающихся горизонтальных и наклонных оконтуривающих линий элементов взрываемого блока в разрезе названа координатной сеткой взрываемого блока (см. рис. 1, а, 2, а). Совокупность взаимно пересекающихся деформированных го-
ризонтальных и наклонных оконтуривающих линий элементов уступа в развале названа координатной сеткой взорванного блока (см. рис. 1, б, 2, б, 2, б).
К геометрическим характеристикам элементов взорванного блока относятся их площади, оконтуривающие линии, коэффициенты разрыхления элементов по горизонтальным и наклонным слоям
VTT'W
TVYv\
б)
J. V-tt:
VT \
n \ \ \ \ 1
Рис. 2. Схема размещения элементов уступа во взрываемом (а) и взорванных блоках при четырехрядном КВЗ на свободную поверхность (б) и с подпорной стенкой (в)
Fig. 2. Arrangement of elements of bench in blast block (a) and blasted block in short-delay blasting in four rows of boreholes with free surface (b) and with retention wall (c)
Площади и коэффициенты разрыхления пород в слоях взорванного блока при однорядном взрывании
Areas and degrees of rock fragmentation in layers of blasted blocks in blasting in single row
S" 11 S" 12 S" 13 S" 14 S" 21 S" 22 S" 23 S" 24
10,59 21,33 32,15 44,85 25,03 50,47 76,94 107,68
S" 31 S" 32 S" 33 S" 34 S" 41 S" 42 S" 43 S" 44
39,75 81,82 125,14 175,04 54,72 113,83 174,96 244,18
k",, oil k ".2 k" 013 k" 14 014 k" k 021 k",, k" k "24 024
1,13 1,14 1,14 1,20 1,33 1,35 1,37 1,44
k" k 031 k" 32 k" 33 k "34 d34 k" K 041 k d42 k "43 043 k "44 044
1,41 1,45 1,48 1,56 1,46 1,52 1,56 1,63
и координаты узловых точек координатной сетки взорванного блока. Они представлены на рис. 1, б, 2, б, 2, в.
На основе работ [2, 3], ниже приведены примеры расчетного определения указанных характеристик элементов взорванного блока для диорит-порфиритов
Сарбайского карьера. Характеристики пород: р0 = 2830 кг/м3, с = 5100 м/с, а = 190 МПа, а = 16 МПа, V = 0,24.
сж р
Параметры расположения зарядов: высота уступа Ь = 15,0 м, линия сопротивления по подошве уступа № = 10,0 м, расстояние между скважинами а = 7,5 м,
S" J ii S" J 12 S" J 13 S" J 14 S" J 15 S" J 16 S "17 S" J 18 S "19 S" J 1,10
14,54 29,86 45,41 61,10 77,03 93,74 106,13 118,71 131,58 145,07
S" 21 S" 22 S" 23 S" 24 S" 25 S" 26 S "27 S" 28 S "29 S" 2,10
29,96 60,96 93,18 126,68 161,39 198,40 225,37 253,17 282,14 312,78
S" J 31 S" J 32 S" J 33 S" J 34 S" J 35 S" J 36 S "37 S" J 38 S "39 S" J 3,10
46,99 94,49 145,35 197,55 252,42 309,88 352,04 396,33 442,74 491,90
S" 41 S" 42 S" 43 S" 44 S" 45 S" 46 S "47 S" 48 S "49 S" 4,10
64,37 130,75 198,92 271,68 347,75 426,57 484,92 545,93 610,02 678,96
k",, oil k k" 013 k" ,A 014 k" d15 k" 016 k" d17 k" K o18 k 019 k 01,10
1,03 1,06 1,08 1,09 1,10 1,11 1,13 1,15 1,17 1,19
k" k 021 k",, 022 k" 023 k" ,4 024 k" „ 025 k" ,6 026 k" 27 027 k" K 028 k 029 k 02,10
1,07 1,08 1,10 1,13 1,15 1,18 1,20 1,23 1,25 1,28
k 031 k" 032 k 033 k" 34 034 k" 035 k" 36 036 k "37 037 k" K 038 k 039 k 03,10
1,11 1,12 1,15 1,17 1,20 1,22 1,25 1,28 1,31 1,35
k' 041 k "n4, d42 k "43 043 k "p44 d44 k" 45 045 k "p46 d46 k "47 047 k" K 048 k 049 k" 04,10
1,14 1,16 1,18 1,21 1,24 1,26 1,29 1,32 1,36 1,39
Таблица 2
Площади и коэффициенты разрыхления пород в слоях взорванного блока при четырехрядном КВЗ на свободную поверхность Areas and degrees of rock fragmentation in layers of blasted blocks in short-delay blasting in four rows with free surface
Площади и коэффициенты разрыхления пород в слоях взорванного блока
при четырехрядном КЗВ с подпорной стенкой
Areas and degrees of rock fragmentation in layers of blasted blocks
in short-delay blasting in four rows with retention wall
S" 11 S" 12 S"i3 S" 14 S" 16 S" 17 S" 18 S" 19 S" 1,10
14,13 28,72 43,40 58,76 74,63 90,88 103,67 116,90 129,23 142,96
S* 21 S" 22 S"23 S" 24 S"25 S" 26 S" 27 S" 28 S" 29 S" 2,10
28,39 57,51 87,48 119,46 152,55 186,29 212,42 240,05 268,85 298,26
S" 31 S" J 32 S"33 S" J 34 S"35 S" J 36 S" J 37 S" J 38 S" J 39 S" J 3,10
42,94 87,43 134,03 184,29 236,26 289,55 329,74 372,81 418,68 465,82
S* 41 S" 42 S"43 S" 44 S45 S" 46 S" 47 S" 48 S" 49 S" 4,10
58,39 118,71 183,41 252,50 324,61 399,44 455,04 515,08 579,55 647,29
k'„ pii k"p12 k" p13 k" p14 k" p15 k" p16 k" 17 p17 k" K p18 k p19 k p1,10
1,00 1,02 1,03 1,04 1,06 1,08 1,11 1,13 1,15 1,17
k" k p21 k" 22 p22 k" 23 p23 k" 24 p24 k" 25 p25 k" 26 p26 k" 27 p27 k" K p28 k p29 k p2,10
1,01 1,02 1,04 1,06 1,08 1,10 1,13 1,16 1,19 1,22
k" k p31 k" 32 p32 k"33 p33 k" p34 k" p35 k" p36 k" 37 p37 k" K p38 k p39 k p3,10
1,02 1,04 1,06 1,09 1,12 1,14 1,17 1,21 1,24 1,27
k" p41 k"p42 p42 k" 43 p43 k" p44 k"^ p45 k" 46 p46 k"p47 p47 k" K p48 k p49 k" p4,10
1,04 1,06 1,09 1,12 1,15 1,18 1,21 1,25 1,29 1,33
между рядами скважин а = 7,5 м, длина перебура 1п = 2,0 м, длина забойки 13 = 7,0 м, масса скважинного заряда ф = 450 кг. Расположение скважин однорядное, четырехрядное. Взрывание короткозамедленное (КЗВ) на свободную плоскость и с подпорной стенкой. Ширина отброшенной части развала 22,0 м. Ширина подпорной стенки по верху Впс = 3,0 м, по низу Впс = 6,0 м. Ширина подпорной стенки по низу входит в состав ширины отброшенной части развала. Удельный расход ВВ при однорядном взрывании q = 0,4 кг/м3, четырехрядном взрывании q = 0,49 кг/м3.
Схемы размещения элементов взорванного уступа при различных условиях взрывания в масштабе 1:468,75 приведены на рис. 1, б, 2, б, 2, б. Площади и коэффициенты разрыхления пород в зонах взорванного блока вычислены
с использованием программы AutoCad. В данном случае слоями взорванного блока, соответствующими искомым точкам координатной сетки, являются зоны массива или развала породы, ограниченные откосом уступа в тыльной части, нижней площадкой уступа и оконтуривающими линиями, которым принадлежат искомые точки. Они сведены в табл. 1-3. Нумерация слоев по вертикали у (у = 1, 2...4) ведется снизу вверх, а по горизонтали 5 (5 = 1, 2.10) слева на право. Этому порядку соответствуют индексы рассматриваемых характеристик элементов взорванного
блока (S*1,1, S*1,2, ..., S*4,10; k *p1,1, kp1,2, ...,
k ). Площади приведены в м2.
Найденные значения S" , k" служат исходными данными для аналитического определения узловых точек координатной сетки взорванного блока уступа.
Аналитическое определение положения узловых точек координатной сетки взорванного блока
Как видно из рис. 1, а, 2, а, координаты узловых точек координатной сетки взорванного блока могут быть определены по следующим формулам:
х.. = кх. (С + т.), у.. = ку.. * /., (1)
I . . р X ' .. . . . ' * '
где х.., у.. — координаты узловых точек координатной сетки, т.е. точек пересечения /-й и у-й оконтуривающих линий элементов взорванного блока уступа, / (/ = 1, 2...5),у у = 1, 2...11); кх. ку. . — коэффициенты пропорциональности изменения координат узловых точек соответственно по горизонтальному и вертикальному направлениям; €х, т., /. — константы уравнения.
Константа С. одинакова для однорядного и многорядного расположения скважин и определяется по зависимости:
т
С = (Л — / .) • ctga.
(2)
Здесь индексы С. обозначают номера узловых точек на линии отрыва пород от массива снизу вверх (X = 1, 2.5), Л — высота уступа, а — угол откоса уступа, /. — ордината /-й окон-туривающей горизонтальной линии взрываемого блока.
/1 = 0; /2 = Л/4; /3 = Л/2;
/4 = 3Л/4; /5 = Л. (3)
Константа т. — расстояние (по горизонтали) между рассматриваемой окон-туривающей наклонной линией и линией отрыва пород от массива во взрываемом блоке.
При однорядном взрывании: т, = 0; т2 = №/4; т3 = №/2;
т4 = 3Щ/4; т5 = Щ (4)
при четырехрядном КЗВ:
0; т2 = ар /2; т3 = а р; т. = 3а /2;
т
т5 = 2ар; тб = 5ар/2; т7 = 3ар;
т
10
3ар + Щ/4; т9 = 3ар + Щ2; 3ар + 3Щ/4; ти = 3ар + Щ. (5)
Численные значения коэффициентов кх. ., ку.. найдены по эмпирическим значениям х.у.. при заданных Сх, т., /..
Аналитическое определение коэффициентов
пропорциональности изменения координат узловых точек координатной сетки взорванного блока
Для теоретического обобщения полученных экспериментальных результатов коэффициенты пропорциональности изменения координат узловых точек необходимо взаимоувязать с интегральной характеристикой состояния развала — с коэффициентом разрыхления пород в слоях развала, соответствующих точкам координатной сетки и с соотношением, учитывающим длину отброшенной части развала Во, константы /., т., т.е. с параметрами:
рх - к 0У - к
т!
где у = 1, 2, 3, 4, 5 = 1, 2.10, i = у + 1, у = 5 + 1.
По экспериментально установленным значениям были найдены аппроксимирующие полиномиальные уравнения для определения искомых коэффициентов пропорциональности точек координатной сетки развала пород [21]:
кх.. = а.х рх 55 + Ь.х рх 54
. ! ~ у5 . ~ у5
+ с.х рх 3
у5
+ ёх рх^2 + е.х рху5 + /х,
ку..
' РУ^2 + Ь.у РУу5 + с.у, (7)
где а х, Ь х, с.х, ёх, е х, /х — константы первого уравнения (7), связанные с /-й горизонтальной оконтуривающей линией; ау, Ьу, су — константы второго уравнения (.7), связанные с у-й наклон-
+
а
ной оконтуривающей линей; также стоит учесть, что индексы слоев и оконту-ривающих линий нумеруются следующим образом: если i = 1, то у = i, если же i > 1, то у = i —1; аналогично, если j = 1, то 5 = j, если же j > 1, то ô = j —1.
По известным значениям pxYÔ, pYÔy в слоях развала были найдены значения констант a.x, bx, c.x, d.x, e.x, fx, ay, by, cy они представлены в табл. 4, 5.
По предложенной методике аналитического определения узловых точек и координатной сетки блока была создана компьютерная программа в среде Microsoft Visual Studio 2019 [22]. Она
позволяет проводить удобный и гибкий расчет координат узловых точек и построить профиль развала взорванных пород. По разработанной программе были смоделированы результаты массовых взрывов при заданных параметрах БВР на 15-метровом уступе. Примеры расчетов и автоматизированного построения профиля развала горных пород представлены на рис. 1, 2. Они подтверждают работоспособность разработанного аналитического метода определения координат узловых точек координатной сетки взорванного блока и автоматизированного прогнозирования
i При однорядном взрывании
ax bx cx dx ex f
1 0,00 0,00 5,64 -7,76 4,42 0,26
2 0,00 0,00 23,56 -35,74 18,90 -1,91
3 0,00 0,00 15,45 -28,75 18,76 -2,34
4 0,00 0,00 19,63 -41,27 29,12 -4,60
5 0,00 0,00 -1,68 6,00 -5,75 4,23
i При четырехрядном КЗВ на свободную поверхность
ax bx cx dx ex f
1 2,5 3 -10,87 18,66 -15,66 6,32 0,03
2 11,17 -48,79 83,02 -68,27 27,02 -3,06
3 6,94 -32,11 57,68 -49,75 20,45 -2,05
4 4,42 -21,94 42,48 -39,50 17,42 -1,62
5 -2,66 14,57 -30,98 32,41 -16,97 5,04
i При четырехрядном КЗВ с подпорной стенкой
ax bx cx dx ex f
1 0,25 -0,89 1,59 -1,53 0,73 0,88
2 0,91 -3,51 5,47 -4,15 1,50 0,82
3 0,33 -1,02 1,23 -0,57 0,04 1,08
4 0,05 0,20 -0,99 1,64 -1,23 1,51
5 -2,24 13,12 -29,70 32,99 -18,17 5,35
Таблица 4
Константы уравнений коэффициентов пропорциональности изменения координат узловых точек в горизонтальном направлении Constants in equations of proportionality factors of horizontal change in coordinates of grid points
Константы уравнений коэффициентов пропорциональности изменения координат узловых точек в вертикальном направлении Constants in equations of proportionality factors of vertical change in coordinates of grid points
j При однорядном взрывании При четырехрядном КЗВ При четырехрядном КЗВ с подпорной стенкой
О.y b. y с.y b. y с.y b. y с.y
1 2,60 -2,15 1,11 -0,17 1,03 -0,32 1,10
2 2,81 -2,36 1,11 -0,09 1,07 -0,12 1,05
3 0,99 -0,98 0,74 -0,04 1,08 -0,0024 1,02
4 1,86 -3,00 1,62 0,03 1,06 0,07 1,01
5 0 4,30 -5,45 0,08 1,01 0,09 1,02
6 - - - 0,11 0,95 0,08 1,02
7 - - - 0,07 0,89 0,02 1,04
8 - - - 0,05 0,80 -0,0025 1,01
9 - - - 0,04 0,68 0,0045 0,92
10 - - - 0,13 0,38 -0,01 0,80
11 - - - 1,40 -1,26 -0,06 0,57
размещения различных элементов уступа во взорванной горной массе в различных условиях взрывания.
Заключение
• Для аналитического определения местоположения отдельных элементов взрываемого блока уступа в развале введены понятия «координатная сетка взрываемого блока», «координатная сетка взорванного блока».
• Совокупность взаимно пересекающихся горизонтальных и наклонных оконтуривающих линий элементов взрываемого блока в разрезе названа координатной сеткой взрываемого блока. Совокупность взаимно пересекающихся деформированных горизонтальных и наклонных оконтуривающих линий элементов уступа в развале названа координатной сеткой взорванного блока. Показаны примеры их построения.
• Разработан аналитический метод определения узловых точек и самой ко-
ординатной сетки взорванного блока. Он включает в свой состав вычисление координат этих точек по установленным зависимостям и автоматизированное построение профиля развала пород.
• Установлены зависимости коэффициентов пропорциональности изменения координат узловых точек координатной сетки взорванного блока от параметров, учитывающих коэффициент разрыхления пород в слоях, длину отброшенной части развала и константы уравнения координат узловых точек во взрываемом массиве.
• Разработана компьютерная программа автоматизированного определения узловых точек и координатной сетки развала горных пород. Полученные результаты служат надежным инструментом для оперативного определения размещения различных элементов уступа в развале: рудных включений различной конфигурации, различных зон дробления и т.д.
СПИСОК ЛИТЕРАТУРЫ
1. Покровский Г. И., Федоров И. С. Действие удара и взрыва в деформируемых средах. - М., 1957. - 276 с.
2. Ракишев Б. Р. Автоматизированное проектирование и производства массовых взрывов на карьерах. - Алматы: Рылым, 2016. - 340 с.
3. Ракишев Б. Р., Шампикова А. Х., Казангапов А. Е. Горно-геологические характеристики взорванных сложноструктурных блоков // Взрывное дело. - 2018. - № 120-77. -С. 82-93.
4. Викторов С.Д., Франтов А. Е., Закалинский В. М. Теория - техника - технология взрывных работ с применением конверсионных ВВ в процессах горного производства. -М.: ИПКОН РАН, 2019. - 384 с.
5. Викторов С.Д. Взрывное разрушение массивов горных пород - основа прогресса в горном деле // Горный информационно-аналитический бюллетень. - 2015. - № S1. -С. 63-75.
6. Казаков Н. Н., Шляпин А. В. Определение тензорного напряженно-деформированного состояния горных пород при взрыве скважинного заряда // Горный информационно-аналитический бюллетень. - 2018. - № S1. - С. 112-126.
7. Парамонов Г. П., Ковалевский В. Н., Мысин А. В. Численное моделирование разрушения блока горной породы взрывом с учетом лабораторных экспериментов // Взрывное дело. - 2019. - № 122-79. - С. 19-33.
8. Яницкий Е. Б., Кобелко С. Г., Дунаев В. А., Рахманов Р. А. Компьютерное моделирование смещения горной массы и оценка разубоживания руды в результате массового взрыва при открытой разработке месторождений // Взрывное дело. - 2018. - № 120-77. -С. 94-108.
9. Тюпин В. Н., Анисимов В. Н. Разработка методов сохранения устойчивости открытых поверхностей трещиноватых горных массивов при проведении массовых взрывов // Горный информационно-аналитический бюллетень. - 2019. - № 4. - С. 53-62. DOI: 10.25018/0236-1493-2019-04-0-53-62.
10. Халкечев Р. К. Нечеткая математическая модель изменения концентрации трещин в минерале под действием внешней нагрузки // Горный информационно-аналитический бюллетень. - 2019. - № 6. - С. 97-105. DOI: 10.25018/0236-1493-2019-06-0-97-105.
11. Сергунин М. П., Еременко В. А. Обучение нейронной сети предсказывать параметры сдвижения горных пород налегающей толщи на основании данных о трещинова-тости массива на примере рудника «Заполярный» // Горный информационно-аналитический бюллетень. - 2019. - № 10. - С. 106 -116. DOI: 10.25018/0236-1493-2019-10-0106-116.
12. Бабелло В. А., Бейдин А. В., Овсейчук В. А., Смолич С. В. Оценка состояния горного массива на основе анализа горно-геологической обстановки и моделирования его напряженности // Горный информационно-аналитический бюллетень. - 2019. - № 12. -С. 41-54. DOI: 10.25018/0236-1493-2019-12-0-41-54.
13. Казаков Н. Н. Разрушение и дробление горных пород в карьерах // Взрывное дело. -2018. - № 119-76. - С. 5-19.
14. An H. M., Hongyuan Liu, Haoyu Han, Xin Zheng, Wang X. G. Hybrid finite-discrete element modelling of dynamic fracture and resultant fragment casting and muck-piling by rock blast // Computers and Geotechnics. 2017. Vol. 81. Pp. 322-345.
15. Onederra I. A., Furtney J. K., Sellers E., Iverson S. Modelling blast induced damage from a fully coupled explosive charge // International Journal of Rock Mechanics and Mining Sciences. 2013. Vol. 58.Pp. 73-84.
16. Bakhshandeh Amnieh Hassan, Moein Bahadori Numerical analysis of the primer location effect on ground vibration caused by blasting // International Journal of Mining and Geo-Engineering. 2017. Vol. 51. No 1. Pp. 53-62. DOI: 10.22059/IJMGE.2017.62153.
17. Branko D, Detournay C., Cundall P.A. Application of particle and lattice codes to simulation of hydraulic fracturing // Computational Particle Mechanics. 2016. Vol. 3. No 2. Pp. 249-261.
18. Furtney J. K., Andrieux P., Hall A. K. Applications for numerical modeling of blast induced rock fracture / 50th US Rock Mechanics/Geomechanics Symposium. American Rock Mechanics Association, 2016.
19. Zirui Mao, Guirong Liu, Yu Huang, Yangjuan Bao A conservative and consistent lagran-gian gradient smoothing method for earthquake-induced landslide simulation // Engineering Geology. 2019. Vol. 260. Article 105226. DOI: 10.1016/j.enggeo.2019.105226.
20. Yingguo Hu, Wenbo Lu, Xinxia Wu, Meishan Liu, Peng Li Numerical and experimental investigation of blasting damage control of a high rock slope in a deep valley // Engineering Geology. 2018. Vol. 237. Pp. 12-20. DOI: 10.1016/j.enggeo.2018.01.003.
21. Лоран П.Ж. Аппроксимация и оптимизация. - М.: Мир, 1975. - С. 496.
22. Powers L., Snell M. Microsoft Visual Studio 2015 Unleashed, 3rd edition. Indianapolis, Imprint Sams, 2015. итш
REFERENCES
1. Pokrovskiy G. I., Fedorov I. S. Deystvie udara i vzryva v deformiruemykh sredakh [The action of shock and explosion in deformable media], Moscow, 1957, 276 p.
2. Rakishev B. R. Avtomatizirovannoe proektirovanie i proizvodstva massovykh vzryvov na karerakh [Automated design and production of mass explosions in quarries], Almaty, Gylym, 2016, 340 p.
3. Rakishev B. R., Shampikova A. Kh., Kazangapov A. E. Mining and geological characteristics of blown up structural blocks. Vzryvnoe delo. 2018, no 120-77, pp. 82-93. [In Russ].
4. Viktorov S. D., Frantov A. E., Zakalinskiy V. M. Teoriya - tekhnika - tekhnologiya vzryvnykh rabot s primeneniem konversionnykh VV v protsessakh gornogo proizvodstva [Theory - technology - blasting technology using conversion explosives in mining processes], Moscow, IPKON RAN, 2019, 384 p.
5. Viktorov S. D. Explosive destruction of rock masses - the basis of progress in mining. MIAB. Mining Inf. Anal. Bull. 2015, no S1, pp. 63-75. [In Russ].
6. Kazakov N. N., Shlyapin A. V. Determination of tensor stress-strain state of rocks during explosion of a borehole charge. MIAB. Mining Inf. Anal. Bull. 2018, no S1, pp. 112-126. [In Russ].
7. Paramonov G. P., Kovalevskiy V. N., Mysin A. V. Numerical simulation of the destruction of a rock block by an explosion taking into account laboratory experiments. Vzryvnoe delo. 2019, no 122-79, pp. 19-33. [In Russ].
8. Yanitskiy E. B., Kobelko S. G., Dunaev V. A., Rakhmanov R. A. Computer simulation of rock mass displacement and evaluation of ore dilution as a result of a mass explosion in open pit mining. Vzryvnoe delo. 2018, no 120-77, pp. 94-108. [In Russ].
9. Tyupin V. N., Anisimov V. N. Methods of stability retention of exposed rock surfaces in fractured rock mass under large-scale blasting. MIAB. Mining Inf. Anal. Bull. 2019;(4):53-62. [In Russ]. DOI: 10.25018/0236-1493-2019-04-0-53-62.
10. Khalkechev R.K. Fuzzy mathematical model of crack density variation in mineral under external loading. MIAB. Mining Inf. Anal. Bull. 2019;(6):97-105. [In Russ]. DOI: 10.25018/02361493-2019-06-0-97-105.
11. Sergunin M. P., Eremenko V. A. Learning of neural network to predict overlying rock mass displacement parameters by the data on jointing in terms of the Zapolyarny Mine. MIAB. Mining Inf. Anal. Bull. 2019;(10):106-116. [In Russ]. DOI: 10.25018/0236-1493-2019-10-0-106-116.
12. Babello V. A., Beydin A. V., Smolich S. V., Ovseychuk V. A. Assessment of rock mass behavior based on its geology analysis and stress state modeling. MIAB. Mining Inf. Anal. Bull. 2019;(12):41-54. [In Russ]. DOI: 10.25018/0236-1493-2019-12-0-41-54.
13. Kazakov N. N. The destruction and crushing of rocks in quarries. Vzryvnoe delo. 2018, no 119-76, pp. 5-19. [In Russ].
14. An H. M., Hongyuan Liu, Haoyu Han, Xin Zheng, Wang X. G. Hybrid finite-discrete element modelling of dynamic fracture and resultant fragment casting and muck-piling by rock blast. Computers and Geotechnics. 2017. Vol. 81. Pp. 322-345.
15. Onederra I. A., Furtney J. K., Sellers E., Iverson S. Modelling blast induced damage from a fully coupled explosive charge. International Journal of Rock Mechanics and Mining Sciences. 2013. Vol. 58.Pp. 73-84.
16. Bakhshandeh Amnieh Hassan, Moein Bahadori Numerical analysis of the primer location effect on ground vibration caused by blasting. International Journal of Mining and Geo-Engineering. 2017. Vol. 51. No 1. Pp. 53-62. DOI: 10.22059/IJMGE.2017.62153.
17. Branko D., Detournay C., Cundall P.A. Application of particle and lattice codes to simulation of hydraulic fracturing. Computational Particle Mechanics. 2016. Vol. 3. No 2. Pp. 249-261.
18. Furtney J. K., Andrieux P., Hall A. K. Applications for numerical modeling of blast induced rock fracture. 50th US Rock Mechanics/Geomechanics Symposium. American Rock Mechanics Association, 2016.
19. Zirui Mao, Guirong Liu, Yu Huang, Yangjuan Bao A conservative and consistent lagran-gian gradient smoothing method for earthquake-induced landslide simulation. Engineering Geology. 2019. Vol. 260. Article 105226. DOI: 10.1016/j.enggeo.2019.105226.
20. Yingguo Hu, Wenbo Lu, Xinxia Wu, Meishan Liu, Peng Li Numerical and experimental investigation of blasting damage control of a high rock slope in a deep valley. Engineering Geology. 2018. Vol. 237. Pp. 12-20. DOI: 10.1016/j.enggeo.2018.01.003.
21. Laurent P. J. Approksimatsiya i optimizatsiya [Approximation and optimization], Moscow, Mir, 1975, pp. 496.
22. Powers L., Snell M. Microsoft Visual Studio 2015 Unleashed, 3rd edition. Indianapolis, Imprint Sams, 2015.
ИНФОРМАЦИЯ ОБ АВТОРАХ
Ракишев Баян Ракишевич1 - академик НАН РК,
д-р техн. наук, профессор, e-mail: b.rakishev@mail.ru,
Ракишева Зауре Баяновна - канд. физ.-мат. наук,
профессор, зав. кафедрой, Казахский национальный университет
имени аль-Фараби, Алматы, Казахстан,
Орынбай Асфандияр Айтказыулы1 - докторант PhD,
старший преподаватель,
1 Казахский национальный исследовательский технический университет имени К.И. Сатпаева (Satbayev Unuversity), Алматы, Казахстан. Для контактов: Ракишев Б.Р., e-mail: b.rakishev@mail.ru.
INFORMATION ABOUT THE AUTHORS
B.R. Rakishev1, Academician of NAS of the Republic of Kazakhstan, Dr. Sci. (Eng.), Professor, e-mail: b.rakishev @mail.ru, Z.B. Rakisheva, Cand. Sci. (Phys. Mathem.), Professor, Head of Chair, Al-Farabi Kazakh National University, Almaty, Kazakhstan, OrynbayAsfandiyarAytkazyuly1, PhD Student, Senior Lecturer, 1 K.I. Satpayev Kazakh National Research Technical University (Satbayev Unuversity), Almaty, Kazakhstan. Corresponding author: B.R.Rakishev, e-mail: b.rakishev@mail.ru.
Получена редакцией 21.02.2020; получена после рецензии 30.04.2020; принята к печати 20.07.2020. Received by the editors 21.02.2020; received after the review 30.04.2020; accepted for printing 20.07.2020.