.jpg)
1.서론
감마 스펙트로메트리는 환경시료 중 천연방사능핵종 혹 은 인공핵종을 분석하기 위하여 자주 이용되며 그 중 정밀분 석을 위해 가장 많이 사용되는 장비는 HPGe 검출기 (High Purity Germanium Detector, 고순도 게르마늄 검출기)이다. 환경시료에 대한 정확한 방사능 농도 분석을 위해서는 우선 기준물질을 이용하여 주어진 기하학적 구조를 바탕으로 효 율 교정이 이루어진다. 각 시료에 대한 정확한 분석을 위해서 는 교정용 표준시료와 실제시료가 물리적 및 화학적으로 동 일해야 하나 이는 현실적으로 어려우며 따라서 이에 대한 보 정이 필요하다. 즉 환경시료는 밀도 및 화학적 구성 등 다양 한 특성을 가지며 이는 분석결과에 영향을 미치게 된다. 특히 교정용 표준시료의 밀도는 한 가지 값으로 고정되어 있으나 0.3 g/cm3에서 1.5 g/cm3까지 다양한 범위의 밀도를 가진 환 경시료에 적용하기 위해서는 밀도에 대한 보정이 필요하다. 저에너지 감마선의 경우 밀도 차이에 의한 self-attenuation 에 의해 분석 결과가 크게 달라질 수 있으므로 밀도 보정은 환경시료를 정확하게 분석하는데 있어 중요한 역할을 하게 된다[1]. 밀도에 의한 효율 보정은 실험적 혹은 이론적으로 접 근 가능하지만 실험적으로 각 밀도에 대한 보정인자를 구하 기 위해서는 다양한 밀도를 가지는 시료들을 준비해야 하는 번거로움이 있다. 이를 대체할 수 있는 방법으로는 주로 몬테 카를로 시뮬레이션을 이용한 코드 사용이 가능하며 이에 대 한 연구는 계속해서 진행되어 왔다[2-4]. 그러나 코드를 사용 한 밀도보정인자 도출의 경우 각 실험실이 가지고 있는 검출 기의 구조 및 구성 물질에 대한 정확한 모사가 필요하며 이 에 대한 검증 또한 필요하다. 즉 계측 효율에 가장 크게 영향 을 미칠 수 있는 크리스털의 지름과 높이를 포함한 형태 특 성 및 크기에 대한 고려가 필요하므로 검출기가 달라질 경우 에는 새로운 모사가 필요할 수 있다. 따라서 본 연구에서는 몬테카를로 코드인 MCNPX 코드를 사용하여 크리스털과 같 은 주요 특성이 다른 세 대의 HPGe 검출기를 모사하고 실제 실험을 통하여 모델링을 검증한 이후 밀도보정인자 F (selfabsorption correction factor)를 도출함으로써 검출기의 특 성에 따른 밀도보정인자의 의존도를 평가하였다.
2.모델링 및 검증
본 연구에서는 Table 1과 같이 세 가지 p-type 검출기에 대한 모델링 및 시뮬레이션을 통하여 효율을 구하였다. 세 검 출기 모두 효율은 30%으로 동일하나 크리스털과 코어 등 주 요 특성이 모두 다르다. 따라서 이를 시뮬레이션하기 위해서 는 우선 제작사에서 제공된 자료를 바탕으로 모델링을 하였 으며 이후 실험값과 비교하여 검증하였다.
각 검출기를 모델링 후 MCNPX 코드를 이용하여 계산한 효율 값과 밀도 1 g/cm3의 마리넬리 1리터 교정용 표준시료 를 이용하여 계측한 효율 값을 비교하였으며 Detector I의 구 조를 모델링한 결과는 Fig. 1과 같다. MCNPX 코드를 이용하 여 계측기의 효율을 비교하기 위해서 F8 tally를 사용하였으며 history number는 108개로 설정하였다. 또한 시뮬레이션 결 과 statistical error는 5% 이하임을 확인하였다. 교정용 표준 시료에는 Am-241 (59.54 keV), Cd-109 (88.03 keV), Co-57 (122.06, 136.47 keV), Ce-139 (165.86 keV), Sn-113 (391.70 keV), Cs-137 (661.66 keV), Y-88 (898.04, 1836.04 keV), Co- 60 (1173.23, 1332.49 keV)등의 핵종이 포함되어 있으며 이 에 해당하는 에너지에 대해 MCNPX를 이용하여 계산한 효율 과 실제 계측한 효율의 비 Esimulation/Eexperiment를 구하였다.
Fig. 2에서 보듯이 높은 에너지 영역에서는 그 비가 1로 수렴하며 계산값과 측정값이 유사함을 확인하였으나 저에너 지 영역에서는 그 차이가 커짐을 확인할 수 있다. 계측 효율 은 게르마늄 크리스털을 둘러싸고 있는 dead layer 영역의 두께에 영향을 많이 받게 된다. 이는 게르마늄에서 광자의 감쇄계수가 크기 때문이며 따라서 낮은 에너지에서는 dead layer의 두께에 매우 민감하게 반응하므로 이에 대한 조정 이 필요하다[5]. 이에 대한 보정이 없을 경우 저에너지 영역 에서의 불확도는 커질 수 있으며 이는 밀도에 대한 시뮬레이 션 결과값에 대한 신뢰도 확보가 어려울 수 있으므로 계측값 과 비교하여 dead layer의 두께를 조정함으로써 가장 유사한 값을 보이는 두께를 구할 필요가 있다. 따라서 dead layer의 두께를 0.1 mm씩 증가시키며 비교 과정을 통해 계측기 I, II 의 경우 dead layer를 0.7 mm에서 1.1 mm로 조정하였으며, 가장 차이를 많이 보였던 계측기 III의 경우에는 1.4 mm로 조정하였다. Dead layer의 두께를 조정한 후, 시뮬레이션 값 과 실험값의 비는 Fig. 3과 같이 에너지 전 영역에서 1에 가 까운 분포를 보임을 확인하였다. 따라서 이러한 검증을 토대 로 각 검출기에 대한 모델링을 확정하였으며, 이를 이용하여 여러 밀도 값에 대한 효율을 계산하였다.
3.결과
검증된 모델을 바탕으로 마리넬리 1리터 비이커에 담긴 샘플의 밀도 변화에 따른 효율을 MCNPX 코드를 이용하여 구하였으며 검출기 I에 대한 결과는 Table 2와 같다. 즉 교정 용 표준시료의 밀도인 1 g/cm3와 비교하였을 때 낮은 밀도에 서는 저에너지 영역에서 효율의 차이가 있음을 확인할 수 있 다. 이를 이용하여 밀도보정인자 F는 밀도 1 g/cm3에 대한 효 율에 대한 비로 나타낼 수 있으며 그 결과는 Table 3과 같이 0.9 - 1.17 값을 가진다. 즉 밀도가 0.3 g/cm3인 샘플에 대해 서 Am-241 (59.54 keV)에 대한 밀도보정인자는 1.17이 되는 것이다. 이와 같은 방법으로 검출기 II, III에 대해서도 밀도 변화에 따른 에너지별 효율과 밀도보정인자 F를 구한 결과는 Table 4, 5와 같다. 밀도보정인자를 비교하였을 때 검출기를 I를 기준으로 각 밀도에 대해 전 에너지 범위에서 그 차이가 1%이하임을 확인할 수 있다. 즉 그 값이 거의 유사함을 확인 할 수 있으며, Fig. 4, 5에 밀도 0.3과 1.5 g/cm3에 대한 밀도 보정인자 F를 비교하였다. Table 5
4.결론
본 연구에서는 크리스털 및 코어와 같은 효율 측정에 중 요한 영향을 미치는 주요 특성이 다른 p-type 고순도 게르마 늄 검출기에 대해 모델링하고 실제 표준시료를 이용하여 실 험값과 시뮬레이션 값을 비교함으로써 모델링을 검증하였 다. 이는 dead layer에 의해 저에너지 영역에서 생길 수 있는 실험값과 계산값의 차이를 줄이기 위함이며 이를 통해 모델 링을 확정하였다. 이 결과를 바탕으로 전 에너지에 걸쳐 0.3, 0.6, 0.9, 1.2, 1.5 g/cm3 밀도에 대한 밀도보정인자를 도출 하였으며, 그 결과 그 차이가 크지 않아 밀도보정인자가 크 리스털과 같은 검출기의 주요 특성에 독립적임을 확인하였 다. 다만 본 연구에서는 일부 p-type의 HPGe 검출기만 다루 었기 때문에 다른 다양한 p-type 및 n-type HPGe 검출기에 대한 추가 검증 및 연구가 필요하다고 판단된다.
