Предварительно преобразуем исходные значения признаков выборочной совокупности к нормированному и центрированному виду (таблица 6).

Таблица 6

Нормированные значения элементов в пробах

№№ проб	Cu	Zn	V	Cорг	Zr	Ti
1	2	3	4	5	6	7
1	-0,6658	-0,8055	-1,0644	-0,6911	3,7887	1,7853
2	-0,6867	-0,8055	-1,1534	-0,6911	0,5186	1,7853
3	2,4376	2,6045	1,3405	3,2703	-0,5559	-1,8371
4	2,4376	2,6045	1,3405	1,0034	-0,5559	-0,8021
5	-0,6867	-0,7707	-0,4409	-0,6571	0,0514	0,7504
6	0,3478	-0,5272	-0,4409	-0,1301	-0,7428	-1,3196
7	0,3478	0,1688	1,3405	1,0034	-0,4158	-0,8021
8	-0,6940	-0,5272	-1,0644	-0,6798	0,5186	1,2678
9	-0,5926	-0,5272	-1,0644	-0,6571	-0,4158	-0,0259
10	-0,5926	-0,5272	0,4498	-0,5268	0,5186	0,7504
11	-0,3836	-0,7011	-0,4409	-0,5268	-0,7428	-1,3196
12	-0,3836	-0,5272	0,4498	-0,1301	-0,4158	-0,5434
13	0,3478	0,1688	1,3405	1,0034	0,0514	-0,0259
14	-0,6658	-0,1792	-1,0644	-0,6401	-0,4158	0,2329
15	1,3927	0,8647	-0,4409	-0,1301	0,5186	0,7504
16	-0,6658	-0,8403	-0,4409	-0,5268	-0,7428	-0,8021

Продолжение таблицы 6

1	2	3	4	5	6	7
17	0,3478	0,1688	1,3405	1,0034	0,0514	0,2329
18	-0,6658	0,1688	-0,8862	-0,6401	-0,6494	-0,2846
19	-0,5926	0,1688	-0,4409	-0,5268	-0,4158	-0,0259
20	-0,3836	-0,1792	1,3405	-0,1301	0,0514	0,2329

Корреляционная матрица нормированных признаков R представлена в таблице 7, элементы которой могут быть также вычислены с помощью компьютерных программ в Exel или Statistica.

Таблица 7

Корреляционная матрица

Признаки	Коэффициенты корреляции
	Cu	Zn	V	Cорг	Zr	Ti
Cu	1,00	0,92	0,60	0,82	-0,20	-0,48
Zn	0,92	1,00	0,56	0,79	-0,25	-0,44
V	0,60	0,56	1,00	0,77	-0,24	-0,44
Coрг	0,82	0,79	0,77	1,00	-0,23	-0,54
Zr	-0,20	-0,25	-0,24	-0,23	1,00	0,73
Ti	-0,48	-0,44	-0,44	-0,54	0,73	1,00

Примечание: коэффициенты корреляции являются значимыми (0,05) при их абсолютном значении не менее 0,44.

Анализ коэффициентов корреляции показывает, что накопление элемента Cu происходило одновременно с Zn (0,92), Cорг (0,82) и V (0,60), но биполярно и не связано с накоплением элементов Zr и Ti. Содержание элементов Zr и Ti коррелирует между собой (0,73). Эту связь между элементами можно проиллюстрировать на корневой диаграмме признаков для выявления структуры признаков и предварительного анализа факторов.

Рисунок 5 Диаграмма связи признаков

Основной математический метод выделения факторов и их нагрузок основан на нахождении собственных чисел и собственных векторов корреляционной матрицы.

Таблица 8

Факторные нагрузки, собственные значения и веса признаков

Признаки	Факторные нагрузки
	F1	F2	F3	F4	F5	F6
Cu	0,893	0,291	0,269	-0,005	0,111	-0,180
Zn	0,876	0,266	0,331	-0,146	0,024	0,173
V	0,783	0,154	-0,574	-0,141	0,118	0,015
Coрг	0,914	0,228	-0,123	0,143	-0,276	-0,018
Zr	-0,478	0,827	-0,056	0,275	0,085	0,045
Ti	-0,724	0,590	-0,006	-0,336	-0,110	-0,049
Собственные значения	3,77	1,26	0,53	0,25	0,12	0,07
Веса факторов, %	62,78	21,06	8,83	4,17	2,04	1,12

Примечание: факторные нагрузки являются значимыми (0,05) при их абсолютном значении не менее 0,44

Величины собственных чисел (факторные нагрузки ) и векторов (веса факторов,%) факторной нагрузки признаков табл.7 показывают, что песчано-глинистые отложения характеризуются на 62,78% фактором F₁ и на 21,06% фактором F₂.

На первом этапе анализа необходимо определить минимальное число факторов, адекватно воспроизводящие наблюдаемые корреляции. Анализ признаковой структуры преимущественного фактора F₁ (62,78%) показывает, что нагрузка этого фактора имеет сложный характер и определяется как содержанием элементов Cu (0,893), Zn (0,876), Cорг (0,914), V (0,783) так и значимой отрицательной связью Ti (-0,724).

Признаковая структура однополярного фактора F₂ (21,06%) характеризует влияние элементов Zr (0,827) и Ti (0,590) на свойства отложения. Влияние элемента V (-0,574) выявляется фактором F₃ (8,83%). Наиболее часто число выделяемых факторов определяется количеством собственных чисел больше единицы. В нашем случае это первые два фактора. Другой критерий определяется через графическое изображение собственных чисел. Выделение заканчивается на том факторе, после которого исследуемая зависимость близка к горизонтальной линии. Как видно на рисунке 6, выделяется не более 4 факторов. Окончательное решение должно базироваться на приемлемости с точки зрения научных представлений в данной области.

Рисунок 6

Следующим шагом с помощью процедуры вращения выявляются наиболее легко интерпретируемые факторы. Существует три подхода к этой проблеме. Геометрический подход используется, когда число факторов не более двух и имеются отдельные скопления (кластеры) признаков. Аналитический – выбирается критерий, на основе которого производится вращение (ортогональное или косоугольное) осей.

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

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

Таблица 9

Факторные нагрузки, собственные значения и веса признаков

после вращения

Признаки	Факторные нагрузки
	F1	F2
Cu	0,929	0,140
Zn	0,903	0,155
V	0,769	0,214
Cорг	0,920	0,206
Zr	-0,056	-0,953
Ti	-0,382	-0,852
Собственные значения	3,263	1,767
Веса факторов, %	54,4	29,5

Рисунок 7

Рисунок 8

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

2. Выделяем три фактора и видим, что элемент V по своим свойствам выделяется из первого фактора в отдельный третий фактор F₃, а элемент Cорг распределилось по факторам F₁ и F₃ . В новых координатах можно провести дополнительные исследования свойств залежи.

Таблица 10

Факторные нагрузки, собственные значения и веса факторов

трех признаков после вращения

Признаки	Факторные нагрузки
	F1	F2	F3
Cu	0,921	0,149	0,290
Zn	0,932	0,168	0,225
V	0,325	0,166	0,912
Cорг	0,698	0,188	0,617
Zr	-0,050	-0,954	-0,044
Ti	-0,299	-0,847	-0,257
Собственные значения	2,403	1,741	1,416
Веса факторов, %	40,0	29,0	23,6

3. Задаем требование – выделить четыре фактора и произвести вращение системы координат.

Таблица 11

Факторные нагрузки, собственные значения и веса факторов

четырех признаков после вращения

Признаки	Факторные нагрузки
	F1	F2	F3	F4
Cu	0,916	0,085	0,276	0,181
Zn	0,942	0,159	0,233	0,052
V	0,336	0,140	0,919	0,091
Cорг	0,677	0,056	0,581	0,354
Zr	-0,080	-0,981	-0,084	-0,124
Ti	-0,263	-0,631	-0,202	-0,691
Собственные значения	2,373	1,415	1,361	0,662
Веса факторов, %	39,5	23,6	22,7	11,0

В новой системе координат четвертый фактор F₄ определяет накопление Ti в отложениях, однако с такой же долей участия этот элемент входит во второй фактор F₂. Видимо, такое положение происходит из-за избыточности четвертого фактора F₄. Таким образом, как следует и из графического представления и численных вычислений, два-три фактора наиболее полно характеризуют признаковую структуру данных отложений (рис.9)

Рисунок 9

Вывод

Анализ признаковой структуры трех выявленных факторов позволяет предполагать, что фактор F₁ может быть интерпретирован как действие сульфидного диагенеза в осадке, вызванного разложением органики в анаэробной среде и развитием процессов бактериальной сульфатредукции. С этим процессом и связано образование сульфидов меди, цинка, накопление Сорг. Кроме того, органика выступает еще и как концентратор ванадия, что описывается фактором F₃. Фактор F₂ может быть интерпретирован как терригенное накопление в осадке аллотигенных минералов титана и циркония.

ЗАКЛЮЧЕНИЕ

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