Файл: Ядернофизические методы анализа и контроля технологических процессов [сборник статей]..pdf

А Ф,О к ьк

к

ехр

( Г2

Гм) РV

CJ -

* “

16

(Ш.39)

где ij — массовый

коэффициент

поглощения ^-излучения для

компонента

у;

г,,

г 10,

г2,

г20 — соответственно

расстояния от

А

до С, от А

до

стенки

скважины

в

направлении г,, от В до

С,

от В до стенки скважины

вдоль г 2;

р — плотность анализируемой

среды.

d v

=

гl sin B2d r2dB2d<?2.

(111.40)

Тогда плотность

потока

-^-квантов

от компонента k в

точке В

со всего объема

суть

Л

= J ~

d v =

dv

=

JSinM62j dr2 r2 ' eXP - \ r * ~

П

К

■r 2 0 ) ^

' j

Cj - { r i

- r u ) Ly N j * j

(III.41)

гю = ( & + rlo - 2 R r20cos 02y/e ’

(1П,42)

где R — расстояние между точками

А и В.

Интенсивность от

компонента k

пропорциональна

величине

/* . Таким образом,

просчитав зависимость f kB от ck,

мы сможем

установить закономерность изменения интенсивности,

создавае­

мой компонентом k,

в зависимости

от концентрации.

между интенсивностью и концентрацией в виде

Ilj = aijCj + bijcj ,

(III.43)

где

аг„

Ьц — искомые константы для компонента j

и

канала I.

Суммарная интенсивность в канале i

А ,

=

2

h i + Ф, =

Фг +

2 (

с +

Ь1} с)

)

(III.44)

7=1

7=

1

является

линейной

системой

относительно

неизвестных

Ф .,

a tj,

bi}.

На эталонах

1 < r < t

получим

А \ ; тогда для

фиксированно­

го канала

I имеем систему

уравнений

'а;= 2 К<7,+а</4)+ф<-

<I1L45>

7=1

для

решения которой относительно atj,

Ф(

необходимо

взять

1 < ^ г < £ > 2 « - |- 1

эталонов

(« — число компонентов).

Повторив

подобную

процедуру

для всех

г (1^ г - ^ т ) ,

мы

вычислим

мат­

рицы (#,,)> ( 6у) и все Фг по

измерениям А[

на

эталонах,

пос­

ле

чего можно

анализировать

пробы, решая

систему

A t = 2

( а ц cj + b u c)

) + ф«

1 < 1 <

m

(m-46)

7 =

1 '

относительно Cj. Как видно,

эта система

нелинейна

относитель­

но Cj и

ее

можно решить приближенными методами.

По

аналогии

с (III.22)

и

с

учетом

нелинейности

функции

с. = / ( / у )

запишем

выражение

т

с ,

= 2

[

(А -

Ф,) + Н

(А -

ФА ]

<ш-47>

7=1

и приведем его

к

виду

т

cJ =

a] +

2 [ d j i A l +

ajlA l2 ),

(Ш.48)

где

т

«, =

2

А

К А - * / ,) •

(ш -49>

7-1

dn =

Ьн -

2ап фс.

(III.50)

Для определения концентрации с. необходимо знать

величи­

ны

a.,

djt,

CLjt, которые

можно

найти

из

измерений

эталонов

1

Подставляя

результаты

измерений эталонов

в

(Ш.48),

получаем

систему

т

с',

=

+

2

К А ,

+ ° Х ) ■

(41.51)

7=1 х

Систему (III.51)

можно решить относительно ар djt,

ан

при

условии, что 1

< г < г ? > 2 / г е + 1 .

Если

t = 2 / n + l ,

то

система

(III.51) решается непосредственно. При

£ > 2 / я + 1

систему мож­

но решить методом

наименьших

квадратов, после

чего

подста­

вив

найденные

a.,

aJt в (III.48), вычислить концентрацию

с;.

по

результатам

At

измерения анализируемой пробы.

в интервале

i\ //г — ин­

тенсивность от компонен­

та j

в интервале

I.

Учитывая

постоянство

спектральных коэффици­

ентов k.tj

[8— 10],

можем

записать

k i, =

-Tl

=

hi

(IV.2)

I

ij

hi

/ г ш6

где штрихи указывают на

опытные

значения, полу­ Рис. 6. Моноэлементные спектры Hg, Са. Fe.

ченные

на

моноэлемент-

ном

(/)

спектре.

Все спектральные коэффициенты можно в

оСз

(IV.3)

&22 —

Са’

«1

где 5 гСа — площадь спектра Са' в

интервале

i.

Элемент

/

Спектральные коэффициенты

/ = 2

i = 3

Fe

1

0 ,4 9 8

0 ,3 3 2

Са

2

0,3 2 1

0 ,0 2 0

Hg

3

0 ,0 1 8

0

Из

(2)

и (1) следует

система

3

7! =

2

* / Л .

(iv -4>

j=i

решая

которую

относительно

неизвестных

/ 1;.,

для интенсивности

Hg в первом интервале

получаем

г

( k 2i —

k m ) ( /,f e 3, —

/ 3 )

—

( &31 —

k 3y)

( / ife;i

—

/> )

/ т у t - \

13

(*21 — ^ 22) (&S1 ■— * 33) — (^31 — ^ 32) (£21 — ^ 23)

Подставляя

сюда

найденные

выше

значения

спектральных

коэффициентов, находим выражения для определения

интенсив­

ности

ртутной

компоненты:

Лз = 1,07/, -

3,42/2 +

1,77/3.

(IV.6)

На спектрограммах

ртутного

месторождения,

интенсивность

в

интервале

3 составляет — 20% от интенсивности

в

интервале

2,

следовательно,

уравнение (1V.6)

примет вид

Лз =

1,07/, -

3,С7/2.

(IV.7)

С учетом

фонового

излучения элементов

вмещающей породы

система (IV.4)

запишется в форме

з

Mi =

/, -

i f

=

2 V

y

С < г' <

3).

(IV.8)

j - 1

Решение

уравнения

(IV.8) для /13

дает

выражение

/« = 1.07ДЛ -

3,42Д/2 +

1.77Д/,.

(IV.9)

Зная

для

интервала

3 его

приближенное

значение

Д/3= 0 ,2 Д /2,

получаем

Лз =

1.07ДЛ - 3,07Д/2,

(IV. 10)

Лз =

1,07/, -

3,42/2 +

1,77/з - 1>07ЛФ+

3,42/ 2Ф-

1,77/* =

=

1,07/, -

3,42/2 + 1,77Л + С.

( IV. 11)

Поскольку

/* ,

/* , /3 — постоянные

величины,

то

С = const,

т. е. учет

фона i f

не влияет на форму кривой

/ 13, а лишь

см е­

щает ее на постоянную величину С.

Величина фона

в энергетическом интервале

i

может

быт ь

найдена из выражения

l f

= ~г~— I*

*

(IV.12)

* *

b — a J i

a

где Л — глубина скважины;

b, a — пределы

интервала по глуби­

не скважины.

Так как L = / ( h.}, то

П=-пЪлН,)’

0V-13)

j - 1

Рис. 7. Представление

результатов по скважине

№

178.

(IV. 10)

при i f и

i f ,

найденных из

(IV. 13) в интервале от 20

до

21

м,

в котором

1и

/ 2 изменяются

достаточно плавно. На кривой

Аз

(№

7) против

ртутных пластов

ртути видны аномальные

от­