книги из ГПНТБ / Чесноков, Н. И. Оптимизация решений при разработке урановых месторождений
.pdfіо
.табл
Продолжение сч •*»
=п
С )
>»
соо О!—Iсосооо05со——сч—_ счо юо о сосо05N. |
_ |
со |
||||
|
'S*—о счСОсоN N 05сч |
N со сч о со |
сосо |
|
|
|
о Nю —СОСОсо—о. юсо—соо сосч—со |
сч_ сососч Осо _ о |
|
|
|||
~ |
СО о Nюо ю О)N о со со05 ( |
>1* |
соо ІОсчо — |
СО |
||
|
о Мю — оо ю |
|||||
о о о о о _ о о о юо о о _ о сч*фсоо соо __ |
Счсч соо |
|
— |
|||
|
|
|
|
|
|
05 |
|
сч— |
счсч |
|
|
|
сч |
с о с о о 05 с ч — w o o f o o — — — o i — с о с о л с о ^і* о о с о — o o n - n - G 5 o o i o c o o o
с о с > . і л с о с 5 с о о > ю ю с о с о — ю с о с о с о ю с о г * - - « ч о с о е ч а о о с о ,« ' с о ю о о е о
NSIsVCCOtOinW^'WOOK'JlOlßiO — СО —iW(OOifM'.tDOOlO*T<ßO<ß
coo)toococo5Oioo^wc)0ortNo)NcO’“"-owioo5Occcoe'iioc'iTn сот*0»гсо —^noirtiflN'j'NCiOOiocoioco^rciaiONioc'WrcowN
ю —о^ |
■«• со ■«■ w о <о ю сил — о N w w o -ia |
——coco |
— |
05 с— —■«■осо - ~ - сч ю |
— і'- |
|
|
|
со |
|
|
|
ю |
-г-і'-іосч —сію —oo-«<cooicocoocioooN.C'j^j‘cocoioioco«oioo)iomeo |
|||
СітгСОСЧСОСіСОС'-СО—СОЮЮСОтгСОіЖООсОЮСО^СОі'-СЧСоС'ІСО—CNN |
|||
. — CO — N c O C O l O 5 ) « C O O ( f l f f ) O M O O ) l O l O O C O C O ' 4 , 0 ) - CO —. « « |
|||
airrOocoococoNcnco—c4O5NC4C0~аоюоспсосчсосоюососою — |
|||
c4N-e‘C4T*ococ4cocoNC4iOG5NC4oo—con-^co сч а «• — сч со со |
|||
— |
■«■ СЧ |
« n v — |
--------со |
N 0 0 5 — C O O N M O C S O — 0 « С О — ' O ' c S O O O ' O f l O — C O T t N ^ Ö C O N i r :
оО -Оооосчооо— О С І — І ' - О -ОсооосчО О Ю т г О -—‘Cсо0 0054 J -ч—'0C кO C--Oсl'OіосчN O O 'l іO,lеоооосчсt ^ O — — C O'Oі
ООООйООООЮ-OOOOOCOx-OOOOOOOQ-QOOOf
О О О О О О О О О о О о О о О О О О О О С І О О О О О О О О О о - С
СГ)-аоЮ(5'Г«''(£іО5)'ДіЛО-ШЩ-<0---ЮОЮ0(3«іЛ'--М |
|||||
О«оосчід««союоолгчоммлсчюочс-сіі©и)—«гѵою«о |
|||||
NOC4’rOCO«C4NU:OoNiO*T>tNNCO(NON'r005)<OC4lOCOCON |
|||||
M i n « ' f l U ) l ß i 0 5 > C 5 ( O C ' H e « C O C O < Ö O N O N 1 N v W « U 3 0 N 0 5 0 0 |
|||||
N N . •«• о о — locoscooifloooocooiPJcoioN^cocsrtw |
|||||
*«>со •«< СЧ — |
СО СЧ N |
СЧ О CO |
с о |
COCO C0 eO CO CO |
— 05 — |
coco |
—I |
— ю со |
— |
„«-..см |
о |
— |
— — |
|
о |
с ч |
|
|
|
|
|
|
|
|
|
|
|
Щ |
|
|
|
|
|
сч |
|
N — O O l O t O M C 4 C 5 « O f O C l O O |
— ■S, 4 , 0 5 l O C O O J < O N ' 4 , < O t O C 4 l O — C 5 N |
||||
C 5 N I O C 4 N . — — — COCON - О О С Ч С О — N * a 0 C O < O l O C 4 4 r N . C O < O O ) — i ß l O C D O
со со со |
сосочг —со —оо—аосчсчсч |
со сч со о со |
|
|
СО СЧ |
— — СО ГГ — |
ю |
C4ON0100O5T#->«><N*«'C04**0>O*'t0<0t0tÖN<0OOlQN.t0NC5C0C0G5O
05OC4N-OC0 — СО — -
со — ечосч — со —
c o d о — C l C O - s M O C O N O O O l O — C 1 C O T M O C O N . O O O O — С Ч С О ^ Ю С О С ^ С О *
с о с о . т ^ , ^ , - ^ ' ' ^ , , г т? ' ' ,С ' ^ ‘ ',* ' 1/31л ю і л ю и э ю и э ю і о с о с о < о с о с о с о с о у э с о М
Для аппроксимации |
было |
выбрано |
уравнение вида |
||
■ |
- |
|
, |
Ь |
|
|
|
У = |
а -\- -XI . |
|
|
Значение |
коэффициента а у х |
определено по формуле |
|||
|
V |
V |
—___ V - L |
|
|
|
j L A |
J l Z j |
£ |
- .V: - |
X[ |
a y x |
— |
|
|
|
|
1 |
Z |
4 |
6 8 10 |
ZO |
4 0 60 |
Lg |
P H c. 10. Эмпирическая линия регрессии в логарифмическом масштабе.
после подстановки численных значений
аѵг = |
758,02 - 4,4205 — 279,3679 - 11,0627 |
, . . |
-------------------------------------------68 - 4,4205— 11,06272 |
1,46. |
|
у х |
|
Значение коэффициента Ьух определено по формуле
іП m l |
V i V |
У і |
|
> |
Xi |
_ j *, и |
|
Ьѵх = |
|
|
|
•2 |
* - c b £ 7 |
|
|
71
после подстановки численных значений
|
_ |
68-279,3679 — 11,0627-758,02 |
59,54. |
|||||
у х ~ |
68-4,4205— 11,06273 |
|
|
|||||
|
|
|
|
|||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
Уі = 1,46 + |
^ і . |
|
|
|
(2.34) |
|
|
|
|
|
Xi |
|
|
|
|
Проведенная проверка |
|
|
|
|
|
|||
|
|
|
V* |
0 |
|
|
|
|
|
|
|
" 4-' Ус |
|
|
|
|
|
|
R |
= ± |
V |
2 |
|
|
2 |
|
|
|
|
п - f Vi - |
|
|
|
|
|
= ± V |
68 -17 737,36 ■38 — 574 594,3204 |
|
± 1/0,8030 == |
|||||
68 • 20 015,7862 — 574 594,3204 |
|
|
||||||
|
|
= |
+ 0,896, |
|
|
|
(2.35) |
|
показала, |
что |
аппроксимирующее |
уравнение хорошо |
|||||
соответствует |
экспериментальным |
|
данным |
выборки. |
||||
Уравнение (2.34) для определенных условий может |
||||||||
служить математической |
моделью |
прямых |
затрат на |
|||||
добычу условной единицы продукции очистными рабо тами. При этом оно в скрытой форме учитывает взаимо связь комплекса определенных технических параметров, находящих отражение в затратах на добычу, с геоло гическим — продуктивностью отрабатываемой жильной площади.
При невозможности определения характера связи между изучаемыми признаками при качественном ана лизе их значений в поле корреляции определяют и коэф
фициент |
корреляции, |
и корреляционное отношение. |
Если они |
оказываются |
равными или близкими друг |
другу, то зависимость линейна, если же корреляцион ное отношение существенно больше значения коэффи циента корреляции, то связь нелинейна. В этих случаях характеристикой тесноты связи служит корреляционное отношение (2.35), интерпретация которого не зависит от вида исследуемой корреляционной зависимости.
Подробнее вопросы аппроксимации эмпирических зависимостей изложены, например, в работах [6, 8, 30].
Понятие о множественной корреляции. Многообра зие существенных факторов, действующих на нзучае мый производственный процесс, предопределяет необхо
72
димость использования методов множественной корре ляции, с помощью которых исследуется совместное влияние двух и более факторов (х, v, z) на результи рующий признак у. Как отмечалось раньше, состояние многих систем горнорудного производства (в данном случае под системой понимается совокупность взаимо связанных признаков, определяющих тот или иной про изводственный процесс или явление) характеризуется
сложной |
статистической |
взаимосвязью их элементов. |
|
В определенных случаях |
глубокие |
внутренние связи, |
|
присущие |
процессам горнорудного |
производства, не |
|
позволяют их расчленить до состояния парных зависи мостей, для анализа которых пригодны изложенные методы корреляционного и регрессионного анализов. В этих случаях для анализа взаимосвязей переменных величин и их множественного количественного воздей ствия на изучаемые признаки применяют методы мно жественной корреляции.
Применение множественного корреляционного и ре грессионного анализов предполагает выполнение пред посылок, сформулированных ранее для парного корре ляционного и регрессионного анализов.
Рассмотрим метод корреляционного анализа влия ния двух переменных факторов на результирующий признак на примере. Предположим, что необходимо изучить действие водовяжущего отношения и модуля крупности песка на прочность твердеющей закладки *. Решение этой задачи имеет существенное значение при разработке урановых месторождений, поскольку на них широко применяют твердеющую закладку, а излагае мый метод может быть с успехом использован для вы бора ее состава и технологии приготовления.
Исходные данные, а также результаты промежуточ ных расчетов удобно представить в форме, приведен ной в табл. 6.
Рассмотрим последовательность действий. В общем виде связь между изучаемыми явлениями может носить как линейный, так и нелинейный характер. В последнем случае применяют обычные методы линеаризации зави симостей— логарифмирование, дифференцирование, а затем проводят анализ влияния переменных факторов на результирующий признак.
* Пример взят из работы [46].
73
Номер опыта (замеса) |
Водовяжущее отношение |
Модуль крупности песка |
Прочность. кГ/см2 |
1 |
|
1 |
1 |
] |
2 |
3 |
4 |
Т а б л и ц а 6
5? |
гз |
Ч |
Ö |
о |
=Г) |
|
гл |
ч |
|
||||
5 |
6 |
7 |
8 |
9 |
10 |
11 |
•
Предположим, что зависимость между указанными в примере факторами линейная и выражается уравне нием регрессии
УXV= «о + ЩХ + а&. |
(2.36) |
Необходимо отыскать значения параметров уравнения «о, Яь Яг, при которых
2 ( у XV — У ? = min.
Используя метод наименьших квадратов, получим для этой цели следующую систему нормальных урав нений:
па0 + а1У х + а2^Іѵ = УІу,
<з0 2S + |
аі 2 |
“Ь |
S |
|
a 0 V V + a Y У х и + a 2 У u z = v y v . |
||||
Следующий |
этап |
расчета — подстановка в систему |
||
рассчитанных |
сумм |
(см. |
табл. |
6). Решив систему ли |
нейных уравнений с тремя неизвестными, получим па раметры ао, аи а.2 . Затем в уравнение регрессии (2.36),
представленное в общем виде, подставим эти значения параметров и найдем эмпирическое уравнение регрессии для данного исследования.
При решении системы нормальных уравнений ис пользуют так называемое правило Крамера, согласно которому система из линейных уравнений с п неизвест ными, содержащая отличный от нуля определитель, всегда имеет единственное решение. Каждый из неиз
вестных |
равен дроби, где знаменатель — определитель |
[ВФО), |
а числитель получен из определителя системы |
74
путем замены столбца коэффициентов при неизвестном столбце свободных членов.
Тесноту связи между факториальными |
признаками |
.г, у и результативным признаком у можно |
определить |
с помощью множественного корреляционного отношения
|
аи ~ |
(XV)% |
(2.37) |
|||
|
RX V |
|
о |
|
||
где о 2у — дисперсия признака |
у |
под |
воздействием |
всех |
||
факторов; |
(.Vl,)a® (читается — сигма |
квадрат |
игрек |
при |
||
исключенном хи) — дисперсия |
признака у за |
счет всех |
||||
факторов, |
за исключением х и ѵ. |
|
|
|
|
|
Ъ { у — Уху)2
(xv)Gy —
N
Расчет тесноты связи между факториальными при знаками X, V и результативным признаком у целесооб разно вести в табличной форме (табл. 7).
|
|
|
|
Т а б л и ц а 7 |
Номер опыта |
У |
« X V |
<«-«XV> |
(«-«х-Х- |
|
||||
• |
|
• |
|
■ |
|
|
|
|
Полезно также для целей анализа рассчитать част ные корреляционные отношения, оценивающие влияние фактора V на у при исключенном х, а также фактора х на у при исключенном ѵ:
|
|
/ о |
2 |
|
<.x)R« = |
± | / |
~ |
(хѵ)Ѵу |
(2.38) |
|
|
|||
(v)R x —■± | / |
(а)0 # |
(хѵ)°у |
(2.39) |
|
|
|
|||
(«Гу
Технику расчета приведенного примера рассмотрим на конкретных цифрах (табл. 8). При этом читатель может выполнить часть расчетов самостоятельно и за полнить пустующие графы. Суммарные значения (пред-
75
Номеропыта )замеса( |
Цементно-водное отношениех |
Модулькрупнос тиV |
Прочностьбетона см{кГ,у 2 |
3 |
Е» |
j!» |
ъ |
|
|
|
|
|
Расчет данных для решения системы нор |
||
1 |
1,953 |
2,38 |
445 |
869,085 |
1059,100 |
3,81420 |
5,6644 |
2 |
2,388 |
2,38 |
399 |
952,812 |
949,62 |
5,70254 |
5,6644 |
3 |
1,910 |
2,49 |
356 |
679,960 |
886,44 |
3,64810 |
6,2001 |
4 |
1,942 |
2,П |
322 |
|
|
|
|
5 |
1,878 |
2,29 |
458 |
|
|
|
|
6 |
1,966 |
2,51 |
448 |
|
|
|
|
7 |
1,870 |
2,60 |
364 |
|
|
|
|
8 |
1,688 |
2,46 |
368 |
|
|
|
|
9 |
1,811 |
2,54 |
425 |
|
|
|
|
10 |
1,763 |
2,60 |
385 |
|
|
|
|
И |
1,938 |
2,13 |
404 |
|
|
|
|
12 |
1,849 |
2,43 |
406 |
|
|
|
|
13 |
1,785 |
2,40 |
396 |
|
|
|
|
14 |
1,836 |
2,53 |
399 |
|
|
|
|
15 |
1,939 |
2,62 |
381 |
|
|
|
|
16 |
2,476 |
2,49 |
463 |
|
|
|
|
17 |
1,769 |
2,36 |
355 |
|
|
|
|
18 |
1,837 |
2,47 |
370 |
|
|
|
|
19 |
1,856 |
2,66 |
317 |
|
|
|
|
20 |
2,121 |
2,46 |
318 |
|
|
|
|
2 |
38,574 48,91 |
7849 |
15 182,640 |
19 163,02 |
75,12808 |
120,0133 |
|
Сред- |
1,9287 2,4455 392,45 |
759,132 |
958,165 |
3,756 |
6,000 |
||
ние зна чения
последняя строчка таблицы) могут служить в этом слу чае в качестве контрольных цифр.
Подставим найденные суммы в систему нормальных уравнений
20а0 + 38,574а! + 48,91а2 = 7849,0,
38,574а0 + 75,12808аі + 94,27844а2 = 15 182,640,
48,91а0 + 94,27844а! -f 120,0133а2 = 19 163,020.
Решив систему уравнений, получим
(h = — 70,7710, öl = 55,4861, а0 = 458,5045.
|
|
|
|
|
Т а б л и ц а 8 |
|
мальпых уравнений |
Расчет «теоретических» значений прочности и |
|||||
|
|
|
корреляционного отношения |
|
||
|
|
н |
Е. |
|
|
|
|
|
о |
|
|
|
|
|
|
□о |
I»- |
|
|
+ |
о |
|
ю |
t - |
О |
|
|
|
О |
=7 |
зГ |
|||
н |
|
ю |
f- |
|||
4,64814 |
198,025 |
108,42 |
168.43 |
398,49 |
46.5 |
2162,25 |
5,68344 |
159,201 |
132,50 |
168.43 |
422,57 |
23.6 |
556,96 |
4,75590 |
129,736 |
105,98 |
176,22 |
388,26 |
32,3 |
1043,29 |
93,27844 |
3115,002 |
27584,30 |
4,71392 |
1379,21 |
|
Уравнение множественной регрессии имеет вид
УXV = 48,5045+ 55,4861 х — 70,771 Оу.
Расчет коэффициентов корреляции (табл. 9):
|
|
|
Т а б л и ц а 9 |
Переменные |
г у |
|
Т |
|
|
|
|
У |
1 |
0,280 |
—0,269 |
X |
0,280 |
1 |
—0,101 |
V |
—0,269 |
—0,101 |
1 |
76 |
77 |
= / |
155,750— 154,017 |
= /Т /7 33 |
= 1,31, |
|||
ax = |
/3 ,7 5 6 — 1,9287я |
= /сЦіЗёГ = 0,19, |
||||
aD= /6 ,0 0 0 —2,445“2= |
/ОДЖГ = |
0,14, |
||||
|
_ 759,132 — 756,918 |
_ |
2,214 |
= + |
0,280, |
|
|
41,64-0,19 |
~ |
7,912 |
|||
|
" |
|
|
|||
r |
958,165 — 959,736 |
_ |
1,571 |
= |
— 0,269, |
|
|
5,830 |
~ |
5,830 |
|
|
|
r |
4,7139 — 4,7166 |
_____0,0027 |
— |
0,101. |
||
0,0266 |
~ |
|
= |
|||
|
0,0266 |
|
|
|||
Расчет корреляционного отношения:
R = |
l / |
1733--_1379121_ = |
A |
V |
1733 |
у о Щ - = 0>453. r
Рассмотренный пример относительно прост и может быть решен с помощью счетной линейки. Однако, ког да надо проанализировать влияние нескольких факто ров на результирующий признак, задача значительно осложняется. Прежде всего надо весьма четко пред ставлять качественные взаимоотношения (взаимосвязь) исследуемых факторов, и только после всестороннего качественного анализа объекта исследования и четкой формулировки целей и задач исследования можно пере ходить к выбору метода и сбору данных применитель но к выбранному методу исследования.
При подготовке анализа процессов горного произ водства с применением методов множественной корре
ляции |
необходимо установить и конкретизировать ос |
|||
новные |
геологические, |
технические, |
организационные |
|
факторы и параметры |
П |
на исследуемый |
||
(х), влияющие |
||||
|
|
і=і |
|
|
процесс или явление (у). |
|
|
||
|
•Чі.) Х\ |
|
Таблица |
10 |
W# л/п |
X, |
х п - 1 |
х п |
|
1 |
|
|
|
|
2 |
|
|
|
|
120 |
|
|
|
|
78
Собранные исходные данные систематизируют и обобщают в таблицах (табл. 10).
По |
систематизированным |
данным определяются: |
1) |
регрессионные уравнения |
(уравнения парной ре |
грессии), выражающие зависимость между у и каждым
из А',- при / = (1, п) в отдельности; |
выражаю |
2) множественное уравнение регрессии, |
|
щее зависимость у от всех ад при і=(1, |
п) одновре |
менно. |
|
Для этого вычисляют средние величины у, х, ух,
у2, а2, ухі и основные отклонения аѵ, аЛ/ при і=(1, п)
по формулам ау= V У1—У2, оХі = У х 2—х2. а также
коэффициенты корреляции и их погрешности по фор
мулам (2.30) |
и (2.31). |
Вычисление коэффициентов корреляции проводят |
|
для каждого |
параметра с каждым из остальных. Ре |
зультаты вычислений сводят в табл. 11 (квадратную матрицу коэффициентов парной корреляции).
|
|
|
|
|
Таблица |
П |
||
|
|
У |
Хі |
Х2 |
Аз |
Х[і—1 |
Хп |
|
У |
\ |
/ |
|
|
|
|
|
|
/ |
\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
*1 |
|
\ |
/ |
|
|
|
|
|
|
/ |
\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
х 2 |
|
|
\ |
/ |
|
|
|
|
|
|
/ |
\ |
|
|
|
|
|
|
|
|
|
|
|
|
||
Хз |
|
|
|
\ |
/ |
|
|
|
|
|
|
|
/ |
\ |
|
|
|
|
|
|
|
|
\ |
/ |
|
|
|
|
|
|
|
/ |
\ |
|
|
х п—1 |
|
|
|
|
|
\ |
/ |
|
|
|
|
|
|
/ |
\ |
|
|
|
|
|
|
|
|
|
||
Хп
Для нахождения множественного уравнения регрес сии можно применить методы множественного регрес
79
смежного анализа только в случае, если регрессионные зависимости і/ и .ѵ,- представляют собой уравнения пря мых (или близких к ним).
Уравнение множественной регрессии изучаемого признака, изменяющегося в зависимости ’от комплекс ного воздействия установленных основных геологиче
ских, технических, организационных |
факторов и пара |
|
метров, имеет следующий вид: |
|
|
У ( х і = х 1—хп) = ± ^ 1Х1 ± |
ЬпХ2 ± |
. . - ± |
± |
± Ьахя. |
(2.40) |
П
X
і~\
чае применения метода на ЭВМ решается система из п алгебраических уравнений первого порядка, матрица которых составлена нз вычисленных коэффициентов корреляции.
Полученные таким образом уравнения множествен ной регрессии позволяют исследовать совокупное влия-
ние многих переменных |
П |
результативный |
факторов .ѵ на |
||
признак у. |
1=1 |
|
При исследовании производственных систем урано- |
||
добывающих предприятий, состояние |
которых может |
|
быть описано системой уравнений первой степени, мо гут быть использованы методы матричного исчисления. Применение методов линейной алгебры (матричного ис числения) для исследования сложных интегрированных систем производства основано на возможности в опре деленных случаях расчленения их структуры на про стейшие линейные зависимости. В этих случаях мате
матическая модель системы может |
быть представлена |
в форме системы матриц, элементы |
которых вычис |
ляются на ЭВМ. |
|
2.ЭЛЕМ ЕНТЫ М А ТРИ Ч Н О ГО И СЧИСЛЕНИЯ
При проектировании и планировании горных работ суммарные затраты часто определяют умноже нием объемов отдельных видов работ на их себестои мость и последующим суммированием полученных про изведений (табл. 12).
80
