Файл: Применение ЦВМ и средств вычислительной техники в геологии и геофизике [сборник]..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 30.10.2024

Просмотров: 51

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

(а), дисперсия (Ь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