- •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
Расчет случайной погрешности химического анализа
Номер пробы ni |
Содержание меди, % |
Разность xi – yi |
Квадрат разности (xi – yi)2 |
||
Основные пробы xi |
Контрольные пробы yi |
||||
1 |
2,74 |
2,70 |
0,04 |
0,0016 |
|
2 |
2,14 |
2,44 |
–0,30 |
0,0900 |
|
3 |
2,33 |
2,19 |
0,14 |
0,0196 |
|
4 |
2,57 |
2,54 |
0,03 |
0,0009 |
|
5 |
2,16 |
2,24 |
–0,08 |
0,0064 |
|
6 |
1,27 |
1,21 |
0,06 |
0,0036 |
|
7 |
1,00 |
1,23 |
–0,23 |
0,0529 |
|
8 |
0,95 |
0,59 |
0,36 |
0,1296 |
|
9 |
1,72 |
1,28 |
0,44 |
0,1936 |
|
10 |
2,06 |
1,76 |
0,30 |
0,0090 |
|
11 |
1,06 |
1,43 |
–0,37 |
0,1369 |
|
12 |
1,83 |
1,83 |
0,00 |
0,0000 |
|
13 |
2,13 |
1,81 |
0,32 |
0,1024 |
|
14 |
3,04 |
3,16 |
–0,12 |
0,0144 |
|
15 |
1,52 |
1,34 |
0,18 |
0,0324 |
|
16 |
1,48 |
1,63 |
–0,15 |
0,0225 |
|
17 |
0,78 |
0,82 |
–0,04 |
0,0016 |
|
18 |
0,92 |
0,60 |
0,32 |
0,1024 |
|
19 |
2,17 |
2,62 |
–0,45 |
0,2025 |
|
20 |
2,96 |
2,56 |
0,40 |
0,1600 |
|
21 |
1,45 |
1,79 |
–0,34 |
0,1156 |
|
22 |
1,82 |
1,83 |
–0,01 |
0,0001 |
|
23 |
2,51 |
2,29 |
0,22 |
0,0484 |
|
24 |
1,70 |
2,18 |
–0,48 |
0,2304 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Сумма |
44,31 |
44,07 |
0,24 |
1,7578 |
|
Среднее |
1,85 |
1,84 |
– |
– |
|
Абсолютная случайная погрешность 0,191 |
|||||
Относительная случайная погрешность 10,4 % |
|||||
3.2.4. Внешний контроль химических анализов
В геологической практике принято регулярно оценивать систематическую погрешность измерений. Наиболее часто определяется погрешность опробования или ее составная часть – погрешность химического анализа, для чего выполняется внешний контроль анализов. Главное требование при изучении систематических погрешностей (или систематических расхождений) – неравноточность основных и контрольных измерений. Для выполнения этого условия основные пробы посылают в одну лабораторию, а контрольные – в другую, где анализ выполняют, как правило, по более совершенной методике. Сравнение анализов основных и контрольных проб позволяет оценить систематическую погрешность анализов.
Обозначим
xi – данные основных проб,
yi – данные контрольных
проб. Для выявления систематической
погрешности применяются графический
и аналитический методы. При графическом
анализе проверяется расположение точек
графика. При отсутствии систематической
погрешности они должны располагаться
вдоль биссектрисы у = х. Из-за
наличия неизбежных случайных погрешностей
точки рассеиваются около биссектрисы,
образуя облако. Если о
блако
точек смещено относительно биссектрисы,
то можно предполагать, что в основных
(а иногда и в контрольных) пробах имеется
систематическая погрешность (рис.3.10).
при
,
(3.34)
где n – число контрольных проб.
Если критерий t окажется больше допустимого tдоп, то систематическая погрешность доказана. Допустимое значение критерия tдоп находят на основе распределения Стьюдента при вероятности = 0,05 и числе степеней свободы k = n – 1 (см. табл.2.10). При увеличении числа контрольных проб распределение Стьюдента приближается к нормальному и в пределе tдоп = 1,960.
Метод выявления систематической погрешности, основанный на сравнении средних значений и , обладает существенным недостатком: равенство средних еще не гарантирует наличия систематических расхождений при низких и высоких содержаниях. Последние могут быть направлены в разные стороны и в среднем компенсируют друг друга.
Более рациональной является методика, основанная на оценке коэффициентов уравнения регрессии y = ax + b. При отсутствии систематической погрешности должны выполняться условия а = 1 и b = 0. Поэтому вычисляют уравнение регрессии между содержаниями в основных и контрольных пробах и находят соответственно два критерия:
