Расчёт поглощённых доз фотонного и нейтронного излучения в эмали и дентине зубов человека методом Монте-Карло

Автор: Хайлов А.М., Иванников А.И., Орленко С.П., Борышева Н.Б., Скворцов В.Г.

Журнал: Радиация и риск (Бюллетень Национального радиационно-эпидемиологического регистра) @radiation-and-risk

Статья в выпуске: 2 т.24, 2015 года.

Бесплатный доступ

С помощью расчётов методом Монте-Карло (программа MCNP-4С) получены значения доз, поглощённых в тканях зуба - эмали и дентине, а также в различных органах тела человека при внешнем облучении нейтронами с энергией от 10 -9 до 20 МэВ и фотонами с энергией от 0,05 до 10 МэВ в шести стандартных геометриях. В качестве фантома тела человека была использована одна из модификаций математического фантома MIRD с добавлением подробной дентальной области, включавшей в себя описание 32 зубов. Верификация путём сравнения значений доз в различных органах с литературными данными показала хорошее соответствие, для быстрых нейтронов и фотонов расхождения составили менее 7%. Приведены формулы для расчёта коэффициентов перехода от доз нейтронов и фотонов, поглощённых в эмали или дентине, к дозе в органах и эффективной дозе во всём теле человека. Полученные значения могут быть использованы для оценки индивидуальных доз облучения, а также дозового распределения в теле пострадавшего по значениям доз в дентине и эмали зубов, определённых методом ЭПР-спектроскопии, с использованием имеющейся информации о радиационном инциденте.

Еще

Эпр-спектроскопия, ретроспективная дозиметрия, дентин зубов, эмаль зубов, метод монте-карло, математический фантом, нейтронное излучение, индивидуальная доза, смешанное излучение, дозовые коэффициенты перехода

Короткий адрес: https://sciup.org/170170209

IDR: 170170209

Текст научной статьи Расчёт поглощённых доз фотонного и нейтронного излучения в эмали и дентине зубов человека методом Монте-Карло

Метод определения доз радиации на основе спектроскопии электронного парамагнитного резонанса образцов твёрдых тканей тела человека (эмаль и дентин зубов, кости, ногти и волосы) и некоторых материалов (стекло, пластик и др.) широко применяется для ретроспективной оценки доз фотонного излучения у пострадавших от атомных бомбардировок и испытаний ядерного оружия, ликвидаторов радиационных аварий, населения загрязнённых радиацией территорий Урала и Чернобыля, а также используется для оценки накопленной дозы от естественной фоновой радиации у населения, проживающего на «чистых» территориях России.

В результате проведённых ранее авторами экспериментальных исследований [1, 2] была подтверждена низкая (порядка 2-3%) относительная радиационная чувствительность эмали зубов к нейтронам [3, 4], а также исследованы дозиметрические характеристики дентина зубов человека к нейтронному и фотонному излучениям, в т.ч. при облучении образцов в фантоме головы in vivo . Также была показана возможность совместного использования эмали и дентина зубов для решения задачи корректной, точной и достоверной оценки поглощённых доз в тканях зуба методом ЭПР-спектроскопии в случае наличия в поле излучения значительного нейтронного компонента.

При этом для реконструкции величины индивидуальной поглощённой дозы нейтронного излучения, а также распределения дозы в теле пострадавшего необходимо осуществить переход от дозы, поглощённой в образцах эмали и дентина, к дозам в различных тканях и органах

Хайлов А.М.* – ст. научн. сотр., к.б.н.; Иванников А.И. – вед. научн. сотр., к.ф.-м.н.; Орленко С.П . – научн. сотр.; Борышева Н.Б. – ст. научн. сотр., к.ф.-м.н.; Скворцов В.Г. – зав. лаб., к.б.н. МРНЦ им. А.Ф. Цыба – филиал ФГБУ «НМИРЦ» Минздрава России.

тела человека. Для расчёта коэффициентов перехода, связывающих дозу в тканях зуба с дозой в органах и всём теле, в свою очередь, необходимо рассчитать значения доз, поглощённых в перечисленных тканях и органах с использованием расчётных программ, основанных на методе Монте-Карло и реалистичных фантомов человека [5, 6]. Поскольку нейтронное излучение на практике всегда сопровождается гамма-излучением, необходимо наличие коэффициентов перехода от поглощённых доз в тканях зубов к дозам в различных тканях и органах тела человека при первичном облучении фотонами.

К настоящему времени в литературных источниках представлены результаты целого ряда работ, посвящённых расчёту доз в тканях зубов, а также расчёту коэффициентов перехода различного типа при внешнем фотонном и нейтронном облучении.

В частности, опубликованы результаты расчётов доз на единичный флюенс в тканях зубов при внешнем фотонном облучении с использованием математического фантома [7] и подробного воксельного фантома человека Голем (Golem) [8]. Несмотря на фундаментальный характер обеих работ и большой объём представленных данных, расчёт дозовых величин при внешнем нейтронном облучении в них выполнен не был.

В работе [9] приведены значения отношений доз, поглощённых в органах и всём теле, а также доз в тканях зубов – дентине и эмали, рассчитанные для случаев внешнего облучения нейтронами, фотонами и облучения радионуклидными источниками (137Cs, 60Co, 252Cf и другие). Однако данные отношения приведены для суммарной поглощённой дозы нейтронов, включающей как непосредственно создаваемую нейтронами в процессе замедления дозу, так и дозу, обусловленную вторичным гамма-излучением. При этом разделение суммарной дозы нейтронов на два перечисленных компонента необходимо для корректной интерпретации результатов измерений спектров ЭПР зубных тканей. Это вызвано большими различиями в значениях радиационной чувствительности эмали и дентина к компоненту суммарной дозы нейтронов, обусловленного непосредственно нейтронными взаимодействиями, и к вторичному гамма-излучению.

Авторами также ранее был проведён расчёт доз нейтронов в зубной эмали в диапазоне энергий 10-9–20 МэВ с использованием упрощённой модели дентальной области [10]. Дозы нейтронов в эмали были рассчитаны с разделением её на дозу непосредственно нейтронов и дозу вторичных фотонов. Однако, к сожалению, расчёт доз нейтронов в дентине зубов в данном исследовании нами не проводился.

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

Материалы и методы

Расчёт доз в тканях зубов и органах тела человека в работе проводился путём моделирования облучения математического фантома тела человека, согласно стандартам Международной комиссии по радиологической защите (МКРЗ) с помощью программы MCNP-4С, основанной на методе Монте-Карло [11].

В качестве базового фантома тела человека была использована одна из модификаций широко применяемого в радионуклидной диагностике и терапии математического фантома MIRD (Medical Internal Radiation Dose) [12]. Для целей нашего исследования данный фантом (рис. 1а) был модифицирован путём добавления в него подробной дентальной области, которая включала в себя 32 зуба (рис. 1б, в). Описание геометрии отдельных зубов было проведено с использованием модели зуба, максимально точно учитывающей взаимное расположение всех его компонентов – эмали, дентина и пульпы. Подробное описание фантома зуба ранее приведено авторами в работе [1].

Рис. 1. Внешний вид фантома человека (а), использованного при расчётах, взаимного расположения костей головы (б), а также дентальной области (в).

Вычисления доз были проведены для случаев общего облучения мононаправленными изотропными полями нейтронов энергией от 10-9 до 20 МэВ и фотонов энергией от 0,05 до 10 МэВ. Вычисления доз нейтронов выполнялись в режиме комбинированного нейтроннофотонного транспорта (NP), где вторичные фотоны производятся в ходе нейтронных взаимодействий. При этом компонент дозы, обусловленный нейтронными взаимодействиями, а также доза вторичных фотонов в дентине и эмали коренных зубов и во всех органах (за исключением кожи) были рассчитаны в керма-приближении, поскольку при облучении соблюдалось равновесие вторичных заряженных частиц (протонов и ядер отдачи, а также вторичных электронов). Доза в коже была рассчитана с учётом транспорта вторичных электронов для случаев облучения как нейтронами, так и фотонами. Для учёта эффекта химического связывания атомов водорода при энергиях нейтронов меньше 1 эВ использовались таблицы сечений теплового рассеяния (S(α, β)). Для определения поглощённых доз в тканях зубов и различных областях фантома применялся тип расчёта F6 (energy deposition tally), для определения поглощённых доз в коже – тип расчёта *F8 (pulse height tally). В ходе моделирования транспорта нейтронов нами была ис- пользована библиотека сечений ENDF60, для транспорта фотонов использовалась библиотека сечений MCPLIB02, для электронного транспорта – библиотека EL03. Время расчётов выбиралось для достижения заданной относительной погрешности значений (приведены далее).

В целях как можно более широкого практического использования результатов, расчёт был проведён для целого ряда стандартных геометрий облучения (рис. 2): передне-задняя (AP), задне-передняя (PA), правая боковая (RLAT), левая боковая (LLAT) ротационная (ROT) и изотропная геометрии (ISO) облучения человека в воздухе, что делает возможным их использование как для случаев разового переоблучения единичными источниками, так и для равномерного облучения рассеянным внешним излучением радионуклидов.

АР РА RLAT LLAT ROT ISO

Рис. 2. Геометрии внешнего облучения фантома (согласно [6]).

Рис. 3. Поперечный срез головы фантома с указанием расположения зубов во рту и их порядковых номеров.

Результаты и обсуждение

Внешнее нейтронное облучение

Поскольку полное представление значений доз в тканях всех зубов челюсти невозможно ввиду ограничений на размер публикации, дозы в эмали и дентине представлены в табл. 1 и 2 для коренного зуба или моляра под номером L1 (согласно рис. 3), который наиболее часто используется для целей ЭПР-дозиметрии. Значения доз нейтронного и фотонного компонентов, а также суммарной дозы нейтронов нормированы на единичный флюенс нейтронов источника и представлены в единицах пГр·см2. Относительная погрешность результатов не превышала 3% для нейтронного компонента и 5% для фотонного компонента суммарной дозы нейтронов. В случаях расчёта доз медленных и промежуточных нейтронов в геометриях, в которых эмаль наиболее удалена от потока нейтронов (PA и RLAT), дозу нейтронов рассчитывали до погрешности 10%.

Таблица 1 Поглощённая доза на единичный флюенс в эмали зуба под номером L1 (согласно рис. 3), в единицах пГр·см² при моноэнергетическом облучении нейтронами различных энергий. В столбиках значений сверху вниз представлены: нейтронный компонент, фотонный компонент и суммарная доза нейтронов соответственно

Энергия (МэВ)	Геометрия облучения
Энергия (МэВ)	AP	PA	ISO	LLAT	RLAT	ROT
1E-9	0,0008 0,587 0,588	0,0002 0,351 0,351	0,0005 0,454 0,454	0,0018 0,819 0,821	0,0002 0,290 0,290	0,0008 0,512 0,513
1E-8	0,0013 0,806 0,807	0,0002 0,396 0,396	0,0009 0,557 0,557	0,0032 1,18 1,19	0,0002 0,357 0,357	0,0013 0,686 0,687
1E-7	0,0026 1,22 1,22	0,0004 0,573 0,574	0,0015 0,724 0,726	0,0065 1,86 1,86	0,0004 0,497 0,497	0,0025 1,04 1,04
1E-6	0,0046 1,71 1,71	0,0007 0,810 0,811	0,0023 1,01 1,01	0,0107 2,58 2,59	0,0007 0,673 0,674	0,0042 1,44 1,45
1E-5	0,0047 1,61 1,61	0,0009 0,830 0,831	0,0025 0,946 0,949	0,0094 2,19 2,20	0,0009 0,706 0,707	0,0040 1,33 1,34
1E-4	0,0050 1,55 1,56	0,0010 0,813 0,814	0,0025 0,974 0,977	0,0097 2,06 2,07	0,0010 0,732 0,733	0,0042 1,29 1,29
1E-3	0,0052 1,52 1,52	0,0011 0,808 0,809	0,0030 0,909 0,912	0,0104 1,87 1,88	0,0011 0,690 0,692	0,0045 1,22 1,23
0,01	0,0070 1,44 1,45	0,0012 0,780 0,782	0,0044 0,851 0,855	0,0169 1,80 1,82	0,0012 0,706 0,707	0,0066 1,18 1,19
0,1	0,041 1,54 1,58	0,003 0,886 0,889	0,024 0,813 0,837	0,117 1,73 1,84	0,003 0,774 0,777	0,041 1,23 1,27
0,5	0,411 1,62 2,03	0,033 1,07 1,10	0,203 1,00 1,20	0,821 1,70 2,52	0,020 1,00 1,02	0,321 1,35 1,67
1	0,903 1,37 2,28	0,111 0,973 1,08	0,476 0,928 1,40	1,77 1,36 3,13	0,087 0,910 0,997	0,717 1,15 1,87
2	1,79 1,43 3,22	0,541 1,13 1,67	0,940 1,02 1,96	2,42 1,36 3,77	0,474 1,09 1,56	1,30 1,25 2,56
4	4,38 1,10 5,47	1,85 0,957 2,81	2,55 0,747 3,30	5,53 0,981 6,51	1,68 0,864 2,54	3,36 0,975 4,33
6	7,46 1,14 8,60	3,88 0,986 4,86	4,60 0,846 5,45	8,93 1,11 10,0	3,50 0,864 4,36	5,94 1,03 6,97
8	9,19 1,86 11,1	5,54 1,50 7,04	6,14 1,48 7,62	10,5 1,87 12,4	5,17 1,43 6,60	7,61 1,67 9,28

Энергия	Геометрия облучения
(МэВ)	AP	PA	ISO	LLAT	RLAT	ROT
	12,1	7,48	8,08	13,9	6,91	10,1
10	2,42	1,88	1,85	2,43	1,78	2,13
	14,5	9,36	9,93	16,3	8,69	12,2
	14,5	9,35	9,91	16,5	8,68	12,2
12	2,44	2,01	1,95	2,46	1,92	2,21
	16,9	11,4	11,9	18,9	10,6	14,5
	16,6	11,2	11,7	18,8	10,4	14,3
14	2,53	2,04	1,89	2,62	2,05	2,31
	19,2	13,2	13,6	21,4	12,5	16,6
	18,1	12,6	12,9	20,4	11,9	15,7
16	2,29	1,98	1,87	2,37	1,83	2,12
	20,4	14,6	14,8	22,8	13,7	17,9
	19,6	14,1	14,1	21,9	13,5	17,3
18	1,92	1,67	1,64	1,94	1,58	1,78
	21,5	15,7	15,7	23,9	15,1	19,0
	22,0	16,3	16,3	24,4	15,6	19,6
20	1,55	1,38	1,42	1,58	1,35	1,47
	23,5	17,7	17,7	26,0	16,9	21,0

Таблица 2

Поглощённая доза на единичный флюенс в дентине зуба под номером L1 (согласно рис. 3), в единицах пГр·см² при моноэнергетическом облучении нейтронами различных энергий. В столбиках значений сверху вниз представлены: нейтронный компонент, фотонный компонент и суммарная доза нейтронов соответственно

Энергия	Геометрия облучения
(МэВ)	AP	PA	ISO	LLAT	RLAT	ROT
	0,0005	0,0001	0,0003	0,0011	0,0001	0,0005
1E-9	0,615	0,370	0,510	0,834	0,297	0,529
	0,615	0,370	0,510	0,835	0,297	0,529
	0,0008	0,0001	0,0005	0,0020	0,0001	0,0008
1E-8	0,798	0,377	0,617	1,25	0,355	0,696
	0,799	0,377	0,617	1,26	0,356	0,697
	0,0016	0,0002	0,0012	0,0040	0,0003	0,0015
1E-7	1,30	0,583	0,979	1,99	0,532	1,10
	1,31	0,584	0,980	2,00	0,532	1,10
	0,0028	0,0004	0,0018	0,0064	0,0004	0,0025
1E-6	1,83	0,845	1,26	2,64	0,685	1,50
	1,84	0,845	1,26	2,65	0,685	1,50
	0,0029	0,0006	0,0019	0,0058	0,0006	0,0025
1E-5	1,68	0,842	1,12	2,28	0,718	1,38
	1,69	0,843	1,12	2,29	0,719	1,38
	0,0031	0,0006	0,0022	0,0065	0,0007	0,0027
1E-4	1,62	0,844	1,10	2,18	0,746	1,35
	1,63	0,844	1,10	2,19	0,747	1,35
	0,0044	0,0008	0,0035	0,0122	0,0008	0,0045
1E-3	1,57	0,848	1,16	1,97	0,713	1,27
	1,58	0,848	1,17	1,98	0,714	1,28
	0,0182	0,0017	0,0200	0,0738	0,0020	0,0239
0,01	1,48	0,815	1,04	1,93	0,717	1,24
	1,50	0,816	1,06	2,01	0,719	1,26
	0,304	0,017	0,287	1,02	0,014	0,338
0,1	1,60	0,903	1,01	1,82	0,829	1,29
	1,90	0,920	1,30	2,83	0,843	1,63
	2,91	0,230	2,18	6,05	0,125	2,33
0,5	1,71	1,10	1,10	1,77	1,02	1,40
	4,62	1,33	3,28	7,83	1,15	3,73

Энергия (МэВ)

Геометрия облучения

ISO

LLAT

RLAT

ROT

5,19

1,41

6,60

0,767

0,991 1,76

3,90 1,01 4,91

9,88

1,39

11,3

0,567 0,932 1,50

4,10

1,18

5,28

12,2

1,49

13,7

3,67

1,18

4,86

8,59 1,01 9,60

16,8

1,36

18,1

3,25

1,10

4,35

8,97

1,28

10,3

20,5

1,13

21,7

9,2 1,00 10,2

15,5 0,830

16,3

26,1 1,00 27,2

8,7 0,914 9,62

16,1

1,01

17,2

26,1

1,18

27,2

14,3

1,02

15,3

20,2 0,854 21,1

31,8

1,14

32,9

13,3 0,913 14,2

21,4

1,06

22,4

29,6

1,87

31,5

18,4

1,50

19,9

23,5

1,50

24,9

35,2

1,84

37,0

17,5

1,56

19,0

25,2

1,69

26,9

33,3

2,39

35,7

21,5

1,89

23,3

26,8

1,89

28,7

38,9

2,43

41,3

20,4

1,78

22,1

28,5

2,12

30,6

36,9

2,44

39,3

24,6

2,03

26,6

30,3

1,74

32,1

42,6

2,48

45,1

23,6

1,91

25,5

31,9

2,21

34,1

39,3

2,50

41,8

27,6

2,03

29,6

32,6

2,03

34,6

45,4

2,61

48,0

26,4

2,05

28,5

34,7

2,30

37,0

42,9

2,28

45,2

30,6

1,94

32,5

35,7

1,99

37,7

49,1

2,38

51,5

29,8

1,87

31,7

38,1

2,12

40,2

44,6

1,94

46,5

33,3

1,65

35,0

38,0

1,68

39,6

51,0

1,90

52,9

32,2

1,60

33,8

40,3

1,77

42,1

47,2

1,53

48,8

36,5

1,35

37,8

41,7

1,35

43,0

53,6

1,57

55,2

35,6

1,31

37,0

43,2

1,44

44,7

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

^Дn ⁼ ^Дnден ^/ ^kден ^,

Д n = Д n эм / k эм , где kден

ли ( K_эм

= K / K , k = K / K – отношения кермы нейтронов в дентине ( K ) или эма-

) и кермы нейтронов в мягкой биологической ткани ( K ).

С целью определения величин k и k , нами были рассчитаны значения кермы нейтронов в дентине, эмали и мягкой биологической ткани путём суммирования соответствующих значений кермы всех элементов [13], входящих в их состав. Рассчитанные по приведённым выше формулам значения величин k и k даны в табл. 3.

Таблица 3

Энергетическая зависимость отношений кермы нейтронов в тканях зуба (эмали и дентине) и кермы нейтронов в мягкой биологической ткани

Энергия (МэВ)	^эм	^ден
1E-9	0,080	0,013
1E-8	0,084	0,022
1E-7	0,125	0,044
1E-6	0,198	0,089
1E-5	0,415	0,274
1E-4	0,334	0,572
1E-3	0,088	0,538
0,01	0,056	0,523
0,1	0,058	0,524
0,5	0,064	0,529
1	0,129	0,531
2	0,073	0,534
4	0,137	0,574
6	0,192	0,595
8	0,207	0,610
10	0,247	0,610
12	0,260	0,622
14	0,276	0,635
16	0,295	0,645
18	0,301	0,660
20	0,305	0,656

Чувствительность тканей зубов к фотонному излучению низкой энергии (до 300 кэВ) превосходит их чувствительность к фотонам калибровочного излучения изотопов 60Со или 137Cs примерно на порядок [14]. Однако при облучении нейтронами доля вторичных фотонов с энергией менее 300 кэВ составляет менее 1% [10]. Следовательно, дозу вторичных фотонов в тканях зуба можно считать численно равной дозе в мягкой биологической ткани в месте расположения зуба.

Внешнее фотонное облучение

В табл. 4 и 5 представлены зависимости дозы фотонов в эмали и дентине моляра под номером L1 (согласно рис. 3) от энергии фотонов. Значения доз нормированы на единичный флюенс фотонов источника и приведены в единицах пГр·см2. Относительная погрешность вычисленных методом Монте-Карло значений доз в эмали и дентине не превышала 3% для энергий фотонов более 0,25 МэВ и 5% для энергий фотонов в диапазоне 0,05-0,25 МэВ.

Таблица 4 Поглощённая доза на единичный флюенс в эмали зуба номер L1 (согласно рис. 3), в единицах пГр^см² при моноэнергетическом облучении фотонами различных энергий

Энергия (МэВ)	Геометрия облучения
Энергия (МэВ)	AP	PA	ISO	LLAT	RLAT	ROT
0,05	1,14	0,436	0,608	2,09	0,309	0,994
0,1	1,22	0,614	0,670	1,65	0,541	1,01
0,14	1,24	0,676	0,684	1,55	0,596	1,02
0,25	1,40	0,902	0,832	1,73	0,844	1,22
0,4	1,97	1,28	1,32	2,33	1,22	1,70
0,67	2,75	2,04	1,90	3,41	1,93	2,53
0,85	3,52	2,61	2,51	4,08	2,45	3,17

Энергия (МэВ)	Геометрия облучения
Энергия (МэВ)	AP	PA	ISO	LLAT	RLAT	ROT
1,25	4,55	3,77	3,28	5,38	3,63	4,33
1,75	5,40	4,43	4,26	6,66	4,80	5,32
2,5	6,98	6,26	5,68	8,72	6,22	7,05
3,5	9,14	8,11	7,29	10,7	8,36	9,09
5	11,3	10,1	9,59	14,1	10,1	11,4
7	14,4	12,9	12,3	17,9	12,8	14,5
10	19,2	17,2	17,2	23,9	17,1	19,3

Таблица 5

Поглощённая доза на единичный флюенс в дентине зуба номер L1 (согласно рис. 3), в единицах пГр·см2 при моноэнергетическом облучении фотонами различных энергий

Энергия (МэВ)	Геометрия облучения
Энергия (МэВ)	AP	PA	ISO	LLAT	RLAT	ROT
0,05	0,621	0,263	0,364	1,25	0,178	0,579
0,1	0,843	0,397	0,457	1,15	0,355	0,685
0,14	0,867	0,505	0,494	1,18	0,422	0,744
0,25	1,26	0,794	0,755	1,56	0,730	1,08
0,4	1,77	1,19	1,23	2,24	1,11	1,58
0,67	2,74	1,97	1,88	3,31	2,00	2,50
0,85	3,66	2,58	2,53	4,12	2,44	3,20
1,25	4,50	3,81	3,43	5,58	3,68	4,39
1,75	5,33	4,90	4,38	7,12	5,11	5,61
2,5	7,96	7,00	6,28	9,26	6,49	7,68
3,5	10,6	9,23	8,08	11,6	8,97	10,1
5	12,9	11,3	10,2	15,0	10,5	12,4
7	16,4	14,4	13,4	19,0	13,3	15,8
10	21,8	19,2	18,7	25,4	17,8	21,1

Расчёт дозы в органах и всём теле человека

В ходе проведённых методом Монте-Карло расчётов нами также были получены значения суммарной дозы нейтронов и доз фотонов на единичный флюенс в различных критических органах тела человека, для которых производится расчёт эффективной дозы. Сравнение рассчитанных поглощённых доз нейтронов и фотонов в различных органах с аналогичными значениями, приведёнными в публикации МКРЗ [6], показало хорошее соответствие. Для случая облучения нейтронами расхождение подавляющего большинства результатов составило менее 8-12% для медленных и промежуточных нейтронов, 5-7% для быстрых нейтронов (в зависимости от размера органа). Для случая облучения фотонами расхождение значений составило менее 5-7% для энергий фотонов в диапазоне 0,05-0,4 МэВ и менее 5% для фотонов с энергиями 0,67 МэВ и более. Поскольку расхождение значений доз в органах, рассчитанных в нашем исследовании, и значений, приведённых в литературе, в большинстве случаев было меньше погрешности вычислений, рассчитанные нами данные не приводятся, а сами значения доз в органах и эффективной дозы в теле человека могут напрямую быть взяты из [6].

Вычисление коэффициентов перехода от дозы в эмали и дентине к дозе в органах и всём теле человека

На предыдущих этапах работы с помощью методов стохастического моделирования транспорта ионизирующего излучения в математическом фантоме человека были рассчитаны дозы, поглощённые в эмали, и дентина коренных зубов, нормированные на основную характеристику поля излучения - флюенс.

Для определения дозового распределения в теле пострадавшего, введём определение коэффициента перехода от дозы нейтронов, поглощённой в эмали или дентине, к дозе в органах и эффективной дозе во всём теле человека. Для случая нейтронного облучения коэффициент перехода ( С _п ) может быть определён по известной геометрии облучения и месту положения зуба как:

C n эм

Д 1Ф n орг , n сум эм

C n ден

Д 1Ф n орг n сум ден

где Дп сум эм / Ф , Д п сум d^ IФ - суммарная поглощённая доза нейтронов в эмали и дентине на единичный флюенс (приведена выше для коренных зубов в табл. 1-2); Дп о/Ф - суммарная поглощённая доза нейтронов в органе на единичный флюенс (приведена в табл. А.26-А.41 Приложения 1 в [6]).

Аналогично, для фотонного компонента поля излучения:

П 1Ф П 1Ф

C y эм

^Д Y орг — ^Д y орг

— , . , C y Ден = — , . ,

Д Y1 эм 1 Ф Д y 1 Ден 1 Ф где Дy1 эм / Ф , Дy1 dен / Ф - поглощённая доза фотонов в эмали или дентине на единичный флюенс (приведена выше для коренных зубов в табл. 4-5); Дy орг /Ф - поглощённая доза фотонов в органе на единичный флюенс.

Поглощённая доза фотонов в органе может быть получена простым умножением соответствующих значений доз, нормированных на керму в воздухе Д_{п ора} / к_а в табл. А.2-А.20 Приложения 1 [6] на значения кермы в воздухе, нормированные на единичный флюенс фотонов к_а / ф в табл. А.1 Приложения 1 [6] согласно формуле:

^Д y орг ^Д y орг K а

-------= -------X --

Ф KФ

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

Заключение

Полученные в результате проведённых расчётов данные могут быть использованы для перехода от поглощённой дозы в местах расположения зубов к дозам в различных тканях и органах тела человека в целях корректной и точной реконструкции методом ЭПР-дозиметрии индивидуальной поглощённой дозы и её распределения в теле пострадавшего от смешанного облучения с присутствием нейтронного компонента.

При этом, учёт ранее полученных значений радиочувствительности эмали и дентина к нейтронам [1, 2], сделает возможным корректную оценку дозовых величин, обусловленных воздействием смешанного излучения на организм человека даже в случае, когда доза каждого из компонентов поля изначально не известна.

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 12-04-31645 мол_а.

Авторы также выражают признательность за финансирование в рамках пилотного проекта номер U19-AI091173, программы Дартмутского Центра физических методов по исследованию медицинских контрмер против радиации (Pilot Project Program of the Dartmouth Physically Based Center for Medical Countermeasures Against Radiation), полученное от Национального института аллергии и инфекционных заболеваний (National Institute of Allergy and Infectious Diseases), входящего в Национальные институты здравоохранения (NIH).

Список литературы Расчёт поглощённых доз фотонного и нейтронного излучения в эмали и дентине зубов человека методом Монте-Карло

Хайлов А.М., Иванников А.И., Скворцов В.Г. Оценка дозы смешанного гамма-нейтронного излучения по результатам спектроскопии ЭПР образцов зубной эмали и дентина при облучении в фантоме головы человека//Радиация и риск. 2013. Т. 22, № 4. С. 30-36.
Хайлов А.М., Орленко С.П., Скворцов В.Г. Определение методом спектроскопии электронного парамагнитного резонанса радиационной чувствительности дентина зубов человека in vitro к нейтронному и фотонному излучениям//Радиация и риск. 2013. Т. 22, № 4. С. 21-29.
Khailov A., Ivannikov A., Tikunov D., Skvortsov V., Stepanenko V., Zhumadilov K., Tanaka K., Endo S., Hoshi M. Tooth enamel EPR dosimetry of neutrons: Enhancement of the apparent sensitivity at irradiation in the human head phantom//Radiat. Meas. 2007. V. 42, N 6-7. P. 1171-1177.
Tikunov D., Trompier F., Ivannikov A., Herve M., Khailov A., Skvortsov V. Relative Sensitivity of Tooth Enamel to Fission Neutrons: Effect of Secondary Protons//Radiat. Meas. 2005. V. 39, N 5. P. 509-514.
International Commission on Radiological Protection. 1990 Recommendations of the International Commission on Radiological Protection. ICRP Publication 60. Oxford: Pergamon press, 1991.
International Commission on Radiological Protection. Conversion Coefficients for use in Radiological Protection against External Radiation. ICRP Publication 74, Annals of the ICRP(3-4). Oxford: Elsevier Science, 1997.
Takahashi F., Yamaguchi Y., Iwasaki M., Miyazawa C., Hamada T. Relations Between Tooth Enamel Dose and Organ Doses for Electron Spin Resonance Dosimetry Against External Photon Exposure//Radiat. Prot. Dosim. 2001. V. 95, N 2. P. 101-108.
Ulanovsky A., Wieser A., Zankl M., Jacob P. Photon dose conversion coefficients for human teeth in standard irradiation geometries//Health Phys. 2005. V. 89, N 6. P. 645-659.
Herve M.L., Clairand I., Trompier F., Tikunov D., Bottollier-Depois J.F. Relation between organ and whole body doses and local doses measured by ESR for standard and realistic neutron and photon external overexposures//Radiat. Prot. Dosim. 2007. V. 125, N 1-4. P. 355-360.
Khailov A., Ivannikov A., Skvortsov V., Stepanenko V., Tsyb A.F., Trompier F., Hoshi M. The neutron dose conversion coefficients calculation in human tooth enamel in an anthropomorphic phantom//Health Phys. 2010. V. 98, N 2. P. 369-377.
Briesmeister J.F. MCNP -A general Monte Carlo N-particle transport code. Version 4C. Report LA-13709-M. Los Alamos National Laboratory. Los Alamos, 2000.
Cristy M., Eckerman F. Specific absorbed fraction of energy at various ages from internal photon sources. Appendix A. Description of the mathematical phantom. ORNL/TM-8381/V1. Oak Ridge: Oak Ridge National Laboratory, 1997.
Caswell R.S., Coyne J.J., Randolph M.L. Kerma factors of elements and compounds for neutron energies below 30 MeV//Appl. Radiat. Isot. 1981. V. 33. P. 1227-1262.
Ivannikov A.I., Tikunov D.D., Borysheva N.B., Skvortsov V.G., Stepanenko V.F., Trompier F., Hoshi M. Calibration of EPR signal dose response of tooth enamel to photons: experiment and Monte Carlo simulation//Radiat. Prot. Dosim. 2004. V. 108, N 4. P. 303-315.

Еще

Статья научная