- •1.1. Геологические объекты и их свойства
- •Характеристика текстуры руды
- •1.2. Понятие о математическом моделировании
- •2.1. Одномерная статистическая модель
- •Содержание никеля в руде х
- •Расчет начальных моментов случайной величины
- •Расчет центральных моментов случайной величины
- •Частота и частость содержания железа в руде
- •Расчет статистических характеристик по сгруппированным данным (по данным гистограммы рис.2.1)
- •2.2. Законы распределения случайных величин
- •Функции нормального распределения
- •Функция распределения Стьюдента Fk (t) в зависимости от числа степеней свободы k
- •Плотность вероятности распределения Стьюдента fk(t) в зависимости от числа степеней свободы k
- •Коэффициенты вероятности t распределения Стьюдента при заданной вероятности и степени свободы k
- •Коэффициенты вероятности распределения 2 при заданных вероятности и числе степеней свободы k
- •Расчет графика плотности вероятности нормального закона
- •Подбор функции преобразования случайной величины х (содержание TiO2 в магнетите)
- •2.3. Геологические приложения
- •Пример выявления аномальных значений
- •Результат вычисления истинных характеристик
- •3.1. Двухмерная статистическая модель
- •Расчет коэффициента корреляции
- •Расчет параболической зависимости, аппроксимирующей изменение мощности рудного тела
- •Сравнение фактической и расчетной (теоретической) мощности
- •Дисперсии отклонений
- •3.2. Геологические приложения двухмерной статистической модели
- •Расчет случайной погрешности химического анализа
- •Расчет систематической погрешности анализов серы
- •Результаты опытов по измерению радиоактивности препарата
- •Содержание железа в руде и плотность руды
- •Расчет частот n
- •4.1. Многомерная статистическая модель
- •Расчет статистических характеристик многофакторной зависимости
- •Сравнение фактических и расчетных содержаний золота
- •Химические анализы проб руды
- •Матрица коэффициентов корреляции
- •Состав магнетита, %
- •Преобразованная матрица коэффициентов корреляции
- •Матрица коэффициентов корреляции
- •Расчет первого собственного числа
- •Расчет матрицы коэффициентов первого фактора
- •Первая остаточная матрица коэффициентов корреляции
- •Факторные нагрузки
- •Нормированный состав магнетита
- •Главные компоненты z
- •Средний состав групп магнетита, %
- •Состав циркона, %
- •Сложные признаки
- •Результаты распознавания объектов
- •Результаты измерения мощности рудного тела
- •Данные по скважинам
- •5.4. Вероятностные модели
- •Результаты измерения мощности рудного тела
- •Значения пространственной переменной
- •Содержание молибдена, %
- •Сглаживание исходных данных
- •Результаты измерений пространственной переменной и некоторые вычисленные данные
- •Результаты расчета всех гармоник
- •Расчет периодической изменчивости
- •Дисперсия различных длин волн
- •5.5. Основы геостатистики
- •Содержание золота в пробах, г/т
- •Данные для кригинга, м
- •Значения вариограммы (h)
- •Реестр разведочных скважин
- •Замеры искривлений скважин
- •Банк данных геологической документации
- •Банк опробования рядовых проб
- •Координаты пунктов искривлений
- •Банк рудных пересечений
- •Банк геологических границ
- •Глава 1. Общие сведения 7
- •Глава 2. Одномерная статистическая модель и ее применение в геологии 23
- •Глава 3. Двухмерная статистическая модель и ее применение в геологии 78
- •Глава 4. Многомерная статистическая модель и ее применение в геологии 115
- •Глава 5. Математическое моделирование пространственных геологических закономерностей 159
- •Глава 6. Основы математического моделирования месторождений 203
Состав магнетита, %
Номер пробы |
TiO2 |
MnO |
V2O5 |
SiO2 |
Al2O3 |
MgO |
CaO |
|
1 |
0,25 |
0,17 |
0,31 |
0,23 |
0,30 |
0,06 |
0,06 |
|
2 |
0,32 |
0,15 |
0,26 |
0,25 |
0,65 |
0,12 |
0,12 |
|
3 |
0,30 |
0,13 |
0,25 |
0,12 |
0,52 |
0,06 |
0,02 |
|
4 |
0,28 |
0,10 |
0,29 |
0,10 |
0,46 |
0,06 |
0,03 |
|
5 |
0,33 |
0,08 |
0,27 |
0,24 |
0,42 |
0,10 |
0,04 |
|
6 |
0,12 |
0,08 |
0,23 |
0,16 |
0,29 |
0,10 |
0,02 |
|
7 |
0,33 |
0,22 |
0,25 |
0,20 |
0,65 |
0,12 |
0,02 |
|
8 |
0,43 |
0,15 |
0,39 |
0,46 |
0,40 |
0,29 |
0,15 |
|
9 |
0,36 |
0,17 |
0,28 |
0,50 |
0,70 |
0,06 |
0,04 |
|
10 |
0,46 |
0,12 |
0,20 |
0,12 |
0,70 |
0,13 |
0,02 |
|
11 |
0,39 |
0,12 |
0,27 |
0,06 |
0,67 |
0,12 |
0,02 |
|
12 |
0,24 |
0,12 |
0,27 |
0,16 |
0,45 |
0,05 |
0,02 |
|
13 |
0,45 |
0,07 |
0,22 |
0,16 |
0,63 |
0,04 |
0,02 |
|
14 |
0,54 |
0,19 |
0,30 |
0,25 |
0,40 |
0,17 |
0,17 |
|
15 |
0,35 |
0,23 |
0,22 |
0,17 |
0,25 |
0,06 |
0,02 |
|
16 |
0,35 |
0,07 |
0,24 |
0,12 |
0,43 |
0,12 |
0,03 |
|
17 |
0,10 |
0,06 |
0,19 |
0,16 |
0,36 |
0,09 |
0,02 |
|
18 |
0,07 |
0,08 |
0,29 |
0,11 |
0,36 |
0,06 |
0,05 |
|
19 |
0,32 |
0,15 |
0,25 |
0,15 |
0,60 |
0,10 |
0,04 |
|
20 |
0,43 |
0,27 |
0,35 |
0,07 |
0,52 |
0,08 |
0,02 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,321 |
0,136 |
0,266 |
0,190 |
0,483 |
0,099 |
0,046 |
|
|
0,119 |
0,057 |
0,047 |
0,111 |
0,138 |
0,055 |
0,044 |
|
Далее
необходимо умножить вектор V или W
на такой множитель, чтобы сумма
квадратов членов вектора была равна
собственному числу .
Нетрудно определить, что к вектору V
нужно применить множитель
а к вектору W множитель
В рассматриваемом примере
Умножая все члены вектора V на
множитель
получим первый собственный вектор Ф1
(первую факторную нагрузку),
Таблица 4.5
Матрица коэффициентов корреляции (rзнач = 0,197)
Компонент |
Si |
Ti |
Al |
Fe |
Mn |
Mg |
Ca |
Na |
K |
P |
Ппп |
CO2 |
S |
Cu |
Zn |
Pb |
Co |
Ni |
Ag |
Si |
1,000 |
0,779 |
0,902 |
–0,900 |
0,604 |
0,417 |
0,596 |
0,781 |
0,582 |
0,300 |
–0,158 |
–0,136 |
–0,612 |
–0,299 |
0,013 |
–0,116 |
–0,715 |
–0,209 |
–0,045 |
Ti |
0,779 |
1,000 |
0,765 |
–0,676 |
0,583 |
0,317 |
0,417 |
0,625 |
0,487 |
0,228 |
–0,250 |
–0,226 |
–0,497 |
–0,145 |
0,080 |
–0,050 |
–0,576 |
–0,250 |
0,021 |
Al |
0,902 |
0,765 |
1,000 |
–0,818 |
0,588 |
0,433 |
0,447 |
0,808 |
0,613 |
0,270 |
–0,255 |
–0,255 |
–0,662 |
–0,245 |
0,014 |
–0,108 |
–0,679 |
–0,255 |
–0,084 |
Fe |
–0,900 |
–0,676 |
–0,818 |
1,000 |
–0,411 |
–0,421 |
–0,798 |
–0,692 |
–0,556 |
–0,266 |
–0,181 |
–0,196 |
0,584 |
0,418 |
0,082 |
0,146 |
0,710 |
0,242 |
0,036 |
Mn |
0,604 |
0,583 |
0,588 |
–0,411 |
1,000 |
0,354 |
0,565 |
0,387 |
0,440 |
0,222 |
–0,062 |
–0,244 |
–0,419 |
–0,201 |
0,136 |
0,057 |
–0,475 |
–0,290 |
0,034 |
Mg |
0,417 |
0,317 |
0,433 |
–0,421 |
0,354 |
1,000 |
0,151 |
0,142 |
0,242 |
0,150 |
–0,112 |
–0,244 |
–0,371 |
–0,185 |
0,030 |
0,013 |
–0,318 |
–0,022 |
–0,039 |
Ca |
0,596 |
0,417 |
0,447 |
–0,798 |
0,565 |
0,151 |
1,000 |
0,356 |
0,344 |
0,174 |
0,488 |
0,544 |
–0,377 |
–0,415 |
–0,104 |
–0,127 |
–0,517 |
–0,241 |
–0,009 |
Na |
0,781 |
0,625 |
0,808 |
–0,692 |
0,387 |
0,142 |
0,356 |
1,000 |
0,319 |
0,204 |
–0,175 |
–0,123 |
–0,542 |
–0,195 |
–0,027 |
–0,121 |
–0,586 |
–0,268 |
0,015 |
K |
0,582 |
0,487 |
0,613 |
–0,556 |
0,440 |
0,242 |
0,344 |
0,319 |
1,000 |
0,272 |
–0,082 |
–0,078 |
–0,398 |
–0,189 |
–0,059 |
–0,091 |
–0,454 |
–0,141 |
–0,020 |
P |
0,300 |
0,228 |
0,270 |
–0,266 |
0,222 |
0,150 |
0,174 |
0,204 |
0,272 |
1,000 |
–0,110 |
–0,117 |
–0,108 |
0,014 |
–0,057 |
–0,002 |
–0,123 |
0,072 |
–0,071 |
Ппп |
–0,158 |
–0,250 |
–0,255 |
–0,181 |
–0,062 |
–0,112 |
0,488 |
–0,175 |
–0,082 |
–0,110 |
1,000 |
0,951 |
0,126 |
–0,267 |
–0,275 |
–0,091 |
0,006 |
–0,015 |
0,049 |
CO2 |
–0,136 |
–0,226 |
–0,255 |
–0,196 |
–0,244 |
–0,244 |
0,544 |
–0,123 |
–0,078 |
–0,117 |
0,951 |
1,000 |
0,084 |
–0,321 |
–0,256 |
–0,099 |
–0,017 |
–0,079 |
0,058 |
S |
–0,612 |
–0,497 |
–0,662 |
0,584 |
–0,419 |
–0,371 |
–0,377 |
–0,542 |
–0,398 |
–0,108 |
0,126 |
0,084 |
1,000 |
0,549 |
0,012 |
0,041 |
0,832 |
0,349 |
0,077 |
Cu |
–0,299 |
–0,145 |
–0,245 |
0,418 |
–0,201 |
–0,185 |
–0,415 |
–0,195 |
–0,189 |
0,014 |
–0,267 |
–0,321 |
0,549 |
1,000 |
0,116 |
0,079 |
0,482 |
0,213 |
0,105 |
Zn |
0,013 |
0,080 |
0,014 |
0,082 |
0,136 |
0,030 |
–0,104 |
–0,027 |
–0,059 |
–0,057 |
–0,257 |
–0,256 |
0,012 |
0,116 |
1,000 |
0,604 |
0,072 |
0,006 |
0,576 |
Pb |
–0,116 |
–0,050 |
–0,108 |
0,146 |
0,057 |
0,013 |
–0,127 |
–0,121 |
–0,091 |
–0,002 |
–0,091 |
–0,099 |
0,041 |
0,079 |
0,604 |
1,000 |
0,128 |
0,114 |
0,369 |
Co |
–0,715 |
–0,576 |
–0,679 |
0,710 |
–0,475 |
–0,318 |
–0,517 |
–0,586 |
–0,454 |
–0,123 |
0,006 |
–0,017 |
0,832 |
0,482 |
0,072 |
0,128 |
1,000 |
0,365 |
0,105 |
Ni |
–0,209 |
–0,250 |
–0,255 |
0,242 |
–0,290 |
–0,022 |
–0,241 |
–0,268 |
–0,141 |
0,072 |
–0,015 |
–0,079 |
0,349 |
0,213 |
0,006 |
0,114 |
0,365 |
1,000 |
–0,081 |
Ag |
–0,045 |
0,021 |
–0,084 |
0,036 |
0,034 |
–0,039 |
–0,009 |
0,015 |
–0,020 |
–0,071 |
0,049 |
0,058 |
0,077 |
0,105 |
0,576 |
0,369 |
0,105 |
–0,081 |
1,000 |
Таблица 4.6
