방사능 토양 감마선 선량 변환 계수 재검토와 시뮬레이션 반경 효과

초록

이 서한은 Gasser 등(2014)의 토양 내 ⁴⁰K, ²³⁶Ra 및 ²³²Th 계열에 대한 감마선‑공기 케르마 변환 계수가 기존 값보다 약 20 % 낮게 보고된 원인을 분석한다. 저자들은 시뮬레이션에서 사용된 원통형 소스의 반경이 제한적이어서 반평면(반무한) 조건을 충분히 재현하지 못했음을 지적하고, 반경을 700 m 이상으로 확대한 PHITS 시뮬레이션 결과를 제시한다. 확대된 반경에서는 변환 계수가 기존 문헌값과 일치함을 보여주며, Gasser 등 결과의 차이는 주로 유한 크기 효과임을 주장한다.

상세 분석

본 서한은 Gasser et al.이 MCNPX를 이용해 1 m 깊이·반경 35 m까지의 원통형 토양 소스를 모델링하고, 지표면 위 1 m에서 공기 케르마를 계산한 방법을 비판한다. ICRU‑53에서 정의한 무한 반평면(half‑space) 조건은 실제 토양이 무한히 넓다고 가정할 때만 적용 가능하므로, 시뮬레이션 반경이 충분히 커야 수렴값을 얻을 수 있다. 저자들은 PHITS 코드를 사용해 반경을 1 m에서 1000 m까지 확대한 일련의 시뮬레이션을 수행했으며, 변환 계수는 반경이 35 m를 넘어 지속적으로 증가한다는 점을 확인했다. 특히 700 m 이상에서야 3자리 유효숫자 수준으로 수렴했으며, 최종 수렴값은 ⁴⁰K = 0.042 nGy h⁻¹ Bq⁻¹ kg, ²³⁸U 계열 = 0.444 nGy h⁻¹ Bq⁻¹ kg, ²³²Th 계열 = 0.592 nGy h⁻¹ Bq⁻¹ kg이다. 이는 Saito & Jacob(1995)과 Askri(2015)의 결과와 일치하고, Gasser 등(2014)의 0.036, 0.357, 0.482와는 현저히 차이 난다.

저자들은 차이의 주요 원인을 “유한 크기 효과”로 규정하고, 남은 소수점 이하 차이는 방사성 붕괴 데이터의 최신화, 에너지 분해능, 통계적 정확도, 그리고 광자 플루언스에서 케르마 변환 과정의 보간 방법 차이 등으로 설명한다. 또한, 시뮬레이션 효율성을 높이기 위해 소스와 검출기 볼륨을 변환하거나, 평면 검출기와 주기적 경계 조건을 적용하는 방안을 제시한다.

결론적으로, 변환 계수를 정확히 산출하려면 반경이 수백 미터에 달하는 가상 토양 소스를 사용해야 하며, 기존 35 m 반경 모델은 충분히 큰 반평면을 대체하지 못한다는 점을 강조한다. 이는 환경 방사선 평가, 특히 자연 방사능에 기반한 외부 피폭 계산에 중요한 교훈을 제공한다.