(57)【特許請求の範囲】
【請求項１】
医療用放射線の照射のために治療対象である目的体系の構成および前記目的体系と前記医療用放射線との相互作用の特性をコンピュータ上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差を、対象を模擬した互いに異なる模擬実験ｉ（ｉ＝１，・・・，ｎ）の結果を用いて推定する照射線量計算誤差推定装置において、
前記シミュレーションのパラメータの入力を外部から受け入れる入力装置と、
前記模擬実験についての前記モデルを用いたシミュレーションで得られた物理量の計算値の当該模擬実験で測定された当該物理量の測定値に対する前記モデルへの入力値に関する相対誤差Ｅpを演算する計算値・測定値相対誤差演算部と、
前記計算値・測定値相対誤差演算部が算出した前記相対誤差Ｅpを記憶する計算値・測定値相対誤差記憶部と、
前記シミュレーションに用いる前記モデルへの前記パラメータの入力値の不確かさの割合の関係を示す入力パラメータ共分散誤差行列Ｗを演算する入力パラメータ共分散誤差行列評価部と、
前記入力パラメータ共分散誤差行列評価部が算出した前記入力パラメータ共分散誤差行列Ｗを記憶する入力パラメータ共分散誤差行列記憶部と、
前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を示す目的体系感度係数ベクトルＳＲを演算する目的体系感度係数ベクトル演算部と、
前記目的体系感度係数ベクトル演算部が算出した前記目的体系感度係数ベクトルＳＲを記憶する目的体系感度係数ベクトル記憶部と、
前記模擬実験の体系について前記モデルを用いてシミュレーションした結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を表す模擬実験体系感度係数ベクトルＳｉを演算する模擬実験体系感度係数ベクトル演算部と、
前記模擬実験体系感度係数ベクトル演算部が算出した前記模擬実験体系感度係数ベクトルＳｉを記憶する模擬実験体系感度係数ベクトル記憶部と、
前記模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲとに基づき算出する線形結合定数演算部と、
前記計算値・測定値相対誤差記憶部に記憶された前記相対誤差ＥＰ、前記入力パラメータ共分散誤差行列記憶部に記憶された前記入力パラメータ共分散誤差行列Ｗ、前記目的体系感度係数ベクトル記憶部に記憶された前記目的体系感度係数ベクトルＳＲおよび前記模擬実験体系感度係数ベクトル記憶部に記憶された前記模擬実験体系感度係数ベクトルＳｉを用いて、前記モデルを用いたシミュレーションの結果得られる前記対象の物理量に含まれる入力パラメータ起因物理量不確かさを算出する目的体系物理量不確かさ演算部と、
を有し、
前記入力パラメータ共分散誤差行列評価部は、前記パラメータの数と前記模擬実験の数とが一致している場合に、前記入力パラメータ共分散誤差行列Ｗを対角行列とし、その成分を次の式により算出することを特徴とする照射線量計算誤差推定装置。
【数１】
0006490535.tif 000027
ｓｉｊは、ｉ番目の模擬実験のｊ番目のパラメータに関する感度係数の要素、ｑｉ２は、測定値と計算値の相対差の２乗、ｎは模擬実験の数、パラメータの数を表す。
【請求項２】
前記入力パラメータ共分散誤差行列評価部は、前記入力パラメータ共分散誤差行列Ｗを対角行列として決定することを特徴とする請求項１に記載の照射線量計算誤差推定装置。
【請求項３】
前記対角行列は単位行列であることを特徴とする請求項２に記載の照射線量計算誤差推定装置。
【請求項４】
医療用放射線の照射のために治療対象である目的体系の構成および前記目的体系と前記医療用放射線との相互作用の特性をコンピュータ上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差を、対象を模擬し、前記物理量を支配するパラメータの個数より少ない回数実施された互いに異なる模擬実験の結果を用いて推定する照射線量計算誤差推定装置において、
前記シミュレーションのパラメータの入力を外部から受け入れる入力装置と、
前記模擬実験についての前記モデルを用いたシミュレーションで得られた物理量の計算値の当該模擬実験で測定された当該物理量の測定値に対する前記モデルへの入力値に関する相対誤差Ｅpを演算する計算値・測定値相対誤差演算部と、
前記計算値・測定値相対誤差演算部が算出した前記相対誤差Ｅpを記憶する計算値・測定値相対誤差記憶部と、
前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を示す目的体系感度係数ベクトルＳＲを演算する目的体系感度係数ベクトル演算部と、
前記目的体系感度係数ベクトル演算部が算出した前記目的体系感度係数ベクトルＳＲを記憶する目的体系感度係数ベクトル記憶部と、
第ｉ番目（ｉは１からｎのそれぞれ）の模擬実験の体系について前記モデルを用いてシミュレーションした結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を表す模擬実験体系感度係数ベクトルＳｉを演算する模擬実験体系感度係数ベクトル演算部と、
前記模擬実験体系感度係数ベクトル演算部が算出した前記模擬実験体系感度係数ベクトルＳｉを記憶する模擬実験体系感度係数ベクトル記憶部と、
前記模擬実験体系感度係数ベクトルＳｉを用いて、同一模擬実験に対するｉ番目の模擬実験の感度係数ベクトル２次形式ＳｉＴＷＳｉを次の式（１）により算出された値として求める感度係数ベクトル２次形式第１種類演算部と、
【数２】
0006490535.tif 000028
前記感度係数ベクトル２次形式第１種類演算部で算出された感度係数ベクトル２次形式ＳｉＴＷＳｉを記憶する感度係数ベクトル２次形式第１種類記憶部と、
第ｉ番目の模擬実験と第ｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊを次の式（２）により算出する感度係数ベクトル２次形式第２種類演算部と、
【数３】
0006490535.tif 000029
前記感度係数ベクトル２次形式第２種類演算部で算出された第ｉ番目の模擬実験と第ｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊを記憶する感度係数ベクトル２次形式第２種類記憶部と、
目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを、次の式（３）の最大の値を与える模擬実験についての前記相対誤差Ｅｐの二乗ｑｍ２を用いて、次の式（４）により、感度係数ベクトル２次形式ＳＲＴＷＳＲを算出する感度係数ベクトル２次形式第３種類演算部と、
【数４】
0006490535.tif 000030
【数５】
0006490535.tif 000031
前記感度係数ベクトル２次形式第３種類演算部で算出された目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを記憶する感度係数ベクトル２次形式第３種類記憶部と、
ｉ番目の模擬実験についての模擬実験体系感度係数ベクトルＳｉと、目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを次の式（５）により算出する感度係数ベクトル２次形式第４種類演算部と、
【数６】
0006490535.tif 000032
前記感度係数ベクトル２次形式第４種類演算部で算出されたｉ番目の模擬実験についての模擬実験体系感度係数ベクトルＳｉと、目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを記憶する感度係数ベクトル２次形式第４種類記憶部と、
模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲとに基づき算出する線形結合定数演算部と、
前記感度係数ベクトル２次形式第１種類記憶部に記憶された第ｉ番目の模擬実験に関する感度係数ベクトル２次形式ＳｉＴＷＳｉと、前記感度係数ベクトル２次形式第２種類記憶部に記憶された第ｉ番目の模擬実験と第ｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊと、前記感度係数ベクトル２次形式第３種類記憶部に記憶された目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲと、前記感度係数ベクトル２次形式第４種類記憶部に記憶されたｉ番目の模擬実験についての模擬実験体系感度係数ベクトルＳｉとを用いて、前記モデルを用いたシミュレーションの結果得られる前記対象の物理量に含まれる入力パラメータ起因物理量不確かさを算出する目的体系物理量不確かさ演算部と、
を有することを特徴とする照射線量計算誤差推定装置。
【請求項５】
前記目的体系感度係数ベクトル演算部による前記目的体系感度係数ベクトルＳＲの演算、および前記模擬実験体系感度係数ベクトル演算部による前記模擬実験体系感度係数ベクトルＳｉの演算は、それぞれの体系において放射線の輸送をモンテカルロ計算で行うことにより算出することを特徴とする請求項１ないし請求項４のいずれか一項に記載の照射線量計算誤差推定装置。
【請求項６】
前記物理量は、前記医療用放射線の到達距離、および患部に照射される放射線量であることを特徴とする請求項１ないし請求項５のいずれか一項に記載の照射線量計算誤差推定装置。
【請求項７】
前記パラメータは、前記医療用放射線の線種、前記医療用放射線のエネルギー、および前記医療用放射線の方向のうち少なくともいずれか一つであることを特徴とする請求項１ないし請求項６のいずれか一項に記載の照射線量計算誤差推定装置。
【請求項８】
医療用放射線の照射のために治療対象である目的体系の構成および前記目的体系と前記医療用放射線との相互作用の特性をコンピュータ上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差をその対象を模擬した模擬実験の結果を用いて推定する照射線量計算誤差推定方法において、
入力装置が、模擬実験体系データおよび測定値、計算のパラメータを入力として受け入れる入力ステップと、
模擬実験体系測定値記憶部が、模擬実験体系データおよび測定値を記憶する実験結果記憶ステップと、
模擬実験体系物理量演算部が、前記模擬実験の体系についてシミュレーションを行い模擬実験体系における物理量を算出する模擬実験体系シミュレーションステップと、
模擬実験体系演算値記憶部が、前記模擬実験体系物理量演算部が算出した模擬実験体系における物理量の計算値を記憶する模擬実験体系計算値記憶ステップと、
前記模擬実験体系計算値記憶ステップの後に、計算値・測定値相対誤差演算部が、前記実験結果記憶ステップで記憶された測定値と、前記模擬実験体系計算値記憶ステップで記憶されたシミュレーション結果の計算値に基づいて両者の相対誤差ＥＰを算出し、前記計算値・測定値相対誤差演算部が算出した前記相対誤差ＥＰを記憶するＥＰ算出ステップと、
計算値・測定値相対誤差記憶部が、前記相対誤差ＥＰを記憶するＥＰ記憶ステップと、
模擬実験体系感度係数ベクトル演算部が、前記模擬実験（ｉ＝１〜ｎ）のそれぞれの体系の感度係数ベクトルである模擬実験体系感度係数ベクトルＳｉを算出するＳｉ算出ステップと、
模擬実験体系感度係数ベクトル記憶部が、前記模擬実験体系感度係数ベクトルＳｉを記憶するＳｉ記憶ステップと、
目的体系物理量演算部が、前記目的体系の物理量の計算結果Ｒ０を算出するＲ０算出ステップと、
目的体系演算値記憶部が、前記計算結果Ｒ０を記憶するＲ０記憶ステップと、
目的体系感度係数ベクトル演算部が、前記目的体系の入力パラメータに対する感度係数ベクトルである目的体系感度係数ベクトルＳＲを算出するＳＲ算出ステップと、
目的体系感度係数ベクトル記憶部が前記目的体系感度係数ベクトルＳＲを記憶するＳＲ記憶ステップと、
感度係数ベクトル２次形式第１種類演算部が、前記相対誤差ＥＰに基づいて第１の感度係数ベクトル２次形式を算出する２次形式第１種類演算ステップと、
感度係数ベクトル２次形第２種類式演算部が、前記相対誤差ＥＰに基づいて第２の感度係数ベクトル２次形式を算出する２次形式第２種類演算ステップと、
感度係数ベクトル２次形式第３種類演算部が、前記相対誤差ＥＰ、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲに基づいて第３の感度係数ベクトル２次形式を算出する２次形式第３種類演算ステップと、
感度係数ベクトル２次形式第４種類演算部が、前記相対誤差ＥＰ、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲに基づいて第４の感度係数ベクトル２次形式を算出する２次形式第４種類演算ステップと、
線形結合定数演算部が、模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲに基づき算出する線形結合定数演算ステップと、
目的体系物理量不確かさ演算部が、前記第１の感度係数ベクトル２次形式、前記第２の感度係数ベクトル２次形式、前記第３の感度係数ベクトル２次形式、前記第４の感度係数ベクトル２次形式、および前記線形結合定数αｉに基づいて、入力パラメータ起因物理量不確かさの二乗（ΔＺ／Ｚ）２を算出する不確かさ算出ステップと、
を有することを特徴とする照射線量計算誤差推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、医療用放射線照射のための照射線量計算誤差推定装置および照射線量計算誤差推定方法に関する。
【背景技術】
【０００２】
放射線は医療分野でこれまで診断に広く活用されてきた。Ｘ線によるレントゲン写真、Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ：コンピュータ断層撮影）、ＰＥＴ（Ｐｏｓｉｔｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇｒａｐｈｙ：陽電子放出断層撮影)が現在広く利用されている。
【０００３】
医学において、放射線を診断のみではなく積極的に利用して治療効果を得ようするものに、放射線治療装置、放射線治療機器、放射線治療システムと呼ばれるものがある。ガンマナイフ、サイバーナイフ、重粒子線がん治療装置、ＢＮＣＴ（Ｂｏｒｏｎ  Ｎｅｕｔｒｏｎ  Ｃａｐｔｕｒｅ  Ｔｈｅｒａｐｙ：ホウ素中性子捕捉療法）と呼ばれるものがその例であり、人体に発生した好ましくない腫瘍（悪性の細胞の塊）等を放射線のエネルギーによって破壊することを目的としている。
【０００４】
放射線は細胞を破壊する能力があるので、破壊するべきではない正常細胞や正常組織に対しても好ましくない影響を与える可能性がある。従って、放射線治療器に要求される性能は、医師が患者を診断・判断して、破壊すべき細胞や組織を特定しそれらのみを正確に破壊すること（患者にとっては治療）であり、目的としない他の生体部分へは、影響を最小限に抑えることが要求される。
【０００５】
放射線を照射したくない場所には照射しないことが原則であり、そのためには放射線治療装置から放出する放射線を、治療計画に従って、放射線の方向およびエネルギーを意図した数値で正確に制御して照射する必要がある。即ち、放射線治療装置では、放出される放射線の位置精度、エネルギー精度が極めて重要になる。
【０００６】
現代の放射線治療器による治療の流れは以下のようなものである。
【０００７】
（１）まず放射線診断器によって生体（患者）の患部あるいは患部を含めた人体の正確な画像を得る。（２）その画像によって治療すべき対象物（患部）の３次元的位置を把握し特定する。（３）医師あるいは医療従事者が、患者の合意の下、治療計画を策定する。（４）医師あるいは医療従事者が、特定のコンピュータとコンピュータ・プログラムによる計算機シミュレーションに基づいて、治療計画を実現させる放射線治療器の使用方法（運用方法）を決める。（場合によっては計算機シミュレーションの結果として放射線治療器の動作を制御するソフトウェア（コンピュータプログラム）への入力が得られる）。
【０００８】
このソフトウェアは患者の治療時、即ち患者の患部に放射線を照射する場面において放射線治療器が自動的に動作することを制御するもので、通常、治療時には放射線治療器は人間によって直接制御されない。理由としては放射線治療器の治療成績を高めるために、機器の極めて微細な位置精度の制御や放射線の照射時間（タイミング）の制御を必要とし、人間の対応能力を越えているからである。
【０００９】
したがって、計算機シミュレーションの技術や精度が治療の品質を大きく左右することになる。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２００８−２１７１３９号公報
【特許文献２】特開２０１１−１０６９７０号公報
【特許文献３】特開２０１１−２５８０５２号公報
【特許文献４】特開２０１３−６５２１３号公報
【特許文献５】特開２０１３−２０３８６号公報
【非特許文献】
【００１１】
【非特許文献１】Journal of NUCLEAR SCIENCE and TECHNOLOGY Vol.44, No.12, Pages 1509-1517 (2007) : Theoretical Study on New Bias Factor Methods to Effectively Use Critical Experiments for Improvement of Prediction Accuracy of Neutronic Characteristics
【非特許文献２】Journal of NUCLEAR SCIENCE and TECHNOLOGY Vol.51, No.5, Pages 608-625 (2014) : Development of a “best representativity” method for experimental data analysis and an application to the critical experiments at the Toshiba NCA facility
【発明の概要】
【発明が解決しようとする課題】
【００１２】
一方、人体には、複雑な幾何形状と物質構成を持つ部分がある。現代の計算機シミュレーション技術をもってしても、計算機上に構築する人体の数値データモデル（計算モデル）は、放射線診断器によって得られた画像そのものを厳密な形で取り込むことは計算規模の面からほぼ不可能である。加えて放射線診断器では、人体の正確な画像データは得られても、その画像データが示す各組織の正確な物質組成は分からない。例えば画像データで骨が示されていてもその骨の正確な組成や質量（あるいは密度）は人体から取り出して測定しなければ精度の高い値は得られない。
【００１３】
現代の計算技術として大きく分けて２つの種類がある。一つは、決定論的手法（Ｄｅｔｅｒｍｉｎｉｓｔｉｃ  Ｍｅｔｈｏｄ）である。もう一つは、確率論的方法（Ｐｒｏｂａｌｉｓｔｉｃ  Ｍｅｔｈｏｄ）である。
【００１４】
放射線（光子や荷電粒子）の空間移動を記述する基本方程式はボルツマンの輸送方程式やＭａｘｗｅｌｌの方程式である。この方程式をディジタル・コンピュータで解く場合、２つの方法がある。すなわち、空間やエネルギーを離散化して解く方法（決定論）と、空間やエネルギーを連続変数で扱って確率論的に解く方法（確率論、モンテカルロ法）がある。一般にモンテカルロ計算では計算対象の体系の幾何形状の記述は非常に自由度が高く、ほぼ人体の幾何形状をそのまま計算機上に取り込んでモデル化して計算を実施することが技術的には可能である。
【００１５】
但し、計算機シミュレーションでは、当然ながら、正しい計算値（計算解）を得るには正しい入力値による必要があり、如何に人体の幾何形状が詳細なものであっても、人体を構成する物質組成が正しくなければ正確な数値解は得られない。
【００１６】
一方、非破壊検査的な手法で人体の組織の物質組成を測定して高い精度で数値化する方法は現時点では存在せず、従って、現在どんなに正確で高級な計算機シミュレーションが可能であっても、その計算機シミュレーションの入力が原理的に不十分なものであるという技術課題が存在する。
【００１７】
すなわちモンテカルロ法と放射線診断で得られた画像を基に、正確な幾何形状をとりこんでモンテカルロ法によって計算機シミュレーションを実施しても、人体の組成の記述の不十分さから得られる計算値は、未だ十分には満足できる品質には達していないと言える。
【００１８】
一方、放射線治療器の運用上、年間に治療できる件数は多い方が好ましい。放射線治療器を運用する運用側の経済的側面からもそうであるし、また１名の患者の治療に極めて長い時間がかかれば次の患者は、長い時間、待たされねばならず、好ましいこととは言えない。そのことは１名の患者に対して費やすべき治療計画の策定時間はある制限があり、従って治療計画、放射線治療器の制御パラメータの決定のために実施するコンピュータ・シミュレーションもある時間的制約の下で実施されることを意味する。
【００１９】
計算機シミュレーションに時間をかければ、放出すべき放射線の位置精度等が良くなるので治療精度は向上する。しかしながら現代の最高の技術水準に則り、人体の幾何形状を詳しく計算機中に取り込むほど、また人体の組成を実際の患者の測定値ではなく仮定の数値であっても詳しく計算機上に入力し取り込むほど、計算に必要な時間は長くなる。
【００２０】
そこで、それぞれの放射線治療器について一人あたりの患者に与えられる標準時間から、計算機シミュレーションの規模も自ずと決まり、常に治療精度を高くすることを志向しながらも、やはりそれぞれの場面に従って治療計画の策定や放射線治療器の制御に対しての標準的な計算機シミュレーションが決められている。
【００２１】
現代の計算機技術の発展に従って、計算速度は常に向上している。よって今後、それぞれの放射線治療器の治療計画に係る標準的な計算機シミュレーションもますます詳細化・高度化していくことは自明である。しかしながらそれでも前述のように、個々の患者に関して、計算対象である全ての組織の正確な組成を得ることはできないことから、組織を形成するたとえば蛋白質等の構造を正確に計算機上に取り込んで計算するに至る時期は相当先の将来のことであると判断される。
【００２２】
現在、計算機シミュレーションでしばしば用いられる近似は、人体の組成の一部あるいは大部分を水（Ｈ２Ｏ）として計算することである。この近似はそれほど悪いものではない。簡単な表現では、人体の約６０％が水でできているといえる。ちなみに、体重に占める水の割合は、胎児は約９０％、１〜２歳児は約７５％、児童は約７０％、成人で約６０％、高齢者では約５０％程度であり、年齢の上昇に従って、体重に占める水の割合は減少する。水分の割合の減少は、加齢によって体内の脂肪分が減少し、その脂肪分に含まれている水が同時に減少することによる。
【００２３】
このように、計算機上に正確に人体の幾何形状モデルを構築しても、人体の組成を全て水（Ｈ２Ｏ）で近似するとすれば、計算機シミュレーションで計算誤差が生じることは明らかである。
【００２４】
計算機によるコンピュータ・シミュレーションの規模は、ケースバイケースで選択されるが、幾何形状、物質構成などは何らかの近似を行うことになる。この様な近似に基づいて計算企画された治療計画、放射線治療器の制御は、やはり厳密な計算による解とは異なるはずであるが、これまで述べたように、その近似は医師や医療従事者、放射線治療器の専門家によって選択・決定される。
【００２５】
放射線治療器から放出される放射線について、到達する位置や与えるエネルギーが治療計画で定めた位置や値と異なることは、治療すべき患部への計画していた治療効果が損なわれるのみならず、放射線を照射すべきではない正常組織に放射線で損傷を引き起こす（副作用）可能性が生ずるので決して好ましいことではない。
以上に述べたように
（１）放射線治療器の治療全体に計算機シミュレーションが利用されている。
（２）治療精度の向上には、より高度のコンピュータ・シミュレーションを実施することが望ましいが、放射線治療器の運用や患者の待ち時間等の面から、現代の最高の技術水準の計算機シミュレーションを実施することは困難であり、それぞれの放射線治療器の標準的な計算機シミュレーションが設定され運用されている。
（３）人体の幾何形状は極めて正確に計算機上にモデル化できるが、人体の組成を正確に数値化できない。人体を多くの部分で水（Ｈ２Ｏ）として近似して計算することがある。
【００２６】
現代の最高の技術水準で、患者本人の幾何形状を詳細に計算機上にモデル化して、加えて組織の物質組成に関しても最も詳しい入力を行えば、きわめて正確なシミュレーションが実現できると期待される。
【００２７】
しかしながら、放射線治療器から放出される放射線の治療対象の患部への位置精度や与えるエネルギーの精度と比較しても、標準的な計画手法で決定された放射線治療器の制御パラメータによる治療精度（放射線の到達位置、与えるエネルギー等）は、未だ完全とは言えず、幾ばくかの誤差あるいは精度の劣化が見られる。
【００２８】
治療計画の策定時に例えば患者の身体を水（Ｈ２Ｏ）として計算して得られる、放射線治療器の制御パラメータに従って患者に放射線を照射した場合、本来、目的とする患部における放射線の到達位置や、放射線の与えるエネルギーなどに大きな差ではなくてもいくらかの差（誤差）が生じるのは当然である。
【００２９】
治療の前に実施された模擬実験ケースからの測定情報に従って、目的とする治療に関して計算機シミュレーションで生じる誤差を評価できれば、その誤差を修正することで更に放射線の治療精度を上げることができる。例えば放射線の到達距離に誤差が生じれば放射線治療器の放射線の放出位置をずらすとか、あるいは放射線の放出エネルギーを少し変化させるとか、改善する方法が存在するからである。
【００３０】
そこで、本発明の実施形態は、放射線治療器の計算機シミュレーションを補完して治療精度を向上させる照射線量計算誤差推定装置、照射線量計算誤差推定方法を提供することを目的とする。
【課題を解決するための手段】
【００３１】
上述の目的を達成するため、本実施形態は、医療用放射線の照射のために治療対象である目的体系の構成および前記目的体系と前記医療用放射線との相互作用の特性をコンピュータ上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差を、対象を模擬した互いに異なる模擬実験ｉ（ｉ＝１，・・・，ｎ）の結果を用いて推定する照射線量計算誤差推定装置において、前記シミュレーションのパラメータの入力を外部から受け入れる入力装置と、前記模擬実験についての前記モデルを用いたシミュレーションで得られた物理量の計算値の当該模擬実験で測定された当該物理量の測定値に対する前記モデルへの入力値に関する相対誤差Ｅpを演算する計算値・測定値相対誤差演算部と、前記計算値・測定値相対誤差演算部が算出した前記相対誤差Ｅpを記憶する計算値・測定値相対誤差記憶部と、前記シミュレーションに用いる前記モデルへの前記パラメータの入力値の不確かさの割合の関係を示す入力パラメータ共分散誤差行列Ｗを演算する入力パラメータ共分散誤差行列評価部と、前記入力パラメータ共分散誤差行列評価部が算出した前記入力パラメータ共分散誤差行列Ｗを記憶する入力パラメータ共分散誤差行列記憶部と、前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を示す目的体系感度係数ベクトルＳＲを演算する目的体系感度係数ベクトル演算部と、前記目的体系感度係数ベクトル演算部が算出した前記目的体系感度係数ベクトルＳＲを記憶する目的体系感度係数ベクトル記憶部と、前記模擬実験の体系について前記モデルを用いてシミュレーションした結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を表す模擬実験体系感度係数ベクトルＳｉを演算する模擬実験体系感度係数ベクトル演算部と、前記模擬実験体系感度係数ベクトル演算部が算出した前記模擬実験体系感度係数ベクトルＳｉを記憶する模擬実験体系感度係数ベクトル記憶部と、前記模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲとに基づき算出する線形結合定数演算部と、前記計算値・測定値相対誤差記憶部に記憶された前記相対誤差ＥＰ、前記入力パラメータ共分散誤差行列記憶部に記憶された前記入力パラメータ共分散誤差行列Ｗ、前記目的体系感度係数ベクトル記憶部に記憶された前記目的体系感度係数ベクトルＳＲおよび前記模擬実験体系感度係数ベクトル記憶部に記憶された前記模擬実験体系感度係数ベクトルＳｉを用いて、前記モデルを用いたシミュレーションの結果得られる前記対象の物理量に含まれる入力パラメータ起因物理量不確かさを算出する目的体系物理量不確かさ演算部と、を有し、前記入力パラメータ共分散誤差行列評価部は、前記パラメータの数と前記模擬実験の数とが一致している場合に、前記入力パラメータ共分散誤差行列Ｗを対角行列とし、その成分を次の式により算出することを特徴とする照射線量計算誤差推定装置。
【数１】
0006490535.tif 000002
ｓｉｊは、ｉ番目の模擬実験のｊ番目のパラメータに関する感度係数の要素、ｑｉ２は、測定値と計算値の相対差の２乗、ｎは模擬実験の数、パラメータの数を表す。
【００３２】
また、本実施形態は、医療用放射線の照射のために治療対象である目的体系の構成および前記目的体系と前記医療用放射線との相互作用の特性をコンピュータ上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差を、対象を模擬し、前記物理量を支配するパラメータの個数より少ない回数実施された互いに異なる模擬実験の結果を用いて推定する照射線量計算誤差推定装置において、前記シミュレーションのパラメータの入力を外部から受け入れる入力装置と、前記模擬実験についての前記モデルを用いたシミュレーションで得られた物理量の計算値の当該模擬実験で測定された当該物理量の測定値に対する前記モデルへの入力値に関する相対誤差ＥＰを演算する計算値・測定値相対誤差演算部と、前記計算値・測定値相対誤差演算部が算出した前記相対誤差ＥＰを記憶する計算値・測定値相対誤差記憶部と、前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を示す目的体系感度係数ベクトルＳＲを演算する目的体系感度係数ベクトル演算部と、前記目的体系感度係数ベクトル演算部が算出した前記目的体系感度係数ベクトルＳＲを記憶する目的体系感度係数ベクトル記憶部と、第ｉ番目（ｉは１からｎのそれぞれ）の模擬実験の体系について前記モデルを用いてシミュレーションした結果のそのモデルへの入力値の単位変化に対する前記物理量の変化量を表す模擬実験体系感度係数ベクトルＳｉを演算する模擬実験体系感度係数ベクトル演算部と、前記模擬実験体系感度係数ベクトル演算部が算出した前記模擬実験体系感度係数ベクトルＳｉを記憶する模擬実験体系感度係数ベクトル記憶部と、前記模擬実験体系感度係数ベクトルＳｉを用いて、同一模擬実験に対するｉ番目の模擬実験の感度係数ベクトル２次形式ＳｉＴＷＳｉを次の式（１）により算出された値として求める感度係数ベクトル２次形式第１種類演算部と、
【００３３】
【数１】
0006490535.tif 000003
【００３４】
前記感度係数ベクトル２次形式第１種類演算部で算出された感度係数ベクトル２次形式ＳｉＴＷＳｉを記憶する感度係数ベクトル２次形式第１種類記憶部と、第ｉ番目の模擬実験と第ｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊを次の式（２）により算出する感度係数ベクトル２次形式第２種類演算部と、
【００３５】
【数２】
0006490535.tif 000004
【００３６】
前記感度係数ベクトル２次形式第２種類演算部で算出された第ｉ番目の模擬実験と第ｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊを記憶する感度係数ベクトル２次形式第２種類記憶部と、目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを、次の式（３）の最大の値を与える模擬実験についての前記相対誤差Ｅｐの二乗ｑｍ２を用いて、次の式（４）により、感度係数ベクトル２次形式ＳＲＴＷＳＲを算出する感度係数ベクトル２次形式第３種類演算部と、
【００３７】
【数３】
0006490535.tif 000005
【００３８】
【数４】
0006490535.tif 000006
【００３９】
前記感度係数ベクトル２次形式第３種類演算部で算出された目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを記憶する感度係数ベクトル２次形式第３種類記憶部と、ｉ番目の模擬実験についての模擬実験体系感度係数ベクトルＳｉと、目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを次の式（５）により算出する感度係数ベクトル２次形式第４種類演算部と、
【００４０】
【数５】
0006490535.tif 000007
【００４１】
前記感度係数ベクトル２次形式第４種類演算部で算出されたｉ番目の模擬実験についての実験体系感度係数ベクトルと、目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲを記憶する感度係数ベクトル２次形式第４種類記憶部と、模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲとに基づき算出する線形結合定数演算部と、前記感度係数ベクトル２次形式第１種類記憶部に記憶された第ｉ番目の模擬実験に関する感度係数ベクトル２次形式ＳｉＴＷＳｉと、前記感度係数ベクトル２次形式第２種類記憶部に記憶された第ｉ番目の模擬実験と第ｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊと、前記感度係数ベクトル２次形式第３種類記憶部に記憶された目的とする患者の体系についての目的体系感度係数ベクトルＳＲの感度係数ベクトル２次形式ＳＲＴＷＳＲと、前記感度係数ベクトル２次形式第４種類記憶部に記憶されたｉ番目の模擬実験についての模擬実験体系感度係数ベクトルＳｉとを用いて、前記モデルを用いたシミュレーションの結果得られる前記対象の物理量に含まれる入力パラメータ起因物理量不確かさを算出する目的体系物理量不確かさ演算部と、を有することを特徴とする。
【００４３】
また、本実施形態は、医療用放射線の照射のために治療対象である目的体系の構成および前記目的体系と前記医療用放射線との相互作用の特性をコンピュータ上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差をその対象を模擬した模擬実験の結果を用いて推定する照射線量計算誤差推定方法において、入力装置が、模擬実験体系データおよび測定値、計算のパラメータを入力として受け入れる入力ステップと、模擬実験体系測定値記憶部が、模擬実験体系データおよび測定値を記憶する実験結果記憶ステップと、模擬実験体系物理量演算部が、前記模擬実験の体系についてシミュレーションを行い模擬実験体系における物理量を算出する模擬実験体系シミュレーションステップと、模擬実験体系演算値記憶部が、前記模擬実験体系物理量演算部が算出した模擬実験体系における物理量の計算値を記憶する模擬実験体系計算値記憶ステップと、前記模擬実験体系計算値記憶ステップの後に、計算値・測定値相対誤差演算部が、前記実験結果記憶ステップで記憶された測定値と、前記模擬実験体系計算値記憶ステップで記憶されたシミュレーション結果の計算値に基づいて両者の相対誤差ＥＰを算出し、前記計算値・測定値相対誤差演算部が算出した前記相対誤差ＥＰを記憶するＥＰ算出ステップと、前記計算値・測定値相対誤差記憶部が、前記相対誤差ＥＰを記憶するＥＰ記憶ステップと、模擬実験体系感度係数ベクトル演算部が、前記模擬実験（ｉ＝１〜ｎ）のそれぞれの体系の感度係数ベクトルである模擬実験体系感度係数ベクトルＳｉを算出するＳｉ算出ステップと、模擬実験体系感度係数ベクトル記憶部が、前記模擬実験体系感度係数ベクトルＳｉを記憶するＳｉ記憶ステップと、目的体系物理量演算部が、前記目的体系の物理量の計算結果Ｒ０を算出するＲ０算出ステップと、目的体系演算値記憶部が、前記計算結果Ｒ０を記憶するＲ０記憶ステップと、目的体系感度係数ベクトル演算部が、前記目的体系の入力パラメータに対する感度係数ベクトルである目的体系感度係数ベクトルＳＲを算出するＳＲ算出ステップと、目的体系感度係数ベクトル記憶部が前記目的体系感度係数ベクトルＳＲを記憶するＳＲ記憶ステップと、感度係数ベクトル２次形式第１種類演算部が、前記相対誤差ＥＰに基づいて第１の感度係数ベクトル２次形式を算出する２次形式第１種類演算ステップと、感度係数ベクトル２次形式第２種類演算部が、前記相対誤差ＥＰに基づいて第２の感度係数ベクトル２次形式を算出する２次形式第２種類演算ステップと、感度係数ベクトル２次形式第３種類演算部が、前記相対誤差ＥＰ、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲに基づいて第３の感度係数ベクトル２次形式を算出する２次形式第３種類演算ステップと、感度係数ベクトル２次形式第４種類演算部が、前記相対誤差ＥＰ、前記模擬実験体系感度係数ベクトルＳｉおよび前記目的体系感度係数ベクトルＳＲに基づいて第４の感度係数ベクトル２次形式を算出する２次形式第４種類演算ステップと、線形結合定数演算部が、模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、前記模擬実験体系感度係数ベク���ルＳｉおよび前記目的体系感度係数ベクトルＳＲに基づき算出する線形結合定数演算ステップと、前記目的体系物理量不確かさ演算部が、前記第１の感度係数ベクトル２次形式、前記第２の感度係数ベクトル２次形式、前記第３の感度係数ベクトル２次形式、前記第４の感度係数ベクトル２次形式、および前記線形結合定数αｉに基づいて、入力パラメータ起因物理量不確かさの二乗（ΔＺ／Ｚ）２を算出する不確かさ算出ステップと、を有することを特徴とする。
【発明の効果】
【００４４】
本発明の実施形態によれば、放射線治療器の計算機シミュレーションを補完して治療精度を向上させる照射線量計算誤差推定装置、照射線量計算誤差推定方法を提供することが可能となる。
【図面の簡単な説明】
【００４５】
【図１】第１の実施形態に係る照射線量計算誤差推定装置の構成を示すブロック図である。
【図２】第１の実施形態に係る照射線量計算誤差推定装置の演算部の構成を示すブロック図である。
【図３】第１の実施形態に係る照射線量計算誤差推定装置の記憶部の構成を示すブロック図である。
【図４】第１の実施形態に係る照射線量計算誤差推定装置の演算部における入力パラメータ共分散誤差行列評価部の構成を示すブロック図である。
【図５】第１の実施形態に係る照射線量計算誤差推定装置の演算部における入力パラメータ共分散誤差行列評価部に関連する照射系の第１の構成例を示すブロック図である。
【図６】第１の実施形態に係る照射線量計算誤差推定装置の演算部における入力パラメータ共分散誤差行列評価部に関連する照射系の第２の構成例を示すブロック図である。
【図７】第１の実施形態に係る照射線量計算誤差推定方法の手順を示すフロー図である。
【図８】第１の実施形態に係る照射線量計算誤差推定方法の準備ステップの手順の詳細を示すフロー図である。
【図９】第１の実施形態に係る照射線量計算誤差推定方法の調整ステップの手順の詳細を示すフロー図である。
【図１０】第２の実施形態に係る照射線量計算誤差推定装置の演算部の構成を示すブロック図である。
【図１１】第３の実施形態に係る照射線量計算誤差推定装置の演算部の構成を示すブロック図である。
【図１２】第３の実施形態に係る照射線量計算誤差推定装置の記憶部の構成を示すブロック図である。
【図１３】第３の実施形態に係る照射線量計算誤差推定方法の準備ステップの手順の詳細を示すフロー図である。
【図１４】第３の実施形態に係る照射線量計算誤差推定方法の調整ステップの手順の詳細を示すフロー図である。
【発明を実施するための形態】
【００４６】
以下、図面を参照して、本発明の実施形態に係る照射線量計算誤差推定装置および照射線量計算誤差推定方法について説明する。ここで、互いに同一または類似の部分には、共通の符号を付して、重複説明は省略する。
【００４７】
［第１の実施形態］
図１は、第１の実施形態に係る照射線量計算誤差推定装置の構成を示すブロック図である。照射線量計算誤差推定装置１０は、医療用放射線の照射のために、治療対象である目的体系の構成および目的体系と医療用放射線との相互作用の特性を計算機２０上に表現したモデルを用いてシミュレーションした照射線量計算の結果で得られる物理量に含まれる誤差をその対象を模擬した模擬実験の結果を用いて推定する装置である。
【００４８】
照射線量計算誤差推定装置１０は、演算部１１０、制御部１３０、記憶部１４０、入力装置１６０および出力装置１７０を有する計算機を用いた装置である。演算部１１０および制御部１３０は、計算機のＣＰＵ（中央演算処理装置）１００内のそれぞれの機能部分である。演算部１１０、制御部１３０および記憶部１４０は、互いにバス３０を介して信号授受可能に構成されている。
【００４９】
記憶部１４０は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）あるいはＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）などの半導体素子を主に用いた主記憶装置の場合、および、磁気ドラム、光ディスクなどの補助記憶装置の場合の両者を含むものとする。制御部１３０は、演算部１１０、記憶部１４０での処理を管理する。制御部１３０の入力制御部１３１は入力装置１６０での処理を、また、出力制御部１３２は出力装置１７０での処理を管理する。
【００５０】
図２は、第１の実施形態に係る照射線量計算誤差推定装置１０の演算部１１０の構成を示すブロック図である。演算部１１０は、模擬実験体系物理量演算部１１１、計算値・測定値相対誤差演算部１１２、模擬実験体系感度係数ベクトル演算部１１３、目的体系感度係数ベクトル演算部１１４、入力パラメータ共分散誤差行列評価部３００、線形結合定数演算部１１６、代表性因子演算部１１７、目的体系物理量不確かさ演算部１１８、目的体系物理量演算部１１９を有する。
【００５１】
図３は、第１の実施形態に係る照射線量計算誤差推定装置１０の記憶部１４０の構成を示すブロック図である。記憶部１４０は、模擬実験体系測定値記憶部１４１、模擬実験体系演算値記憶部１４２、計算値・測定値相対誤差記憶部１４３、模擬実験体系感度係数ベクトル記憶部１４４、目的体系感度係数ベクトル記憶部１４５、入力パラメータ共分散誤差行列記憶部１４６、線形結合定数記憶部１４７、代表性因子記憶部１４８、目的体系演算値記憶部１４９を有する。
【００５２】
以下、図２に示す演算部１１０、および図３に示す記憶部１４０の各要素について説明する。
【００５３】
模擬実験体系測定値記憶部１４１は、実施したそれぞれの模擬実験で得られた着目する物理量Ｒの測定値を、入力装置１６０が受け入れた後に、これを記憶する。この模擬実験ｉで得られた物理量Ｒの測定値をＥｉ（ｉ＝１，２，・・・，ｎ）とする。なお、模擬実験は、同一のものはなく、互いに少なくとも一部が異なるものとする。
【００５４】
模擬実験体系物理量演算部１１１は、評価の目的とする目的体系を模擬するために行った単数または複数の模擬実験体系について、物理現象を記述したモデル方程式を用いて対象とする物理量Ｒを算出する。模擬実験体系演算値記憶部１４２は、模擬実験体系物理量演算部１１１によって算出された物理量Ｒの解析値Ｃｉを記憶する。
【００５５】
計算値・測定値相対誤差演算部１１２は、単数または複数の模擬実験のそれぞれで得られた所定現象の物理量Ｒの測定値Ｅｉから測定誤差を除いた値と、模擬実験体系物理量演算部１１１でそれぞれの模擬実験をモデル方程式で確認計算を行い得られた物理量Ｒの解析値Ｃｉとの相対誤差Ｅｐｉ（ΔＥｉ／Ｅｉ＝（Ｃｉ―Ｅｉ）／Ｅｉ）を算出する。計算値・測定値相対差記憶部１４３は、計算値・測定値相対誤差演算部１１２によって算出された相対誤差Ｅｐｉ（ｉ＝１，２，・・・，ｎ）を記憶する。
【００５６】
模擬実験体系感度係数ベクトル演算部１１３は、模擬実験体系について物理現象を記述したモデル方程式を用いて模擬実験体系の物理量Ｒを算出する場合のそれぞれの入力パラメータの値の単位変化に対して計算により得られる物理量Ｒの値が変化する割合を示す模擬実験体系感度係数ベクトルＳｉ（ｉ＝１，２，・・・，ｎ）を算出する。なお、模擬実験体系感度係数ベクトルＳｉを総称して、模擬実験体系感度係数ベクトルＳＥと呼ぶことがある。
【００５７】
すなわち、模擬実験体系について物理現象を記述したモデル方程式において用いるパラメータをｘｊ（ｊ＝１，２、・・・、ｍ）とすると、ｘｊに関する物理量Ｒの模擬実験体系感度係数ベクトルのｊ番目の要素ｓｊは、次の式（１）で与えられる。
【００５８】
【数６】
0006490535.tif 000008
【００５９】
模擬実験体系ｊ（ｊ＝１，２、・・・、ｎ）について模擬実験体系感度係数ベクトルＳｊが求められる。模擬実験体系感度係数ベクトルＳｊの要素ｓｊ（ｊ＝１，２、・・・、ｍ）は、ｍ個でありしたがって模擬実験体系感度係数ベクトルＳｉは、ｍ次元ベクトルである。
【００６０】
模擬実験体系感度係数ベクトル記憶部１４４は、模擬実験体系感度係数ベクトル演算部１１３によって算出された模擬実験体系感度係数ベクトルＳｉ（ｉ＝１，２、・・・、ｎ）を記憶する。
【００６１】
目的体系感度係数ベクトル演算部１１４は、目的体系について物理現象を記述したモデル方程式を用いて目的体系の物理量Ｒを算出する場合のそれぞれの入力パラメータの値の単位変化に対して計算により得られる物理量Ｒの値が変化する割合を示す目的体系感度係数ベクトルＳＲを算出する。
【００６２】
すなわち、目的体系について物理現象を記述したモデル方程式において用いるパラメータをｘｊ（ｊ＝１，２、・・・、ｍ）とすると、ｘｊに関する物理量Ｒの目的体系感度係数ベクトルＳＲの成分ＳＲＪ（ｊ＝１，２、・・・、ｍ）は、次の式（２）で与えられる。目的体系感度係数ベクトルＳＲは、ｍ次元ベクトルである。
【００６３】
【数７】
0006490535.tif 000009
【００６４】
目的体系感度係数ベクトル記憶部１４５は、目的体系感度係数ベクトル演算部１１４によって算出された目的体系感度係数ベクトルＳＲを記憶する。
【００６５】
入力パラメータ共分散誤差行列評価部３００は、モデル方程式を解く際の各入力パラメータの誤差、すなわち入力パラメータ共分散誤差行列Ｗの各要素が設定された結果を評価し、設定する。設定される入力パラメータ共分散誤差行列Ｗは、目的体系と模擬実験体系で使用する同じ構造の入力値の不確かさを表す共分散行列であり、シミュレーションに用いるモデルへのパラメータの入力値の不確かさの割合の関係を示している。
【００６６】
感度係数ベクトルは、パラメータの数ｍを次元とするｍ次元ベクトルであるため、入力パラメータ共分散誤差行列Ｗは、ｍ×ｍの行列である。入力パラメータ共分散誤差行列Ｗの各要素も単位の無い無次元値である。
【００６７】
入力パラメータ共分散誤差行列Ｗの対角成分は、それぞれのパラメータの相対誤差を２乗した値が収納されている。また、非対角成分（ｉ，ｊ）はパラメータｉの相対誤差とパラメータｊの相対誤差の積が収納されている。本実施形態におけるアルゴリズムで用いる入力パラメータ共分散誤差行列Ｗの各要素の大きさは、相互の大きさの関係が正しければよく、絶対値は問題ではない。また場合によって相関係数（Ｃｏｒｒｅｌａｔｉｏｎ  Ｃｏｅｆｆｉｃｉｅｎｔ）の考え方に沿って対角線成分を１に規格化して非対角成分の値を分かりやすくする場合がある。
【００６８】
一般にｉ番目の模擬実験体系の計算値に関する不確かさは、通常２次以上の微小量を無視してサンドイッチルール（Ｓａｎｄｗｉｃｈ  ｒｕｌｅ）と呼ばれる次の式（３）のような感度係数ベクトルの２次形式（以下、感度係数ベクトル２次形式）で表されるものとする。本実施形態は、このサンドイッチルールが成立することを前提として構成されている。
（ΔＲｉ／Ｒｉ）２＝ＳｉＴＷＳｉ                    …（３）
【００６９】
ここで、（ΔＲｉ／Ｒｉ）は、入力パラメータに含まれる不確かさに起因して模擬実験体系ｉの着目している物理量Ｒｉの計算値に含まれる不確かさの相対値、Ｓｉは模擬実験体系ｉの感度係数ベクトル、Ｗは入力パラメータ共分散誤差行列、Ｔは、ベクトルの転置を示す。
【００７０】
図４は、第１の実施形態に係る照射線量計算誤差推定装置の演算部における入力パラメータ共分散誤差行列評価部の構成を示すブロック図である。入力パラメータ共分散誤差行列評価部３００は、集積データ平均値演算部３０１、および入力パラメータ共分散誤差行列設定部３０２を有する。
【００７１】
図５は、入力パラメータ共分散誤差行列評価部に関連する照射系の第１の構成例を示すブロック図である。図５に示す、高速スキャニング照射方式照射装置５１０は、スキャニング電磁石５１１、およびレンジフィルタ５１２を有する。
【００７２】
図示しないサイクロトロンにより所定のエネルギーまで加速された荷電粒子のビームは、スキャニング電磁石５１１により水平方向の方角および鉛直方向の方角を調整される。スキャニング電磁石５１１により方角を調整された荷電粒子のビームは、レンジフィルタ５１２により飛程を微調整された後に、がん病巣６００に照射される。
【００７３】
図６は、入力パラメータ共分散誤差行列評価部に関連する照射系の第２の構成例を示すブロック図である。図６に示す、ブロードビーム方式照射装置５２０は、ワブラー電磁石５２１、散乱体５２２、多葉コリメータ５２３、補償フィルタ５２４、患者コリメータ５２５を有する。
【００７４】
図示しないサイクロトロンにより所定のエネルギーまで加速された荷電粒子のビームは、ワブラー電磁石５２１により、照射量の平坦化のために、所定の角度範囲で所定の順序に従った角度方向に変化するように制御される。次に、散乱体５２２により、さらに平坦化される。
【００７５】
次に多葉コリメータ５２３により、がんの形状に対応した分布に調整され、補償フィルタ５２４を経て、患者コリメータ５２５によりがん病巣の目的部分に収束するように調整された後に、がん病巣６００に照射される。
【００７６】
以上のように、がん病巣６００に荷電粒子のビームを照射する際の、荷電粒子のビームを正確にがん病巣６００の目的位置に照射する上で、照射の実際の結果と目的とする結果との誤差に影響を与える重要なパラメータが調整され、その結果が、全体の精度となって現れる。
【００７７】
重要なパラメータとしては、荷電粒子のビームのエネルギー、荷電粒子のビームの放射線量、荷電粒子のビームの方向、荷電粒子のビームの収束具合などがある。なお、以上のパラメータは、例として示したものであり、これに限定されるものではない。荷電粒子のビームを正確にがん病巣６００の目的位置に照射する上で、その誤差に影響を与える重要なパラメータであれば考慮すべきである。
【００７８】
また、これらのパラメータに直接影響する要因としては、荷電粒子のビームの照射のＯＮ／ＯＦＦのタイミング、スキャニング電磁石５１１またはワブラー電磁石５２１の電流値、スキャニング電磁石５１１またはワブラー電磁石５２１のＯＮ／ＯＦＦのタイミング、各要素の設置位置、レンジフィルタ５１２などにより微調整される飛程などがある。各要素の設置位置とは、電磁石、フィルタ、あるいは照射の目的であるがん病巣６００の位置を意味する。
【００７９】
これらの要因が前記のパラメータにどのように影響するかについては、それぞれ評価が可能である。たとえば、スキャニング電磁石５１１またはワブラー電磁石５２１の電流値と、荷電粒子のビームの方向との関係は、電磁気学的に計算可能である。したがって、これらの要因の誤差が分かれば、前記の重要なパラメータの誤差が求められる。
【００８０】
図４の集積データ平均値演算部３０１は、これらの要因についての誤差を集積して、平均化する。
【００８１】
たとえば、荷電粒子のビームの照射のＯＮ／ＯＦＦのタイミングの誤差は、ＯＮ／ＯＦＦに係る処理系、信号系およびスイッチング部分のそれぞれの誤差の総合誤差である。したがって、処理系、信号系、およびスイッチング部分のそれぞれの誤差についてフィールドデータを集積することが可能である。
【００８２】
また、たとえば、がん病巣６００の位置の誤差については、照射時の患者の位置の誤差、事前のＣＴ検査の精度なども含めたがん病巣６００の患者の体内における位置の評価の誤差などの総合的な誤差である。たとえば、がん病巣６００の患者の体内における位置の評価の誤差については、治療後の画像診断等の確認によって、経験的に誤差がどの程度あるかデータが集積され、あるいはデータの集積が可能である。
【００８３】
図４の、入力パラメータ共分散誤差行列設定部３０２は、集積データ平均値演算部３０１で求められた各要因の誤差の平均値に基づいて、前記の重要なパラメータの誤差を算出する。
【００８４】
ここで、前記の重要なパラメータ同士の関係を見てみる。重要なパラメータとして、荷電粒子のビームのエネルギー、荷電粒子のビームの放射線量、荷電粒子のビームの方向、荷電粒子のビームの収束具合を挙げた。
【００８５】
荷電粒子のビームのエネルギーは、シンクロトロンによる加速と、照射系における微調整の結果による。荷電粒子のビームの放射線量は、シンクロトロンにおける荷電粒子の量と照射系への取り出し部の設定による。荷電粒子のビームの方向は、電磁石の電流およびオンオフのタイミングによる。荷電粒子のビームの収束具合は、取出し系の性能による。
【００８６】
このように、重要なパラメータの誤差に影響する要因は、それぞれのパラメータごとに異なっている。この結果、重要なパラメータの誤差は、相互に関連せずに独立と考えても誤差評価上は問題がないと考えられる。
【００８７】
この結果、入力パラメータ共分散誤差行列Ｗの成分は、対角成分であるそれぞれのパラメータ自身の誤差の二乗のみとなる。すなわち、入力パラメータ共分散誤差行列Ｗは対角行列となる。なお、仮に、相互に関連する部分があったとしても、対角行列と近似することにより、演算を簡素化することができる。
【００８８】
入力パラメータ共分散誤差行列記憶部１４６は、入力パラメータ共分散誤差行列評価部３００により設定された入力パラメータ共分散誤差行列Ｗを記憶する。
【００８９】
線形結合定数演算部１１６は、模擬実験体系感度係数ベクトルＳｉを線形結合した線形結合ベクトルＳｋの線形結合定数αｉを、感度係数ベクトル２次形式ＳｋＴＷＳｋが感度係数ベクトル２次形式ＳＲＴＷＳＲと等しくかつ目的体系感度係数ベクトルＳＲと線形結合ベクトルＳｋとのなす角が最小になるように求める。
【００９０】
まず、ｉ番目の模擬実験体系の感度係数ベクトルＳｉを用いて、式（４）を満足する感度係数ベクトルの線形結合（感度係数線形結合ベクトル）Ｓｋを、次のように定義し、算出する。
【数８】
0006490535.tif 000010
【００９１】
ここで、αｉは任意の定数であり、これを線形結合定数と呼ぶ。また、αｉ（ｉ＝１，２，…，ｎ）で構成されるベクトルを以下で線形結合定数ベクトルαと呼ぶこととする。
【００９２】
次に、線形結合定数αｉを、次の式（５）のように、線形結合ベクトルＳｋの２次形式の値が、目的体系感度係数ベクトルＳＲの２次形式の値と等しくなり、目的体系感度係数ベクトルＳＲと線形結合ベクトルＳｋとのなす角θが最小になるように求める。
ＳｋＴＷＳｋ＝ＳＲＴＷＳＲ      …（５）
ここで、Ｔはベクトルや行列の転置を示す記号である。この式（５）は、相対誤差の総量が同じであるという条件を表している。
【００９３】
目的体系感度係数ベクトルＳＲと線形結合ベクトルＳｋとのなす角θが最小になるように求める方法はいくつか想定される。
【００９４】
たとえば、（ＳＲ−Ｓｋ）ＴＷ（ＳＲ−Ｓｋ）の絶対値を最小とすることにより、目的体系感度係数ベクトルＳＲと線形結合ベクトルＳｋとのなす角θを最小にすることができる。（ＳＲ−Ｓｋ）ＴＷ（ＳＲ−Ｓｋ）の値が零（０）となった場合、すなわち、次の式（６）を満足する場合は、自動的に式（７）が成立する。
（ＳＲ−Ｓｋ）ＴＷ（ＳＲ−Ｓｋ）＝０      …（６）
ＳｋＴＷＳＲ／｛（ＳｋＴＷＳｋ）１／２（ＳＲＴＷＳＲ）１／２｝＝１
…（７）
この場合、式（７）は、後述する代表性因子ＲＦが１となっていることを示している。
【００９５】
式（７）の左辺の分子は、入力パラメータ共分散誤差行列Ｗを考慮した線形結合ベクトルＳｋと目的体系感度係数ベクトルＳＲとの内積を示している。式（７）の左辺の分母は、入力パラメータ共分散誤差行列Ｗを考慮した線形結合ベクトルＳｋと目的体系感度係数ベクトルＳＲのベクトルの大きさになっている。
【００９６】
（ＳＲ−Ｓｋ）ＴＷ（ＳＲ−Ｓｋ）の絶対値を最小とする線形結合定数αｉは、たとえばＬａｇｒａｎｇｅの未定定数法を用いて求めることができる。すなわち、以下のとおりである。
ｄを、ｄ＝ＳｋＴＷＳｋ        …（８）
とし、この式を束縛条件とし、Ｌａｇｒａｎｇｅの未定定数λを用いて、以下のＬを導出する式（９）を作る。
Ｌ＝｛（ＳＲ−Ｓｋ）ＴＷ（ＳＲ−Ｓｋ）｝＋λ｛ＳｋＴＷＳｋ−ｄ｝
Ｌ＝｛２ｄ−２（ＳｋＴＷＳＲ）｝＋λ｛ＳｋＴＷＳｋ−ｄ｝    …（９）
この式（９）に関して極値条件から、
【数９】
0006490535.tif 000011
として、αｉの条件式（連立方程式）を作成する。この条件式から、λとαｉ（ｉ＝１，２，…，ｎ）を決定する。
【００９７】
今、ｒｉｊ＝ＳｉＴＷＳｊ、ｒＲｉ＝ＳｉＴＷＳｊとすれば、Ｌａｇｒａｎｇｅ未定定数λを用いて、λとαｉ（ｉ＝１，２，…，ｎ）とを解く方程式は次の式（１１）のように書ける。
Ｒ・λ・α＝ｒ      …（１１）
ここで、Ｒはｒｉｊで構成される行列、αはαｉで構成される列ベクトル、ｒはｒＲｉで構成される列ベクトルである。
【００９８】
Ｒは、ｎを実験の数としてｎ行ｎ列の行列であるが、例として、Ｒが４行４列の行列である場合について説明する。以下、各実施形態を含めて同様である。この場合、上述の方程式は、次の式（１２）となる。
【００９９】
【数１０】
0006490535.tif 000012
【０１００】
この式を解いて、λαｉ（ｉ＝１，２，…，ｎ）を求める。
すなわち、
【０１０１】
【数１１】
0006490535.tif 000013
により、λαｉ（ｉ＝１，２，…，ｎ）を算出する。
加えて、αｉは、
【０１０２】
【数１２】
0006490535.tif 000014
が成り立つように規格化される。
このようにして、λとαｉ（ｉ＝１，２，…，ｎ）が決定される。
【０１０３】
なお、λとαｉ（ｉ＝１，２，…，ｎ）は、正負で与えられる。原理的には、２ｄ−２・（ＳｋＴＷＳＲ）が小さくなるほうの符号を選ぶことになる。しかし、後述する代表性因子ＲＦがλに等しいことがわかっており、代表性因子は１に近い値であるべきなので、λを負の値とすることは目的と合致しない。そこで、λとしては常に正の値を選択すればよいことになる。
【０１０４】
線形結合定数記憶部１４７は、線形結合定数演算部１１６で算出された線形結合定数αｉ（ｉ＝１，２，…，ｎ）を記憶する。
【０１０５】
代表性因子演算部１１７は、次の式（１５）により代表性因子ＲＦ（Ｒｅｐｒｅｎｔａｔｉｖｅ  Ｆａｃｔｏｒ）を算出する。
ＲＦ＝ＳｋＴＷＳＲ／｛（ＳｋＴＷＳｋ）１／２（ＳＲＴＷＳＲ）１／２｝
…（１５）
【０１０６】
代表性因子ＲＦは、目的体系感度係数ベクトルＳＲと線形結合ベクトルＳｋとのなす角をθとしたときの、ｃｏｓθを表している。式（１５）の右辺の分子は、入力パラメータ共分散誤差行列Ｗを考慮した線形結合ベクトルＳｋと目的体系感度係数ベクトルＳＲとの内積を示している。式（１５）の右辺の分母は、入力パラメータ共分散誤差行列Ｗを考慮した線形結合ベクトルＳｋと目的体系感度係数ベクトルＳＲのベクトルの大きさになっている。
【０１０７】
代表性因子ＲＦの値が１のとき、すなわちｃｏｓθ＝１のときに、θ＝０となり、線形結合ベクトルＳｋと目的体系感度係数ベクトルＳＲとが重なった状態、すなわち、数理的に完全に一致していることになる。逆に何の関係もない場合は、代表性因子は０となる。
【０１０８】
代表性因子記憶部１４８は、代表性因子演算部１１７で算出された代表性因子ＲＦを記憶する。
【０１０９】
目的体系物理量不確かさ演算部１１８は、次の式（１６）により模擬実験ｉにより求められる物理量と模擬実験ｊで求められる物理量の間の相関を表わすＣＯＲｉｊを算出し、このＣＯＲｉｊを用いて、次の式（１７）により相対誤差Ｅｐｉ（すなわち、入力パラメータ起因物理量不確かさ）の二乗（ΔＺ／Ｚ）２を算出する。
【０１１０】
【数１３】
0006490535.tif 000015
【０１１１】
【数１４】
0006490535.tif 000016
【０１１２】
目的体系物理量不確かさ演算部１１８は、相対誤差Ｅｐｉの二乗（ΔＺ／Ｚ）２を最小とするように、線形結合定数αｉ（ｉ＝１，２，…，ｎ）を調整する。
【０１１３】
目的体系物理量演算部１１９は、目的体系における物理量Ｒを算出する。目的体系演算値記憶部１４９は、目的体系物理量演算部１１９によって算出された目的体系における物理量Ｒを記憶する。
【０１１４】
図７は、第１の実施形態に係る照射線量計算誤差推定方法の手順を示すフロー図である。最初のステップは、準備ステップＳ１００である。準備ステップＳ１００の後に、調整ステップＳ２００となる。準備ステップＳ１００および調整ステップＳ２００の詳細を、以下に説明する。
【０１１５】
図８は、第１の実施形態に係る照射線量計算誤差推定方法の準備ステップＳ１００の手順の詳細を示すフロー図である。まず、実施する、あるいは実施した模擬実験のそれぞれについて、模擬実験体系物理量演算部１１１が、その体系、条件のもとでモデル方程式に基づいて解析を実施する（ステップＳ１０１）。この解析により、それぞれの模擬実験の場合の、着目する物理量Ｒについての解析値Ｃｉを算出する。ただし、ｉ＝１，２，・・，ｎ（ｎは模擬実験数）。これらの解析値Ｃｉは、模擬実験体系演算値記憶部１４２により記憶される。
【０１１６】
次に、模擬実験の結果に基づき入力される着目する物理量Ｒについての測定値Ｅｉを、入力装置１６０が受け入れる（ステップＳ１０２）。これらの測定値Ｅｉは、模擬実験体系測定値記憶部１４１が記憶する。
【０１１７】
計算値・測定値相対誤差演算部１１２は、模擬実験体系演算値記憶部１４２に記憶された解析値Ｃｉと、模擬実験体系測定値記憶部１４１に記憶された測定値Ｅｉとに基づき、相対誤差Ｅｐｉ（ΔＥｉ／Ｅｉ＝（Ｃｉ―Ｅｉ）／Ｅｉ＝（Ｃｉ／Ｅｉ）−１）を算出する。また、計算値・測定値相対誤差記憶部１４３が、相対誤差の二乗（Ｅｐｉ）２を記憶する（ステップＳ１０３）。
【０１１８】
また、目的体系感度係数ベクトル演算部１１４が、目的体系感度係数ベクトルＳＲを算出する。また、模擬実験体系感度係数ベクトル演算部１１３が、模擬実験体系感度係数ベクトルＳｉ（ｉ＝１，２，・・，ｎ）を算出する（ステップＳ１０４）。算出された目的体系感度係数ベクトルＳＲは、目的体系感度係数ベクトル記憶部１４５に記憶される。また、算出された模擬実験体系感度係数ベクトルＳｉは、模擬実験体系感度係数ベクトル記憶部１４４に記憶される。
【０１１９】
また、入力パラメータ共分散誤差行列評価部３００は、入力パラメータ共分散誤差行列Ｗを設定する（ステップＳ１０５）。設定された入力パラメータ共分散誤差行列Ｗは、入力パラメータ共分散誤差行列記憶部１４６に記憶される。なお、ステップＳ１０５は、ステップＳ１０１に先立って、行われることでもよい。
【０１２０】
次に、線形結合定数演算部１１６は、前述の式（４）に基づいて、模擬実験体系感度係数ベクトルＳｉの線形結合としての線形結合ベクトルＳｋを設定する（ステップＳ１０６）。
【０１２１】
以上が、準備ステップＳ１００の手順の詳細である。
【０１２２】
図９は、第１の実施形態に係る照射線量計算誤差推定方法の調整ステップＳ２００の手順の詳細を示すフロー図である。
【０１２３】
まず、線形結合定数演算部１１６が、線形結合ベクトルＳｋの線形結合定数αｉを算出する（ステップＳ２１１）。具体的には、前述の式（１３）を用いて、行列方程式を解き、線形結合定数ベクトルα、すなわち、αｉ（ｉ＝１，２，…，ｎ）を求める。なお、この段階では、αｉ（ｉ＝１，２，…，ｎ）の相対値および、Ｌａｇｒａｎｇｅ未定定数λが求められる。
【０１２４】
次に、線形結合定数演算部１１６が、前述の式（５）を満足し、かつ目的体系感度係数ベクトルＳＲと線形結合ベクトルＳｋとのなす角θが最小になるように、線形結合定数αｉ（ｉ＝１，２，…，ｎ）を算出する（ステップＳ２１２）。線形結合定数演算部１１６で算出された線形行列ベクトル、およびこれに基づく線形結合ベクトルＳｋは、線形結合定数記憶部１４７に記憶される。
【０１２５】
次に、代表性因子演算部１１７が、線形結合定数記憶部１４７に記憶された線形結合ベクトルＳｋ、目的体系感度係数ベクトル記憶部１４５に記憶された目的体系感度係数ベクトルＳＲ、および入力パラメータ共分散誤差行列記憶部１４６に記憶された入力パラメータ共分散誤差行列Ｗに基づいて、前述の式（１５）を用いて代表性因子ＲＦを算出する（ステップＳ２１３）。この結果、目的体系感度係数ベクトルＳＲを線形結合ベクトルＳｋが、どの程度模擬しているかを判断することができる。算出された代表性因子ＲＦは、代表性因子記憶部１４８に記憶される。
【０１２６】
次に、目的体系物理量不確かさ演算部１１８は、前述の式（１６）および式（１７）を用いて、相対誤差の二乗（ΔＺ／Ｚ）２を算出する（ステップＳ２１４）。次に、目的体系物理量不確かさ演算部１１８は、相対誤差の二乗（ΔＺ／Ｚ）２を求める（ステップＳ２１５）。詳細には、相対誤差の二乗（ΔＺ／Ｚ）２から正負の判断を行い（ΔＺ／Ｚ）あるいは−（ΔＺ／Ｚ）を求める。その後、この値を０に近付けるように単数、あるいは複数の制御パラメータを調整する。
【０１２７】
以上のように、本実施形態によれば、放射線治療器の計算機シミュレーションにおける誤差を推定し、かつ、誤差を最小することにより、治療精度を向上させることができる。
【０１２８】
［第２の実施形態］
図１０は、第２の実施形態に係る照射線量計算誤差推定装置の演算部の構成を示すブロック図である。本実施形態は、第１の実施形態の変形である。第１の実施形態における入力パラメータ共分散誤差行列評価部３００に代えて、本第２の実施形態における演算部２１０は、入力パラメータ共分散誤差行列評価部３０５を有する。以下、入力パラメータ共分散誤差行列評価部３０５について説明する。
【０１２９】
いまパラメータの数をｎとする。また、模擬実験の回数が同じｎ回である場合を考える。以下、ｎをたとえば４回とした場合を示す。
【０１３０】
模擬実験体系感度係数ベクトルＳｉを、Ｓｉ＝（ｓｉ１，ｓｉ２，ｓｉ３，ｓｉ４）（ただし、ｉ＝１，２，３，４）として、入力パラメータ共分散誤差行列Ｗを、次の式（１８）のように仮定する。
【０１３１】
【数１５】
0006490535.tif 000017
【０１３２】
また測定値と計算値の相対誤差の二乗ｑｉ２を、次の式（１９）で定義すれば、式（２０）が成立し、具体的には、式（２１）が得られる。ただし、ｒｉＥは基準値、Δｒｉは基準値からのずれを表す。
【０１３３】
【数１６】
0006490535.tif 000018
【０１３４】
【数１７】
0006490535.tif 000019
【０１３５】
この式から、次の式（２１）が得られる。
【０１３６】
【数１８】
0006490535.tif 000020
【０１３７】
式（２１）において、右辺の行列部分の逆行列を両辺に乗じると、次の式（２２）により入力パラメータ共分散誤差行列Ｗの各対角要素が求められる。
【０１３８】
【数１９】
0006490535.tif 000021
【０１３９】
入力パラメータ共分散誤差行列評価部３０５によって算出された入力パラメータ共分散誤差行列Ｗは、入力パラメータ共分散誤差行列記憶部１４６により記憶される。
【０１４０】
［第３の実施形態］
図１１は、第３の実施形態に係る照射線量計算誤差推定装置の演算部の構成を示すブロック図である。また、図１２は、記憶部の構成を示すブロック図である。
【０１４１】
本実施形態は、第１の実施形態の変形である。第２の実施形態は、模擬実験の回数がパラメータの数と同じ場合であるが、本第３の実施形態は、模擬実験回数がパラメータの数より小さい場合である。
【０１４２】
第１および第２の実施形態においては、入力パラメータ共分散誤差行列Ｗを求め、求められた入力パラメータ共分散誤差行列Ｗと、それぞれの感度係数から、ＳｉＴＷＳｉ、ＳｉＴＷＳｊ、ＳＲＴＷＳＲ、およびＳｉＴＷＳＲを求めている。一方、本第３の実施形態においては、ＳｉＴＷＳｉ、ＳｉＴＷＳｊ、ＳＲＴＷＳＲ、およびＳｉＴＷＳＲに対応する値を直接求める点で異なる。
【０１４３】
本第３の実施形態においては、演算部３１０は、第１の実施形態の演算部１１０における入力パラメータ共分散誤差行列評価部３００に相当する部分を有さない。これに代えて、本第３の実施形態における演算部３１０は、感度係数ベクトル２次形式第１種類演算部３１１、感度係数ベクトル２次形式第２種類演算部３１２、感度係数ベクトル２次形式第３種類演算部３１３、および感度係数ベクトル２次形式第４種類演算部３１４を有する。
【０１４４】
これに対応して、本第３の実施形態における記憶部３４０は、入力パラメータ共分散誤差行列記憶部１４６に相当する部分を有さず、感度係数ベクトル２次形式第１種類記憶部３４１、感度係数ベクトル２次形式第２種類記憶部３４２、感度係数ベクトル２次形式第３種類記憶部３４３、および感度係数ベクトル２次形式第４種類記憶部３４４を有する。
【０１４５】
感度係数ベクトル２次形式第１種類演算部３１１は、同一模擬実験に対するｉ番目の模擬実験の感度係数ベクトル２次形式ＳｉＴＷＳｉを算出する。
【０１４６】
まず、ｈを模擬実験回数として、模擬実験体系の測定値Ｅｉと模擬実験体系に関する解析値Ｃｉとの相対誤差の二乗ｑｉ２を次の式（２３）により算出する。
【０１４７】
【数２０】
0006490535.tif 000022
【０１４８】
この結果、まず、ＳｉＴＷＳｉが、ＳｉＴＷＳｉ＝ｑｉ２により得られる。感度係数ベクトル２次形式第１種類記憶部３４１は、この結果を記憶する。
【０１４９】
感度係数ベクトル２次形式第２種類演算部３１２は、ｉ番目の模擬実験とｊ番目の模擬実験に対する感度係数ベクトル２次形式ＳｉＴＷＳｊを次の式（２４）により算出する。
【０１５０】
【数２１】
0006490535.tif 000023
【０１５１】
感度係数ベクトル２次形式第２種類記憶部３４２は、この結果を記憶する。
【０１５２】
感度係数ベクトル２次形式第３種類演算部３１３は、目的とする患者の体系についての感度係数ベクトルＳＲの２次形式ＳＲＴＷＳＲを算出する。
【０１５３】
まず、次の式（２５）の値の絶対値が最大となる模擬実験の番号ｉを選ぶ。その番号をｍとする。
【０１５４】
【数２２】
0006490535.tif 000024
【０１５５】
これにより選定されたｍに基づき、ＳｍＴＷＳｍ＝ｑｍ２、ここで、ｑｍ２＝（Ｅｐｍ）２＝（（Ｃｍ−Ｅｍ）／Ｅｍ）２から、次の式（２６）により、ＳＲＴＷＳＲを算出する。
【０１５６】
【数２３】
0006490535.tif 000025
【０１５７】
感度係数ベクトル２次形式第３種類記憶部３４３は、この結果を記憶する。
【０１５８】
感度係数ベクトル２次形式第４種類演算部３１４は、ｉ番目の模擬実験についての感度係数ベクトルと、目的とする患者の体系についての感度係数ベクトルＳＲの２次形式ＳＲＴＷＳＲを算出する。具体的には、式（２７）の右辺をｑＲ２と置けば、ＳｉＴＷＳＲは次の式（２７）により求められる。
【０１５９】
【数２４】
0006490535.tif 000026
【０１６０】
感度係数ベクトル２次形式第４種類記憶部３４４は、この結果を記憶する。
【０１６１】
以上のように、第１の実施形態のステップＳ１０５を実施せずに、直接、ステップＳ２１３での代表性因子ＲＦ、ステップＳ２１４での相対誤差の二乗（ΔＺ／Ｚ）２の算出に必要なＳｉＴＷＳｉ、ＳｉＴＷＳｊ、ＳＲＴＷＳＲ、およびＳｉＴＷＳＲに対応する値を直接求めることができる。
【０１６２】
図１３は、本実施形態の準備ステップの手順の詳細を示すフロー図である。第１の実施形態における準備ステップＳ１００においては、ステップＳ１０５で、目的体系に係る入力パラメータの共分散誤差行列Ｗを設定する。一方、本第３の実施形態における準備ステップＳ１３０においては、ステップＳ１３１として感度係数ベクトル２次形式第１種類ないし感度係数ベクトル２次形式第４種類を算出する。
【０１６３】
図１４は、調整ステップの手順の詳細を示すフロー図である。図１４のステップＳ２３１は、第１の実施形態の図９で示すステップＳ２１１に対応する。図１４のステップＳ２３２は、第１の実施形態の図９で示すステップＳ２１２に対応する。図１４のステップＳ２３３は、第１の実施形態の図９で示すステップＳ２１３に対応する。同様に、図１４のステップＳ２３４は、第１の実施形態の図９で示すステップＳ２１４に対応する。
【０１６４】
いずれも、第１の実施形態を示す図９と、本第３の実施形態を示す図１４とでは、同様の内容となっているが、それぞれに用いられている感度係数ベクトルの２次形式は、入力パラメータ共分散誤差行列Ｗを用いて算出するのではなく、前述のように、感度係数ベクトル２次形式第１種類演算部３１１ないし感度係数ベクトル２次形式第４種類演算部３１４で直接算出されるものである。
【０１６５】
以上のように、本実施形態によれば、入力パラメータ共分散誤差行列Ｗの設定が難しく、またパラメータの数に比べて模擬実験の回数が少ない場合でも、必要な感度係数ベクトル２次形式をそれぞれ算出して、目的とする物理量についての誤差を評価することができる。
【０１６６】
［その他の実施形態］
以上、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。また、各実施形態の特徴を組み合わせてもよい。さらに、これらの実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１６７】
１０…照射線量計算誤差推定装置、２０…計算機、３０…バス、１００…中央演算処理装置（ＣＰＵ）、１１０…演算部、１１１…模擬実験体系物理量演算部、１１２…計算値・測定値相対誤差演算部、１１３…模擬実験体系感度係数ベクトル演算部、１１４…目的体系感度係数ベクトル演算部、１１６…線形結合定数演算部、１１７…代表性因子演算部、１１８…目的体系物理量不確かさ演算部、１１９…目的体系物理量演算部、１３０…制御部、１３１…入力制御部、１３２…出力制御部、１４０…記憶部、１４１…模擬実験体系測定値記憶部、１４２…模擬実験体系演算値記憶部、１４３…計算値・測定値相対誤差記憶部、１４４…模擬実験体系感度係数ベクトル記憶部、１４５…目的体系感度係数ベクトル記憶部、１４６…入力パラメータ共分散誤差行列記憶部、１４７…線形結合定数記憶部、１４８…代表性因子記憶部、１４９…目的体系演算値記憶部、１６０…入力装置、１７０…出力装置、２１０…演算部、３００…入力パラメータ共分散誤差行列評価部、３０１…集積データ平均値演算部、３０２…入力パラメータ共分散誤差行列設定部、３０５…入力パラメータ共分散誤差行列評価部、３１０…演算部、３１１…感度係数ベクトル２次形式第１種類演算部、３１２…感度係数ベクトル２次形式第２種類演算部、３１３…感度係数ベクトル２次形式第３種類演算部、３１４…感度係数ベクトル２次形式第４種類演算部、３４１…感度係数ベクトル２次形式第１種類記憶部、３４２…感度係数ベクトル２次形式第２種類記憶部、３４３…感度係数ベクトル２次形式第３種類記憶部、３４４…感度係数ベクトル２次形式第４種類記憶部、５１０…高速スキャニング照射方式照射装置、５１１…スキャニング電磁石、５１２…レンジフィルタ、５２０…ブロードビーム方式照射装置、５２１…ワブラー電磁石、５２２…散乱体、５２３…多葉コリメータ、５２４…補償フィルタ、５２５…患者コリメータ、６００…がん病巣
【図１】
0006490535.tif 000033
【図２】
0006490535.tif 000034
【図３】
0006490535.tif 000035
【図４】
0006490535.tif 000036
【図５】
0006490535.tif 000037
【図６】
0006490535.tif 000038
【図７】
0006490535.tif 000039
【図８】
0006490535.tif 000040
【図９】
0006490535.tif 000041
【図１０】
0006490535.tif 000042
【図１１】
0006490535.tif 000043
【図１２】
0006490535.tif 000044
【図１３】
0006490535.tif 000045
【図１４】
0006490535.tif 000046