В данной статьи приведены зависимости приборных интегральных и дифференциальных массовых коэффициентов ослабления рентгеновского излучения для стали Ст45 от толщины и аналитические выражения для их аппроксимации. Выведено уравнение рентгеновского измерителя толщины.
Геометрическая схема рентгеновского измерителя толщины
На рис. 1 приведена геометрическая схема рентгеновского измерителя толщины стального проката.
К выбираемым параметрам конических коллиматоров источника рентгеновского излучения и радиометрического детектора относятся их толщины hK1, hK2 и диаметры выходных окон dK1, dK2. Одним из наиболее распространённых критериев выбора толщин коллиматоров является ограничение на k – кратность ослабления рентгеновского излучения, например, кратность ослабления k не менее 1000.
Для выбора фокусного расстояния – F=1000 мм, чувствительный объём радиометрического детектора – сцинтиллятор CsI имеет размеры: диаметр dсц=25 мм; толщину hсц=45 мм. В результате расчетов получили , dK1=0,8 мм, dK2=23,9 мм.
Расчёт зависимостей приборных интегральных и дифференциальных массовых коэффициентов ослабления рентгеновского излучения от толщины объекта контроля и от максимальной энергии рентгеновского излучения mint(ρH,Emax) и mdiff(ρH,Emax)
Параметры, характеризующие процесс ослабления и регистрации рентгеновского излучения, зависят от толщины ослабляющего барьера – ρH, максимальной энергииEmax рентгеновского излучения и энергетического спектра – f(E,Emax), материала сцинтиллятора и его толщины – hсц,. Приборный интегральный массовый коэффициент ослабления (МКО) рентгеновского излучения вводится следующим образом
, (1)
где I0 – величина поглощенной в детекторе энергии рентгеновского излучения с максимальной энергией в спектре Emax без ОК при той же геометрии; Eab(E) – среднее значение поглощенной энергии зарегистрированного кванта с энергией E;μсц(E),m(E)– энергетические зависимости линейного и массового коэффициентов ослабления фотонного излучения материалом сцинтиллятора и ослабляющего материала.
Формула для вычисления приборного интегрального массового коэффициента ослабления рентгеновского излучения mint(ρH,Emax) выводится из (1)
. (2)
Выражение для вычисления приборного дифференциального МКО mdiff(ρH,Emax)имеет вид
. (3)
Следует отметить, что интегральные и дифференциальные линейные коэффициенты ослабления связаны друг с другом. Несложно вывести выражение, связывающее mint(ρH,Emax) и mdiff(ρH,Emax). Уравнение упомянутой связи имеет вид
. (4)
Была проведена серия расчётов по формулам (2), (3) интегральных mint(ρH,Emax) и дифференциальных mdiff(ρH,Emax) МКО рентгеновского излучения для стали Ст45 в диапазоне максимальных энергий Emax от 100 до 450 кэВ. Толщины стального проката варьировали от 1 до 20 г/см2. Для описания энергетического спектра применяли выражение, близкое к формуле Крамерса, . Данные по коэффициентам ослабления гамма-излучения с веществом заимствованы из 127 групповой библиотеки данных [4], по зависимости Eab(E) – [5]. Расчёты проводились в системе MathCad. Результаты расчётов mint(ρH,Emax) и mdiff(ρH,Emax) сведены в табл. 1 и 2.
Таблица 1. Приборные интегральные массовые коэффициенты ослабления рентгеновского излучения – mint(ρH,Emax)
Emax, кэВ |
ρH, г/см2 |
|||||||||
1 |
2 |
3 |
4 |
5 |
6 |
8 |
10 |
15 |
20 |
|
100 |
0,536 |
0,518 |
0,506 |
0,497 |
0,49 |
0,483 |
0,473 |
0.465 |
0,451 |
0,443 |
150 |
0,273 |
0,258 |
0,247 |
0,236 |
0,227 |
0,218 |
0,203 |
0,19 |
0,161 |
0,145 |
200 |
0.216 |
0,208 |
0,201 |
0,196 |
0,191 |
0,186 |
0,179 |
0,172 |
0,158 |
0,149 |
250 |
0,188 |
0,182 |
0,177 |
0,173 |
0,17 |
0,167 |
0,162 |
0,157 |
0,148 |
0,142 |
300 |
0,171 |
0,166 |
0,162 |
0,158 |
0,156 |
0,153 |
0,149 |
0,146 |
0,139 |
0,134 |
350 |
0,159 |
0,154 |
0,151 |
0,148 |
0,145 |
0,143 |
0,14 |
0,137 |
0,131 |
0,127 |
400 |
0,149 |
0,145 |
0,142 |
0,14 |
0,137 |
0,135 |
0,132 |
0,129 |
0,124 |
0,121 |
450 |
0,142 |
0,138 |
0,135 |
0,133 |
0,131 |
0,129 |
0,126 |
0,124 |
0,119 |
0,116 |
Таблица 2. Приборные дифференциальные массовые коэффициенты ослабления рентгеновского излучения – mdiff(ρH,Emax)
Emax, кэВ |
ρH, г/см2 |
|||||||||
1 |
2 |
3 |
4 |
5 |
6 |
8 |
10 |
15 |
20 |
|
100 |
0,495 |
0,478 |
0,466 |
0,456 |
0,448 |
0,442 |
0,436 |
0,428 |
0,411 |
0,402 |
150 |
0,23 |
0,211 |
0,194 |
0,18 |
0,167 |
0,156 |
0,137 |
0,123 |
0,1 |
0,088 |
200 |
0,187 |
0,178 |
0,171 |
0,166 |
0,16 |
0,156 |
0,148 |
0,137 |
0,129 |
0,118 |
250 |
0,163 |
0,157 |
0,153 |
0,149 |
0,147 |
0,144 |
0,14 |
0,137 |
0,13 |
0,124 |
300 |
0,146 |
0,142 |
0,139 |
0,137 |
0,134 |
0,133 |
0,13 |
0,127 |
0,123 |
0,12 |
350 |
0,135 |
0,131 |
0,129 |
0,127 |
0,125 |
0,124 |
0,121 |
0,119 |
0,116 |
0,114 |
400 |
0,126 |
0,123 |
0,121 |
0,119 |
0,118 |
0,116 |
0,114 |
0,113 |
0,11 |
0,108 |
450 |
0,119 |
0,116 |
0,114 |
0,113 |
0,111 |
0,11 |
0,109 |
0,107 |
0,105 |
0,103 |
Для иллюстрации на рис. 3 приведены зависимости mint(ρH) и mdiff(ρH) для нескольких уровней энергий Emax.
Рис. 3 – Зависимости mint(ρH) и mdiff(ρH)
Уравнение измерителя толщины
Уравнение измерителя толщины должно связывать оцениваемый параметр объекта контроля – ρH с измеряемой физической величиной – I(ρH,Emax).
, (5)
здесь I(ρH,Emax) – величина поглощенной в детекторе энергии рентгеновского излучения с максимальной энергией в спектре Emax за барьером толщиной ρH при фиксированной геометрии контроля; I0 – величина поглощенной в детекторе энергии рентгеновского излучения с максимальной энергией в спектре Emax без объекта контроля при той же геометрии. Величины I(ρH,Emax) и I0 на момент начала этапа оценки толщины являются известными, поэтому (7) может быть переписано в следующем виде
. (6)
Уравнение (6) является нелинейным относительно неизвестной ρH.
В качестве аппроксимации mint(ρH) при фиксированном значении Emax было решено использовать функцию, зависящую от четырёх параметров. Указанная функция является естественным усовершенствованием и выглядит следующим образом
. (7)
В качестве аппроксимации mint(ρH) при фиксированном значении Emax было решено использовать функцию, зависящую от четырёх параметров. Указанная функция является естественным усовершенствованием [17] и выглядит следующим образом
. (8)
Подставим (8) в уравнение (6)
. (9)
Уравнение (9) и может быть интерпретировано как уравнение рентгеновского измерителя толщины. Коэффициенты a, b, c,d определяются на стадии калибровки для конкретной максимальной энергии рентгеновского излучения Emax и конкретного материала объекта контроля
Список литературыАртемьев, Б.В., Шубочкин, А.Е. Рентгеновская толщинометрия // Контроль. Диагностика. – 2014. – № 2. – С. 24–31.
Артемьев, Б.В. Рентгеновские толщиномеры // Контроль. Диагностика. – 2009. – № 4. – С. 22–25.
Артемьев, Б.В., Маслов, А.И., Потапов, В.Н., Ведерников, М.Б. Использование рентгеновских толщиномеров в производстве проката цветных металлов // Дефектоскопия. – 2003. – № 6. – С. 55–62.
http://www.ippe.ru/podr/abbn/libr/groupkon.php – 127 групповая библиотека данных о взаимодействии гамма- квантов с веществом.
Будаи, Б.Т., Касаткин, Н.В. Измерение параметров листового горячего проката // Инженерно-физический журнал. – 2014. – Т. 87. – № 1. – С. 225–228.