Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Рыжов, П. А. Математическая статистика в горном деле учеб. пособие

.pdf
Скачиваний:
25
Добавлен:
19.10.2023
Размер:
7.54 Mб
Скачать

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

Выдано за отчетный период

Из блока Ai Из блока А<>

Из блока A3

Общая числен­ ность

 

 

 

 

 

 

Т а б л и ц а 15

Число

вагонеток ( ч а с т о т а )

 

с

содержанием

 

 

металла

п %

(С^)

 

 

 

1-3

3-5

 

 

 

1 7 - 9

В с"е г о

5 - 7

 

Ci

с в

С

 

j

С,

 

o

 

 

 

3

 

 

 

О tо

 

 

 

 

 

 

 

ОО ссо

80

60

 

20

340

 

 

 

140

80

 

20

330

 

140

80

 

50

330

330

360

220

 

90

1000

Прямоугольная таблица, подобная табл. 15, называ­ ется матрицей. Порядок этой матрицы /If-С,-, в нашем случае 3X4 (3 — число строк, 4 — число столбцов).

Выпишем на билеты все номера вагонеток, с указа­ нием того, из какого блока она выдана и с каким содер­ жанием металла.

Свернем билеты в трубочки и вложим в патроны оди­ накового цвета, размера п веса. Поместим эти патроны в урну и перемешаем. Наудачу достанем один билет (этим достигается равновозможность для каждого биле­ та быть вынутым):

Общее число вагонеток, выданное из очистных бло­ ков, называют безусловной абсолютной частотой появле­ ния вагонеток из данного блока.

Эти частоты: для

блока

Лі = 340; для блока

Л 2 = 3 3 0

и для блока Л 3

= 330.

 

частота

появления

вагоне­

Абсолютная

безусловная

ток с содержанием металла 7—9%

(Ci)

составит 90, с со­

держанием С3 — 220;

С2 — 360 и С\ — 330.

 

Если все элементы

матрицы,

представляющие абсо­

лютные безусловные частоты, разделить на общую чис­ ленность вагонеток (N=1000), получим результаты, при­ веденные в табл. 16.

Каждый элемент матрицы табл. 16 выражает стати­ стическую вероятность события, заключающуюся в одно­ временном выполнении условий А І и С], и имеет следую-

щии вид

 

 

 

 

 

 

 

 

P{AiCj)

= ^j-

= Pij,

 

 

 

 

 

 

 

Т а б л и ц а 16

Блок

у /

 

Содержание

металла в %

 

 

 

 

 

 

 

- ^ С т а т н с т н -

 

3 - 5

 

 

Сумма

 

ческне,

1- 3

5 - 7

7 - 9

 

s'

вероят-

С,

с.

С,

с,

 

^

НОСПІ

 

 

 

 

 

 

 

0,180

0,080

0,060

'0,020

0,340

 

 

0,090

0,140

0,080

0,020

0,330

As

 

0,060

0,140

0,080

0,050

0,330

Общая

числен­

0,330

0,360

0,220

0,090

1,000

ность

 

 

 

 

 

 

Следовательно, вероятность того, что взятая наудачу вагонетка окажется с содержанием С2 (3—5%) из бло­

ка Л 3

Р ( Л 3 С 2 ) = 140:1000 = 0,140.

Величина Р(Л 3 С 2 ) называется совместной

вероятно­

стью (или вероятностью совмещения) Ai и Cj.

 

В случае,

когда

рассматривается значение только од­

ного признака А (или только С), будем

иметь

безуслов­

ную вероятность Р ( А І )

или

P(Cj).

 

 

Безусловная вероятность события С2 (вагонетка с

содержанием

металла 3—5%), выразится как

 

Р(С 2 ) =

Р (А{С2)

+

Р 2 С2 ) + Р (А3С2)

= 0,080 +

 

+ 0,140 + 0,140 =

0,360 и т. д.

 

Какова вероятность того, что взятая наудачу вагонет­

ка скажется

с содержанием

металла

С2 при условии,

что эта вагонетка

будет из выемочного

блока

Л3 ? Это

записывается выражением Р(С 2 3 ) .

 

 

Вероятность С2

при условии выполнения Л 3

эта вели­

чина будет условной вероятности С2 при данном А 3 .

Из 330 вагонеток, выданных из блока Л3 , вагонеток с

содержанием С 2 — 140.

 

того, что вагонетка ока­

Следовательно,

вероятность

жется с содержанием С2 , будет

140:330 = 0,424. Это мо-

— 61 —

жет быть выражено формулой

Р(С23)

= -

(2.9)

Условная вероятность события Со при условии выпол­ нения А 3 равна вероятности совместного наступления событий Со и /1з, деленной на безусловную вероятность события А 3 . В нашем случае

Р(Сг3)

0,140

0,424.

 

0,330

Рис. 11

На рис. 11 события / 1 2 , А 3 , С2 представлены в виде подмножеств. Вероятность события заключается в том,

что наступит или А 2 ,

или С2 :

 

 

Р(А2 или С2 )

= Р (А2) + Р(С2)-Р

2А2).

(I)

Вероятность того, что наудачу отобранная

вагонетка

ока­

жется с содержанием Со (3—5%), или она окажется из

очистного

блока А 2 , равна сумме безусловной

вероятно­

сти того,

что она окажется с содержанием

С2Р(С2)

=

= 0,360 и безусловной вероятности того, что она

выдана

из блока А2Р(А2)

=0,330, минус вероятность

одновремен­

ного появления этих событий Р(С2А2)

=0,140:

 

Р(А2

или d) = 0,360 + 0,330 — 0,140 =

0,550.

Из рис. 11 видно, что при суммировании

безусловных

вероятностен величина Р(А2С2),

заштрихованная на

рисунке

.площадь

S, — приплюсовывается дважды: один

раз — в виде слагаемого в Р(А2),

а второй — в

Р(С2).

Чтобы избежать

этого, вычитают

Р(А2С2).

Следователь­

но, выражение (I) есть формула сложения вероятностей.

Какова вероятность события 2 или Л3 )? Из табл. 16 следует, что Р(А2) =0,330; Р(А3) =0,330. События А 2 и

А3 — исключающие друг друга события, следовательно,

событие (Л2 Л3 )невозможное, его вероятность равна нулю. По формуле (I) имеем

Р 2 или А І ) = 0,330 + 0,330 - 0 = 0,660.

— 62

Событие ( А 2 ) противоположно в статистическом смысле событию (Ai или Л 3 ) .

Событие (Лі или /13) обозначим через (Л2 ), тогда

Р(А2) + Р(А,) = 1.

Аналогично

Я(С2 ) + Р ( С 2 ) = 1.

Условная вероятность события Ai (при условии, что осуществлялось событие Cj) равна условной вероятности события А при условии, что совершилось событие Cj,

P(At/C})

= P (At/Ci).

(II)

События Ci и Aj будут независимыми в статистическом смысле. Если верно выражение ( I I ) , то

Р ( Л , / ^ ) = Р ( Л , ) .

(III)

Д о к а з а т е л ь с т в о .

 

P(Ai) =

P(AiCj) +

 

P(AiCj).

 

Следовательно,

 

 

 

 

 

 

Р (АІ)

=

Р (Ai/Cj) Р (Cj) +

Р (Ai/Cj)

Р (СІ) .

 

Из выражения

(III) имеем

 

 

 

 

 

P(Ai)

= P(Ai/Cj)

[Р(С,)

 

+

Р(П,п

 

Так как

 

Р(АІ)

= Р(АІ/СЛ.

 

 

 

 

 

 

 

 

 

 

 

 

 

Р(АІ/С)

 

Р

ІМ/Ci)

 

 

 

 

 

=

 

 

 

 

 

 

 

 

Р(СІ)

'

 

 

то

 

 

 

 

 

 

 

 

 

 

Р(А

СІ)

= Р(АІ)Р(СІ).

 

 

(IV)

Если два

события

независимы

в статистическом

(ве­

роятностном)

смысле, то вероятность их совместного

по­

явления равна произведению соответствующих безуслов­ ных вероятностей.

Равенство (IV) выражает статистическую независи­ мость. Какова вероятность, что наудачу отобранные две вагонетки окажутся обе из блока Л2 ? Вероятность того, что первая вагонетка окажется из блока Ля, равна 0,330.

— 63 —

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

Вероятность того, что обе вагонетки окажутся из бло­ ка Л2 ,

Р(А2,Ап)

= Р 2 ) • Р 2 ) =

0,330 • 0,330 = 0,109.

Теорема

гипотез (формула

Бейеса). Теорема гипотез

является следствием теоремы умножения и формулы пол­

ной

вероятности.

 

 

 

 

 

 

 

Кі,

К.2,

Имеется полная

группа

несовместных

гипотез

 

Кп, вероятности которых известны и равны соот­

ветственно Р{Кі),

Р2),

Р(Кп).

 

 

 

Произведен опыт, в результате которого отмечено по­

явление события В . Как следует

изменить

вероятности

гипотез в связи с появлением

события В ? Задача сводит­

ся к нахождению условной вероятности Р /В)

для каж­

дой из гипотез.

 

 

 

 

 

 

 

 

Из теоремы умножения вероятностей имеем:

 

Р(В)-Р

(КІ/В)

=

Р

(КІ) Р

(ВІ/КІ)

(і = 1,2,3,.... п),

откудаР{Кі/В)=щ^т

 

 

„_,,„

„).

 

 

 

 

Р(В)

 

 

 

 

 

Выражая Р ( В )

по формуле полной вероятности

(2.8),

получим формулу

 

Бейеса:

 

 

 

 

 

"^ЩШт- <'-'•** »>•

Вероятность события КІ, вычисленную при условии, что имело место другое событие В , обозначим через P(Uj) и введем символ умножения П. Тогда формула Бейеса примет вид

P(U,) =

(2.11)

m

2 Pill Pu

Если априорные вероятности гипотез неизвестны, то гипотезы принимают равновероятными и формула Бене-

64 —

са будет иметь вид

р т = — 11— •

,(2.і2)

m

 

3 = 1

Так как событие U может происходить при осущест­ влении одной или нескольких несовместных гипотез, то

m

2 ^ з ) = 1 -

(2-13).

3 = 1

 

 

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

Пример 2.11. Прогнозирование смещений кровли под­ готовительной выработки. При проектировании проход­ ки подготовительных выработок важно правильно вы­ брать крепь, определить необходимый типо-размер крепи с заданным интервалом ее податливости.

Величина податливости крепи определяется смещени­ ем кровли выработки.

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

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

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

факторов, от которых зависит устойчивость

выработки.

В приложении X I I * приведены данные

результатов

многолетних наблюдений за смещением кровли в 69 вы­ работках шахт комбината Донецуголь, охраняемых мас­ сивом пород и бутовой полосой. Время, за которое про-

* Материалы таблицы и пример предрасчета смещении кровли выработок взяты из докторскоіі диссертации И. Л. Черняка.

3—1924 — 65 —

изошло зафиксированное смещение, определяется деле­ нием / (графа 12) на ѵ (графа 9).

Из приложения X I I видно, что каждый фактор, от ко­ торого зависит устойчивость выработки, варьирует в оп­ ределенных пределах и может иметь количественное (глубина заложения выработки, мощность угольного пласта и т. д.) или качественное выражение (тип пород непосредственной кровли). Однако каждый конкретный участок выработки, где наблюдалось смещение кровли, характеризуется для каждого фактора только одним его значением. Данные таблицы показывают, что глубина заложения выработки варьирует в пределах от 84 до 1038 м; мощность угольного пласта — от 0,45 до 1,7 м; крепость пород непосредственной кровли при одноосном сжатии 170—1086 кГ/см2; расстояние от места наблюде­ ния в выработке до очистного забоя — 0—250 м; сечение выработки в свету—5,0—17,4 м2; ширина бутовой поло­ сы — 5—37 м; длина очистного забоя— 100—305 м\ ско­ рость подвпгания очистного забоя—10—70 м/мес; мощ­ ность пород непосредственной кровли от 1,5 до 16,7 м.

Наблюдаемое смещение пород кровли также варьи­ рует в пределах от 12 до 912 мм.

Разделим весь диапазон смещений на следующие три интервала:

0 <

Ui <

100,

 

 

100 <

Uz<

300,

 

 

 

 

U3 >

300.

 

 

Рассортируем выработки по смещениям

кровли, в ре­

зультате найдем, что в интервал

11\ попало

34 выработ­

ки, в интервал U2 — 22 выработки и в интервал

U% —

13 выработок (всего выработок 34 + 22+ 13 = 69).

 

Найдем также, что из 34 выработок на глубине

зале­

гания Hi до 500 м оказалось 22 выработки,

11 выработок

с глубиной залегания Н2

от 500 до 700 м и одна выработ­

ка на глубине # з более 700 м.

 

 

 

Вычислим вероятности Рц:

 

 

 

РІЛ = 22:34 =

0,65;

 

 

Ріл=

11:34 =

0,32;

 

 

P j . 3 =

1:34 =

0,03;

 

 

— 66 —

Результаты 'вычислении запишем соответственно в графах 4, 5 и 6 на первой строке, соответствующей смещению пород кровли до 100 мм.

Аналогично вычисляем вероятности Рц и Р^у.

Ри

= 6:22

=

0,27

Рг.2 = 9:22 =

0,41

Рг.з — 7:22

=

0,32

2 ; == 1.

 

 

Л и

= 2:13

=

0,15

Р3

= 8:13 =

0,62

 

 

 

 

 

 

Рз.а = 3:13 =

0,23

 

 

 

 

 

 

 

 

 

 

 

3 ; == 1.

 

 

 

 

 

 

Результаты

 

вычислений

записывают

соответственно

на второй и третьей строке в графы 4, 5, 6.

 

 

Рц

Аналогично

вычисляют

вероятности

совмещения

по всем

другим

признакам — факторам. Результаты вы­

числений записывают в табл. 16а.

 

 

 

Т а б л и ц а

16а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р а с с т о я н и е '

 

 

 

Ширина

буто ­

Прочность

 

 

от

очистного

Глубина

И, м

вой

полосы

породы о,

 

 

забоя

 

Л м

 

 

 

 

/„. •*

 

кГІсм

 

Интервалы

 

100—

 

 

 

 

 

15—

 

 

 

 

смещений,

мм

100

>200

<500

5 0 0 -

>700

<15

>25

<500

600 -

>800

 

 

200

 

700

25

800

 

 

1

о

 

 

3

4

5

6

7

8

э

10

11

12

<100

 

0,7б|о,21 0,Оз|о,65 0,32 0,Оз|о,59 0,26 0,15 0,38 0,36 0,26

100—300 0,4з|о.41 0,14 0,27 0,4і|о,32 0,68 0,27 0,05 0,32 0,45 0,23

>300 0.08 0,38 0,54 0,15 0,62 0.23jo. 69 0,16 0,15 0,38 0,47 0,15

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

Питерит е р- налы с м е щ е ­

нии, ММ

<100

 

 

 

 

 

 

 

 

Продолжение

табл.

Іба

 

 

 

 

 

 

Скорость

Мощность

Наименование

Сомснне

Мощность

Длина

непосред ­

иодннга-

 

пород

 

ственной

 

 

5",

/л -

ігіаста m

лавы Z,

ШПІ V,

 

 

 

 

 

кровли

 

 

 

 

 

 

 

 

 

м'.мес

'"н.к

глинисі сланец

песчанист, сланец

песчані

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<8

>8

<1

M

•;• ISO M 80

<30

>зо

<5

>5

 

 

 

ІЯ

м

15

111

17

18

19

20

21

22 23

24

25

0,44 0,56

0,29 0,71

0,56 0.44 0,59 0,41

0,39 0,63 0,53 0,44 0,03

100-300 0,50 0,50 0 , 2 з | о , 7 7 0,41 0,59 0,55 0,45 0,32 0,68 0,64 0,26 0,10

>30 0 0,46 0,54 0,15 0.S5 0,31 0.69 0.60 0,31 0,08 0,92 0,85 0,05 0,10

На основании этоіі таблицы, по формуле (2.12) де­ лается предрасчет величины смещения пород кровли в проектируемой выработке. Например, по данным геоло­ гической разведки и проекту плана горных работ выра­ ботка должна проводиться вслед за очистным забоем и охраняться по схеме «бутовая полоса — массив». Глуби­ на залегания выработки #=1000 м, сечение выработки 5—14 м2. В кровле выработки залегает глинистый сла­ нец мощностью / н , і . к = 10 м и прочностью при одноосном сжатии а,; = 400 кГ/см2. Очистной забой длиной Z=250 м разрабатывает пласт мощностью /»„=1,3 м со скоростью

подвигаиия ѵ = 20 м/мес.

Ширина бутовой полосы /я =

= 8/«.

 

Определить интервал

с/,-, в который может попасть

выработка, по величине смещения кровли иа расстоянии от очистного забоя / = 250 м.

В соответствии с постановкой задачи выдвигаются ги­ потезы попадания места выработки (250 м от очистного за­

боя)

в каждый из трех интервалов по смещениям

U\, U2,

Uz. Из табл. 16а выписываются значения Рц

в

соответ­

ствии

с заданными признаками-факторами,

влияющими

на величину смещения кровли выработки, из граф 3, 6, 7, 10, 14, 16, 18, 19, 22, 23 и сводятся в табл. 166.

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

— 68 —

Интерпалы смещен ніі

1

и

 

а

s

m у

z

•J, мм

<100

0,03

0,03

0,59

0,38

0.56

0,71

0,44

100—300

0,14

0,32

0,68

0,32

0,50

0,77

0,59

>300

0,54

0,23

0,69

0,38

0,54

0,85

0,69

По данным предпоследней колонки по формуле Бейе-

са (формула 2.12) вычисляют вероятности

Р\, Рч и Ръ-

Из приведенного в табл. 166 решения

видно,

что в

проектируемой выработке смещение кровли будет

свыше

300 мм, отсюда необходимо задаться податливостью кре­ пи и определить ее необходимый типоразмер с заданным интервалом податливости.

§ 3. СЛУЧАЙНЫЕ ВЕЛИЧИНЫ И ИХ РАСПРЕДЕЛЕНИЯ

Случайной называют величину, которая в результате опыта может принять любое значение, причем неизвест­ но, какое именно.

Величины,

значения

которых возможно

заранее

указать, в отличие от случайных, называют

определен­

ными.

 

 

 

Случайные величины могут быть прерывными, или

дискретными, и непрерывными.

 

Примером

дискретной

случайной величины может

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

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

В дальнейшем случайные величины будем обозначать большими буквами (X, Y и т. д.), а их возможные зна-

— 69 —

Соседние файлы в папке книги из ГПНТБ