Математическое моделирование средней естественной отдельности Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

© Г.М. Рсдькин, A.B. Коновалов, В.Ю. Алсксанов, 2013

УДК 621.32

Г.М. Редькин, A.B. Коновалов, В.Ю. Алексанов

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СРЕДНЕЙ ЕСТЕСТВЕННОЙ ОТДЕЛЬНОСТИ

На основе вектора частоты системы трещин разработана методика определения формы и ориентировки в пространстве средней естественной отдельности и блочности массива горных пород. Приведены примеры определения этих параметров по участку Коробковского месторождения.

Ключевые слова: системы трещин, трещиноватость, вектор системы трещин, средняя естественная отдельность, блочность.

Горные породы всех генетических типов (магматические, метаморфические, осадочные) характеризуются ориентированными закономерно разрывами сплошности, которые называют трещинами. Трещины, имеющие одинаковую (параллельную) или близкую ориентировку, объединяют в системы трещин, совокупность которых называют трещиноватостью [1, 2].

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

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

Пусть в массиве горных пород развито три системы трещин. Если систем трещин более трёх, то их можно с некоторой погрешностью аппроксимировать тремя системами. Итак исходим из того, что массив горных пород характеризуется тремя системами.

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

Рис. 1. Средняя естественная отдельность, образованная тремя системами трещин

= « (- sin Ai sin 5г; cos Ai sin 5г; cos 5г), (1) i = 1,2,3,

где - частота трещин i-ой системы; AjSj - соответственно азимут линии простирания и угол линии падения плоскости i-ой системы трещин;

Щ = 1.

Тогда, из векторного представления (1) следует, что обусловленная этими системами трещин интенсивность тре-щиноватости массива горных пород адекватно выражается суммой модулей проекций векторов (1) систем трещин на направление e = (cosa;cosP;cosу); \e\ = 1 [3].

3

= ^юг |—sin Ai sin 5г cos a + cos Ai sin 5г cosp + cos5i cos у |. (2)

i=i

Далее среднюю естественную отдельность и блочность моделируемого массива горных пород определим следующим образом.

Обозначим через H - вектор перпендикулярной плоскости системы трещин <, модуль которого равен среднему расстоянию между трещинами системы. r

Тогда векторам систем трещин (1) будут соответствовать векторы Hi =—,,

«г

i=1,2,3, которые можно выразить через векторы (1)

йг = ÍTL«г, i = 1,2,3. (3)

Ы

Кроме того, векторам систем (1) и индуцированной ими трещиноватости (2) отвечает средняя естественная отдельность в форме параллелепипеда, высотами которого являются векторы (3) (см. рис. 1).

Требуется по данным высотам (3) найти образующие параллелепипеда, которые определяют его форму и ориентировку в пространстве, и объём параллелепипеда — блочность массива горных пород.

Векторные произведения векторов (3) определяют направления в пространстве рёбер, образующих параллелепипед

г, =

г, =

Гя =

щ

к

и\

<И=1<

= г

ини высот (3)

.Щ= 1.

(4)

Единичные векторы (4) рёбер параллелепипеда образуют углы с векторами его

которые определяются по формулам

п-т

еоза,- = ■ 1 Ш

, I = 1,2,3-

______(5)

Обозначим через Я.1, Я.2'^-1- векторы, совпадающие с рёбрами параллелепипеда и образующие его. Тогда модули этих векторов можно выразить через модули векторов высот параллелепипеда (3)

(6)

Сами векторы С — получим, умножая их модули (6) на соответст-

вующие направления (4)

(7)

где единичные векторы Г определяются равенствами (4). Векторы (7) определяют форму и ориентировку в пространстве средней естественной отдельности.

Модуль же смешанного произведения векторов (7) определяет блочность моделируемого массива горных пород

_ _ _ Н^Р-Н^Р-и^Р

(8)

Выразим в равенствах (4), (7), (8) векторы И (3) через векторы систем трещин (1) и получим в зависимости от элементов залегания систем трещин: векторы, образующие среднюю естественную отдельность

п2 X щ

С0± ■ ¡П! х п2 ■ щ

й3=-

■ (Л, У.Щ-Щ

О

■ [И, У.Щ-Щ

и объём средней естественной отдельности

V =

- " «V К^Г X ■ ПрГ

где

(10)

(11)

причём индексы 1 и т принимают значения согласно формуле (10).

Лля наглядности рассмотрим условный пример. Пусть есть 3 системы трещин, которые параллельны координатным плоскостям и расстояние между трещинами равно единице. Данным системам трещин отвечают: векторы систем трещин

= = -■' с}, (.'.'а = £0^0! 1}; интенсивность трещиноватости в

направлении единичного вектора е по формуле (2) равна

■'■ 1 .'■ 1 - (рис. 2, а). Рассчитаем направляющие вектора и объём блока средней естественной отдельности: максимальные интенсивности трещиноватости систем = = = 1, направления векторов

= у П; о], = {П; 1; о} ??.3 = {П; П; 1], их смешанное произведение

!;□]-[0;0я])= 1. Направляющие вектора блока ^средней отдельности подформуле (10) равны

ём блока по формуле (11) равен V =■ 1

(рис. 2, б), объ-

1 -1 • 1 • |1|

= 1.

Приведём пример расчёта ориентировки и размера средней естественной отдельности для участка №402 Коробковского месторождения, где есть 3 системы трещин с векторами

(8,53; 8,29; 0,7г}, &£= {2,34; -4,43; 0,Зь}ГШ7= {1,47; 0Д5; 4,6?]

а) б)

Рис. 2. Визуализация модели трещииоватости условного примера: а) интенсивность по направлениям; б) блок средней естественной отдельности

а) б)

Рис. 3. Визуализация модели трещииоватости участка №402 Коробковского месторождения: а) интенсивность по направлениям; б) блок средней естественной отдельности

Разложим каждый вектор системы на его направление и интенсивность:

- 11,917, Т£ = {0,716;0,69&; 0,Оби}. о>2= 5,023, ¡0,466; -0,882: 0,07гЬ

о>з= 4,898, т£= {0,300; 0,031; 0,95*Ь

Интенсивность трещиноватости в направлении единичного вектора е по формуле (2) равна (рис. 3, а)

¿(3)= 11,916 |0,716 ё£+0,696 0,060 Щ \ +5,023 ■ |0,466-(£- 0,882 ■ 0,072 + 898■ |0,300 ё^Н-0,031 0,953-Щ.

Рассчитаем по формулам (10), (11) направляющие векторы и объём V блока средней естественной отдельности (рис. 36):

_{0,466; -0,882; 0,072] х {0,300; 0,031; 0,95з]__

11,917 ■ ({0,716; 0,696; 0,Обо] х {0,466; -0,882; 0,072] ■ {0,300; 0,031; 0,95з])

Е3 = [ОД 50; -0,150; -0,04зЬ

Вш - {-0.024; 0,005; 0,221],

1

V =

11,917 - 5,023 -1-0.8811

= 0,00387.

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

1. Нейштадт Л. И. Методы инженерно-геологического изучения трещиноватости горных пород / Л. Н. Нейштадт, И. А. Пирогов. - М.: Энергия, 1969. - 248с.

2. Рац М. В. Трещиноватость и свойства трещиноватости горных пород / М. В. Рац, С. Н. Чернышев. - М.: Недра, 1970. - 160с.

3. Редькин Г. М. Нестационарное анизотропное математическое моделирование неод-нородностей систем минерального сырья / Г. М. Редькин. - М.: Издательство Ассоциации строительных вузов, 2007. - 500 с.

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

Редькин Г. М. - доктор технических наук, доцент, профессор кафедры прикладной математики, Белгородский государственный технологический университет им. В. Г. Шухова, redking@emai1.ru Коновалов А. В. - аспирант кафедры прикладной математики, Белгородский государственный технологический университет им. В. Г. Шухова, akon32@ramb1er.ru .

Алексанов В. Ю. - аспирант кафедры прикладной математики, Белгородский государственный технологический университет им. В. Г. Шухова, V a1exanov@mai1.ru .