Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 2288-3509(Print)
ISSN : 2384-1168(Online)
Journal of Radiological Science and Technology Vol.39 No.1 pp.89-97
DOI : https://doi.org/10.17946/JRST.2016.39.1.11

Characterization of the Neutron for Linear Accelerator Shielding Wall using a Monte Carlo Simulation

Dong Yeon Lee1), Eun Tae Park2), Jung Hoon Kim3)
1)Dept. of Radiation Oncology, Dongnam Inst. of Radiological & Medical Scienc
2)Dept. of Radiation Oncology, Inje University Busan Paik Hospital
3)Dept. of Radiologic Science, College of Health Sciences, Catholic University of Pusan
Corresponding author: Kim Jung Hoon, 57 Oryundae-ro, Geumjeong-gu, Busan, Rep. of KOREA. Catholic University of Pusan. Dept. of Radiologic Science, College of Health Sciences. / 82-51-510-0583
January 18, 2016 March 8, 2016 March 15, 2016

Abstract

As previous studies to proceed with the evaluation of the radioactive at linear accelerator's shielding concrete wall. And the shielding wall was evaluated the characteristics for the incoming neutron. As a result, the shielding wall is the average amount of incoming neutrons 10 MV 4.63E-7%, 15 MV 9.69E-6%, showed the occurrence of 20 MV 2.18E-5%. The proportion of thermal neutrons of which are found to be approximately 18-33%.

The neutron generation rate can be seen as a slight numerical order. However, in consideration of the linear accelerator operating time we can not ignore the effects of neutrons. Accordingly radioactive problem of the radiation shield wall of the treatment room will be this should be considered.


몬테칼로시뮬레이션을 이용한 선형가속기 차폐벽에 대한 중성자 특성 평가

이 동연1), 박 은태2), 김 정훈3)
1)동남권원자력의학원 방사선종양학과
2)인제대학교 부산백병원 방사선종양학과
3)부산가톨릭대학교 보건과학대학 방사선학과

초록

선형가속기를 차폐하고 있는 콘크리트 차폐벽의 방사화 평가를 진행하기 위한 선행 연구로서 차폐벽이 받는 중 성자에 대한 특성을 평가하였다. 그 결과 차폐벽이 받는 중성자 평균 양은 생성된 광자 대비 10 MV 4.63E-7%, 15 MV 9.69E-6%, 20 MV 2.18E-5%의 발생비율을 보였으며, 이 중 반응단면적이 큰 열중성자의 비율은 대략 18 ~ 33%로 나타났다. 의료용 선형가속기의 경우 수년간 지속적으로 가동되기 때문에 중성자 생성 비율에 관계없이 누 적으로 생성되는 중성자의 양은 무시할 수 없을 것이며, 이에 따른 방사선 치료실 차폐벽의 방사화 문제는 반드시 고려되어야 할 것이다.


    I서 론

    방사선치료에 이용되는 의료용선형가속기는 전자를 가속 시켜 광자선을 얻는 대표적인 치료 장비이며, 임상에서 주 로 사용되는 광자에너지는 6-24 MV 영역이다1).

    현재에는 방사선치료의 목표인 종양만을 선택적으로 사 멸시키고 정상조직의 피해는 최소화하기 위한 노력의 일환으 로 10 MV 이상의 높은 에너지를 선호하고 있다2). 그러나 미 국방사선방호측정심의회 (National Council on Radiation Protection and Measurement; NCRP)에서는 NCRP Report No.79를 통해 높은 에너지의 광자선이 물질과 반응할 때, 그 물질 고유의 임계에너지보다 높을 경우에는 광핵반응으 로 인한 중성자가 발생한다고 보고하였으며3,4), NCRP Report No.151을 통해 높은 에너지의 광자선을 이용하는 선형가속기에서는 중성자의 영향을 고려하여 차폐를 할 것 을 권고하였다5). 위와 같이 높은 에너지의 광자선 사용 시 생성되는 중성자는 치료를 받는 환자와 방사선작업종사자 에게 예상치 못한 피폭을 가져올 수 있다. 이에 방사선방호 를 위하여 치료실내의 중성자 분포, 선량 및 특성에 관한 연 구가 활발히 진행되고 있다6-9). 하지만 발생된 중성자로 인 한 2차적인 방사화와 관련한 연구는 아직까지 미흡한 실정 이다. 최근의 연구에 의하면 치료실내에서 생성된 중성자 의 스펙트럼을 분석한 결과, 전체 중성자에서 핵반응을 일 으킬 확률이 높은 열중성자가 차지하는 비율이 약 20% 정 도 차지한다고 보고하였다1). 이에 따라 선형가속기 내의 부 속품(전자총, 가속관, 편향전자석, 타깃, 조리개, 선속평탄 여과판 등)과 선형가속기를 차폐하고 있는 콘크리트 벽이 방사화가 될 가능성이 충분히 있으며, 그로 인한 오염을 고려하여야 할 것이다.

    일반적으로 선형가속기가 교체되더라도 방사선치료실은 20년 이상 계속해서 가동하는 경우가 많기 때문에 방사화 로 인한 오염은 더욱 위험할 수 있다. 특히 방사화 된 물질 은 지속적으로 방사선이 방출되기 때문에 장비의 유지, 보 수, 해체를 진행 할 경우에는 전체적인 안전방침에 대한 논 의가 있어야 할 것이다. 이와 관련하여 국제원자력기구 (International Atomic Energy Agency; IAEA)에서는 의 료용 소형 연구장비 해체에 대한 권고안을 마련하였지만, 의료용 선형가속기와 방사선치료실의 방사화에 한정된 더 실제적인 연구가 요구되고 있다10).

    이에 본 연구는 치료실내의 방사화에 대한 평가를 진행하 기 위하여 필요한 중성자를 측정하였으며, 특히 치료실을 차폐하고 있는 콘크리트에 대하여 중성자 스펙트럼과 선속 을 측정함으로서 추후 방사화 평가를 진행함에 있어 기초가 되는 중성자 선속에 대한 자료를 제시하고자 한다.

    Ⅱ대상 및 방법

    본 연구에서는 몬테카를로 방법을 사용하여 모의실험을 진행하였다. 현재 상용되는 모의실험 프로그램으로는 Fluka (FLUktuierende KAskade), Geant4 (GEometry ANd Tracking), MCNPX (Monte Carlo N-Patricle Extended) 등이 있으며, 이 가운데 윈도우운영체제를 기반으로 함으로 써 접근성이 용이하고 입자 수송측면에 있어서 장점을 가지 고 있는 MCNPX (Ver.2.5.0)을 사용하여 모의실험을 수행 하였다.

    1선형가속기 모의모사

    본 모의실험에서 모사한 선형가속기는 일반적으로 치료 에 사용되고 있는 선형가속기 장비를 단순화시켜 모사한 것 으로 가속관을 제외한 헤드부분을 중심으로 모사하였으며, 그 모습은 Fig. 1과 같다. 1차 조리개 (Primary collimator), 타깃 (Target), 선속평탄여과판 (Flattening filter), 2차 조리개 (Secondary collimator), 다엽콜리메이터 (Multileaf collimator)를 세부적으로 나타냈으며, 사용한 재질과 밀도 그리고 기하학적 구조에 대한 세부사항은 Table 1.으로 설 명하였다.

    2차폐벽 모사

    선형가속기를 둘러싸고 있는 방사선치료실의 차폐벽은 IAEA11)에서 권고하는 형태를 기본으로 모사하였으며, 간략 한 구성도는 Fig. 2와 같다. 세부적으로 살펴보면 선형가속 기 헤드를 중심으로 좌측과 우측으로 각각 300 ㎝, 바닥 170 ㎝, 천장 100 ㎝, 머리방향 150 ㎝, 출입구 방향 650 ㎝, 통로 (maze) 550 ㎝의 거리를 설정함으로서 전신방사선 조사 (Total Body Irradiation, TBI)와 수술 중 방사선치료 (Intra-opera tive Procedures, IORT)를 시행할 수 있는 최소한의 크기를 확보하였다. 또한 차폐벽의 두께는 100 ㎝ 로 설정하였으며, 선형가속기의 회전으로 인하여 1차선속이 조사될 수 있는 좌측과 우측 방향은 50 ㎝ 두께를 추가하였 다11). 차폐벽을 이루고 있는 콘크리트 재질은 한국원자력연 구소에서 제공하는 재질 중 상용화 되어있는 Ordinary (NBS 03)을 바탕으로 설정하였다12).

    3모의실험 방법

    본 연구에 사용된 전자선의 에너지는 현재 치료에 사용하 는 에너지 중 광핵반응을 일으킬 확률이 높은 10, 15, 20 MeV의 에너지를 대상으로 모의실험을 진행하였다.

    먼저 선형가속기 헤드에서 제동복사로 인해 발생되는 광 자를 측정하였다. 측정 위치는 타깃 바로 아래 5 ㎝ 위치에 반경 3 ㎝인 가상의 구를 설정하여 전자 (e) 1개당 생성된 광자 플럭스 (photon number/㎠·sec·e)를 측정하였으며, 모사한 선형가속기의 신뢰성을 검증하기 위해 10 keV 간격 으로 광자스펙트럼을 측정하였다. 이 때 출력인자 (Tally)는 F5를 사용하였다13).

    다음으로 콘크리트 벽이 받는 중성자를 측정하였다. 세부 적으로 선형가속기 헤드를 중심으로 9개 구역으로 나누어 알파벳 소문자 a~i로 나타내었으며(Fig. 2), Tally F4를 사 용하여 9개 구역의 콘크리트 벽에 대한 중성자의 플럭스 (neutron number/㎠·sec·e)를 개별적으로 측정 하였다. 또 한 중성자의 측정은 열중성자 (Thermal, 0~1 eV), 열외중 성자 (Epithermal, 1~10 keV), 속중성자 (Fast, 10 keV~ 10 MeV)로 구분하여 측정하였다13). 이는 중성자의 경우 에 너지에 따라 핵반응 단면적이 다르기 때문이며, 특히 열중 성자에서 반응단면적이 크기 때문에 열중성자에 대한 에너 지 스펙트럼을 보기 위하여 세 가지로 나누어 스펙트럼을 측정하였다. 또한 열중성자의 경우 범위가 정해져 있는 것 이 아닌 핵반응이 잘 일어나는 구간을 말하고 있으나, 일반 적으로 기존의 이론에서는 수소의 열중성자 1 eV를 기준으 로 기술하고 있기 때문에 이를 바탕으로 설정하였다14).

    이상의 모든 실험은 결과값의 불확도를 3%로 이내로 적 용하여 실험의 신뢰도를 확보하였다.

    Ⅲ실험결과

    1광자선 스펙트럼

    Fig. 3은 전자선을 가속하여 측정한 10, 15, 20 MV 광자 선에 대한 스펙트럼 결과이며, 10 keV 간격으로 플럭스 (photon number/㎠·sec·e)를 측정하였다.

    생성된 광자의 입자수는 10, 15, 20 MV 각각 4.01E-03 개/㎠·sec·e, 1.08E-2 개/㎠·sec·e, 2.05E-2 개/㎠·sec·e였 으며 전자 1개당 광자 발생율이 약 0.4 ~ 2%로 나타났다. 또한 측정된 광자의 특성에너지는 512 keV으로 나타났으 며, 평균에너지는 각각 2.24 MeV, 3.04 MeV, 3.75 MeV로 측정되었다.

    2중성자 스펙트럼과 에너지 영역별 비율

    Fig. 4는 차폐벽 9개 구역에 대한 중성자 스펙트럼 결과 이며 열중성자의 스펙트럼을 대표적으로 나타낸 것이다.

    Table 2.은 생성된 광자에 대한 중성자의 생성 비율이며, 검출된 전체 중성자에 대한 에너지 영역별(열, 열외, 속중성 자) 결과를 Table 3-5.으로 나타냈다. Table 4

    이를 세부적으로 살펴보면, 9개 구역에서 받은 중성자 플 럭스에 대한 평균값은 10 MV 1.86E-11 개/㎠·sec·e, 15 MV 1.05E-09 개/㎠·sec·e, 20 MV 4.47E-09 개/㎠·sec·e으로 산정되었으며, 15 MV는 10 MV에 비해 약 56배 높았으며, 20 MV에 비해 약 240배 높은 결과를 보였다. 생성된 광자 에 대한 비율은 10 MV 4.63E-7%, 15 MV 9.69E-6%, 20 MV 2.18E-5%의 발생비율을 나타냈다.

    다음으로 에너지 영역별 중성자 비율은 열중성자 17.5 ~ 33.4%, 열외중성자 50.3 ~ 54.6%, 속중성자 12.2 ~ 31.7% 를 나타냈다. 특히 열중성자의 평균 플럭스 결과는 10 MV, 3.72E-12 개/㎠·sec·e, 15 MV 2.11E-10 개/㎠·sec·e, 20 MV 8.88E-10으로 나타났으며, 이것은 총 중성자 대비 10 MV 20.0%, 15 MV 20.0%, 20 MV 19.8%의 비율을 차지하였다.

    마지막으로 구역별 중성자 플럭스의 양을 보면, 실험대상 의 모든 에너지 영역에서 선형가속기 헤드와 근접된 부분이 높게 나타났다.

    Ⅳ고 찰

    본 연구는 선형가속기에서 높은 에너지의 광자들로 인하 여 발생된 중성자에 관한 연구로서 선형가속기를 차폐하고 있는 콘크리트가 받는 중성자의 양을 측정하였다.

    선형가속기의 광자선 스펙트럼과 관련하여 Geant3를 이 용한 Mesbahi15) 등의 연구와 PENELOPE-2006을 이용한 Baumgatner16) 등의 연구에서 6 MV의 평균에너지는 각각 1.67, 1.65 MeV로 산정되었다. 또한 MCNPX를 이용한 이 정욱17-19) 등의 연구에서는 24 MV 광자선의 평균에너지는 약 4.90 MeV로 보고하였다. 위의 연구와 본 연구를 비교하 면 평균에너지가 약 4~8% 높은 결과이다. 이는 광자스펙 트럼을 측정한 위치가 본 연구에서는 타깃 바로 아래 지점 에서 측정을 하였기 때문에 상대적으로 저에너지 광자의 수가 증가함으로써 평균에너지가 낮게 측정된 것으로 사료 된다.

    고 에너지의 광자로 인하여 생성된 중성자의 발생 비율을 보면 10 MV 대비 15 MV에서는 약 20배, 15 MV 대비 20 MV에서는 약 2.25배로 중성자의 발생비율이 증가했다. 이 것은 측정된 광자수에 따라 비례적으로 중성자 발생확률이 증가하는 것이 아닌 다른 확률적인 변수가 있다는 점을 의 미하는 것이며, 10 MV를 초과할 경우 중성자 발생확률이 급격히 높아지는 현상을 확인할 수 있었다. 또한 이와 같은 결과는 기존의 박철수 등의 연구8)에서 광자 에너지에 커질 수록 중성자의 발생이 증가하지만 비율적이지 않다는 것과 일정한 광자를 조사하더라도 중성자의 발생은 일정하게 나 타나지 않는다는 결과와 유사한 결과를 얻었다.

    마지막으로 차폐벽이 받은 중성자 양은 치료기 헤드를 중 심으로 거리가 이격될수록 감소하는 경향을 보였다. 하지만 이에 반해 열중성자의 비율은 높게 나타났다. 이것은 발생 된 중성자가 특정한 방향성을 가지고 방사하는 것이 아닌 전체적으로 균등하게 방사하는 것으로 판단되며, 열외중성 자와 속중성자의 경우에는 치료실 차폐벽과 산란으로 인하 여 속도를 잃고 열중성자로 변하게 됨으로써, 중성자가 발 생된 선형가속기 헤드에서 거리가 이격될수록 열중성자의 비율이 높아진 것으로 사료된다.

    Ⅴ결 론

    본 연구는 의료용 선형가속기의 유지, 보수, 해체에 대한 연구의 선행연구로서 치료실 내에서 발생하는 중성자에 대 한 특성을 평가하고자 하였다. 특히 본 논문에서는 차폐벽 이 받는 중성자에 대한 영향을 구역을 설정하여 분석하였 다. 그 결과 높은 에너지에서 콘크리트 차폐벽이 받는 중 성자의 양은 제동복사로 인해 생성된 광자 대비 최대 5.0E-5% 정도로 미미하게 나타났다. 그러나 본 연구에서 사용한 MCNPX 프로그램의 경우 초당 전자 1개에 대하여 측정된 양이므로 실제 치료 시 가속하는 전자의 개수를 감 안한다면, 실질적인 양은 무시할 수 없을 것으로 판단된 다. 또한 치료실의 가동시간과 해체 시점까지의 운영기간 을 역산하여 고려한다면, 수십 년 후 방사화가 진행될 가 능성은 충분이 높을 것으로 생각된다. 또한 본 실험을 통 해 발생한 중성자의 특성은 특정한 방향성 없이 방사하며, 거리에 따라 중성자의 양과 열중성자의 비율이 변화됨을 알 수 있었다.

    Figure

    JRST-39-89_F1.gif

    The linear accelerator head in MCNPX

    JRST-39-89_F2.gif

    A schematic design for shielding wall in MCNPX

    JRST-39-89_F3.gif

    The spectra of thermal neutron measured in the simulation

    JRST-39-89_F4.gif

    The spectra of thermal neutron measured in the simulation

    Table

    The accessories for linear accelerator head in MCNPX

    Neutron Production rate of the photon

    Ratio of the neutron in shielding wall (unit : #/㎠·sec·e)

    Ratio of the neutron in shielding wall (unit : #/㎠·sec·e)

    Ratio of the neutron in shielding wall (unit : #/㎠·sec·e)

    Reference

    1. Dept. of Radiological Science Graduate School Catholic University of PusanEun Tae Park (2015) Study on the Characteristics of Photoneutron Produced in Medical Linear Accelerators ,
    2. Dong Wook Kim , Myong Geun Yoon , Sung Yong Park , Se Byeong Lee (2008) Comparison of Helical TomoTherapy with Linear Accelerator Base Intensity-modulated Radiotherapy for Head & Neck Cases , Korean Society of Medical Physics, Vol.19 (2) ; pp.89-94
    3. NCRP Report No 79 (1984) Neutron Contamination From Medical Electron Accelerators,
    4. Huang WL , Li QF , Lin YZ , Su Q , Luo YS (2006) Measurements of phononeutrons produced by a 15 MeV electron linac for radiography applications , Nucl. Instr. and Meth. in Phy. Res. B, Vol.251 ; pp.361-366
    5. NCRP Report No 151 : Radiation (2005) Structural Shielding Design and Evaluation for Megavoltage Xand Gamma-Ray Radiotherapy Facilities,
    6. Chao JH , Liu WS , Chen CY (2007) Estimation of Argon-41 concentrations in the vicinity of a high energy medical accelerator , Radiation Measurements, Vol.42 ; pp.1538-1544
    7. Sang Koo Kang , Sung Hwan Ahn , Chong Yeal Kim (2011) A Study on Photon Dose Calculation in 6 MV Linear Accelerator Based on Monte Carlo Method , Korean Society of Radiological Science, Vol.34 (1) ; pp.43-50
    8. Cheol Soo Park , Cheong Hwan Lim , Hong Ryang Jung , Seong Soo Shin (2008) A Study on the Neutron Dose Distribution in Case of 10 MV X-rays Radiotherapy , Korean Society of Radiological Science, Vol.34 (4) ; pp.415-427
    9. Dept. of Radiological Science The Graduate School of Health promotion Hanseo UniversityOh Nam Yang (2013) Study on the Photoneutron Produced in Medical Linear Accelerators : Comparison of 3DCRT and IMRT ,
    10. IAEAIndustrial and Research Facilities (1999) Safety Standard Series No. WS-G-2.2, Decommissioning of Medical,
    11. IAEA (2006) Safety Report Series No.47 Radiation Protection in the Design of Radoitherapy Facilities,
    12. Korea Atomic Energy Research Institute (2004) Shielding Technology for High Energy Radiation Production Facility,
    13. Los Alamos National Laboratory (2005) User's manual version 2.5.0, LANL,
    14. Herman Cember , Thomas E Johnson (2006) Introduction to Health Physics, Mv Graw Hill,
    15. Baumgartner A , Steurer A , Maringer F (2009) Simulation of photon energy spectra from Varian 2100C and 2300C/D Linacs: Simplified estimates with PENELOPE Monte Carlo models , Applied Radiation and Isotope, Vol.67 (11) ; pp.2007-2012
    16. Mesbahi A , Fix M , Allahverdi M (2005) Monte Carlo calculation of Varian 2300C/D Linac photon beam Characteristics: a comparison between MCNP4C , GEANT3 and measurements Applied Radiation and Isotopes, Vol.62 (3) ; pp.469-477
    17. Jeong Ok Lee , Dong Hyeok Jeong , Kang Jeong Ku (2005) Neutron Generation from a 24 MV Medical Linac , Korean Society of Medical Physics, Vol.16 (2) ; pp.97-103
    18. Dept. of Radiation Science and Technology The Graduate School Chonbuk National UniversityKang Sang Koo (2012) A study on 6 MV Photon Dose Distributions from Medical Linear Accelerators in phantom based on CT data using GEANT4 ,
    19. Lee Jeong Ok , Jeong Dong Hyeok (2006) Calculation of Energy Spectra for 6 MeV Electron Beam of LINAC Using MCNPX , Korean Society of Medical Physics, Vol.17 (4) ; pp.224-231