УДК 620.179.152.1
РАСЧЕТ ПАРАМЕТРОВ ПИНХОЛА КОЛЛИМАЦИОННОЙ СИСТЕМЫ
АЛЬБЕДНОЙ КОМПЬЮТЕРНОЙ ТОМОГРАФИИ
СТАЛЬНЫХ ИЗДЕЛИЙ
© 2020 г. Е.Е. Журавский1,*, Б.И. Капранов1,**, Д.С. Белкин1,***, С.В. Чахлов1,****,
А.М. Штейн1,*****
1Томский политехнический университет, Россия 634050 Томск, пр-т Ленинa, 30
E-mail: * zhuravskiy@tpu.ru; ** introbob@mail.ru; ***Belkin@tpu.ru; ****chakhlov@tpu.ru;
***** shteynam@gmail.com
Поступила в редакцию 25.10.2019; после доработки 12.12.2019
Принята к публикации 20.12.2019
Рассмотрены особенности коллимационной системы томографа для контроля стальных изделий на комптоновском
обратном рассеянии. Описано устройство и проведен расчет параметров первичного коллиматора «пинхола» коллима-
ционной системы, позволяющее достоверно проводить сканирование стального объекта контроля простым перемеще-
нием пинхола перпендикулярно оси рентгеновского пучка от анода трубки в пределах 50 мм при фокусном расстоянии
163 мм.
Ключевые слова: комптоновская томография, коллимационная система, пинхол, неразрушающий контроль, радиа-
ционный контроль, дефектоскопия, приборостроение.
DOI: 10.31857/S0130308220020050
ВВЕДЕНИЕ
Для внутритрубного диагностирования при одностороннем доступе стальных нефтегазопрово-
дов, применяемых и эксплуатируемых на опасных производственных объектах, используются
следующие виды неразрушающего контроля (НК): ультразвуковой, магнитный, вихретоковый.
Повысить эксплуатационную безопасность нефтегазопроводов возможно путем повышения
достоверности НК, в частности, применением различной комбинации видов (методов) НК.
Повышение достоверности и производительности применения радиационных методов для контро-
ля изделий из углеводородов, алюминия, стали при одностороннем доступе стало возможно благо-
даря разработкам рентгеновской вычислительной томографии (РВТ) на комптоновском обратном
рассеянии (КОР)[1—6]. Известна установка для альбедной компьютерной томографии (АКТ)
Гамбургской лаборатории ComScan фирмы «Philips». Разработанная установка применяется для
изделий с низким атомным номером (Z < 26) и позволяет контролировать при одностороннем
доступе изделия из стали с глубиной контроля до 3 мм, алюминия — с глубиной контроля до
20 мм при напряжении на рентгеновской трубке Emax 160 кВ. Ограничение толщины просвечивае-
мого слоя определяется энергией первичного и рассеянного излучения
E
0
E
(E
;θ)
=
,
s
0
(1)
E
0
1+
·(1cosθ)
0,511
2
где E0 — эффективная энергия рентгеновского излучения
E
, кэВ; Es — энергия рассеянного
max
3
излучения, кэВ; θ — угол рассеяния.
При углах рассеяния θ ≈ 135° величина Es (для E0 = 106 кэВ) составляет 78 кэВ, при заданной
энергии линейный коэффициент ослабления для железа составляет 4,671 см-1. Если принять в
качестве толщины просвечиваемого слоя 3 длины свободного пробега рассеянного излучения, то
для Es = 78 кэВ толщина для стали составит примерно 5 мм. Повысить предельно возможную
контролируемую толщину для стальных изделий возможно путем повышения напряжения на
рентгеновской трубке. При повышении напряжении на рентгеновской трубке параметры КС
также изменятся, вследствие чего необходимо выбрать и оценить параметры разрабатываемой
КС системы АКТ.
ОПИСАНИЕ КОЛЛИМАЦИОННОЙ СИСТЕМЫ
Компания Mitsubishi Heavy Industries Ltd описывает детектор РВТ на КОР в своем патенте,
предназначенном для обнаружения несплошностей (таких как внутренние трещины) в бетоне тун-
нелей и мостов [7].
Расчет параметров пинхола коллимационной системы альбедной компьютерной томографии...
43
Одна из основных составляющих детектирующих систем, использующих обратнорассеянное
рентгеновское излучение (ОРРИ) — коллимационная система, обеспечивающая пространственное
формирование первичного и рассеянного излучения (первичного и вторичного коллиматоров). Ее
основная характеристика — информационная способность, которая полностью определяется фор-
мой и размерами рассеивающего объема (РО), из которого собираются однократно рассеянные
кванты [8—12].
Для описания информационной способности коллимационной системы введено понятие
пространственно-статистической неопределенности (ПСН) по направлению сканирования
(x, y, z) [10].
Если учесть, что поток рассеянных квантов и квантов, попавших в детектор, описывается по
закону Пуассона, то относительная статистическая погрешность измерения его интенсивности
1
равна
δ
=
,
где N — количество квантов. Так как N пропорционально объему V с учетом
ст
N
того, что размеры воксела малы, то можно представить δст по следующей формуле:
K
δ
=
,
(2)
ст
V
где K — коэффициент пропорциональности. Погрешность измерения координат или простран-
ственная разрешающая способность в первом приближении определяется половиной ширины
апертуры на полувысоте по заданному направлению
l
δ =
(3)
р
2
Результирующая пространственно-статистическая неопределенность апертуры может быть
описана произведением δст·δр:
l·K
∆=δ
·
ст р
δ
=
(4)
2 V
Форма и размеры РО определяются формой и размером первичного и вторичного коллимато-
ров. Пространственно-статистическая неопределенность главным образом зависит от соотноше-
ния размеров первичного и вторичного коллиматоров. В качестве первичного коллиматора чаще
всего используют «пинхол» [10], формирующий узкий «карандашный» пучок первичного (зонди-
рующего) излучения. Вторичный коллиматор представляет собой «пинхол», растянутый вдоль
одной оси по направлению сканирования (щель).
Взаимное пересечение первичного пучка и зоны чувствительности вторичного коллиматора
формируют РО, геометрия РО представлена на рис. 1.
Объект контроля
РО
Размер пинхола
вторичного коллиматора
Размер пинхола
b
первичного коллиматора
a
Рис. 1. Вариации формы РО в зависимости от размеров пинхола первичного коллиматора (а) и вторичного колли-
матора (b).
Дефектоскопия
№ 2
2020
44
Е.Е. Журавский, Б.И. Капранов, Д.С. Белкин и др.
1
А
2
А
3
4
АА
5
6
Рис. 2. Конструкция системы коллиматоров:
1 — объект контроля; 2 — вторичный коллиматор; 3 — первичный коллиматор; 4 — пинхол; 5 — рентгеновское излучение;
6 — коллимационная система.
Практическая конструкция, разработанная для контроля стальных изделий при одностороннем
доступе толщиной до 10 мм, представлена на рис. 2.
Рентгеновское излучения 5, проходя через первичный коллиматор 3 преобразуется в узкий
«карандашный» пучок пинхолом 4, далее комптоновски обратнорассеянное рентгеновское излуче-
ние коллимируется пинхолом вторичного коллиматора 2. Разработанная конструкция предполагает
сканирование ОК. Существуют два способа сканирования ОК:
первый, при котором пинхол первичного коллиматора перемещается относительно ОК;
второй, при котором сама КС вместе с источником рентгеновского излучения перемещается
относительно ОК, в данном случае пинхол неподвижен.
Первый способ оказывается намного проще второго по техническому исполнению. При реализа-
ции первого способа к пинхолу предъявляется ряд требований. Конструкция пинхола должна обе-
спечивать постоянное значение поперечного сечения для всех углов наклона, то есть поперечное
сечение рабочего пучка рентгеновского излучения для каждого угла наклона должно быть одинако-
вым. При данном условии пространственно-статистическая неопределенность (в зависимости от угла
наклона рабочего пучка рентгеновского излучения) будет постоянной. Также требованием к пинхолу
является 100-процентное ослабление интенсивности рабочего пучка за стенками пинхола. Рассмотрим
возможность реализации первого способа сканирования применительно к разработанной системе.
МЕТОДИКА РАСЧЕТА ПАРАМЕТРОВ ПИНХОЛА
В данной работе проведен анализ геометрии формирования РО и пространственно-статистиче-
ской неопределенности сканирующей системы в целом. Рассмотрена система на основе рентгенов-
ского излучения с напряжением на трубке 300 кВ. В качестве эффективной энергии рентгеновско-
го излучения принято значение 200 кэВ. Первичный коллиматор (пинхол) выполнен из вольфрама
толщиной 5 мм. Размер пинхола первичного коллиматора (см. рис.1) принят 1 мм. Для данной
энергии рентгеновского излучения ослабление интенсивности первичного пучка в стенках пинхо-
ла составляет более 1000 раз. Сканирование проводилось вдоль продольной оси щели вторичного
с длиной 50 мм. Расстояние от анода рентгеновской трубки до поверхности рассеивателя равно
163 мм (см. рис. 2).
Угол наклона рабочего пучка рентгеновского излучения для крайних точек перемещения α
следующий:
163
α=
arctg
8,72°.
(5)
25
Рассмотрим возможную конструкцию пинхола, выполняющую требование постоянного попе-
речного значения рабочего пучка рентгеновского излучения. Постоянное значение поперечного
сечения во всех углах наклона может быть выполнено для схемы пинхола, приведенной на рис. 3,
где α = 8,72°; a = 1 мм.
Дефектоскопия
№ 2
2020
Расчет параметров пинхола коллимационной системы альбедной компьютерной томографии...
45
α
a
C
A
B
a/cos(α)
Рис. 3. Схема пинхола.
Как видно из рис. 3, схема обеспечивает максимальное отклонение диаметра рабочего пучка не
более 0,2 %.
Геометрия на рис. 3 не обеспечивает требование 100-процентного ослабления интенсивности
рабочего пучка в стенках пинхола, т.к. за границами А и В не происходит необходимого ослабления
рабочего пучка. Физически это требование невыполнимо, поэтому при разработке конструкции
пинхола было принято условие десятикратного ослабления за границами А и В.
Линейный коэффициент ослабления µW вольфрама при эффективной энергии 200 кэВ равен
14,417 см-1. Минимальная толщина вольфрама h, при которой обеспечивается ослабление интен-
сивности пучка рентгеновского излучения в 10 раз, составляет
ln10
L(E
,µ
)
=
=1,6мм.
(6)
0
w
µ
(E
)
w
0
Учтем вычисленную толщину вольфрама для получившейся конструкции на рис. 4.
На левой части рис. 4 изображена схема пинхола с учетом требования по десятикратному
ослаблению. Учесть два требования, предъявляемые к пинхолу, не удается. Поэтому исследуем
А
АА
B2
H1
X
∆X
e
L
c
a0
a
L
Рис. 4. Пинхол.
Дефектоскопия
№ 2
2020
46
Е.Е. Журавский, Б.И. Капранов, Д.С. Белкин и др.
схему пинхола на рис. 4. Выделив часть А левого рис. 4, можно видеть (правая часть рис. 4), что
требуемое ослабление на границах пучка обеспечивается при углах 8,72 и 0°. Для углов, лежа-
щих между этими значениями, требуемого ослабления излучения не будет. Но при увеличении
поперечного сечения рабочего пучка для углов, при которых не обеспечивается достаточного
ослабления интенсивности пучка, радиационная толщина достигнет необходимого значения при
XХ на рис. 4 (правая часть). Здесь X — поперечное сечение рабочего пучка; ΔХ — дополни-
тельная часть рабочего пучка, которая обеспечивает необходимое ослабление. Рассмотрим гео-
метрию рис. 4 с целью определения поперечных размеров пучка XХ для каждого угла α, а
также для определения интенсивности излучения I рабочего пучка для каждого угла. Исходя из
начальных условий, поперечное сечение пучка а (см. рис. 3) равно 1 мм, тогда ΔX будет опреде-
ляться для данного угла только геометрией конструкции при угле α, равном 8,72°:
a
a
0
(8,72)
=1
=
5,874мкм.
(7)
cos(8,72)
Как видно из уравнения (7), поперечное сечение рабочего пучка при угле 0° увеличивается за
счет геометрии конструкции всего на 2a0, то есть на 10 мкм. Величина фаски B2 пинхола
L
B α)
=
tgα
(8)
2
Величина фаски для угла 8,72° равна 0,122 мм. При выполнении требуемых условий изменение
сечения рабочего пучка для углов, отличных от 0° и 8,72° без учета ослабления, будет
L
X(α)
=
2
+∆a
(8,72)cosα
(9)
0
2
Анализ уравнения (9) показывает, что максимальным значение рабочего пучка будет для угла
0° и минимальным — для угла 8,72°.
Для анализа увеличения поперечных размеров пучка при углах, не равных 0° и 8,72°, учитыва-
ющих необходимое ослабление пучка, нужно рассмотреть два случая. Рассмотрим первый случай
B2(8,72)
зависимости поперечного значения пучка для углов, меньших
arcsin
(угол e на рис. 4).
L
Значение ΔХ составляет
L L
X
∆ α)
=
- αcosα
(10)
sinα
2
Анализ уравнения (10) показывает, что сечение рабочего пучка увеличивается при увеличении
B2(8,72)
B2(8,72)
угла до
arcsin
Второй случай возникает для углов, больших
arcsin
(угол c
L
L
на рис. 4). Отсюда значение ΔХ для второго случая
Ltg8,72
L
X
∆ α)
=
-
sinα
(11)
tgα
2
1,15
1,1
1,05
1
0,95
0,9
0
0,8
1,6
2,4
3,2
4
4,8
5,6
6,4
7,2
8
8,72
α, град
Рис. 5. Зависимость диаметра рабочего пучка излучения от угла.
Дефектоскопия
№ 2
2020
Расчет параметров пинхола коллимационной системы альбедной компьютерной томографии...
47
Анализ уравнения (11) показал, что сечение пучка излучения уменьшается, начиная с угла
B2(8,72)
arcsin
,
и минимально при угле 8,72°. На рис. 5 изображена зависимость сечения рабо-
L
чего пучка излучения D для разных углов, обеспечивающая уменьшение интенсивности излучения
в 10 раз.
Для определения зависимости полной интенсивности рабочего пучка от угла α рассмотрим
рис. 4. Величина L должна составлять не менее 1,6 мм. Рабочий пучок излучения, для α отличных
от 0°, будет частично задевать стенки пинхола, тем самым на его границах будет ослабление пучка.
Для вычисления полной интенсивности излучения рабочего пучка проанализируем изменение L
при перемещении вдоль отрезка H1
H1(α)
L(α)
=
(12)
sinα⋅cosα
Суммарная интенсивность излучения получается при интегрировании следующих уравнений:
µw
(
x-H1(α)
)
sinα·cosα
e
0
x< H
1(α)
(13)
I x)
=
1
H
→ α)
x
<
2
(
B
2+
X
)
H1(α)
−µ
x(2
B2+
X
H
1(α)
w
(
(
)
)
sinα·cosα
e
2
(
B2+
X
)
H
− α)
x
<
2
(
B2+
X
)
Проинтегрировав уравнения (13) по углам от 0 до 8,72° включительно, получим следующее
распределение нормированной интенсивности излучения In(x) рабочего пучка для заданных углов
(см. рис. 6).
1,02
1
0,98
0,96
0,94
0,92
0,9
0,88
0,86
0,84
0,82
0
1
2
3
4
4,4
5
6
7
8
8,72
α, град
Рис. 6. Зависимость интенсивности излучения рабочего пучка от угла.
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ
Расчеты показывают, что при выбранной геометрии пинхола отклонение оси пучка на ±8,72°
вызывает отклонение интенсивности излучения не более 11 %. Это обеспечивает возможность
аппаратной или программной коррекции уровней сигналов в системе. На рис. 7 изображена кон-
струкция пинхола со всеми требуемыми геометрическими размерами.
Разработанный макет коллимационной системы для контроля стальных изделий представлен
на рис. 8.
Результаты сканирования стальной пластинки 1 рентгеновским излучением с эффективной
энергией 200 кэВ представлены на рис. 9.
Как видно из рис. 9, разница между максимальным уровнем сигнала In(α) и минимальным
составляет 20 %. Форма профиля экспериментального уровня на рис. 9 соответствует форме уров-
ня расчетного сигнала на рис. 6.
Дефектоскопия
№ 2
2020
48
Е.Е. Журавский, Б.И. Капранов, Д.С. Белкин и др.
1,24
1,01
Рис. 7. Конструкция пинхола.
1
2
3
Рис. 8. Общий вид коллимационной системы:
1 — стальная пластинка; 2 — первичный коллиматор (пинхол); 3 — коллимационная система.
1,2
1
0,8
0,6
0,4
0,2
0
-8,72
0
8,72
α, град
Рис. 9. Экспериментальные данные профиля стальной пластинки.
Дефектоскопия
№ 2
2020
Расчет параметров пинхола коллимационной системы альбедной компьютерной томографии...
49
ВЫВОДЫ
1. Величина поперечного размера рабочего пучка для углов от 0 до 8,72° изменяется не более
чем на 15 %, что незначительно влияет на пространственно-статистическую неопределенность.
2. Интенсивность излучения максимальна при угле 0° и минимальна при угле 8,72°. Разница
между максимальным уровнем интенсивности и минимальным составляет не более 11 %, что
легко компенсируется аппаратно.
3. Разработанная геометрия пинхола позволяет осуществлять сканирование объекта контроля
простым перемещением пинхола перпендикулярно оси рентгеновского пучка от анода трубки в
пределах 50 мм при фокусном расстоянии 163 мм.
4. Формы расчетного и экспериментального уровней сигналов сопоставляются между собой.
Разница между минимальным и максимальным уровнями расчетного и экспериментального сиг-
налов составляет 9 %.
5. Разработанная методика позволяет проводить расчет параметров пинхола первичного колли-
матора КС АКТ для заданной энергии и геометрии контроля.
СПИСОК ЛИТЕРАТУРЫ
1. Berecz T., Jenei P., Csóré A., Lábár J., Gubicza J., Szabó P.J. Determination of dislocation density by
electron backscatter diffraction and X-ray line profile analysis in ferrous lath martensite // Materials
Characterization. 2016. V. 113. P. 117—124.
2. Kolkoori S., Wrobel N., Zscherpel U., Ewert U. A new X-ray backscatter imaging technique for non-
destructive testing of aerospace materials // Ndt & E Internationa. 2015. V. 70. P. 41—52.
3. Margret M., Subramanian V., Baskaran R., Venkatraman B. Detection of scales and its thickness
determination in industrial pipes using Compton backscattering system // Review of Scientific Instruments.
2018. V. 89. No. 11. P. 113117.
4. O’Flynn D., Crews C., Fox N., Allen B.P., Sammons M., Speller R.D. X-ray backscatter sensing of
defects in carbon fibre composite materials / Advanced Photon Counting Techniques XI // International
Society for Optics and Photonics. 2017. V. 10212. P. 102120R.
5. Margret M., Menaka M., Venkatraman B., Chandrasekaran S. Compton back scatter imaging for mild
steel rebar detection and depth characterization embedded in concrete // Nuclear Instruments and Methods in
Physics Research Section B: Beam Interactions with Materials and Atoms. 2015. V. 343. P. 77—82.
6. Margret M., Menaka M., Subramanian V., Baskaran R., Venkatraman B. Non-destructive inspection of
hidden corrosion through Compton backscattering technique // Radiation Physics and Chemistry. 2018. V. 152.
P. 158—164.
7. Shinji Nomura, Kazunori Tejima, Ikuo Wakamoto. Scattered x-ray type defect detector and x-ray
detector. JP. Patent No. 2001208705A. 03 August 2001.
8. Kapranov B.I., Mudrov M.A. Contribution of Compton scattering in problems associated with measuring
the surface density of radiation protection coatings // J. of Physics: Conference Series. IOP Publishing. 2016.
V. 671. No. 1. P. 012038.
9. Kapranov B.I., Vavilova G.V., Volchkova A.V., Kuznetsova I.S. Mathematical modeling of tomographic
scanning of cylindrically shaped test objects // IOP Conference Series: Materials Science and Engineering. IOP
Publishing. 2018. V. 363. No. 1. P. 012015.
10. Капранов Б.И., Короткова И.А., Чахлов В.Л., Филинов В.Н., Маклашевский В.Я. Анализ систем
сканирования в комптоновской томографии. Возможные области применения // Изв. Томского политех.
ун-та. 2003. Т. 306. № 1. С. 122—127.
11. Balamesh A., Salloum M., Abdul-Majid S. Feasibility of a New Moving Collimator for Industrial
Backscatter Imaging // Research in Nondestructive Evaluation. 2018. V. 29. No. 3. P. 143—155.
12. Cui S., Nimmagadda J.K., Baciak J.E. Backscatter Radiography as a Non-Destructive Examination
Tool for Concrete Structures / 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference
(NSS/MIC). IEEE, 2017. P. 1—6.
Дефектоскопия
№ 2
2020