Файл: Ядернофизические методы анализа и контроля технологических процессов [сборник статей]..pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 55
Скачиваний: 0
путем дифференцирования по неизвестным Jlk ,ФГ Здесь 'Л [ —
экспериментальные результаты, соответствующие к \ , которые предполагаются нормированными на показания на некотором эталоне.
Теперь можно вычислить коэффициенты пропорциональности d x.(см. ШЛО), составив сумму (метод наименьших квадратов)
t
S = |
2 |
U |
, - - C; 4 j |
(HI-17) |
|
|
|
r = 1 |
|
|
|
и продифференцировав ее |
по d Vj\ |
|
|||
t |
0 С7 |
|
t |
|
|
г=1 |
|
гS7)=1 |
|
||
V |
1 * СТ. |
|
|
||
t |
( |
сг |
„ ~~ |
* |
(III.18) |
v |
у |
с1 V |
|
||
-J |
1 |
LJ |
) |
|
|
Г—1 |
|
|
|
Г=1 |
|
где l\j определены из системы (III.13). Подставив / | ; в (III.12),.
получим |
при |
фиксированном канале |
i |
(2 |
|
i ■< т — число |
кана |
|||||||||||||
лов) |
систему |
линейных уравнений (если 1< г |
|
— п |
|
1) |
отно |
|||||||||||||
сительно |
Фг, |
ku , которая для данного |
i имеет |
единственное |
ре |
|||||||||||||||
шение. Повторив подобную процедуру |
для |
всех |
каналов |
т — I |
||||||||||||||||
раз, |
найдем |
все |
п ( т — 1) |
элементы |
^.(зам етим , |
что |
|
= |
1) |
|||||||||||
матрицы [ktj) и все значения Ф .. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Д ля |
нахождения матричных |
элементов |
ki} |
и Ф. |
для |
данного |
||||||||||||||
канала |
i |
взамен |
системы |
(III. 12) |
может |
быть решена |
система |
|||||||||||||
п -р 1 |
уравнений |
с п + 1 |
неизвестными |
ktj, |
Фг; |
|
|
|
|
|
|
|||||||||
dS |
= 2 2* ( ^ г |
А \ ) Л** “ Т ~ = |
|
(1 < k < п, 2 < i < т) |
|
|
||||||||||||||
dkik |
|
|
|
|||||||||||||||||
|
|
/•=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Д1Н.19) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
dS |
22(л;-'лг) = °, (2 </<от)‘ |
|
|
|
|
|
|
||||||||||
|
|
|
дФ. |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
Г=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Повторяя |
подобную процедуру |
при |
всех |
i |
(2 < г '< /? г ), |
полу |
||||||||||||||
чаем |
все |
k.. |
и Ф.. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
1} |
|
I |
|
можно |
определить |
иным |
образом. |
||||||||||
Элементы |
матрицы |
|
||||||||||||||||||
Из уравнения |
(III.5) для |
пробы |
„г“ |
согласно |
(111.11) |
можно |
за |
|||||||||||||
писать |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л: |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
CJ |
|
|
ф .. |
|
|
|
|
|
|
|
( 111.20) |
||
|
|
|
|
|
|
А л ' |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
! = 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. Д л я фиксированного |
канала |
i при |
1 |
<. г |
|
|
— t -j- |
1 |
и замене |
Л ; на экспериментальные значения 'А ' имеем систему линейный
30
уравнений |
для непосредственного определения |
Ф.. Повторив |
||
подобную |
процедуру |
для всех |
каналов 1 < г < « , |
получим все |
матричные |
элементы |
Тогда |
из соотношений |
|
|
|
к‘>= - 7 Г - Т |
(HI-21). |
найдем все матричные элементы ktj. С помощью этой матрицы машина как бы запоминает специфические особенности индиви дуальных моноэлементных спектров. При введении величин 34.
нового спектра исследуемой смеси машина решает задачу (ре шая систему III.8), как получить новый спектр из индивидуаль ных спектров тех компонентов, на основе которых определены матричные элементы k.tj и элементы Ф .. Если считать, что вели чины Фг в процессе последующих измерений на анализируемых
пробах остаются постоянными, то систему (III.8) можно решить всего один раз (поскольку коэффициенты kl} постоянны):
т
А; = 2 ^ ( Л - ф ,-)- |
( П1.2 2 ) |
г = 1
Здесь матрица ( ^ г)суть обратная матрице [ktjy Иными словами,,,
стоит лишь найти элементы b /г обратной матрицы, и концентра ция с. вычислится из (III.22), (III.18), (ШЛО).
Для повышения точности необходимо оценить изменения ве личин Ф.. Представим Ф." — новое значение фона в канале i как функцию этого канала:
Фн |
1 |
Ф. |
V |
(III.23) |
I |
I |
|
Т*0
Тогда (III.8) принимает вид
пр
|
■ 2 * , ' „ |
+ |
2 a q |
( i + 1 ) + Фг > |
(III.24) |
|
|
j = 1 |
q |
С |
|
|
|
где р — степень |
полинома. Система |
уравнений |
с п 4- р |
1 не |
||
известными /, . |
и ад должна |
иметь |
т ^ п - \ - р ~ |
1 число |
кана |
лов i, чтобы решаться относительно /, . и aq . Затем можно найти
концентрацию |
по /];. с помощью |
(ШЛО), |
(Ш.18) и величину |
||
Ф" по (111.23). |
|
|
|
|
|
Как видно, |
при |
больших т (1 < / < |
т) для определения коэф |
||
фициентов ktj, a q |
необходимо |
привлекать |
вычислительную ма |
||
шину. |
|
|
|
|
|
, До сих пор мы предполагали, |
что интенсивность радиационного- |
излучения при захвате нейтронов ядрами компонента k пропор циональна концентрации этого компонента. Однако плотность
ч
31
потока уменьшается вследствие процессов взаимодействия нейтро нов с веществом рудного пласта, пересеченного скважиной (рис. 4)
иудаления потока от источника.
Для достаточно быстрых нейтронов, можно, пренебрегая диф
фузией, представить плотность потока от точечного источника в виде
f = _Фо_ |
Г |
к |
|
ехр |
(III.25) |
||
J |
4-г-' |
|
|
где Ф0 — число моноэнергетических найтронов, испускаемых ис точником в единицу времени; N . — число ядер компонента у в
Рис. 4. Схема расположения источника и детектсра
вскважине.
единичном объеме; о — полное сечение ядер компонента у; г, г0 — соответственно расстояния от источника до рассматриваемой точ ки и до стенки скважины по направлению г;
|
|
|
|
aj |
= |
° j + aj’, |
|
|
|
(Ш.26) |
||
of — сечение |
радиационного |
захвата |
ядрами |
компонента |
у; |
|||||||
a,j — сечение, |
обусловленное |
упругим |
и |
неупругим |
рассея |
|||||||
нием, ядерным |
расщеплением, делением |
ядер, |
поглощением |
|||||||||
нетронов и т. д. |
|
|
|
|
|
|
|
|
|
|
||
|
Следовательно, плотность распределения актов радиационного |
|||||||||||
захвата нейтронов |
ядрами компонента k определится |
как |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
(III.27) |
|
где |
N k — плотность |
распределения |
ядер компонента k |
(т. е. чис |
||||||||
ло |
этих ядер |
в |
единице |
объема). |
Переходя к концентрации |
ск |
||||||
с помощью соотношения |
|
|
|
|
|
|
|
|
||||
|
|
|
|
ck = A - 1 N k , |
|
|
|
(П1-28) |
32
лолучаем |
|
|
dnk |
|
|
|
|
|
|
|
= A H |
ck |
|
(III.29) |
|
|
|
|
dv |
|
|||
|
|
|
|
|
|||
Из (Ш.28) следует |
dnk |
|
|
|
|||
|
|
|
п и |
d v = |
|
|
|
|
|
|
dv |
|
|
||
|
|
2я |
x |
|
|
|
|
= Ла! ck -fe- |
о |
J sin bdb J exp |
- (r - r0) 2 |
N j a. |
dr. (Ш.30) |
||
|
|
0 |
|
|
|
|
|
Поскольку во |
внутреннем |
интеграле r 0 не зависит от |
изменения |
||||
г, то число актов |
захвата нейтронов ядрами |
компонента k в за |
|||||
висимости |
от ск |
равно |
|
|
|
|
|
A°l ск ф0 |
Я |
|
|
|
|
|
|
Jsin Mb exp |
|
|
|
|
|||
п к |
2 |
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
A°l |
Ч фо |
j* sin0d0 exp |
|
|
|
||
|
2 |
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
_ |
H ск фо |
|
|
(III.31) |
|
|
|
|
|
|
|
Но^число ^-квантов радиационного захвата пропорционально nk хотя бы в первом приближении, —
ЬФ0 Ск |
(III.32) |
п \ = Ьпк = |
|
ск + ( АЛ ) 1 2 |
Nj °j |
Следовательно, и число ^-квантов, зарегистрированных де тектором, и интенсивность от компонента k будут пропорцио
нальны |
п \ : |
|
|
|
|
|
|
|
|
!„ = В*1 = -------------- Ьл^ |
г Ь -----------• |
(ш -33> |
|||
|
|
|
С* + ( Л31) |
H l NJ |
aJ |
|
|
|
|
|
|
|
1 |
|
|
где В, |
Ьх — коэффициенты |
пропорциональности. |
|
||||
Допустим далее для простоты, что увеличение концентрации |
|||||||
компонента k |
происходит |
за счет |
изменения |
концентрации ком |
|||
понента j по |
закону |
|
|
|
|
|
|
|
|
Nj = Nj ~ |
aj N*> |
|
(ш -34) |
3—99 |
33 |
гд е |
dj |
к он стан та д л я |
к ом п он ен та |
j . |
Т о гд а |
|
|||||
|
|
h |
= |
|
Ь1фо cft |
|
|
|
|
||
|
|
Ч |
|
k~\ |
|
N° - |
aj Nk ) |
||||
|
|
|
Cft + |
~1 2 |
|
||||||
|
|
|
|
|
bl Ф0 Cft |
|
|
|
|
||
|
|
|
’ч |
|
|
ft-1 |
|
|
|
k |
- \ |
|
|
|
|
Nft- sjU |
V |
-j |
|
|
1 |
||
|
|
|
|
фо cft |
|
|
|
|
|
||
|
|
|
|
Ak‘‘ft |
c'ftk + |
Bk |
|
|
|
|
|
где |
A* , |
5 , |
константы |
для |
компонента |
|
k : |
|
|||
|
|
|
= bi 1 [ 1 - ( al ) 12 |
|
°y а/ |
— const, |
|||||
|
|
|
|
|
|
||||||
|
|
B,= N (A°l г |
k-1 |
|
|
|
|
|
|||
|
|
i |
■з. N j |
= const. |
(111.35)
(111.36)
(111.37}
|
|
|
|
|
График |
функции |
(III.35) при |
всех |
||||||||
|
|
|
|
|
положительных |
константах |
имеет |
|||||||||
|
|
|
|
|
вид, показанный |
на |
рис. 5. |
|
Если |
|||||||
|
|
|
|
|
учесть поглощение, то зависимость |
|||||||||||
|
|
|
|
|
интенсивности |
Jk |
от |
концентрации |
||||||||
|
|
|
|
|
|
ск |
нельзя |
представить в элемен |
||||||||
|
|
|
|
|
тарных функциях. Рассмотрим этот |
|||||||||||
|
|
|
|
|
вопрос на примере скважины. |
|
||||||||||
|
|
|
|
|
|
|
Пусть |
в |
точке |
А |
на оси |
|
сква |
|||
Рис. 5. Зависимость |
/* от ск- |
жины |
(см. |
рис. |
4) |
находится |
ис |
|||||||||
|
|
|
|
|
точник |
нейтронов, а в точке Б —де |
||||||||||
тектор -[-квантов. |
В точке |
С среды |
с элементарным объемом d v |
|||||||||||||
происходит |
d n k (см. |
III.28) |
актов захвата нейтронов, |
следователь- |
||||||||||||
но, |
возникает dm h |
-[-квантов: |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
d m k = |
bdnk |
|
|
|
|
|
|
|
(III.38) |
|||
где |
b — коэффициент |
пропорциональности. |
Однако |
точки |
В до |
|||||||||||
стигают не |
все эти |
-[--кванты. |
Плотность |
их |
потока из объема |
|||||||||||
d v |
(точки С) в точке В равна |
|
|
|
|
|
* |
|
|
|
|
|
||||
|
|
|
dmk |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
d f B = |
|
НГ ехр |
- |
( |
г2 - |
г2о) р 2 |
Ci |
|
|
|
||||
|
|
|
|
Г 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
bdih |
|
( Г2 - Г 20) р*^ |
|
т |
с |
|
|
|
|
|||||
|
|
= Т |
^ |
ехР |
|
|
|
|
|
|||||||
|
|
|
j |
CJ |
|
|
|
|
||||||||
|
|
4 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
34