(57)【要約】【課題】放射線照射室の床下などの室外空間の線量を精度良く計算することのできる線量計算方法、線量計算プログラムを提供する。【解決手段】放射線を照射する線源Ｓを有する放射線照射室の周囲に設けられた放射線遮蔽用の遮蔽体ａの外側の室外空間ｂの線量を計算する方法であって、線源Ｓから照射されて遮蔽体ａを透過する放射線のうち少なくとも一部が遮蔽体ａ内で散乱して、室外空間ｂ内に設定した所定の評価点Ｅに到達すると仮定して、線源Ｓから照射された放射線が遮蔽体ａ内で散乱して評価点Ｅに到達する寄与の値と、線源Ｓから照射された放射線が遮蔽体ａを透過した後に、遮蔽体ａの外側の室外空間ｂで多重反射して評価点Ｅに到達する寄与の値とを合算して評価点Ｅの線量を求めるようにする。【選択図】図１
2020110367.tif 000002
【特許請求の範囲】
【請求項１】
放射線を照射する線源を有する放射線照射室の周囲に設けられた放射線遮蔽用の遮蔽体の外側の室外空間の線量を計算する方法であって、
線源から照射されて遮蔽体を透過する放射線のうち少なくとも一部が遮蔽体内で散乱して、室外空間内に設定した所定の評価点に到達すると仮定して、
線源から照射された放射線が遮蔽体内で散乱して評価点に到達する寄与の値と、線源から照射された放射線が遮蔽体を透過した後に、遮蔽体の外側の室外空間で多重反射して評価点に到達する寄与の値とを合算して評価点の線量を求めることを特徴とする線量計算方法。
【請求項２】
室外空間での放射線の多重反射の反射率を、エネルギーの減衰によらず一定値と仮定して、室外空間で多重反射して評価点に到達する寄与の値を求めることを特徴とする請求項１に記載の線量計算方法。
【請求項３】
線源から照射された放射線のエネルギーと、遮蔽体の材質と、室外空間を挟んで遮蔽体の表面に対向する物体の材質と、線源から照射された放射線が遮蔽体を透過した後に、物体の表面で反射して評価点に向かう放射線による線量率と、物体の表面で反射した後、遮蔽体の表面で反射して評価点に向かう放射線による線量率とに基づいて、室外空間での放射線の多重反射の反射率を求め、求めた反射率を用いて、室外空間で多重反射して評価点に到達する寄与の値を求めることを特徴とする請求項１または請求項２に記載の線量計算方法。
【請求項４】
請求項１から請求項３のうちいずれか１項に記載の線量計算方法をコンピュータに実行させることを特徴とする線量計算プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば放射線施設の線量計算方法、線量計算プログラムに関するものである。
【背景技術】
【０００２】
従来、電子リニアック（直線加速器）等の放射線発生装置が設置された放射線照射室では、管理区域外への放射線の漏洩を抑え、管理区域境界の実効線量を法令で定める値以下にするため、コンクリートや鉄を使った厚い遮蔽壁が設けられる。
【０００３】
国内では、放射線施設の管理区域境界の線量評価は、非特許文献１に記載の簡易計算式で計算するのが一般的である。また、より複雑な遮蔽の場合は、対象を３次元モデル化した上で、モンテカルロ法による放射線輸送計算により、実効線量を計算することができる。
【０００４】
上記の他に、放射線施設の線量評価の簡易計算式は非特許文献２〜４にまとめられている。特に、非特許文献２には、Ｘ線の利用線錐が照射室の迷路に向けられ、Ｘ線が迷路壁を透過したのちに、迷路の奥側の壁で１回散乱し、照射室出入り口に到達するＸ線の実効線量の簡易計算式について、図５を用いて以下のように記載されている。
【０００５】
【数１】
2020110367.tif 000003
【０００６】
ただし、ＤwTは照射室入り口Ｅでの線量、Ｗは照射線量、Ｕmは方向利用率、Ｂprは迷路壁Ｍの透過率、αpは壁Ｐにおける反射率、Ａpは壁Ｐにおける最大照射範囲、ｄpは線源Ｓから壁Ｐの中心までの距離、ｄ’’は壁Ｐの中心から照射室入り口Ｅまでの距離である。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】「放射線施設のしゃへい計算実務マニュアル２０１５」、公益財団法人原子力安全技術センター、放射線障害防止法出版物編集委員会発行、平成２７年３月
【非特許文献２】International  Atomic  Energy  Agency,  “Radiation  protection  in  the  design  of  radiotherapy  facilities.  Safety  Reports  Series  No.  47”,  Vienna  (Austria):  International  Atomic  Energy  Agency,  2006,  ISBN-10:  92-0-100505-9
【非特許文献３】National  Council  on  Radiation  Protection  and  Measurements,  “Structural  shielding  design  and  evaluation  for  megavoltage  x-  and  gamma-ray  radiotherapy  facilities.  NCRP  Report  No.  151.”  Bethesda,  MD  (USA):  National  Council  on  Radiation  Protection  and  Measurements,  2005,  ISBN-13:  978-0-929600-87-1
【非特許文献４】McGinley  PH,  “Shielding  techniques  for  radiation  oncology  facilities.”  Madison,  WI  (USA):  Medical  Physics,    2002,  ISBN-13:  978-1930524071
【発明の概要】
【発明が解決しようとする課題】
【０００８】
ところで、近年の免震構造の普及により、照射室直下の免震層の線量評価のニーズが増えている。
【０００９】
上記の従来の簡易計算法では、免震層のある床に向けて照射したＸ線のように、Ｘ線が床を透過し、免震層内の天井の低い空間で反射するような経路を通る評価点については、精度良く評価することが難しい。
【００１０】
また、従来のモンテカルロ法による線量計算は、簡易計算式に比べ計算時間が必要であり、計算体系によっては数日から１週間程度必要である。
【００１１】
本発明は、上記に鑑みてなされたものであって、放射線照射室の床下などの室外空間の線量を精度良く計算することのできる線量計算方法、線量計算プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１２】
上記した課題を解決し、目的を達成するために、本発明に係る線量計算方法は、放射線を照射する線源を有する放射線照射室の周囲に設けられた放射線遮蔽用の遮蔽体の外側の室外空間の線量を計算する方法であって、線源から照射されて遮蔽体を透過する放射線のうち少なくとも一部が遮蔽体内で散乱して、室外空間内に設定した所定の評価点に到達すると仮定して、線源から照射された放射線が遮蔽体内で散乱して評価点に到達する寄与の値と、線源から照射された放射線が遮蔽体を透過した後に、遮蔽体の外側の室外空間で多重反射して評価点に到達する寄与の値とを合算して評価点の線量を求めることを特徴とする。
【００１３】
また、本発明に係る他の線量計算方法は、上述した発明において、室外空間での放射線の多重反射の反射率を、エネルギーの減衰によらず一定値と仮定して、室外空間で多重反射して評価点に到達する寄与の値を求めることを特徴とする。
【００１４】
また、本発明に係る他の線量計算方法は、上述した発明において、線源から照射された放射線のエネルギーと、遮蔽体の材質と、室外空間を挟んで遮蔽体の表面に対向する物体の材質と、線源から照射された放射線が遮蔽体を透過した後に、物体の表面で反射して評価点に向かう放射線による線量率と、物体の表面で反射した後、遮蔽体の表面で反射して評価点に向かう放射線による線量率とに基づいて、室外空間での放射線の多重反射の反射率を求め、求めた反射率を用いて、室外空間で多重反射して評価点に到達する寄与の値を求めることを特徴とする。
【００１５】
また、本発明に係る線量計算プログラムは、上述した線量計算方法をコンピュータに実行させることを特徴とする。
【発明の効果】
【００１６】
本発明に係る線量計算方法によれば、放射線を照射する線源を有する放射線照射室の周囲に設けられた放射線遮蔽用の遮蔽体の外側の室外空間の線量を計算する方法であって、線源から照射されて遮蔽体を透過する放射線のうち少なくとも一部が遮蔽体内で散乱して、室外空間内に設定した所定の評価点に到達すると仮定して、線源から照射された放射線が遮蔽体内で散乱して評価点に到達する寄与の値と、線源から照射された放射線が遮蔽体を透過した後に、遮蔽体の外側の室外空間で多重反射して評価点に到達する寄与の値とを合算して評価点の線量を求めるので、例えば放射線照射室の床下などの室外空間の線量を、モンテカルロ法によることなく迅速かつ精度良く計算することができるという効果を奏する。
【００１７】
また、本発明に係る他の線量計算方法によれば、室外空間での放射線の多重反射の反射率を、エネルギーの減衰によらず一定値と仮定して、室外空間で多重反射して評価点に到達する寄与の値を求めるので、室外空間の線量を迅速に計算することができるという効果を奏する。
【００１８】
また、本発明に係る他の線量計算方法によれば、線源から照射された放射線のエネルギーと、遮蔽体の材質と、室外空間を挟んで遮蔽体の表面に対向する物体の材質と、線源から照射された放射線が遮蔽体を透過した後に、物体の表面で反射して評価点に向かう放射線による線量率と、物体の表面で反射した後、遮蔽体の表面で反射して評価点に向かう放射線による線量率とに基づいて、室外空間での放射線の多重反射の反射率を求め、求めた反射率を用いて、室外空間で多重反射して評価点に到達する寄与の値を求めるので、着目する経路に多重反射を含む放射線の線量を迅速かつ精度良く計算することができるという効果を奏する。
【００１９】
また、本発明に係る線量計算プログラムによれば、上述した線量計算方法をコンピュータに実行させるので、室外空間の線量をモンテカルロ法によることなく迅速かつ精度良く計算することができるという効果を奏する。
【図面の簡単な説明】
【００２０】
【図１】図１は、本発明に係る線量計算方法、線量計算プログラムの実施の形態の説明図である。
【図２】図２は、本発明に係る線量計算方法、線量計算プログラムの実施の形態を示すフローチャート図である。
【図３】図３は、モンテカルロ法による計算結果と、本発明による計算結果の比較図である。
【図４】図４は、反射率の計算フローチャートの一例を示す図である。
【図５】図５は、従来の線量評価方法の説明図である。
【発明を実施するための形態】
【００２１】
以下に、本発明に係る線量計算方法、線量計算プログラムの実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
【００２２】
（線量計算方法）
まず、本発明に係る線量計算方法の実施の形態について説明する。
【００２３】
図１に示すように、線源Ｓが放射線照射室に設置されており、線源ＳからＸ線（放射線）がビーム状（X-ray  beam）に照射されるものする。放射線照射室の床は放射線遮蔽用の遮蔽体ａで構成され、遮蔽体ａの直下には空間ｂ（室外空間）が配置されているものとする。本実施の形態では、放射線照射室として医療用リニアック室を想定し、空間ｂとして免震層を想定している。
【００２４】
線源Ｓからビーム状に照射したＸ線について、遮蔽体ａ内で散乱し空間ｂ内の評価点Ｅに到達するＸ線（経路Ｉ）、遮蔽体ａを透過した後に空間ｂで多重反射を起こし、最終的に空間ｂの下側の床（または奥側の壁）ｅの表面ｄで反射し評価点Ｅに到達するＸ線（経路ＩＩ）、遮蔽体ａを透過した後に空間ｂの中で多重反射を経て、最終的に遮蔽体ａの表面ｃで散乱し評価点Ｅに到達するＸ線（経路ＩＩＩ）の３経路に分けて考える。図１に、Ｘ線の経路を示す。
【００２５】
Ｘ線の経路Ｉの実効線量率ＤＳ（Ｓｖ／ｈ）は、遮蔽体ａと空間ｂの境界の壁表面ｃでの散乱とみなし、以下の式（２）で計算する。ここで、ＷはＸ線の線源Ｓから真下に１ｍ離れた位置での最大線量率（Ｇｙ／ｈ）、Ｂは遮蔽体ａの透過率、ＦはＧｙをＳｖに変える換算係数（Ｓｖ／Ｇｙ）、Ｕは方向利用率、ＡＣはＸ線利用線錐の表面ｃにおける照射面積、αＳは表面ｃでの反射点（Ｘ線利用線錐の中心軸と表面ｃの交点）における角度θＳでの反射率、ｄＣ０は線源Ｓから表面ｃでの反射点までの距離、ｄＣは表面ｃ上の反射点から評価点Ｅまでの距離である。
【００２６】
【数２】
2020110367.tif 000004
【００２７】
Ｘ線の経路ＩＩの実効線量率ＤＦ（Ｓｖ／ｈ）については、次のように考える。Ｘ線の一部は遮蔽体ａを透過した後に遮蔽体表面ｃと床表面ｄを往復するように多重反射を繰り返す。また、床表面ｄでの反射の際に、一部が反射率αＦで評価点Ｅに向かう。多重反射の反射率を一定値ｒとし、公比ｒ２の無限等比級数の和として扱い、Ｘ線の経路ＩＩの実効線量率ＤＦを以下の式（３）で計算する。ここで、ＡＦはＸ線利用線錐の床表面ｄにおける照射面積、αＦは床表面ｄでの反射点（Ｘ線利用線錐の中心軸と床表面ｄの交点）における角度θＦでの反射率、ｄＦ０は線源Ｓから床表面ｄでの反射点までの距離、ｄＦは床表面ｄ上の反射点から評価点Ｅまでの距離である。
【００２８】
【数３】
2020110367.tif 000005
【００２９】
Ｘ線の経路ＩＩＩの実効線量率ＤＣ（Ｓｖ／ｈ）は、遮蔽体ａを透過し、床表面ｄから表面ｃに反射率ｒで反射したＸ線が、空間ｂ内で多重反射を起こすと考え、経路ＩＩと同様に無限等比級数の和の式を使い、以下の式（４）で計算できる。ここで、αＣは表面ｃでの反射点（Ｘ線利用線錐の中心軸と表面ｃの交点）における角度θＣでの反射率である。
【００３０】
【数４】
2020110367.tif 000006
【００３１】
上記の３式（式（２）〜（４））をまとめたものが空間ｂの評価点Ｅにおける実効線量率Ｄ（Ｓｖ／ｈ）である。特に、線源Ｓから真下１ｍの位置での照射面積をＡＩＣとすると、ＡＩＣ＝ＡＣ／ｄＣ０２＝ＡＦ／ｄＦ０２であるため、以下の式（５）で表すことができる。
【００３２】
【数５】
2020110367.tif 000007
【００３３】
評価点Ｅが空間ｂの中間の高さの場合、反射率がαＦ＝αＣとなる。また、反射点から評価点Ｅまでの距離は、Ｘ線利用線錐の中心軸から評価点Ｅまでの距離ｄに近似（ｄ≒ｄＣ≒ｄＦ）することができる。
【００３４】
評価点Ｅが空間ｂ内のＸ線利用線錐が直接届く範囲とその近傍にある場合は、本発明の方法で得られる実効線量率が過小評価になるため注意が必要である。単純なＸ線の透過線が支配的になるからである。
【００３５】
透過率Ｂの計算は周知の手法を使うことができる。一例として、上記の非特許文献１に記載の手法により、Ｘ線のエネルギー、遮蔽体の材質（コンクリート、鉄、鉛など）と厚さから求めることができる。
【００３６】
反射率αＳ，αＦ，αＣの計算は周知の手法を使うことができる。一例として、ＣｈｉｌｔｏｎとＨｕｄｄｌｅｓｔｏｎの半経験式により、Ｘ線のエネルギーと反射角から求めることができる。
【００３７】
多重反射の反射率ｒは、例えば、後述の方法により計算することが可能であるが、後述の方法に限定されない。なお、反射率ｒは、加速器エネルギーや表面ｃ、床表面ｄの材質に応じた値がデータベース化され、計算時に必要な値が得られることが望ましい。
【００３８】
したがって、実際の計算では、Ｘ線の照射条件（線源から１ｍ離れた位置での最大線量率、照射面積、方向利用率）、施設のレイアウト、遮蔽体ａの遮蔽率、評価点Ｅの位置、空間ｂ内の多重反射の反射率ｒ、表面ｃ、床表面ｄでの反射角に応じた反射率が必要である。換算係数Ｆは定数であり、１を使用する。
【００３９】
（線量計算プログラム）
次に、本発明に係る線量計算プログラムの実施の形態について説明する。本実施の形態に係る線量計算プログラムは、上記の線形計算方法をコンピュータに実行させるものである。
【００４０】
図２は、本実施の形態の線量計算プログラムのフローチャートの一例である。この計算プログラムでは、空間ｂ内の任意の点の実効線量を計算する。
【００４１】
まず、Ｘ線の照射条件（線源Ｓから１ｍ離れた位置での最大線量率、照射面積、方向利用率、加速器エネルギー）、遮蔽体ａの構造（材質、厚さ）、評価点Ｅの位置をパラメーターとして入力する（ステップＳ１）。透過率や反射率を入力パラメーターとしてもよい。
【００４２】
加速器エネルギーと遮蔽体ａの構造から、遮蔽体ａの透過率Ｂが得られる（ステップＳ４）。評価点Ｅの位置に応じて、反射点から評価点Ｅの距離と反射角を求め（ステップＳ２）、エネルギーと反射角に応じた反射率αＳ，αＦ，αＣが得られる（ステップＳ３）。予め準備したデータベースから加速器エネルギーや表面ｃ、床表面ｄの材質に応じた多重反射の反射率ｒが得られる（ステップＳ５）。以上の値を使うことで、評価点Ｅの実効線量が計算される（ステップＳ６）。
【００４３】
本実施の形態によれば、線源Ｓから照射されて遮蔽体ａを透過するＸ線のうち少なくとも一部が遮蔽体ａ内で散乱して評価点Ｅに到達すると仮定して、線源Ｓから照射されたＸ線が遮蔽体ａ内で散乱して評価点Ｅに到達する寄与の値と、線源Ｓから照射されたＸ線が遮蔽体ａを透過した後に、遮蔽体ａの外側の空間ｂで多重反射して評価点Ｅに到達する寄与の値とを合算して評価点Ｅの実効線量を求める。このため、従来、適切な簡易計算式がなかった医療用リニアック室直下の免震層の実効線量を、モンテカルロ法によることなく迅速かつ精度良く計算することができる。
【００４４】
また、空間ｂのような層間でのＸ線の多重反射がある場合に、エネルギーの減衰によらず多重反射の反射率ｒを一定値と仮定して、無限等比級数の和の考え方を使って定式化し、空間ｂで多重反射して評価点Ｅに到達する寄与の値を求めるので、空間ｂの実効線量を迅速に計算することができる。
【００４５】
なお、上記の実施の形態は、ビーム状にＸ線が照射される施設であればいかなる施設にも適用可能であり、医療用リニアックに限定されない。また、線源からのＸ線の照射方向は床方向に限定されず、例えば照射室の天井方向や、側壁方向に遮蔽体と同様の構造がある場合に、適用可能である。
【００４６】
（線量計算装置）
次に、本発明に係る線量計算装置の実施の形態について説明する。
【００４７】
本実施の形態に係る線量計算装置は、上述した線量計算方法を装置として具現化したものであり、例えば入力部、記憶部、演算部、出力部とからなる。この線量計算装置は、例えばＣＰＵを有するコンピュータ、メモリ、ディスプレイ、キーボード等のハードウェア、これらハードウェアを用いて実行されるコンピュータプログラム等のソフトウェアにより構成することができる。
【００４８】
入力部は、上記のステップＳ１の入力パラメータを入力するためのものであり、例えばキーボードなどで構成することができる。記憶部は、加速器エネルギーや表面ｃ、床表面ｄの材質に応じた値などの計算用データを記憶するデータベースなどで構成することができる。
【００４９】
演算部は、上記のステップＳ２〜Ｓ６の各計算を実行するものであり、例えばコンピュータと演算ソフトウェアなどで構成することができる。出力部は、上記の演算部による演算処理結果を出力するものであり、例えばディスプレイやプリンタなどで構成することができる。
【００５０】
このように構成した線量計算装置によれば、従来、適切な簡易計算式がなかった医療用リニアック室直下の免震層の実効線量を、モンテカルロ法によることなく迅速かつ精度良く計算することができる。
【００５１】
（本発明の効果の検証）
次に、本発明の効果を検証するために行った計算結果について説明する。
【００５２】
本検証では、上記の線量計算方法による計算結果（実施例）と、モンテカルロ法による計算結果（比較例）を比較することによって、本発明の効果を検証した。計算条件は、Ｘ線が遮蔽体ａから２．２９５ｍ離れた位置から真下に円錐状に照射されるものとし、照射条件はＷ＝３６０Ｇｙ／ｈ、ＡＩＣ＝０．１６ｍ２、Ｕ＝１とした。遮蔽体ａと床ｅは普通コンクリートとし、遮蔽体ａの厚さは１．８ｍ、透過率は上記の非特許文献１に記載の手法により求めてＢ＝１．１８×１０−４２．１／１８０とした。空間ｂの高さは１．０ｍとし、評価点ＥはＸ線利用線錐の中心軸から水平方向に３〜９ｍ離れ、空間ｂの中間の高さに位置するものとした。また、空間ｂ内の多重反射の反射率ｒは０．４とした。比較例のモンテカルロ計算コードには、ＭＣＮＰ５を使用し、利用線錐の中心軸からの水平距離１ｍ毎の点を評価した。
【００５３】
図３に比較結果を示す。この図に示すように、本実施例と比較例の結果の差は最大で１９％であった。したがって、本発明によれば、従来、適切な簡易計算式がなかった医療用リニアック室直下の免震層の実効線量を比較的精度良く計算できることがわかる。
【００５４】
（多重反射の反射率の計算方法）
次に、上記の多重反射の反射率の計算方法の一例について説明する。
【００５５】
上記の実施の形態では、空間ｂ内でＸ線が多重反射を生じ、その時の反射率（多重散乱率）をエネルギーの減衰によらず一定値として扱うことで実効線量を簡易に計算できる方法を示したが、Ｘ線のエネルギーや周囲の材質に応じた多重反射の反射率が必要である。そこで、反射率を以下のようにして計算する。
【００５６】
まず、評価点Ｅが空間ｂの中間の高さの場合の実効線量率ＤＦ、ＤＣをそれぞれＤＦ’、ＤＣ’とする。反射率がαＦ＝αＣに、反射点から評価点Ｅまでの距離がｄＦ＝ｄＦになるため、比を取ると以下の式（６）が成り立ち、反射率ｒを求めるにはＤＣ’とＤＦ’を求めればよいことがわかる。
【００５７】
【数６】
2020110367.tif 000008
【００５８】
図４は、多重反射の反射率ｒを求める計算プログラムのフローチャートの一例である。この図に示すように、まず、加速器エネルギーと遮蔽体ａおよび床ｅの材質をパラメーターとして入力する（ステップＳ１１）。床表面ｄおよび床ｅでのＸ線の相互作用を無視させたモンテカルロ計算により、空間ｂの中間高さに位置する評価点Ｅで実効線量ＤＳが算出される（ステップＳ１３）。同じ位置にある評価点Ｅに対し、遮蔽体ａ側から評価点Ｅに到達するＸ線のみの実効線量を求めるモンテカルロ計算により、実効線量ＤＳ＋ＤＣ’を計算する（ステップＳ１２）。さらに、床ｅ側から評価点Ｅに到達するＸ線のみの実効線量を求めるモンテカルロ計算により、実効線量ＤＦ’を計算する（ステップＳ１５）。ＤＳとＤＳ＋ＤＣ’から、ＤＣ’を求め（ステップＳ１４）、ＤＦ’とＤＣ’から反射率ｒを求める（ステップＳ１６）。
【００５９】
上記の反射率ｒの計算方法によれば、１つの照射条件から、プログラム内の実効線量の計算方法を変えるだけで遮蔽体ａ透過後の空間ｂにおける多重反射の反射率ｒを計算により求めることができる。着目する経路に多重反射を含むＸ線の実効線量計算が可能になる。
【００６０】
以上説明したように、本発明に係る線量計算方法によれば、放射線を照射する線源を有する放射線照射室の周囲に設けられた放射線遮蔽用の遮蔽体の外側の室外空間の線量を計算する方法であって、線源から照射されて遮蔽体を透過する放射線のうち少なくとも一部が遮蔽体内で散乱して、室外空間内に設定した所定の評価点に到達すると仮定して、線源から照射された放射線が遮蔽体内で散乱して評価点に到達する寄与の値と、線源から照射された放射線が遮蔽体を透過した後に、遮蔽体の外側の室外空間で多重反射して評価点に到達する寄与の値とを合算して評価点の線量を求めるので、例えば放射線照射室の床下などの室外空間の線量を、モンテカルロ法によることなく迅速かつ精度良く計算することができる。
【００６１】
また、本発明に係る他の線量計算方法によれば、室外空間での放射線の多重反射の反射率を、エネルギーの減衰によらず一定値と仮定して、室外空間で多重反射して評価点に到達する寄与の値を求めるので、室外空間の線量を迅速に計算することができる。
【００６２】
また、本発明に係る他の線量計算方法によれば、線源から照射された放射線のエネルギーと、遮蔽体の材質と、室外空間を挟んで遮蔽体の表面に対向する物体の材質と、線源から照射された放射線が遮蔽体を透過した後に、物体の表面で反射して評価点に向かう放射線による線量率と、物体の表面で反射した後、遮蔽体の表面で反射して評価点に向かう放射線による線量率とに基づいて、室外空間での放射線の多重反射の反射率を求め、求めた反射率を用いて、室外空間で多重反射して評価点に到達する寄与の値を求めるので、着目する経路に多重反射を含む放射線の線量を迅速かつ精度良く計算することができる。
【００６３】
また、本発明に係る線量計算プログラムによれば、上述した線量計算方法をコンピュータに実行させるので、室外空間の線量をモンテカルロ法によることなく迅速かつ精度良く計算することができる。
【産業上の利用可能性】
【００６４】
以上のように、本発明に係る線量計算方法、線量計算プログラムは、医療用リニアック室などの放射線照射室の直下にある免震層などの室外空間の光子線による放射線量を計算するのに有用であり、特に、実効線量を迅速かつ精度良く計算するのに適している。
【符号の説明】
【００６５】
Ｅ  評価点
Ｓ  線源
ａ  遮蔽体
ｂ  空間
ｃ  表面
ｄ  表面
ｅ  床（または壁）
【図１】
2020110367.tif 000009
【図２】
2020110367.tif 000010
【図３】
2020110367.tif 000011
【図４】
2020110367.tif 000012
【図５】
2020110367.tif 000013