125I粒子源剂量计算参数模拟研究

Monte Carlo Simulation of Dose Calculation Parameter of Radioactive 125I Brachytherapy Source

  • 摘要: 国产125I粒子源的银棒末端结构为直角型,与典型的6711型粒子源结构略有不同,结构不同会对剂量计算参数产生一定影响。本文针对国产粒子源结构,利用蒙特卡罗方法计算美国医学物理学家协会(AAPM)在TG43-U1报告中推荐的剂量计算参数,并分析研究银棒末端结构对剂量计算参数的影响。模拟得到国产125I粒子源剂量率常数为0.955 cGy·h-1·U-1(空气比释动能强度基于点探测器计算得到),与TG43-U1推荐值较接近,两者仅相差1.03%,更加精细地计算了在源中垂线0.05~10 cm(1 cm间隔)范围内的径向剂量函数,拟合得到较好的经验公式,得到在0°~90°(5°间隔)、距源中心0.25~7 cm(2 cm间隔)范围内的二维各向异性函数,通过对比分析得到银棒末端为直角型结构时的二维各向异性函数在r=0.25 cm处会引起驼峰区。

     

    Abstract: The end of silver rod of the domestic 125I brachytherapy source is right angle type, which is slightly different from the typical model of 6711 125I brachytherapy source. And it can have influence on the dose calculation parameters. Based on the structure of domestic 125I brachytherapy source, dose calculation parameters which are recommended by AAPM TG43-U1 were calculated by Monte Carlo method. The influence of the end of silver rod on the dose calculation parameters was studied. The simulation result of dose rate constant is 0.955 cGy·h-1·U-1 when the air kerma strength was calculated by the point detector, and the difference with the result of the TG43-U1 is within 1.03%. The radial dose function g(r) in the range of 0.05-10 cm at the transverse axis was calculated precisely. Then empiric equation was acquired by curve fitting. 2D anisotropy function F(r,θ) was calculated in 0°-90° and 0.25-7 cm. The source of the right angle structure of the end of the silver rod would cause a hump area of 2D anisotropy function when r equals 0.25 cm.

     

/

返回文章
返回