
книги из ГПНТБ / Применение ЦВМ и средств вычислительной техники в геологии и геофизике [сборник]
..pdf(а), дисперсия (Ь2) и стандартное отклонение >(6) при помощи максимально правдоподобной оценки [2]
а = \ 0 ^ х е1{\ |
t(t + 1) I |
+ 22t + 21) |
I |
|
n |
ört2 |
/ |
||
|
На печать результаты выдаются в последовательности:
и» lg X, S]g , Sij,, -Aig, 3(^4ig), |
-^ig • ^(•Aig), |
5(*^ig)i |
Э1е:а{ЭІЯ), hu, |
V, a, b\ b. |
|
Условиями нормальности распределения будет соблюдение неравенств:
А |
< 3, |
Э |
< 3 . |
|
а(А) |
Ф ) |
|||
|
|
Аналогичные условия необходимы для соблюдения логариф мической нормальности.
Наряду с анализом статистических оценок параметров распределения элементов в природных объектах возникает необходимость выяснения зависимости элементов друг от дру га и построения парагенетических ассоциаций.
В качестве оценки силы парагенетических связей между элементами и минералами в геохимии используются различ ные модификации коэффициентов корреляции [5].
Подготовка материала для обработки на ЭВМ аналогич на выше описанной. Если используется перфорационный ма териал после получения статистических оценок функции рас пределения, то из него убираются признаки окончания масси вов 777 7777 7777 7777 и КЪ.
Дополнительно пробивается карточка информации следую
щего вида: |
|
|
|
|
0 |
0500 |
0 |
0 |
КА |
0 |
0 |
я |
0 |
|
0 |
0 |
т |
0 |
|
0 |
0 |
1 |
0, |
|
где я — число анализов, т — число признаков, I — число мас сивов.
Подготовительный материал обрабатывается по специаль ной программе, которая из данного массива исключает те пробы, в которых отсутствует хотя бы один анализ, т. е. произ водится построение новой матрицы без —1.
Преобразованная таким образом матрица используется для подсчетов коэффициентов корреляции по следующим фор мулам [6, 7].
130
1. П арны й коэф ф ициент корреляции
В качестве статистического критерия для проверки гипоте зы Н0: Q= 0 вычисляется величина
которая в условиях нулевой гипотезы распределена по закону Стыодента с п-2 степенями свободы. При t> t табл, при за
данном уровне значимости |
(q) гипотеза Н0: Q— 0 отвергает |
ся и зависимость считается |
установленной. |
Если п мало, |
то и в качестве оценки Н0: q = 0 |
используют |
|
величину |
|
|
|
где z — преобразование Фишера |
|
||
Если |
|«|>1,96, |
то гипотеза H0:q = 0 может быть |
отклонена |
при |
q = 0,05. |
|
|
Результаты подсчета парного коэффициента корреляции с критериями t и и выдаются в виде таблиц на широкую печать.
12. Частный коэффициент корреляции
где Pljt Pih Pjj — алгебраические дополнения корреляци онной матрицы Р.
Необходимость подсчета частного коэффициента корреляции часто возникает при анализе различных парагенетических ассоциаций элементов с целью выявления ложных гі;- , причем, Рі / дает возможность выделить действительную корреляцию между парой элементов, исключив влияние всех остальных, а P,j (k) позволяет оценить последовательно влияние каждого элемента на данную пару.
3. Множественный коэффициент корреляции.
9* |
131 |
гДІ, 2. . . . і — 1, і + 1 . . . т)~ j/" 1
где
Этот коэффициент дает оценку связи одного элемента одно временно со всей совокупностью.
В качестве примера реализации выше описанных программ при решении геохимических задач приведем выделение в од нородных геологических объектах геохимического фона по степени битуминизации (ß). Предполагается, что содержание битумоидов в осадочных породах обусловлено наличием орга нического вещества и вторичного аллохтонного битумоида. Так как при решении этой задачи необходимо учитывать все факторы, обусловливающие процесс автохтонного битумообразования, исследуемые образцы с помощью перфокартотеки были разбиты на однородные по условиям битумообразования совокупности: по литологии, степени карбонатное™, составу органического вещества, глубине залегания и окислительновосстановительным условиям.
Подсчет парных коэффициентов корреляций между орга ническим углеродом (Сорг ) и битумоидом (Б) в песчаниках, алевролитах и глинах в большинстве случаев оказался очень низким (т<0,5), что указывает на слабую взаимосвязь иссле
дуемых компонентов. Однако |
t> t табл, дают возможность |
|
предположить, что в данном |
случае хотя связь и низка, |
но |
реальна. Затушевывание взаимосвязи (соответственно и |
по |
нижение коэффициента корреляции) между органическим ве ществом и битумоидом обусловлено в первую очередь нали чием миграционного битумоида. Для выявления пород, содер жащих аллохтонный (аномальный, не обусловленный содержа нием в породе органическим веществом) битумоид, необходи мо определить для данной однородной геологической совокуп ности геохимический фон по степени битуминозное™. Предварительно проверялась гипотеза о непротиворечивости распределения коэффициента битуминозное™ пород (ß) нор мальному и логнормальному закону.
Анализ функции распределения ß показал, что в большин стве случаев закон распределения ß не противоречит нормаль ному. В случае логнормального распределения исследовались максимально правдоподобные оценки среднего значения и
132
стандартного отклонения в предложении, что в исследуемых случаях возикиовение логнормального распределения скорее обусловлено аналитическими данными, чем распределением этой величины в природном объекте.
Геохимический фон в пределах одной выделенной совокуп ности определялся последовательно по формуле:
фон |
:<ß» |
= X ± t jS |
i — 1, 2, 3, 4. |
*і = 1, |
*2 = |
0,75, /з = |
0,5, *4=0,25. |
Подобный выбор *; обусловлен тем, что стандартное отклоне ние в выделенных совокупностях в большинстве случаев со размерно со средним значением и даже превосходит его. Для каждой выделенной таким образом новой совокупности под считывался гсорг_ £ - Анализ материалов показал, что им так
же определялся фон, т. е., чем меньше было t г, тем выше коэф фициент корреляции. В различных однородных совокупностях гСорг_Б достигал высоких значений при различных tt. Для разделения пород с автохтонным и аллохтонным битумоидом выбирался такой размах фона, при котором і'сорг—Б достигал 0,7—0,9 при t > i табл. Такой уровень Г£орг_£, подсчитанный по
Данным А. П. Лисицына [8], характеризует взаимосвязь орга нического вещества и битумоида в современных осадках Ку рило-Камчатской впадины.
Выделенные подобным образом автохтонные и аллохтон ные битумоиды различаются и по качественному составу: аллохтонные, в отличие от автохтонных, имеют более высокую степень восстановленности, нейтральности, содержат больше масел, а в элементарном составе хлороформенной фракции— больше углерода, водорода и меньше гетероатомов (A+0 + S).
В качестве другого примера использования коэффициен тов корреляции приведем построение парагенетических ас социаций химических элементов в песчано-алевритовых гли нах апта междуречья Урал—Волга. Известно, что распределе ние химических элементов в породе зависит от форм их ми грации, которые в свою очередь обусловлены физико-геогра фическими условиями седиментации [9]. Накоплен некоторый материал по построению парагенетических ассоциаций хими ческих элементов в различных типах пород (магматических, метаморфических, осадочных) на основе корреляционного ана лиза, который показывает, что в условиях преобладания на водоразделах химической денудации над механической в оса дочных породах намечается распад группы железа, обособле-
133
а
3
Значения парных коэффициентов корреляции и статистического критерия значимости (и)
см |
|
+ 0 .0 96 |
1 0.858 |
+ 0 .0 22 |
0.195 |
+ 0 .3 03 |
2.78 |
I[-0 .0 80 |
0.719 |
—0.128 |
1.14 |
—0.220 |
1,99 |
|
|
||||||||||||
- |
'S, |
+0.112 |
|
іО |
to |
+0.066 |
0.587 |
—0.124 |
|
- 0 .2 0 2 |
|
—0.209 |
|
& |
1.00 |
ю |
© |
1.10 |
1.82 |
1.88 |
|||||||
|
О |
тг |
|||||||||||
|
и |
О о |
|||||||||||
|
|
+ |
|
||||||||||
о |
»0 |
О |
о |
+0.183 |
1.64 |
+0.166 |
1.49 |
— 0.146 |
1.31 |
— 0.171 |
1.54 |
—0.224 |
2.03 |
|
О; |
іО |
to |
||||||||||
|
U- |
О |
Th |
||||||||||
|
|
О |
О |
||||||||||
|
|
і |
|
||||||||||
© |
<u |
1660 - ог и |
О о |
© р |
05 |
ю |
|
|
-0 .0 4 6 |
0.410 |
|||
“ |
оо |
О |
-0 .1 51 |
1.35 |
|||||||||
|
U . |
о |
— |
|
|
© |
о |
||||||
|
|
Ö Ö |
|
|
©* |
||||||||
|
|
-1- |
■ + |
~ |
Ч- |
|
|||||||
со |
|
+0.312 |
|
<Г5 |
|
+0.171 |
1 53 |
+0.041 |
0.371 |
ГГ СО |
+0.208 |
1.88 |
|
І) |
2.87 |
іС |
|
О |
© |
||||||||
t-ч со |
© |
||||||||||||
|
U-. |
О |
~ |
© |
© |
||||||||
«О |
|
-і- |
|
1 |
|
||||||||
|
+0.125 |
|
—0.004 |
0.037 |
, 2000+ |
0 990 |
+0.029 |
0.263 |
1+0.367 |
|
—01.401 |
|
|
<D |
|
1.12 |
3.42 |
3.78 |
|||||||||
S |
|
||||||||||||
О) |
|
||||||||||||
О) |
|
0.1U |
|
-0.051 |
0.455 |
-0 .072 |
0.649 |
-0 .1 4 9 |
|
СО |
|
|
|
о |
cs |
1.02 |
1.34 |
Th 1- |
|
|
|||||||
іО |
|
|
|
||||||||||
|
CO |
© іо |
|
|
|||||||||
|
|
4- |
|
|
|
||||||||
|
|
СО© |
-0 .1 1 7 |
1.05 |
-0 .1 3 9 |
1.24 |
-0.0Э 2 |
( 0.824 |
|
|
|
|
|
ю |
со |
О |
-** |
|
|
|
|
||||||
|
о |
о |
|
|
|
|
|||||||
|
|
1 |
|
|
|
|
|
||||||
|
|
-0 .0 3 5 |
0.760 |
-0.061 |
0.649 |
О |
|
|
|
|
|
|
|
|
|
см со |
|
|
|
|
|
|
|||||
|
** |
см © |
|
|
|
|
|
|
|||||
|
© ^ |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
||||||
|
|
+ |
|
|
|
|
|
|
|
||||
СО |
Pf. |
+0.269 |
2.45 |
+ 0 118 |
1.05 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
||||
см |
Си |
0.270 |
2.47 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
- |
>■ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
|
СМ |
ео |
|
Tf |
|
lO |
© |
|
.049 |
-443 |
—0 |
0 |
|
I |
—0.032 |
0.288 |
.1 5 7 |
1.41 |
-0 |
|
+0.012 |
I 0.108 |
0.127 |
1.13 |
+ |
|
и-
|
|
|
|
|
. |
|
—0.109 |
0.973 |
—0.052 |
0.468 |
+ 0 .4 4 2 |
! 4.23 + 0 .6 9 8 |
1 7-68 |
+ 0 .082 |
0.735 |
—0,113 |
1.00 |
+ 0.567 |
5.71 |
|
.0 9 4 |
0.843 |
.1 7 0 |
1.52 |
|
|
|
-0 |
|
-0 |
|
|
|
|
0+ Г 60 |
8260 |
|
|
|
|
|
©
00 ©
134
Значения множественных и частных коэффициентов корреляций
с ч
—
о
и
о
t u
05 CD
Cu
00 Cu
Л
н
я t-ч n-H
о -J
S
<т> о cd CQ
юCO
со 2
СЧ 3
и
^ ч >
1
cfe |
S - |
о |
|
т-ч |
со |
г о |
со |
сО |
Г5 1> Тч- |
||
CO |
«o |
со |
с ч |
t o |
г—< |
< . |
05 |
N |
* |
О |
|
о |
о |
•*+ |
г*4 |
о |
f-ч |
о |
о |
со СО |
|||
о |
о |
о |
о |
|
о |
о |
о |
|
О |
О |
О |
1 |
1 |
4 - |
1 |
+ |
1 |
1 1 + |
1 |
+ |
- Г |
||
СЧ |
Ю |
ОТ) |
|
г - |
ю |
сч |
т#* |
СЧ |
ю |
с ч |
|
СО Г"- |
по |
|
ь |
СЧ |
1-0 |
СО |
Ю ю |
|
|||
О» |
о |
со г -“ г - |
^Ч |
1— |
о |
тг р |
|
||||
о о о о о о* о о о о о |
|
||||||||||
+ |
1 1 |
+ |
1 |
1 |
+ |
+ |
1 |
+ |
+ |
|
|
Ю О) |
00 |
ю |
05 |
о |
|ч |
05 |
Г ч |
1*4 |
|
|
|
со С-» |
г - |
|
0 5 |
СО |
о |
о |
о |
|
|
||
<м »-* |
СЧ |
о о |
о |
о |
о |
t-ч |
|
|
|||
о |
о |
о |
о |
о |
о |
о |
о |
|
|
||
1 |
+ |
4 |
+ |
+ |
1 |
1 |
1 |
1 |
+ |
|
|
Г- -Ч г—ч h> |
о |
со со |
оо |
о |
|
|
|
||||
05 |
СЧ |
о |
05 |
со |
С55 |
|
Ч.) |
со |
|
|
|
о |
о |
С.5 |
ю |
CN о |
о |
со |
|
|
|
||
о |
о |
о |
о |
о |
о |
о |
|
|
|
||
1 |
-ь |
+ |
+ |
1 |
+ |
+ |
+ |
+ |
|
|
|
СЧ |
г ч |
ft* |
»П о |
со 1-4 |
__, |
|
|
|
|
||
Г}* |
СО |
|
(ГЛ со |
ю |
05 |
со |
|
|
|
|
|
СЧ |
О |
|
о |
о |
о |
о |
|
|
|
|
|
О |
О о |
о |
о |
о |
о |
|
|
|
|
||
+ |
4 |
+ |
+ |
1 |
+ |
+ |
t- |
|
|
|
|
СОо |
іО |
о |
со |
05 |
со |
|
|
|
|
|
|
Юrf |
и.5 |
«) |
Гч со |
0 5 |
|
|
|
|
|
||
о |
— о |
о |
о |
сч ю |
|
|
|
|
|
||
о |
о |
о |
о |
о |
о |
|
|
|
|
|
|
+ |
+ |
1 |
+ |
+ |
+ |
|
|
|
|
|
|
05 |
*#• |
о |
05 |
СО |
со |
|
|
|
|
|
|
05 |
Гч |
|
с - |
со |
со |
|
|
|
|
|
|
о |
•— |
г-ч |
’-и ю |
| ч |
|
|
|
|
|
|
|
о |
о |
о |
о |
о |
о |
|
|
|
|
|
|
+• |
1 |
+ |
1 |
+ |
+ |
|
|
|
|
|
|
СЧ |
t-и |
ю |
со |
со |
|
|
|
|
|
|
|
СО |
О |
о |
0 5 |
|
|
|
|
|
|
|
|
о |
о |
о |
о |
05 |
|
|
|
|
|
|
|
о |
о |
о |
|
|
|
|
|
|
|
||
1 |
1 |
1 |
+ |
+ |
|
|
|
|
|
|
|
со |
со |
|
00 |
|
|
|
|
|
|
|
|
о |
|
|
СО |
|
|
|
|
|
|
|
|
—< |
<N |
со |
|
|
|
|
|
|
|
|
|
о |
о |
о |
о |
|
|
|
|
|
|
|
|
1 |
1 |
+ |
+ |
|
|
|
|
|
|
|
|
T f |
ю |
f-H |
|
|
|
|
|
|
|
|
|
О |
Ю |
Г ч |
|
|
|
|
|
|
|
|
|
со |
О |
со |
|
|
|
|
|
|
|
|
|
о |
о |
о |
|
|
|
|
|
|
|
|
|
++ +
СЧ СО
Tft о
СЧ Ю
О о
++
СЧ
со
о
4
і-н СЧ со ю со г- СО 05 О ГН СѵІ
135
нйе ассоциации Sr—Ba— (еще в континентальных условиях) и, наконец, в условиях морской седиментации происходит рас пад и этой ассоциации и возникают новые, например Sr —
Mn [10].
Результаты подсчета ги и и между микроэлементами, аутигенно-минералогическиіми формами железа, органическим углеродом и битумоидом представлены в табл. 1.
Для наглядности нанесем на «корреляционное кольцо» [11] только те связи, которые соответствуют 5-процентному уровню значимости (см. рис.). Анализ рисунка позволяет выделить три группы ассоциации элементов:
V, Cu, Ni, Mn, F e -, Fe"
Схема выделения парагенетических ассоциаций химических элементов в аптских песчано-алевритовых глинах междуречья Урал — Волга.
Связь Б—Ni мы относим к межгрупповой. Истинные взаимо отношения элементов между собой внутри группы проконтро лируем р;у (табл. 2). Первое число в каждой строке табл. 2 является множественным коэффициентом корреляции, отра жающим связь одного элемента со всей совокупностью.
136
Контроль rt] по |
Pij осуществляем по следующему |
прин |
||||
ципу: если |
Pij~fij |
или pij > ги , то элементы действительно |
||||
связаны между собой; если Рц существенно |
меньше гі} |
или |
||||
равны нулю, |
то следует признать, что рассматриваемые |
эле |
||||
менты между собой не связаны и высокое значение rtj |
ложно. |
|||||
Тогда по Pu (k), где исключается последовательно влияние |
||||||
каждого элемента, можно выделить |
причину |
возникновения |
||||
ложной корреляции. |
показывает, |
что выделенная |
первая |
|||
Анализ материала |
группа не претерпевает никаких изменений, вторая состоит из двух подгрупп — стронций-'барий и стронций-барий-литий, в третьей-группе исчезает связь битумоид-пиритное железо. Та ким образом, хотя образование песчано-алевритовых глин апта, отложившихся в основном в западной части междуречья Урал—Волга [12], происходило в условиях начавшегося хими ческого выветривания (отдельные ассоциации в первой группе: V—(Cu, V—Ni; Ni—Мп), однако осадки несут на себе черты континентального характера (резко обособившаяся, но не рас павшаяся ассоциация Sr—Ba).
Корреляция Мп—Fe" позволяет говорить о том, что марга
нец, очевидно, тоже присутствует в .двухвалентной |
форме. В |
|
свою очередь гидраты окисей и закисей |
железа и |
марганца |
сорбируют на своей поверхности другие |
элементы |
(V—Ni), |
являясь для них элементами-хозяевами. |
Органическое веще |
ство в данном случае не оказывает непосредственного влияния на распределение микроэлементов. Ипритное железо фикси рует благоприятную обстановку для битумообразования, в процессе которого часть никеля, очевидно, входит в состав металлоорганических соединений, связывая свою судьбу с би тумной частью породы.
ЛИ Т Е Р А Т У Р А
1.Руководство по предварительной математической обработке гео
химической информации при поисковых работах. |
М., «Недра», 1965. |
2. Р о д и о н о в Д. А. Функции распределения |
содержаний элементов |
и минералов в изверженных горных породах. М., |
«Наука», 1964. |
3. В о р о б ь е в В. Я., Л и п с к а я В. А. О типе закона распределения |
элементов в горных породах, почвах и растениях . Геология и геофизика
№ 7, |
1967. |
|
|
|
|
4. |
К р а м е р Г. |
Математические |
методы |
статистики. М., ИЛ, 1948. |
|
5. |
В и с т е л и у с |
А. Б., З у л ь ф у г а р л ы |
Д. И. Естественные па |
||
рагенезисы |
некоторых компонентов |
нефтей |
Азербайджана. Изв АН |
||
СССР, |
№ |
2, 1962. |
|
|
|
6. |
Ван |
д е р В а р д е н . Математическая статистика. М.( ИЛ, 1960. |
Ю. Заказ 1928 |
137 |
|
7. |
А н д е р с о н Т. Введение |
в многомерный статистический |
анализ. |
||||
М., |
Физматгиз, 1963. |
|
современного |
осадкообразования в |
||||
|
8. |
Л и с и ц ы н |
А. П. Процессы |
|||||
Беринговом |
море. |
М., «Наука», |
1966. |
|
|
|||
|
9. С т р а х о в |
Н. М. Основы теории литогенеза. М., Изд-во АН СССР, |
||||||
1962. |
Б у р к о в |
Ю. К. Линейные парагенезисы малых элементов в оса |
||||||
|
10. |
|||||||
дочных |
толщах как индикаторы |
условий седименгогенеза. К кн. Физиче |
||||||
ские и химические процессы и фации. М., «Наука», 1968. |
|
|||||||
|
II |
Т е р е н т ь е в П. В. Метод корреляционных |
плеяд. Вестник ЛГУ, |
|||||
№ |
9, |
1959. |
|
М. П., Ч а р ы г и н |
М. М. и др. Тектоническое строение |
|||
|
12. К а з а к о в |
|||||||
и история развития Прикаспийской впадины и смежных областей |
в связи |
|||||||
с |
вопросами |
нефтегазоносности. |
М., |
Гостоптехиздат, 1958. |
|
О. к. НАВРОЦКИЙ, Е. В. СОЛОДКОВА
ОБ ИСПОЛЬЗОВАНИИ
tr КРИТЕРИЯ ДЛЯ ИЗУЧЕНИЯ РАСПРЕДЕЛЕНИЯ
ХИМИЧЕСКИХ ЭЛЕМЕНТОВ В ЮРСКИХ И НИЖНЕМЕЛОВЫХ ОТЛОЖЕНИЯХ
ЗАПАДНОЙ ЧАСТИ ПРИКАСПИЙСКОЙ ВПАДИНЫ
Химические элементы, обладающие различными свойства ми и минералогическими характеристиками, по-разному реа гируют на физико-химические условия накопления их в осадке, что обуславливает их разнообразные формы миграции и распределение в породах.
По миграционной способности химических элементов, |
ко |
|
торая оценивается распределением |
средних содержаний |
на |
идеальном литологическом профиле, |
восстанавливаются |
па |
леогеографические условия осадконакопления [1]. Изменение во времени физико-географических условий (размеров бассей на, состава размываемых пород в областях денудации, тектони ческой активности исследуемой территории, климата) обу славливает различное накопление химических элементов в однотипных породах разного возраста. Сравнение средних со держаний элементов в разновозрастных петрографических профилях пород позволит выявить особенности накопления элементов во времени.
Обычное, визуальное, сравнение средних содержаний мо жет привести к ошибкам, так как при незначительной разнице между ними расхождение может быть существенным и, наобо рот, при значительной — несущественным.
Чтобы избежать подобных ошибок, рекомендуется приме
нять ряд статистических критериев: критерий |
Стьюдента — |
ю* |
139 |