Файл: Ядернофизические методы анализа и контроля технологических процессов [сборник статей]..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