• 한국의학물리학회지:의학물리
• /
• 제16권3호
• /
• pp.138-147
• /
• 2005

몬테칼로 선량계산 시 적절한 정확도를 얻기 위해서는 계산입자수를 늘려야 하고 그로 인해 계산 시간이 오래 걸리게되므로 일상적 치료계획의 선량계산 방법으로 이용되지 못했다. 본 연구에서는 몬테칼로 모의실험 시 계산입자 수를 줄여서 선량계산을 수행한 후 잡음 감소 필터를 적용하여 선량계산 결과를 개선하고자 하였다. 이를 위해 국소 최소자승 잡음 감소 필터를 제작하고 광자선 및 전자선 균질/비균질 팬텀 내 선량계산 결과에 대하여 적용하여 선택적 여과를 수행하였으며 그 유효성을 등선량 곡선 비교 및 감마시험을 통하여 검증하였다. 연구 결과 통계적 불확실도를 2$\%$ 이내로 유지하기 위해 필요한 계산입자수의 10$\%$ 이하의 계산입자 수를 이용하여 몬테칼로 선량계산 뒤 후처리한 결과가 기준계산 입자수를 이용하여 얻은 몬테칼로 선량계산 결과와 유사해질 수 있음을 확인하였다.

• Journal of Radiation Protection and Research
• /
• 제39권1호
• /
• pp.30-37
• /
• 2014

최근 표준한국인 성인남녀 복셀 팬텀 한 쌍인 High-Definition Reference Korean-Man (HDRK-Man)과 High-Definition Reference Korean-Woman (HDRK-Woman)을 개발한 바 있다. HDRK-Man과 HDRK-Woman은 각각 한국인 성인 남성과 여성을 대표하므로 이들을 이용하면 한국인 고유의 유효선량을 계산할 수 있다. 이러한 한국인 고유의 유효선량을 계산하기 위해서는 표준한국인 팬텀을 이용한 몬테칼로 선량계산을 수행하여야 한다. 하지만 이러한 몬테칼로 선량계산은 전산인체모델에 대한 이해와 고도의 몬테칼로 전산모사 기술이 요구되기 때문에 몬테칼로 전산모사 경험이 적은 연구자들에게는 접근성이 용이하지 않다. 최근 UFPE (Federal University of Pernambuco) 연구팀에서는 몬테칼로 선량계산 온라인 시스템인 CALDOSE_X를 운영하고 있다. CALDOSE_X는 사용자가 온라인상에서 직접 입력한 선원조건에 따라 곧바로 몬테칼로 선량계산을 수행하고 계산된 결과를 제공해주는 온라인 서비스이다. 따라서 CALDOSE_X를 이용하면 몬테칼로 전산모사 경험 없이도 누구나 쉽게 원하는 선원조건에 대한 선량계산을 수행할 수 있다. 하지만 CALDOSE_X에서 제공하는 선량값은 서양인 전산팬텀을 이용하기 때문에, 이를 한국인 종사자에 대한 방사선 방호목적으로 사용하기에는 여전히 한계가 있다. 이에 본 연구에서는 온라인상에서 사용자가 입력한 선원조건에 따라 실시간으로 표준 한국인에 대한 유효선량을 계산할 수 있는 웹 기반 표준한국인 선량평가 시스템을 구축하였다.

• 한국의학물리학회지:의학물리
• /
• 제14권4호
• /
• pp.240-248
• /
• 2003

몬테칼로 계산은 다른 어떤 알고리즘보다 정확한 선량 계산 결과를 주지만 계산 시간이 오래 걸리는 단점이 있다. 본 연구에서는 Varian 600 C/D 선형가속기로부터지 6 MV 광자선에 대해 몬테칼로 계산을 사용하여 얻은 선량 분포가 측정에 의해 얻은 선량 분포와 2％ 이내에서 서로 잘 일치하도록 하며 분산 감소 기법을 사용하여 계산 시간 단축 정도를 평가하였다. 그리고 연산 능력을 높여 계산 시간 단축 정도를 평가하여 분산 감소 기법을 사용한 경우와 연산 능력을 높인 경우 간에 계산 시간 단축 정도를 비교하였다. 몬테칼로 계산 코드로는 빔 모사를 위해 BEAMnrc 코드, 선량 계산을 위해 DOSXYZnrc 코트를 각각 사용하였는데 분산 감소 기법은 이 코드들에서 지원하는 방법들을 사용하였고 연산 능력을 높이는 방법으로는 컴퓨터 클러스터를 이용한 병렬 처리를 사용하였다. 비교 결과, 분산 감소 기법을 사용하여 계산 시간을 최대 1/25 이상 단축시킬 수 있었고 9대의 컴퓨터를 이용한 병렬 처리 결과 계산 시간을 1/9로 단축시킬 수 있었다. 계산 곁과의 정확성을 만족할 만한 수준으로 유지할 수 있다면 분산감소 기법을 포함한 간략화된 물리의 적용은 현 시점에서 몬테칼로 선량 계산 시간을 획기적으로 단축시킬 대안이 될 수 있다.

• Communications for Statistical Applications and Methods
• /
• 제3권2호
• /
• pp.291-298
• /
• 1996

통계학분야 그리고 또 다른 많은 분야에서 수치적 계산을 다루는 문제가 자주 발생한다. 적당한 컴퓨터 컴퓨터 시간안에 상당한 정도의 정확성을 줄 수 있고 또한 보다 광범위하게 사용 가능한 유용한 알고리즘의 필요성을 느낀다. 이러한 문제에 가능한 하나의 몬테칼로 알고리즘인 주표본 알고리즘을 소개하였다. 그리고 특히 본 눈문에서는 축차확률비검정의 오차확률을 계산하는 곳에 주표본 알고리즘을 적용하고 결과를 비교분석하였다.

• 한국의학물리학회지:의학물리
• /
• 제10권1호
• /
• pp.55-62
• /
• 1999

방사선치료에 이용되는 X-선의 특성을 고찰하기 위하여, 몬테칼로 방법을 이용하여 텅스텐 표적에 입사한 6, 10, 15 MeV 의 전자선에 의해 발생되는 X-선의 에너지분포와 평균에너지를 계산하였다. 빔의 반경의 함수로서 계산된 광자의 평균에너지는 6, 10, 15 MV에 대하여 각각, 1.4-1.6, 2.1-2.5, 2.8-3.3 MeV 범위로서 반경에 크게 의존하지 않고 거의 일정하였다. 표적과의 수직거리 100 cm에서 구해진 6, 10, 15 MV X-선의 에너지분포를 이용하여, 몬테칼로 계산으로 깊이선량율을 계산하였다. 이 중 10 MV에 대한 것을 측정값과 비교하였다. 계산된 10 MV X-선의 깊이선량율은 표면영역을 제외하고 깊이의 증가에 따라 측정값보다 낮게 나타났다. 그 이유는 실제 X-선의 에너지분포에는 편평화여과기에 의한 빔경화효과가 포함되어 있는 반면, 본 연구에서 수행한 몬테칼로 계산결과에는 이 효과가 포함되어 있지 않기 때문으로 생각된다.

• 대한방사선기술학회지:방사선기술과학
• /
• 제31권1호
• /
• pp.89-95
• /
• 2008

국제원자력기구의 TRS-398 측정 프로토콜을 임상에 적용하기 위해서는 사용하는 빔과 전리함에 대한 선질 보정인자가 필요하다. 본 연구에서는 몬테칼로 계산코드(DOSRZnrc/EGSnrc)를 사용하여 상용의 평판형전리함에 대한 고에너지 전자선($4{\sim}20\;MeV$)에서의 선질보정인자를 계산하였다. 계산결과를 프로토콜에서 제시하는 값과 비교한 결과 $5{\sim}20\;MeV$에서 약 1% 이내로 일치하였으며 4 MeV의 경우에는 약 1.9% 차이를 보였다. 본 연구 방법은 선질보정인자를 독립적으로 결정하는 방법의 하나로서 프로토콜에서 주어진 값들의 확인이 필요하거나 또는 새로운 모델의 전리함을 사용하는 경우에 응용될 수 있다.

• 한국의학물리학회지:의학물리
• /
• 제21권3호
• /
• pp.253-260
• /
• 2010

대다수의 근접치료용 방사선치료계획장치는 AAPM TG-43의 계산식에 기반을 둔 선량계산 알고리듬을 적용하고 있으나 이는 조직의 비균질성을 적절히 고려하지 못한다. 본 연구에서는 몬테칼로 방법을 이용하여 강내고선량근접치료계획을 검증하는 체계를 구축하고자 하였으며, 특히 환자의 CT 영상을 이용하여 물질정보로 변환한 후 직접 몬테칼로 계산을 수행하는 방법의 타당성에 초점을 맞추었다. 판형 팬텀 및 자궁경부암 환자의 CT 영상을 Plato (Nucletron, Netherlands) 치료계획장치를 이용하여 근접치료계획을 수행한 후 여기서 얻어진 인자들을 이용하여 EGSnrc 기반의 DOSXYZnrc 코드로 몬테칼로 계산을 수행하였으며, EBT 필름측정 결과와 비교하였다. DOSXYZnrc 코드의 선원 모델링 특성 상 후장전 장치의 $^{192}Ir$ 선원들을 직육면체 형태로 근사화하여 모델링하였으며 계산 시 체적소의 크기는 $2{\times}2{\times}2\;mm^3$로 하였다. 균질 매질 내에서는 TG-43 기반의 선량계산 결과와 몬테칼로 선량계산 결과가 잘 일치함을 확인할 수 있었으나 고밀도 물질이 포함된 비균질 매질 내에서는 오차가 커졌다. 환자의 경우 A점 및 B점의 오차는 3% 이내, 평균선량 오차는 5% 정도였다. 그러나 기존 선량계산 알고리듬의 경우 고밀도 물질의 영향을 적절히 고려하지 못하여 표적의 선량을 과대평가하여 실제로는 더 적은 선량이 들어갈 우려가 있다. 본 연구에서 제안된 선량계산 검증체계는 타당하며 선량 계산 결과도 실제와 잘 일치함을 확인할 수 있었다. 또한 기존의 선량계산 알고리듬으로 계산된 치료계획결과를 확인할 경우에는 주의가 필요하며, 몬테칼로 방법과 같은 독립적인 검증 시스템이 유용할 것이다.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1997년도 추계학술발표회논문집(1)
• /
• pp.115-120
• /
• 1997

영광3호기 방사선관리구역에 대한 중성자선량률을 정확히 평가하기 위하여 MCNP4A 전산코드를 이용, 방사선관리구역에서의 중성자 스펙트럼 계산을 수행하였다. 영광3호기에 대한 보다 정확하고 정밀한 3차원 몬테칼로 모델을 구축하기 위하여 핵연료집합체 구성요소 및 원자로심을 둘러싸고 있는 baffle, barrel,압력용기 등을 정확하게 묘사하였으며, 특히 방사선관리구역 주위의 구조물에 대해서도 3자원 MCNP 모델을 구축함으로써 원자로심부터 방사선관리구역까지 완전한 몬테칼로 모사(full-scope Monte Carlo simulation)를 이용한 계산을 수행하였다. 계산결과는 에너지 구간에 따른 중성자속 스펙트럼으로 나타내었으며 이 결과를 바탕으로 중성자속에 대한 선량률 환산인자를 고려하여 중성자선량률을 계산할 수 있다.

• 한국의학물리학회지:의학물리
• /
• 제14권4호
• /
• pp.234-239
• /
• 2003

본 연구에서는 몬테칼로 계산을 이용하여 외부 가로 자기장에 의한 깊이선량율(PDD)의 변화를 고찰하였다. 몬테칼로 계산은 자기장에서 전자의 편향을 고려하도록 수정한 EGS4 몬테칼로 코드를 사용하였다. 자기장에서 깊이선량율의 변화를 기술하기 위하여, 선량증가(DI； dose improvement)와 선량감소(DR； dose reduction)를 정의하였다. 10 MV 광자선에 대하여 1-5 T 자기장 범위에서 계산한 결과, 자기장의 세기에 따라 DI와 DR은 거의 선형으로 각각 증가, 감소하였다. 자기장 3 T의 경우에 조사면 10${\times}$10 $\textrm{cm}^2$와 자기장 인가깊이 5-10 cm에서 DI는 1.56 (56% 증가), DR은 0.68 (32% 감소)로 나타났다. 깊이선량율 변화의 원리는 로렌츠 법칙과 전자평형 개념으로부터 설명하였으며, 이러한 특성을 이용하여 방사선치료의 최적화를 달성할 수 있음을 제안하였다.

• 한국통계학회:학술대회논문집
• /
• 한국통계학회 2000년도 추계학술발표회 논문집
• /
• pp.135-140
• /
• 2000

베이지안 신경망 모형(Bayesian Neural Networks Models)에서 주어진 입력값(input)은 블랙 박스(Black-Box)와 같은 신경망 구조의 각 층(layer)을 거쳐서 출력값(output)으로 계산된다. 새로운 입력 데이터에 대한 예측값은 사후분포(posterior distribution)의 기대값(mean)에 의해 계산된다. 주어진 사전분포(prior distribution)와 학습데이터에 의한 가능도함수(likelihood functions)를 통해 계산되어진 사후분포는 매우 복잡한 구조를 갖게 됨으로서 기대값의 적분계산에 대한 어려움이 발생한다. 이때 확률적 추정에 의한 근사 방법인 몬테칼로 적분을 이용한다. 이러한 방법으로서 Hybrid Monte Carlo 알고리즘은 우수한 결과를 제공하여준다(Neal 1996). 본 논문에서는 Hybrid Monte Carlo 알고리즘과 기존에 많이 사용되고 있는 Gibbs sampling, Metropolis algorithm, 그리고 Slice Sampling등의 몬테칼로 방법들을 비교한다.