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

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

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

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

Добавлен: 29.10.2024

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

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

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

где А = а d — новая константа. Если т]0 измерить на эталонах с нулевой концентрацией, то в (11.57) останется только две неиз­ вестные, A, k, которые можно вычислить по экспериментальным данным для эталонов из системы (полученной методом наимень­ ших квадратов)

 

 

 

 

 

Ct Дт)' — k \ r H ) ct =

 

 

 

 

 

2 (Ас1— С1“Н"

 

 

 

 

 

(11.58)

 

 

 

) Л7>1 = °

 

 

тле

п — число измерений эталонов.

 

Однако

■»)„

можно

вычислить

совместно с константами

A, k, для

чего (II 57)

представим в виде

 

 

 

 

 

 

Дт]'

 

 

 

 

(11.59)

 

 

 

 

 

С ~

В — D t{

 

 

 

 

 

 

 

 

 

 

 

 

где

Д = ^Л — \ ^k

\

D — k

1— новые неизвестные

величины.

Для неизвестных В , D, т}0 методом

наименьших

квадратов сос­

тавим систему

 

 

 

 

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

 

V

О ( в

~

D r ii

j -

\

\

ct—О

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

V

С1

-

Dri'i) -

ги + у

СЛ

=

°}.

(11.60)

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

V

С1 [В -

 

) - \ + гн

= О

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив найденные из этой системы значения В, D, т]0 в (11.59),

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

Ш . Математический аппарат интерпретации у-спектрограмм радиационного захвата с помощью ЭЦВМ

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

26


Так, важнейшая систематическая ошибка в методике интерпре­ таций спектральных диаграмм многокомпонентной смеси возника­ ет часто из-за того, что исследуемая смесь предполагается би­ компонентной (в лучшем случае трехкомпонентной) из элементов, доминирующих по тем или иным ядернофизическим свойствам над остальными. Подобная методика позволяет лишь грубо оценивать процентное содержание основного или интересующего нас элемен­ та, однако, она подкупает простотой и достаточно высокой экспрессностью.

Более достоверная оценка концентрации элемента стала воз­ можной благодаря высокому современному уровню автоматиза­ ции вычислений с использованием ЭЦВМ.

Ниже рассматривается возможность матричной идентификации на ЭЦВМ «немых» спектров захватного у_излУчения-

Интенсивность у-излучения радиационного захвата нейтронов пропорциональна концентрации данного элемента при малой' его концентрации. В этом можно убедиться из следующих сообра­ жений.

Число актов радиационного захвата нейтронов в единицу времени количеством ядер Nj компонента у в единице объема равно

 

 

 

 

=

 

 

 

 

(Ш.1)

^ / — плотность

потока нейтронов;

aj — сечение

захвата для ядер

компонента yj,

а

число

квантов

nt . в канале i

энергией

е,)

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

dNj

т. е.

 

 

 

 

 

 

at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nH = aijfajNr

 

 

 

(Ш-2)

где

atj — коэффициент

пропорциональности (доля выхода у-кван-

тов в канал i от компонента у). Интенсивность

I t

в

 

канале i

от

компонента у равна

 

 

 

 

 

 

 

 

 

 

 

In = bi ni,*l'

 

 

 

(Ш.З)

где

— энергия

у-кванта, попадающего в канал

i;

 

bt — аппара­

турный коэффициент пропорциональности для канала г.

 

Следовательно

 

 

 

 

 

 

 

 

 

1,: =

 

 

С ш

 

с

 

(III.4)

 

Ц av fz\ Nj = Ь.ач г/ а ] - j — = d „

 

 

ч

 

 

 

J

ч

У’

 

 

 

 

 

 

 

 

 

 

 

 

где т — удельный вес образца в граммах; Aj — вес атома компо­ нента у в граммах; d tj — константа, включающая все указанные константы; с . — концентрация компонента у.

и

27


Величина / у пропорциональна

концентрации компонента у .

Однако рассматриваемые выкладки справедливы при

постоянст­

ве плотности /

нейтронного

потока по

всей

глубине

проник­

новения

этого

потока в

рассматриваемый

объем.

В реаль­

ных

же

условиях плотность

потока

по

мере

проникнове­

ния

в вещество

ослабляется

тем

больше,

чем больше

концен­

трация компонента.

По принципу суперпозиции суммарная интенсивность в кана­

ле i от всех компонентов /

исследуемого объекта может

быть

записана в виде

П

 

 

 

Л =

2 7Ц + Ф0

(Ш.5)

 

J- 1

 

где п — число компонентов; Ф; — фон в канале i.

Из формулы (III.4) для одного и того же компонента у следует постоянство отношения интенсивностей в двух каналах г и 1.

- r

= T T l - r = klj = const

<IIL6>

чj

°\ “lj E1

 

Как видно, спектральное отношение ktj зависит от выбора кана­ ла г и компонента у и не зависит от абсолютных значений / у и /j., т. е. спектры одного и того же компонента у при различных его концентрациях „подобны". Из |III.4) так же следует, что в одном и том же канале i интенсивности от различных компо­ нентов пропорциональны их концентрациям, т. е.

Л!

= ац Ч

bi f ° ) ci Ak

(III.7)

hk

aik si

bi

Ck A j

 

где a.jk — постоянная величина для компонентов j, k и канала i. Из (III.5) с учетом (III.6) получим

Л = 2 Х - А / + ф г

(Ш-8)

i= 1

Грубо коэффициенты k tj могут быть определены промером ин­

дивидуальных моноэлементных спектров у-лучей захватного излу­ чения компонента / из отношения

k t, = - j -

( HI 9 )

У

 

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

28


Более надежный по точности способ определения коэффициен­

тов k tj может

быть рекомендован из следующих соображений.

Проведем в

тождественных условиях, соблюдая

геометрию

опыта, измерения на эталонных объектах с номерами

число эталонов), где известны концентрации всех / компонентов.

На основании (III.4) для всех г

эталонов и канала

1 можем запи­

сать

 

 

 

 

 

 

т. е.

 

 

 

 

 

 

(ШЛО)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(III. 11)

Тогда уравнение (III.8) для эталона

г принимает

вид

 

 

 

 

 

 

 

 

(III. 12)

Для

первого канала получаем

 

 

 

 

 

 

 

 

 

 

 

 

(III.13)

так

как

 

 

 

 

 

 

 

 

 

 

 

 

 

(III. 14)

где j\j,

— неизвестные величины,

которые необходимо

опре­

делить

в дальнейшем; ст} , с! — известные

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

компо­

нента у в эталонах г, 1.

 

 

 

 

 

Для

решения системы (III. 13) с л +1

неизвестной величиной

необходимо исследовать t ^ n + l

эталонов.

Если t = n 4-1, то полу­

чим обычную систему линейных уравнений, решение которой даст величины l \j, Фь Однако для повышения точности целесообразно взять i > n + 1, чтобы иметь систему линейных условных уравнений, решение которой зависит от выбора метода решения. Если тако­

вым выбрать

метод наименьших

квадратов, то (III. 13) сводится к

системе п + 1

линейных уравнений

 

 

dS_

 

 

 

' A \ ) 4 - = 0 , (1 < £ < « )

 

 

r=1

Cb

 

 

k

(111,15)

 

 

,

которая следует из минимизации

выражения

 

 

 

 

(III.16)

29